Advertisement
Guest User

diamantcode

a guest
Feb 23rd, 2018
75
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.96 KB | None | 0 0
  1. def canonisch(codon):
  2. codon2 = codon[::-1]
  3. if codon[1] == "T":
  4. codon1 = codon[0] + "A" + codon[2]
  5. elif codon[1] == "A":
  6. codon1 = codon[0] + "T" + codon[2]
  7. elif codon[1] == "C":
  8. codon1 = codon[0] + "G" + codon[2]
  9. elif codon[1] == "G":
  10. codon1 = codon[0] + "C" + codon[2]
  11. codon3 = codon1[::-1]
  12.  
  13. L = [codon,codon1,codon2,codon3]
  14. L.sort()
  15. return L[0]
  16.  
  17. def codon2aa(codon):
  18. #stap1
  19. c = canonisch(codon)
  20.  
  21. #stap2
  22. n = c.replace("G","0").replace("T","1").replace("C","2").replace("A","3")
  23. w1 = n[0]
  24. w2 = n[1]
  25. w3 = n[2]
  26. p = (int(w1) + 4 * int(w2) + 16 * int(w3)) % 25
  27.  
  28. #stap3
  29. alfabet = "ABCDEFGHIJKLMNOPQRSTUVWXYZ"
  30. return alfabet[p]
  31.  
  32. def dna2eiwit(DNA):
  33. # codon bepalen
  34. str=""
  35. l = int(len(DNA))
  36. for i in range(0,l- 2):
  37. codon = DNA[i:i + 3]
  38. letter = codon2aa(codon)
  39. str += letter
  40. return str
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement