python
from Bio import SeqIO
sequence = SeqIO.read("sequence.fasta", "fasta")
base_counts = dict(sequence.seq.count(x) for x in "ACGT")
for base, count in base_counts.items():
print(f"{base}: {count}")
python
from Bio import pairwise2
from Bio import Align
alignments = pairwise2.align.globalxx("AGTACACTGG", "ACTGACTG")
for align in alignments:
print(align)
tree = Align.TreeConstructor()
evolutionary_tree = tree.build_tree(alignments)
print(evolutionary_tree)
