pip install bccb
import bccb
python
from bccb import seq
dna_seq = seq.read_dna('dna_sequence.fasta')
print(dna_seq)
python
from bccb import align
seq1 = 'AACCTTGG'
seq2 = 'ACCTTG'
alignment = align.SmithWaterman(seq1, seq2)
print(alignment)
python
from bccb import seq
dna_seq = 'ACGTACGTACGT'
base_count = seq.count_bases(dna_seq)
print(base_count)
