samtools-rmdup - removes duplicate reads (obsolete)
samtools rmdup [
-sS]
input.srt.bam out.bam
This command is obsolete. Use markdup instead.
Remove potential PCR duplicates: if multiple read pairs have identical external
coordinates, only retain the pair with highest mapping quality. In the
paired-end mode, this command
ONLY works with FR orientation and
requires ISIZE is correctly set. It does not work for unpaired reads (e.g. two
ends mapped to different chromosomes or orphan reads).
- -s
- Remove duplicates for single-end reads. By default, the
command works for paired-end reads only.
- -S
- Treat paired-end reads and single-end reads.
- o
- Samtools paired-end rmdup does not work for unpaired reads
(e.g. orphan reads or ends mapped to different chromosomes). If this is a
concern, please use Picard's MarkDuplicates which correctly handles these
cases, although a little slower.
Written by Heng Li from the Sanger Institute.
samtools(1)
Samtools website: <
http://www.htslib.org/>