samtools`markdup`操作的正确顺序

具体实例在5 一步法找基因变异流程

samtools**markdup**操作的正确顺序

1. The first sort can be omitted if the file is already name ordered

代码语言：javascript复制

samtools sort -n -o namesort.bam example.bam

必须是-n，不能省略

1. Add ms and MC tags for markdup to use later

代码语言：javascript复制

samtools fixmate -m namesort.bam fixmate.bam

-m不能少

1. Markdup needs position order

代码语言：javascript复制

samtools sort -o positionsort.bam fixmate.bam

1. Finally mark duplicates

代码语言：javascript复制

samtools markdup positionsort.bam markdup.bam

add duplicates file position tags

1 人点赞