pip install biopython
python
from Bio import SeqIO
from Bio.SeqUtils import GC
python
for record in SeqIO.parse("genome.fasta", "fasta"):
sequence = record.seq
python
gc_content = GC(sequence)
python
pattern = "ATGC"
if pattern in sequence:
print("Pattern found!")
else:
print("Pattern not found.")
python
from Bio.Seq import Seq
from Bio import pairwise2
sequence1 = Seq("ATCG")
sequence2 = Seq("AGTC")
alignments = pairwise2.align.globalxx(sequence1, sequence2)
for alignment in alignments:
print(alignment)
python genomics_analysis.py
