Estou tentando resolver um exercício no qual foi solicitado uma função com um start e um stop, onde o start deve ocorrer quando encontrar "ATG"
e o stop quando encontrar o "TAA"
, "TAG"
, "TGA"
.
Com ajuda dos comentários, consegui ajeitar um pouco o código, porém o stop está limitado a somente ao que eu colocar no s.find()
. Gostaria de ajuda para fazer ele procurar por um dos três stop codons proposto no exercício.
Segue o enunciado do exercício que gerou a pergunta:
Faça uma função
contacodon
em Python que receba uma sequência de letras representando nucleotídeos (e.g.,A
,C
,G
,T
), verifique se a sequência é válida (i.e. contém somenteA
s,C
s,G
s, eT
s) e retorne com o número de ocorrências de cada codon contando a partir do start codon"ATG"
até um dos stop codons,"TAA"
,"TAG"
,"TGA"
Exemplo de entrada:
contacodon("AGCGATCGAGATGAGCATCGCATCGCGGACTACCGCGCGCGCGCGCGGGAGATGAGCATCGACGACTCGACTAG")
Saída para a entrada acima:
{
'ATG': 1,
'AGC': 1,
'ATC': 1,
'GCA': 1,
'TCG': 1,
'CGG': 1,
'ACT': 1,
'ACC': 1,
'GCG': 2,
'CGC': 2,
'GGG': 1,
'AGA': 1,
'TGA': 1
}
Meu código:
def verificar (s):
s = s.upper()
for ent in s:
if not ent in "ACGT":
return False
return True
while True:
s = str(input("Entre com a seq: \n")).upper()
if verificar(s):
break
print("Seq inválida")
count={}
for i in range(s.find('ATG'),s.rfind('TAA')+1,3):
codon = s[i:i+3]
if codon in count:
count[codon] += 1
else:
count[codon] = 1
print('\n', count, '\n',)
range
) é completamente inválida. Você provavelmente terá que pensar em outra coisa (e, nesse caso, não faremos isso para você aqui, já que não é o objetivo do site).str.find()
para localizar o indice da primeira ocorrência deATG
.range
começar em zero, deve começar do índice do primeiro codon ATG (que dá pra encontrar usando o métodos.find()
do seu string de entrada, como o @AugustoVasques comentou). Outra coisa, norange
não há necessidade de calcular o final exato da sequência comlen(s)-len(s)%3
, já que você pega 3 caracteres com slices, e slices nunca dãoIndexError
.s[s.find("ATG"):]
e a fragmentar em porções de três caracteres com essa função pt.stackoverflow.com/a/496160/137387