Entering edit mode
4.2 years ago
pingu77
▴
20
Dear All,
I am trying to find the PCR duplicate reads from the bam/sam. I have used picard to mark duplicates in the bam/sam file and then samtools flagstats. I obtained the following output:
196499952 + 0 in total (QC-passed reads + QC-failed reads)
125920476 + 0 secondary
0 + 0 supplementary
11613716 + 0 duplicates
195906701 + 0 mapped (99.70% : N/A)
70579476 + 0 paired in sequencing
35289738 + 0 read1
35289738 + 0 read2
68711404 + 0 properly paired (97.35% : N/A)
69739378 + 0 with itself and mate mapped
246847 + 0 singletons (0.35% : N/A)
127264 + 0 with mate mapped to a different chr
31498 + 0 with mate mapped to a different chr (mapQ>=5)
Is the number of duplicate reads calculated with respect to the number of the QC-passed reads?
Thank you in advance