# Lista de stop codons:
stops = ["TAA", "TAG", "TGA"]
count = {}
for i in range(s.find('ATG'), len(s) + 1, 3):
codon = s[i: i + 3]
if codon in count:
count[codon] += 1
else:
count[codon] = 1
# Se tiver encontrado um stop codon, pare de iterar:
if codon in stops:
break
print(count)
OmitiOmiti, por brevidade, o restante do código. Veja funcionando no Ideone.
Vale lembrar que, embora funcione com um dicionário comum, o Python oferece outras abordagens melhores para realizar a contagem, tal como a classe Counter
, disponível no módulo collections
. Veja:
from collections import Counter
stops = ["TAA", "TAG", "TGA"]
count = Counter()
for i in range(seq.find('ATG'), len(seq) + 1, 3):
codon = seq[i: i + 3]
count.update([codon])
if codon in stops:
break
print(dict(count))
Veja funcionando no Ideone.
Uma outra alternativa é utilizar o defaultdict(int)
. Embora menos poderoso que o Counter
, já torna implícita a inicialização de valores 0
no dicionário, o que é feito manualmente no primeiro exemplo de código. Veja mais sobre eles no SOen.
Vale observar, por brevidadefim, que os dois códigos apresentados nesta resposta imprimirão o resultado da contagem mesmo que um stop codon não seja encontrado. Ademais, se o restanteúltimo "códon" dessa string não terminada tiver um ou dois caracteres (e não três, conforme esperado), este também será contabilizado. Para prevenir isso, você pode implementar algum tipo de verificação. Um exemplo é utilizar uma expressão regular, tal como excelentemente demonstrado pela resposta do Augusto Vasques, de modo a invalidar sequências que fogem do códigopadrão esperado.