Are bam2fastq --no-aligned and samtools fastq -f 4 same?
0
0
Entering edit mode
5 months ago
Yogesh • 0

hi i am trying to explore the unmapped reads of my genome assembly. I straight forward tried samtools fastq -f 4 bla.bam > unmapped.fq but it gave me single file which was unusable for contig generation using spades. So i shifted to bam2fastq --no-aligned. Although i got things working but i am not pretty sure if both approach are same and are right to extract unmapped reads only.

samtools assembly spades fastq • 282 views
ADD COMMENT
0
Entering edit mode

1) samtools sort -n or samtools collate must be called before samtools fastq

2)

it gave me single file wh

see option -1 and -2

ADD REPLY

Login before adding your answer.

Traffic: 2495 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6