python
from Bio import SeqIO
record = SeqIO.read("genome.fasta", "fasta")
python
sequence_length = len(record.seq)
reverse_seq = record.seq.reverse_complement()
search_seq = "ATG"
if search_seq in record.seq:
else:
python
record = SeqIO.read("genome.gb", "genbank")
gene_name = record.name
for feature in record.features:
if feature.type == "CDS":
cds_seq = feature.qualifiers["translation"][0]
protein_seq = Seq(cds_seq)
python
from Bio import Align
aligner = Align.PairwiseAligner()
alignments = aligner.align("ACGT", "ACT")
for alignment in alignments:
print(alignment)
