Advertisement
Guest User

Untitled

a guest
Mar 2nd, 2015
185
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.71 KB | None | 0 0
  1.  
  2. 2.
  3. #!/usr/bin
  4. file = open("seq.txt")
  5. seq = file.readline().strip()
  6. dic = {"A":"T" , "T":"A" , "G":"C" , "C":"G"}
  7. lista = []
  8. for n in seq:
  9. n=dic[n]
  10. lista.append(n)
  11. print "".join(lista)
  12.  
  13. 3. zliczanie wystapien
  14. #!/usr/bin
  15. file = open("seq.txt")
  16. seq = file.readline().strip()
  17. #dic = {"A":"T" , "T":"A" , "G":"C" , "C":"G"}
  18. lista = []
  19. liczA=0
  20. liczT=0
  21. liczG=0
  22. liczC=0
  23. hexamer = []
  24. for n in seq:
  25. se = seq[
  26. if n=='A':
  27. liczA = liczA+1
  28. if n=='T':
  29. liczT = liczT+1
  30. if n=='G':
  31. liczG = liczG+1
  32. if n=='C':
  33. liczC = liczC+1
  34. #print "".join(lista)
  35. print "Tyle jest A:"
  36. print liczA
  37. print "Tyle jest T:"
  38. print liczT
  39. print "Tyle jest C:"
  40. print liczC
  41. print "Tyle jest G:"
  42. print liczG
  43.  
  44. 4. hexamery
  45. #!/usr/bin
  46. file = open("seq.txt")
  47. seq = file.readline().strip()
  48. lista = []
  49. dl = len(seq)
  50. dic = {}
  51. for n in range(0, (dl-5)):
  52. a=seq[n:n+6]
  53. if a not in dic:
  54. dic[a]=0
  55. dic[a]+=1
  56. print dic
  57.  
  58. 5.tlumaczenie na bialka
  59. #!/usr/bin
  60. file = open("seq.txt")
  61. seq = file.readline().strip()
  62. lista = []
  63. dl = len(seq)
  64. dic = {}
  65. f = open("ramki.txt",'r')
  66. line1 = f.readlines()
  67. for line in line1:
  68. string=line[0:3]
  69. klucz=line[4]
  70. dic[string]=klucz
  71.  
  72. lizda = []
  73. print "ramka odczytu od 0"
  74. for n in range(0, dl,3):
  75. a=seq[n:n+3]
  76. if a in dic:
  77. lizda.append(dic[a])
  78. print lizda
  79. lizda = []
  80.  
  81. print "ramka odczytu od 1"
  82. for n in range(1, dl-1,3):
  83. a=seq[n:n+3]
  84. if a in dic:
  85. lizda.append(dic[a])
  86. print lizda
  87. lizda = []
  88.  
  89. print "ramka odczytu od 2"
  90. for n in range(2, dl-2,3):
  91. a=seq[n:n+3]
  92. if a in dic:
  93. lizda.append(dic[a])
  94. print lizda
  95. lizda = []
  96.  
  97. print "ramka odczytu od N"
  98. for n in range(dl, 0,-3):
  99. a=seq[n:n+3]
  100. if a in dic:
  101. lizda.append(dic[a])
  102. print lizda
  103. lizda = []
  104.  
  105. print "ramka odczytu od N-1"
  106. for n in range(dl-1, 0,-3):
  107. a=seq[n:n+3]
  108. if a in dic:
  109. lizda.append(dic[a])
  110. print lizda
  111. lizda = []
  112.  
  113. print "ramka odczytu od N-2"
  114. for n in range(dl-2, 2,-3):
  115. a=seq[n:n+3]
  116. if a in dic:
  117. lizda.append(dic[a])
  118. print lizda
  119. lizda = []
  120.  
  121. 6.generowanie
  122. #!/usr/bin
  123. import random
  124.  
  125. dlugosc = raw_input('Podaj dlugosc sekwencji :')
  126. dlugosc=float(dlugosc)
  127. prawdC = dlugosc * 0.25
  128. prawdC=int(prawdC)
  129. prawdT = dlugosc * 0.25
  130. prawdT = int(prawdT)
  131. prawdG = dlugosc * 0.25
  132. prawdG=int(prawdG)
  133. prawdA = dlugosc * 0.25
  134. prawdA=int(prawdA)
  135.  
  136. literki = ["A", "T", "C", "G"]
  137. sekwencja = []
  138.  
  139. for i in range(0, prawdC-1):
  140. sekwencja.append(literki[0])
  141. for i in range(0, prawdT-1):
  142. sekwencja.append(literki[1])
  143. for i in range(0, prawdG-1):
  144. sekwencja.append(literki[2])
  145. for i in range(0, prawdA-1):
  146. sekwencja.append(literki[3])
  147.  
  148. from random import shuffle
  149. shuffle(sekwencja)
  150.  
  151. print sekwencja
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement