Advertisement
Guest User

Untitled

a guest
Dec 18th, 2017
113
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 2.12 KB | None | 0 0
  1. #kodony stop UAA UAG UGA
  2. #kodon start AUG
  3.  
  4. gb = 'sequence2.gb'
  5. START = 'AUG'
  6. STOP1 = 'UAA'
  7. STOP2 = 'UAG'
  8. STOP3 = 'UGA'
  9. origin_flag=False
  10. sequence = ''
  11. seq=''
  12. pozycja = ''
  13. slownik = {"UUU":"F", "UUC":"F", "UUA":"L", "UUG":"L",
  14.     "UCU":"S", "UCC":"s", "UCA":"S", "UCG":"S",
  15.     "UAU":"Y", "UAC":"Y", "UAA":"STOP", "UAG":"STOP",
  16.     "UGU":"C", "UGC":"C", "UGA":"STOP", "UGG":"W",
  17.     "CUU":"L", "CUC":"L", "CUA":"L", "CUG":"L",
  18.     "CCU":"P", "CCC":"P", "CCA":"P", "CCG":"P",
  19.     "CAU":"H", "CAC":"H", "CAA":"Q", "CAG":"Q",
  20.     "CGU":"R", "CGC":"R", "CGA":"R", "CGG":"R",
  21.     "AUU":"I", "AUC":"I", "AUA":"I", "AUG":"M",
  22.     "ACU":"T", "ACC":"T", "ACA":"T", "ACG":"T",
  23.     "AAU":"N", "AAC":"N", "AAA":"K", "AAG":"K",
  24.     "AGU":"S", "AGC":"S", "AGA":"R", "AGG":"R",
  25.     "GUU":"V", "GUC":"V", "GUA":"V", "GUG":"V",
  26.     "GCU":"A", "GCC":"A", "GCA":"A", "GCG":"A",
  27.     "GAU":"D", "GAC":"D", "GAA":"E", "GAG":"E",
  28.     "GGU":"G", "GGC":"G", "GGA":"G", "GGG":"G",}
  29.  
  30.  
  31. with open(gb, 'r') as file:
  32.     for line in file:
  33.         line.strip('\n')
  34.         if line.startswith('ORIGIN'):
  35.             origin_flag = True
  36.             continue
  37.         else:
  38.             pass
  39.         if origin_flag:
  40.             sequence_line = line.strip()
  41.             sequence += ' '.join(line.replace('t', 'u').replace(' ', '').upper().split())
  42.             sequence = ''.join([i for i in sequence if not i.isdigit() and i.isalnum()])
  43.  
  44.  
  45.  
  46. chain_index = 0;
  47. starts = []
  48. stops = []
  49.  
  50. #sequence = 'AUGAAAAAAAAAAUGUUUUUUUUUUAACCCCCCUAA'
  51.  
  52. for i in range(0,len(sequence),3):
  53.     str = sequence[i:i+3]
  54.     if(str == START):
  55.         print('START {}'.format(i))
  56.         starts.append(i)
  57.     elif(str == STOP1 or str == STOP2 or str == STOP3):
  58.         print('STOP')
  59.         stops.append(i)
  60.        
  61. print("------")
  62. print(starts)
  63. print(stops)
  64.  
  65. seq = []
  66.  
  67.  
  68. for i in range(0,len(stops)):
  69.     for j in range(0,len(starts)):
  70.         seq.append(sequence[starts[j]:stops[i]+3])
  71.        
  72.        
  73. #print(seq)
  74. seq.sort(key=len)
  75. seq.reverse()
  76. #print(seq)
  77.  
  78. translated = []
  79. for str in seq:
  80.     tmp = ''
  81.     for x in range(0,len(str),3):
  82.         tmp += slownik[str[x:x+3]]
  83.     translated.append(tmp)
  84.    
  85.    
  86. print(translated)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement