Can you lose unmapped reads in BAM file?
1
0
Entering edit mode
3 months ago
Taktak31 • 0

If you convert a BAM file to FASTQ that you know has unmapped reads, will it contain all of the original raw sequence data including the unmapped reads? Is there a chance you can lose some of the unmapped reads that was originally present in the FASTQ file? If so is dependent on whoever processed it and how can you find out that this might have occurred?

unmapped reads BAM • 291 views
ADD COMMENT
0
Entering edit mode

To verify if reads were lost in the extraction process, I think you should be able to use samtools view -c to count the number of reads in the BAM and wc -l to count the number of lines in FASTQ. The latter should be 4 * the former, as FASTQ has 4 lines per read.

ADD REPLY
0
Entering edit mode

Yes/No/Maybe. It simply depends on the parameters used by the software converting it.

ADD REPLY
0
Entering edit mode
3 months ago

If you convert a BAM file to FASTQ that you know has unmapped reads, will it contain all of the original raw sequence data including the unmapped reads?

yes for example , if 'samtools fastq' is used with option -F

-F INT only include reads with none of the FLAGS in INT present [0x900]

ADD COMMENT

Login before adding your answer.

Traffic: 1116 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