2016年11月17日
mismatch位置（MD tag）- sam/bam格式解读进阶
摘要： 这算是第二讲了，前面一讲是：Edit Distance编辑距离（NM tag）- sam/bam格式解读进阶 MD是mismatch位置的字符串的表示形式，貌似在call SNP和indel的时候会用到。 当然我这里要说的只是利用它来计算mismatch的个数 MD = line.get_tag('MD') pat = "[0-9]+[ATGC]+" MD_list = re.find... 阅读全文
posted @ 2016-11-17 17:59 Life·Intelligence 阅读(2962) 评论(0) 推荐(0) 编辑
TOP