To count nucleotides in a Fasta file:
from Bio import SeqIO
for seq in SeqIO.parse('chrA01.fasta', 'fasta'):
s = str(seq.seq)
print s.count('T') + s.count('G') + s.count('C') + s.count('A')
-------------------------------------------------------
from Bio import SeqIO
for seq in SeqIO.parse('chrA01.fasta', 'fasta'):
print s.count('T') + s.count('G') + s.count('C') + s.count('A') + s.count('N')
-------------------------------------------------------
from Bio import SeqIO
for seq in SeqIO.parse('chrA01.fasta', 'fasta'):
print len(seq.seq)
No comments:
Post a Comment