Bug Bulletin: The recent 3.2 release fixes many issues. If you run into a problem, please try the latest version before posting a bug report, as your problem may already have been solved.

FlagStat

SystemSystem Posts: 226Administrator admin
edited July 2012 in Tool Bulletin

A new tool has been released!

Check out the documentation at FlagStat.

Comments

  • mgujralmgujral Posts: 2Member

    Hi there,

    I have a question related to the output of FlagStat. I used BWA mem with version 0.7.4 to align the reads and following is the output of FlagStat

    1290814928 in total 0 QC failure 0 duplicates 1269372151 mapped (98.34%) 1290814928 paired in sequencing 645318681 read1 645496247 read2 0 properly paired (0.00%) 1264066894 with itself and mate mapped 5305257 singletons (0.41%) 54566084 with mate mapped to a different chr 13502927 with mate mapped to a different chr (mapQ>=5)

    In output,1290814928 paired in sequencing, but 0% are properly paired. What does that mean? What could be reason for 0% properly mapped? I will greatly appreciate your input.

    best, --Madhu

  • ebanksebanks Posts: 678GATK Developer mod

    Hi Madhu,

    The proper pair status is read directly from the SAM record - it's the 0x2 bit in the FLAG field (the second token in a SAM record). See the SAM specification for more details.

    Eric Banks, PhD -- Senior Group Leader, MPG Analysis, Broad Institute of Harvard and MIT

Sign In or Register to comment.