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)