def agreement(xs, ys, scoring, alignment): score, x, y = compute_global_alignment(xs, ys, scoring, alignment) similar = [1. for (a, b) in zip(x, ys) if a == b] return 100. * len(similar) / len(x)
def edit_dist(xs, ys): alphabet = ascii_lowercase scoring = build_scoring_matrix(alphabet, 2, 1, 0) align = compute_alignment_matrix(xs, ys, scoring, True) score, x, y = compute_global_alignment(xs, ys, scoring, align) return len(xs) + len(ys) - score