fork download
  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.  
Success #stdin #stdout 0.02s 9052KB
stdin
Standard input is empty
stdout
{'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}