去除Bam文件中的PCR 重复

1.使用samtools 去除重复

samtools sort  -n -@ 20  file.bam |samtools fixmate -m -@20 - -  |samtools sort -@ 20 |samtools markdup -r -@ 20 - - >file.rmdup.bam

2.使用Picard去除PCR重复

picard MarkDuplicates REMOVE_DUPLICATES=True I=file.bam  O=file.rmdup.bam M=file_marked_dup_metrics.txt

 

posted @ 2023-02-08 23:21  bio调参侠  阅读(342)  评论(0)  编辑  收藏  举报