Remove Duplicates from Illumina AmpliSeq data
Hi Everyone,
Could anyone please recommend any tools to remove duplicates from AmpliSeq paired-end data (not MarkDuplicates)? I am having a difficult time to find such tool. Any help is much appreciated. Please let me know if further details are required.
Thanks.
Best wishes,
Preeti
• 2,837 views
•
link
written
by
preeti.singh •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Tools for strand direction detection RNA-Seq
written by pubsurfted •
-
featureCounts parameters
written by Lucy •
-
Remove duplicate reads for bisulfite sequencing reads
written by qwzhang0601
-
How can I filter RNA-seq reads to remove reads where more than 50% of the read is low quality (Q<20…
written by mar77 •
-
CIRCexplorer2 for paired end data
written by archana.bioinfo87
-
Read Name in IGV and duplicates
written by ww22runner •
-
About clustering genes isoforms in RNA-Seq
written by Lucas Peres
-
Breakdancer + Squaredancer
written by Hiroyuki Mishima •
-
Getting Confused With The Flagstat After Pcr Duplicates Removed
written by KJ Lim •
-
About Paired-End Sequencing
written by Pierre Lindenbaum
why ?
Hi Pierre,
In theory, MarkDuplicates and/or rmdup will remove all reads. In Ampliseq, DNA is fragmented after the target is captured so most of the reads could possibly be from same DNA fragment, also ion Torent suggest not to remove duplicates. So was wondering if any tool is specific for AmpliSeq reads.
Preeti
Are it just fasta files?
No on aligned bam files (after conversion of aligned SAM files)
AmpliSeq is an amplicon based (pcr based) library method, isn't it? Unless you have used UMI, you have no chance to say which read is a pcr duplicate. So don't remove "duplicates", because as you said, you will end up having only 1 read per position.
Thanks finswimmer.
As you suggested I will give it a go without removing duplicates and see how the VCFs look on IGV.
Preeti