Guest User

Untitled

a guest
Jun 18th, 2018
77
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.39 KB | None | 0 0
  1. with open ("dna.txt", "r") as myfile:
  2. data=myfile.readlines()
  3. mRNA=data[1]
  4.  
  5. def next_transcript(mRNA, cur_pos):
  6. initial=mRNA.find("AUG", cur_pos)
  7. for i in range(initial, len(mRNA),3):
  8. if mRNA[i:i+3] == "UAG":
  9. return initial, i
  10.  
  11. map = {"UUU":"F", "UUC":"F", "UUA":"L", "UUG":"L",
  12. "UCU":"S", "UCC":"s", "UCA":"S", "UCG":"S",
  13. "UAU":"Y", "UAC":"Y", "UAA":"STOP", "UAG":"STOP",
  14. "UGU":"C", "UGC":"C", "UGA":"STOP", "UGG":"W",
  15. "CUU":"L", "CUC":"L", "CUA":"L", "CUG":"L",
  16. "CCU":"P", "CCC":"P", "CCA":"P", "CCG":"P",
  17. "CAU":"H", "CAC":"H", "CAA":"Q", "CAG":"Q",
  18. "CGU":"R", "CGC":"R", "CGA":"R", "CGG":"R",
  19. "AUU":"I", "AUC":"I", "AUA":"I", "AUG":"M",
  20. "ACU":"T", "ACC":"T", "ACA":"T", "ACG":"T",
  21. "AAU":"N", "AAC":"N", "AAA":"K", "AAG":"K",
  22. "AGU":"S", "AGC":"S", "AGA":"R", "AGG":"R",
  23. "GUU":"V", "GUC":"V", "GUA":"V", "GUG":"V",
  24. "GCU":"A", "GCC":"A", "GCA":"A", "GCG":"A",
  25. "GAU":"D", "GAC":"D", "GAA":"E", "GAG":"E",
  26. "GGU":"G", "GGC":"G", "GGA":"G", "GGG":"G",}
  27.  
  28. def print_translation(mRNA, start, end):
  29. newseq = ""
  30. for i in range(start, end, 3):
  31. key = mRNA[i:i+3]
  32. newseq = newseq + map[key]
  33. print(newseq)
  34.  
  35. cur_pos = 0
  36. while 1:
  37. start, end = next_transcript(mRNA, cur_pos)
  38. if start == -1:
  39. break
  40. print_translation(mRNA, start, end)
  41. cur_pos = end
Add Comment
Please, Sign In to add comment