Se as trincas estão separadas por espaço, uma alternativa é separá-las usando split
e contar quantas são iguais a "ATC":
seq = 'ATC CAA GTC AGC TAG CGT ATC ATC GTC ATG CTC AAA CAC TAC GAT GCT AAT'
atc = 0 # quantidade de ATC
for s in seq.split():
if s == 'ATC':
atc += 1
Já para contar as letras individualmente, você pode usar um Counter
:
from collections import Counter
c = Counter(seq)
E para obter os totais basta pegá-los do Counter
:
print(f'A quantidade de nucleotideos "A" é {c["A"]}, de "T" é {c["T"]}, de "C" é {c["C"]}, de "G" é {c["G"]}, de "ATC" é {atc}')
Quanto ao total geral, bastaria somar os valores do Counter
, ignorando os espaços:
total = sum(qtd for s, qtd in c.items() if s in ('A', 'T', 'C', 'G'))
Detalhe: quando você itera por uma string (como em for x in seq
), a cada iteração o x
será um dos caracteres da string, e portanto ele nunca será igual a "ATC". Sendo assim, seu código não funcionará.
Agora se não tiver espaços, uma alternativa é iterar pela string e ir pegando pedaços de 3 em 3:
seq = 'ATCGATCTA'
atc = 0
for i in range(0, len(seq), 3):
if seq[i:i + 3] == 'ATC':
atc += 1
Eu uso a sintaxe de slicing para pegar determinado trecho da string: seq[i:i + 3]
pega 3 caracteres a partir da posição i
, e o i
vai de zero até o tamanho da string, pulando de 3 em 3. Ou seja, no exemplo acima eu pego primeiro "ATC", depois "GAT" e depois "CTA", e portanto a quantidade de "ATC" será 1.
Nos comentários você disse que usou count
, mas cuidado que pode dar diferença. Para o exemplo acima, seq.count('ATC')
resulta em 2, pois "ATC" ocorre 2 vezes na string:
ATCGATCTA
^^^ ^^^
| |
| \_ segunda ocorrência
\_____ primeira ocorrência
Só que a segunda ocorrência eu entendo que está errada, pois o "AT" faz parte da trinca "GAT" e o "C" faz parte do "CTA", então na verdade não é uma ocorrência de "ATC".