def compute_perc_id(aln):
	""" Compute percent identity of aligned region on read """
	length = len(aln.query_alignment_sequence)
	edit = dict(aln.tags)['NM']
	return 100 * (length - edit)/float(length)