What Does Qc-Passed/Failed Reads Mean In Samtools Flagstat?
2
1
Entering edit mode
10.9 years ago
KJ Lim ▴ 140
8075928 + 0 in total (QC-passed reads + QC-failed reads)
0 + 0 duplicates
4723128 + 0 mapped (58.48%:nan%)
8075928 + 0 paired in sequencing
4037964 + 0 read1
4037964 + 0 read2
1370328 + 0 properly paired (16.97%:nan%)
1370328 + 0 with itself and mate mapped
3352800 + 0 singletons (41.52%:nan%)
0 + 0 with mate mapped to a different chr
0 + 0 with mate mapped to a different chr (mapQ>=5)

This is the flagstat result I obtained from paired-end mapped BAM file with SRHiMP2.The total paired-end reads are 114 100 948 but only 8 075 928 reads shown in the flagstat.

Could anyone kindly lights me how the QC-passed/failed reads come to?

Thanks a lot. Have a nice day.

samtools qc • 11k views
ADD COMMENT
2
Entering edit mode
10.9 years ago

In the SAM specification flag 0x200 (512) stands for "not passing quality controls". This is a platform/aligner specific flag that does not have a unique definition.

I believe that the flagstat output counts these flags. Most aligners that I have seen do not actually make use of this flag so most of the time this value is not relevant.

ADD COMMENT
0
Entering edit mode

Thanks for your prompt replied, Istvan.

ADD REPLY
0
Entering edit mode

There are several mysterious flags in BAMs

ADD REPLY
0
Entering edit mode
10.9 years ago
Lim • 0

I'm confuse with the flagstat result of BAM file.

The total paired-end reads are 114 100 948 but only 8 075 928 reads were reported mapped. How could I find the related information for the other 106 025 020 reads or what happen to these reads?

Thanks for the time and help.

ADD COMMENT

Login before adding your answer.

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