Решала задачу с сайта rosalind.info. Дано: FASTA-файл. Нужно: получить процент, который составляют G
и C
от всех символов в каждой из последовательностей.
Пришла к следующему решению:
from Bio import SeqIO
filename = "rosalind_gc.txt"
count = 0
for record in SeqIO.parse(filename, "fasta"):
seqs = record
a = seqs[count].count('A')
t = seqs[count].count('T')
g = seqs[count].count('G')
c = seqs[count].count('C')
Q = (g+c)/(a+g+t+c)
print("Record " + record.id + ", GC content " + str(Q))
count = count + 1
Проблема в том, что оно выдаёт вместо необходимого процентного содержания G
и C
нули.
В чём может быть ошибка?