Does one need to sort and index after performing a rmdup in Samtools?
1
2
Entering edit mode
7.6 years ago
Dataman ▴ 350

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!

sequencing alignment samtools variant • 4.2k views
ADD COMMENT
9
Entering edit mode
7.6 years ago

"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

ADD COMMENT

Login before adding your answer.

Traffic: 611 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6