Skip to main content
adicionou 1739 caracteres ao conteúdo
Fonte Link
Luiz Felipe
  • 36,5mil
  • 10
  • 74
  • 154

Vale observar, por fim, 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 ú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 padrão esperado. Deixo um outro exemplo no Ideone.

Vale observar, por fim, 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 ú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 padrão esperado.

Vale observar, por fim, 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 ú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 padrão esperado. Deixo um outro exemplo no Ideone.

adicionou 1739 caracteres ao conteúdo
Fonte Link
Luiz Felipe
  • 36,5mil
  • 10
  • 74
  • 154
# 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.

# 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

Omiti, por brevidade, o restante do código.

# 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)

Omiti, 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 fim, 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 ú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 padrão esperado.

Fonte Link
Luiz Felipe
  • 36,5mil
  • 10
  • 74
  • 154

O código funciona, de modo que irá realizar a contagem a partir do códon ATG até o TAA.

O problema é que existem três stop codons diferentes, e o programa só reconhece o TAA. Veja, o problema está onde você cria o intervalo para as iterações:

range(s.find('ATG'), s.find('TAA') + 1, 3)

Você está basicamente criando um intervalo, que incrementa de três em três, do índice da primeira ocorrência da substring "ATG" até o índice da primeira ocorrência da substring "TAA".

Vejo dois problemas ao especificar um índice final tal como resultado de s.find('TAA'):

  1. Você corre o risco de criar um range que percorre sobre uma sequência com número inválido de códons (pode haver apenas 5 caracteres entre ATG e TAA, por exemplo).
  2. Você não reconhece os outros dois códons de finalização, "TAG" e "TGA".

Para corrigir esses dois problemas, simplesmente não defina um limite superior no range. Simplesmente itere até o comprimento total da string.

Nesse sentido, você pode fazer a verificação do stop codon dentro do próprio laço de repetição, de modo a utilizar a declaração break para interromper as repetições se tiver encontrado um códon de finalização.

Mais ou menos assim:

# 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

Omiti, por brevidade, o restante do código.