removing PCR duplicates from MedipSeq Data
2
0
Entering edit mode
5.6 years ago

Dear all

i have paired end data of medip seq from illumina platform. i have done mapping using bowtie2 and after mapping is it necessary to remove PCR Duplicates from bam file. if it is not removed whether it will effect in DMR or peak calling analysis.

Need help

Thank you

removePCRduplicates Medipseq • 1.8k views
ADD COMMENT
1
Entering edit mode
5.6 years ago
Medhat 9.4k

PCR duplicates are sequence reads that result from sequencing two or more copies of the exact same DNA fragment, which, at worst, may contain erroneous mutations introduced during PCR amplification, or, at the very least, make the occurrence of the allele(s) sequenced in duplicates appear proportionately more often than it should compared to the other allele

It is always recommended to remove PCR duplication and you can use this tool

https://broadinstitute.github.io/picard/command-line-overview.html#MarkDuplicates

you can read more about the topic here:

Evaluating the necessity of PCR duplicate removal from next-generation sequencing data and a comparison of approaches

ADD COMMENT
1
Entering edit mode
5.6 years ago

Whether or not you should remove PCR duplicates depends on many factors.

1) Did you PCR-amplify? If not, then there are no PCR duplicates.

2) Are your reads paired? If not, most of your duplicates will be false-positives (coincidental rather than PCR duplicates).

3) Is your coverage locally super-high (RNA-seq, Chip-seq, etc)? If so, most of your duplicates will be false-positives.

ADD COMMENT

Login before adding your answer.

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