fork download
  1. from collections import Counter
  2.  
  3.  
  4. def contacodon(seq):
  5. stops = ["TAA", "TAG", "TGA"]
  6. count = Counter()
  7. for i in range(seq.find('ATG'), len(seq) + 1, 3):
  8. codon = seq[i: i + 3]
  9. count.update([codon])
  10. if codon in stops:
  11. break
  12. return dict(count)
  13.  
  14.  
  15. seq = "GCGATCGAGATGAGCATCGCATCGCGGACTACCGCGCGCGCGCGCGGGAGATGAGCATCGACGACTCGACTAG"
  16. print(contacodon(seq))
  17.  
Success #stdin #stdout 0.04s 9092KB
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}