Does one need to sort and index after performing a rmdup in Samtools?
Hi,
I would like to know whether one needs to sort and to index the .bam files again after removing the duplicates in Samtools using samtools rmdup -S.
I am under the impression that since we have already sorted and indexed the .bam files, and we are just removing the duplicates, there would be no need to do the sorting and indexing again. Am I wrong?
Thanks in advance for your thoughts!
1 answer
"sort"
no: your bam is already sorted.
"index"
yes if you need to access some specific genomic regions in your bam. The previous index of the bam cannot be used with the new rmdup.bam
• 0 views
•
link
updated
by
Ram
•
written
by
Pierre Lindenbaum
•
written
by
Pierre Lindenbaum
Log in to answer this question.
More posts like this
-
Why did the QV value decrease significantly after using YAHS for HiC scaffolding?
written by xinguok794 •
-
samtools index error - "Unsorted positions on sequence" even though file is sorted
written by Geoffrey •
-
How can I insert a "remove duplicates" command in a samtools script?
written by José •
-
Input and output of picard markduplicate for deepvariant
written by hello_bioinf •
-
samtools mpileup fail to create bcf
written by mthm •
-
samtools index error
written by evelyn
-
should FASTA files be sorted before indexed with SAMtools?
written by Makplus T
-
Does bamtools split preserve sort and rmdup?
written by ccnn •
-
samtools sorting and indexing
written by ggman
-
Remove Duplicates Workflow
written by Anjan
