Azkiin

DNA reading program

Jun 17th, 2019
52
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.22 KB | None | 0 0
  1. DNA = input("Plak de DNA sequentie")
  2. A = DNA.count("A")
  3. T = DNA.count("T")
  4. C = DNA.count("C")
  5. G = DNA.count("G")
  6.  
  7.  
  8. print("De volledige sequentie is :", (A+C+T+G),"bp lang.")
  9. print("Het aantal A bp in de sequentie is :",A,"bp lang.")
  10. print("Het aantal T bp in de sequentie is :",T,"bp lang.")
  11. print("Het aantal C bp in de sequentie is :",C,"bp lang.")
  12. print("Het aantal G bp in de sequentie is :",G,"bp lang.")
  13. print("De GC ratio van de sequence is :" ,(G+C)/(G+C+T+A),'\n')
  14.  
  15.  
  16. def transcribe(sequence):
  17. return sequence.replace('T', 'U')
  18.  
  19. def translate_rna(s):
  20. codon2aa = {"AAA":"K", "AAC":"N", "AAG":"K", "AAU":"N",
  21. "ACA":"T", "ACC":"T", "ACG":"T", "ACU":"T",
  22. "AGA":"R", "AGC":"S", "AGG":"R", "AGU":"S",
  23. "AUA":"I", "AUC":"I", "AUG":"M", "AUU":"I",
  24.  
  25. "CAA":"Q", "CAC":"H", "CAG":"Q", "CAU":"H",
  26. "CCA":"P", "CCC":"P", "CCG":"P", "CCU":"P",
  27. "CGA":"R", "CGC":"R", "CGG":"R", "CGU":"R",
  28. "CUA":"L", "CUC":"L", "CUG":"L", "CUU":"L",
  29.  
  30. "GAA":"E", "GAC":"D", "GAG":"E", "GAU":"D",
  31. "GCA":"A", "GCC":"A", "GCG":"A", "GCU":"A",
  32. "GGA":"G", "GGC":"G", "GGG":"G", "GGU":"G",
  33. "GUA":"V", "GUC":"V", "GUG":"V", "GUU":"V",
  34.  
  35. "UAA":"_", "UAC":"Y", "UAG":"_", "UAU":"T",
  36. "UCA":"S", "UCC":"S", "UCG":"S", "UCU":"S",
  37. "UGA":"_", "UGC":"C", "UGG":"W", "UGU":"C",
  38. "UUA":"L", "UUC":"F", "UUG":"L", "UUU":"F"}
  39.  
  40. l = [codon2aa.get(s[n:n+3], 'X') for n in range(0, len(s), 3)]
  41. return "".join(l)
  42.  
  43.  
  44. RNA=transcribe(DNA)
  45.  
  46. print("De RNA sequentie is:")
  47.  
  48. print(RNA,'\n')
  49.  
  50. print("De Eiwitsequentie is:")
  51. protein=translate_rna(RNA)
  52. print(protein, '\n')
  53. print ("Het aantal aminozuren van het eiwit is:")
  54. print (len(protein)-1,'\n\n\n')
  55. print ("Deze software is het eigendom van Loeke van der Linden ,")
  56. print ("Deze Software kan daarom niet instaan voor de juistheid of (on)volledigheid van de gepubliceerde informatie.")
  57. print ("Loeke aanvaardt geen enkele verantwoordelijkheid of aansprakelijkheid voor eventuele schade, van welke aard dan ook, die op welke wijze dan ook voortvloeit uit het gebruik")
Add Comment
Please, Sign In to add comment