python pip install biopython python from Bio import SeqIO from Bio.Seq import Seq from Bio.Alphabet import IUPAC python file_path = "sequence.fasta" sequences = SeqIO.parse(file_path, "fasta") for sequence in sequences: print("ID:", sequence.id) print("Sequence:", sequence.seq) python protein_seq = Seq("MKWVTFISLLFLFSSAYSRGVFRRDAHKSEVAHRFKDLGEENFKALVLIAFAQYLQQCPFEDHVKLVNEVTEFAKTCVADESAENCDKSLHTLFGDKLCTVATLRETYGEMADCCAKQEPERNECFLSHKDDSPDLPKLKPDPKTNKQCQVTCAGTMDGKNHGLVLGKFDSLDTVKNNDKDWQDQLLentingSLGQYLALQkDPDPNTNCQISLENLDRNDFQDEN", IUPAC.protein) python length = len(sequence) python base_count = sequence.count("A") aa_count = protein_seq.count("L") python reverse_comp = sequence.reverse_complement() python protein_seq = sequence.translate() python pattern = "ACGT" matches = sequence.find(pattern)


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