python from Bio.Seq import Seq from Bio import pairwise2 seq = Seq("ATGGCCATTGTAATGGGCCGCTGAAAGGGTGCCCGATAG") protein_seq = seq.translate() print("Translated protein sequence:", protein_seq) seq1 = Seq("ATCGTACGATC") seq2 = Seq("ATCGTAGCTAC") alignments = pairwise2.align.globalxx(seq1, seq2) for alignment in alignments: print("Alignment score:", alignment[2]) print("Sequence 1 alignment:", alignment[0]) print("Sequence 2 alignment:", alignment[1])


上一篇:
下一篇:
切换中文