SEQanswers

Go Back   SEQanswers > General



Similar Threads
Thread Thread Starter Forum Replies Last Post
Dindel is seeing too many reads, potential bug? Heisman Bioinformatics 2 01-10-2013 04:11 AM
How to differentiate between PCR duplicates and real data? frymor Bioinformatics 5 09-15-2011 12:05 PM
PCR duplicates questions slny Bioinformatics 8 06-07-2011 03:06 AM
PCR duplicates increase when excess of beads tdm SOLiD 10 03-31-2011 08:48 AM
Picard - MakeDuplicates (remove pcr duplicates) dmb Bioinformatics 2 03-16-2011 07:56 AM

Reply
 
Thread Tools
Old 01-18-2010, 10:53 AM   #1
julien
Junior Member
 
Location: Raleigh NC

Join Date: Nov 2009
Posts: 3
Default how critical is the filtering of potential PCR duplicates?

Hello,

I don't really have a good appreciation for how critical it is to remove potential PCR duplicates from the alignments (either via samtools rmdup or picard). It would appear (based on posts in this forum and pipelines that I have found online) that many ??? do not worry about this issue, am I wrong ? It did not seem to affect my data very much but I would like to know how important an issue this is.

Thank you !!!!
julien is offline   Reply With Quote
Old 01-18-2010, 12:56 PM   #2
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

Quote:
Originally Posted by julien View Post
Hello,

I don't really have a good appreciation for how critical it is to remove potential PCR duplicates from the alignments (either via samtools rmdup or picard). It would appear (based on posts in this forum and pipelines that I have found online) that many ??? do not worry about this issue, am I wrong ? It did not seem to affect my data very much but I would like to know how important an issue this is.

Thank you !!!!
Most whole-genome human resequencing papers try to detect PCR duplicates and I find that it is extremely important to remove them (or flag them), especially when dealing with low-complexity libraries. The latter can happen with small amounts of input DNA or re-using the same library. When searching for variants, it is assumed that each read is conditionally independent given the underlying DNA sequence (i.e. not the same DNA fragment read twice).
nilshomer is offline   Reply With Quote
Old 03-26-2010, 05:27 AM   #3
dnusol
Senior Member
 
Location: Spain

Join Date: Jul 2009
Posts: 133
Default

Hi nilshomer, sorry for getting back to this thread, but do you know if Picard is able to remove duplicates from single-reads? Samtools manual does not recommend it but I am not sure if it is a matter of implementation or a real technical issue on the single read part which I don't get.

best
dnusol is offline   Reply With Quote
Old 03-26-2010, 09:24 AM   #4
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

Quote:
Originally Posted by dnusol View Post
Hi nilshomer, sorry for getting back to this thread, but do you know if Picard is able to remove duplicates from single-reads? Samtools manual does not recommend it but I am not sure if it is a matter of implementation or a real technical issue on the single read part which I don't get.

best
If you have high coverage (say 50x for 50bp reads), then it becomes very likely that multiple single end reads will start at the same coordinate. You will then wont be able to tell coverage from PCR duplicates. For PE or LMP reads, both ends can be used to determine duplication, since the insert size has its own distribution.
nilshomer is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 09:30 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO