how samtools rmdup remove cut reads made more reads been drop( cut or uncut reads mapping ratio are similar)?
compare PE cut reads with uncut reads samtools rmdup bam,how samtools rmdup remove cut reads made more reads been drop ?( cut or uncut reads mapping ratio are similar?)
• 1,526 views
•
link
written
by
siqizhao2018 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
samtools markdup in Rsamtools
written by Learner •
-
how to identify where multi-mappers, or "non primary hits" are mapping to
written by bertb •
-
samtools markdup vs samtools rmdup
written by a.abnousi •
-
Does bamtools split preserve sort and rmdup?
written by ccnn •
-
Removal of PCR dupliactes from trimmed reads
written by deepti1rao •
-
Picard vs samtools rmdup
written by xd_d
-
Remove some reads from a BAM file
written by lisadavic66 •
-
Samtools: removing PCR duplicates
written by devinliao0918 •
-
Samtools Rmdup And Picard Mark Duplicates
written by Kssr
-
How Reads Have Been Mapped In My Bam File?
written by Pascal
eg:I had a PE reads,mapping then samtools rmdup,the result remove 20% reads,but when I cut the reads as 36bp(sequencing 150bp PE,36bp no any adapter),then mapping ,the mapping ratio more then 85%,but when I rmdup,it reduce 75%,only 25% reads are keep,I don't know what happend...
Hello,
the estimation of duplicates is based on mapping information. If you cut your reads to that small size of 36bp it becomes more likely that two reads map to the same position. So it is not that surprising that your duplicate number increases.
What program with which parameters are using für trimming? Why do you think this is necessary?
fin swimmer
By the way,
rmdupis deprecated. For single-end reads I still use it because it is fast and can read from STDIN, but for paired-end data, you may considermarkdup(introduced in SAMtools 1.6),SAMBLASTERorMarkDuplicatesfrom Picard.