from Bio import Align aligner = Align.PairwiseAligner() seq1 = "GAACT" seq2 = "GAT" score = aligner.score(seq1, seq2) score 3.0 alignments = aligner.align(seq1, seq2) for alignment in alignments: print(alignment)