pip install biopython
python
from Bio.Seq import Seq
from Bio import SeqIO
from Bio.Align import MultipleSeqAlignment
from Bio.SeqUtils.ProtParam import ProteinAnalysis
python
record = SeqIO.read("protein.fasta", "fasta")
protein_sequence = record.seq
python
protein_analysis = ProteinAnalysis(str(protein_sequence))
amino_acid_composition = protein_analysis.get_amino_acids_percent()
for amino_acid, percentage in amino_acid_composition.items():
print(amino_acid, ":", percentage)
molecular_weight = protein_analysis.molecular_weight()
isoelectric_point = protein_analysis.isoelectric_point()
gravy_score = protein_analysis.gravy()
python
seq1 = Seq("ATTGCATAG")
seq2 = Seq("ATAGCGTAA")
alignment = MultipleSeqAlignment([seq1, seq2])
print(alignment)
