双序列对比：动态规划

from Bio import SeqIO
import swalign
xPA_seq = SeqIO.read("...","genbank")
xHA_seq = SeqIO.read("...","genbank")
xPA = xPA_seq.seq
xHA = xHA_seq.seq
cNP_seq = SeqIO.read("...","genbank")
cHA_seq = SeqIO.read("...","genbank")
cNP = cNP_seq.seq
cHA = cHA_seq.seq
match_score = 2
mismatch_score = -1
matrix = swalign.NucleotideScoringMatrix(match_score, mismatch_score)
lalignment_object = swalign.LocalAlignment(matrix)
alignment_object = lalignment_object.align(xPA, cNP)
alignment_object.dump()