python
from Bio.Seq import Seq
dna_seq = Seq("ATCGTTAGCTGACTG")
print(dna_seq)
python
from Bio import SeqIO
for seq_record in SeqIO.parse("sequence.fasta", "fasta"):
print("ID:", seq_record.id)
print("Seq:", seq_record.seq)
python
from Bio import SeqIO
seq = Seq("ATCGTTAGCTGACTG")
seq_record = SeqIO.SeqRecord(seq, id="seq1", description="Example sequence")
SeqIO.write(seq_record, "output.fasta", "fasta")
python
from Bio import pairwise2
seq1 = "ATCGTTAGCTGACTG"
seq2 = "ATCGCTAGCTAGCTG"
alignments = pairwise2.align.globalxx(seq1, seq2)
for alignment in alignments:
print("Score:", alignment.score)
print("Seq1:", alignment.seqA)
print("Seq2:", alignment.seqB)
python
from Bio import Phylo
alignment = [
Seq("ATCGTTAGCTGACTG"),
Seq("ATCGCTAGCAGACTG"),
Seq("ATCGCTAGCTATCTG"),
Seq("ATCGCTAGCTGACTA")
]
tree = Phylo.Phyml().build_tree(alignment)
Phylo.draw(tree)
