We are currently using V4 Illumina Sequencers. I am trying to estimate optical duplicates and to do so I have recently upgraded PicardTools to V 1.122 to enable readname abilty. The generated metrics file now outputs a non-zero value for READ_PAIR_OPTICAL_DUPLICATES which seems reasonable.
When I change the option OPTICAL_DUPLICATE_PIXEL_DISTANCE in the PicardTool MarkDuplicate I get differing duplicate numbers of course. For example if the option is left at default (100pixels) the tool returns 1111 optical duplicates for a particular BAM. When I change the option to 10pixels the same BAM is shown to have 60 optical duplicates.
My question, what is the ideal pixel distance to set OPTICAL_DUPLICATE_PIXEL_DISTANCE to for this V4 generation of Illumina sequencer? Is the default good enough?
•
written
by
sbilobram •
0 answers
No answers yet.
Log in to answer this question.
-
GATK HaplotypeCaller missed variants in amplicon gene panel data
written by Sofya •
-
samtools depth and sambamba depth produce different results
written by Johannes •
-
Tutorial: Yes .. BBMap can do that! - Part III clumpify (mark (and dedupe) duplicates without align…
written by GenoMax
-
STAR aligner maps 0 reads to multiple loci, but 3.7% to too many loci - meaning?
written by Anand Rao
-
AddOrReplace gp in Pickard Sam validation error
written by KG •
-
GATK - gene panel VCF - too many variants - how to filter out
written by stanedav •
-
Tutorial: Tutorial (How to analyze) on Whole Exome sequencing. Common Errors. Best Practices.
written by Chirag Nepal
-
Should quality scores be considered when discriminating between optical and library duplicates?
written by Dan D
-
Tutorial: Reference Assembly - Mapping Reads To A Reference Genome
written by Joseph Hughes
-
Library Duplicates Vs. Optical Duplicates (Picard Markduplicates)
written by toni
