MD-editDistance-fromBam_2.py

import pysam
import re

def extractMDFromBam(bam):
	'''
	loop over a bam file and get the edit distance to the reference genome
 	mismatches are stored in the MD tag
	scale by aligned read length
	'''
	samfile = pysam.AlignmentFile(bam, "rb")
	return [sum([len(item) for item in re.split('[0-9^]', read.get_tag("MD"))]) / read.query_alignment_length for read in samfile.fetch()]