MarkDuplicates Error: Exception in thread "main" htsjdk.samtools.SAMException: Value was put into...

问题:

    WES 分析中,在标记重复reads这一步出现如下错误:


Exception in thread "main" htsjdk.samtools.SAMException: Value was put into PairInfoMap more than once


解决:

    比对完毕的bam文件中存在 unpaired reads,需要对这一步进行过滤后重新标记,bwa比对完毕以后通过samtools进行过滤,代码如下:

代码:samtools view -f 0x2 -b in.bam > out.bam

随后重新进行排序和MarkDuplicates,无报错。

参考:

    https://bioinformatics.stackexchange.com/questions/19271/how-to-remove-the-unpaired-reads-in-sam-bam-files

    https://gatk.broadinstitute.org/hc/en-us/community/posts/4408717387803-SAMException-Value-was-put-into-PairInfoMap-more-than-once    

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容