• Source
    1. def contacodon(seq):
    2. # Lista de stop codons:
    3. stops = ["TAA", "TAG", "TGA"]
    4.  
    5. count = {}
    6. for i in range(seq.find('ATG'), len(seq) + 1, 3):
    7. codon = seq[i: i + 3]
    8. if codon in count:
    9. count[codon] += 1
    10. else:
    11. count[codon] = 1
    12.  
    13. # Se tiver encontrado um stop codon, pare de iterar:
    14. if codon in stops:
    15. break
    16. return count
    17.  
    18.  
    19. seq = "GCGATCGAGATGAGCATCGCATCGCGGACTACCGCGCGCGCGCGCGGGAGATGAGCATCGACGACTCGACTAG"
    20. print(contacodon(seq))
    21.