Deleted:PCR duplicate reads - samtools
1
0
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

samtools duplicated flagstats picard • 1.5k views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 2996 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