SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   PCR-Free Libraries & Bioinformatic PCR Duplicate detection? (http://seqanswers.com/forums/showthread.php?t=34690)

aeonsim 10-22-2013 08:52 AM

PCR-Free Libraries & Bioinformatic PCR Duplicate detection?
 
Does any one have any thoughts or experience with regards to the Bioinformatics pipeline for Illumina PCR-Free libraries?

Especially with regards to the need to run Picards MarkDuplicates or some other PCR duplicate detection tool?

As we've started a large sequencing project using the PCR-Free libraries and having run Markduplicates over the initial test lane of sequence it's reporting a PCR-duplicate rate of between 0.004-0.006. Some of which it claims are Optical Duplicates and the rest (2.5x the Optical numbers) it claims are Library duplicates. However considering the kit & protocol was one of the PCR-Free ones I'm wondering if these are actual duplicates or misidentification by Markduplicates (potentially reads from repetitive regions).

Any thoughts would be greatly appreciated.

atcghelix 10-22-2013 02:44 PM

What are you sequencing? And are you sequencing paired-ends?

As I understand it tools find "PCR duplicates" by comparing mapping coordinates of reads. However, if in your library prep you have actual fragments that shear in the same place then the resulting reads will also map to the same location. So they'll look like PCR duplicates in that they map to the same spot on the reference, but they actually represent good, unique reads.

The probability of this happening will increase with smaller references (smaller genomes, RNA-seq, targeted enrichment sequencing). The probability should decrease greatly by using paired-end reads as well, as the fragments in the library then have to shear at the same location at each end of the fragment.

I don't know if there's a standard protocol, but I'd be inclined to filter out the optical duplicates and leave in the "PCR" duplicates in your case.

aeonsim 10-22-2013 03:33 PM

Data is for WGS Bos Taurus, ~3GB Genome, using 101bp PE or 109bp PE PCR-Free Libraries.

atcghelix 10-23-2013 01:36 PM

If you pull out the reads that have been marked as duplicates, do they look like they are highly repetitive?

And just to make sure, the reads that are showing up as duplicates are reads that map to the Bos taurus genome, yes?


All times are GMT -8. The time now is 06:48 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.