Advertisement
chilikla

Untitled

Jan 20th, 2020
118
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 2.52 KB | None | 0 0
  1. WT="ATGCTGCTATCCGTGCCGCTGCTGCTCGGCCTCCTCGGCCTGGCCGTCGCCGAGCCTGCCGTCTACTTCAAGGAGCAGTTTCTGGACGGAGACGGGTGGACTTCCCGCTGGATCGAATCCAAACACAAGTCAGATTTTGGCAAATTCGTTCTCAGTTCCGGCAAGTTCTACGGTGACGAGGAGAAAGATAAAGGTTTGCAGACAAGCCAGGATGCACGCTTTTATGCTCTGTCGGCCAGTTTCGAGCCTTTCAGCAACAAAGGCCAGACGCTGGTGGTGCAGTTCACGGTGAAACATGAGCAGAACATCGACTGTGGGGGCGGCTATGTGAAGCTGTTTCCTAATAGTTTGGACCAGACAGACATGCACGGAGACTCAGAATACAACATCATGTTTGGTCCCGACATCTGTGGCCCTGGCACCAAGAAGGTTCATGTCATCTTCAACTACAAGGGCAAGAACGTGCTGATCAACAAGGACATCCGTTGCAAGGATGATGAGTTTACACACCTGTACACACTGATTGTGCGGCCAGACAACACCTATGAGGTGAAGATTGACAACAGCCAGGTGGAGTCCGGCTCCTTGGAAGACGATTGGGACTTCCTGCCACCCAAGAAGATAAAGGATCCTGATGCTTCAAAACCGGAAGACTGGGATGAGCGGGCCAAGATCGATGATCCCACAGACTCCAAGCCTGAGGACTGGGACAAGCCCGAGCATATCCCTGACCCTGATGCTAAGAAGCCCGAGGACTGGGATGAAGAGATGGACGGAGAGTGGGAACCCCCAGTGATTCAGAACCCTGAGTACAAGGGTGAGTGGAAGCCCCGGCAGATCGACAACCCAGATTACAAGGGCACTTGGATCCACCCAGAAATTGACAACCCCGAGTATTCTCCCGATCCCAGTATCTATGCCTATGATAACTTTGGCGTGCTGGGCCTGGACCTCTGGCAGGTCAAGTCTGGCACCATCTTTGACAACTTCCTCATCACCAACGATGAGGCATACGCTGAGGAGTTTGGCAACGAGACGTGGGGCGTAACAAAGGCAGCAGAGAAACAAATGAAGGACAAACAGGACGAGGAGCAGAGGCTTAAGGAGGAGGAAGAAGACAAGAAACGCAAAGAGGAGGAGGAGGCAGAGGACAAGGAGGATGATGAGGACAAAGATGAGGATGAGGAGGATGAGGAGGACAAGGAGGAAGATGAGGAGGAAGATGTCCCCGGCCAGGCCAAGGACGAGCTGTAG"
  2.  
  3. def delete(index1,index2):
  4.   index1=int(index1)
  5.   index2=int(index2)
  6.   length=index2-index1+1
  7.   begin=WT[:index1]
  8.   end=WT[index1+length:]
  9.   print(begin+end)
  10.   return begin+end
  11.  
  12. def addition(index1, nucleotides):
  13.   index1=int(index1)
  14.   first_list = WT[:index1+1]
  15.   second_list = WT[index1+1:]
  16.   print(first_list+nucleotides+second_list)
  17.   return first_list+nucleotides+second_list
  18.  
  19. def main():
  20.   my_string="c.1100_1145insTTAAGGAGGAGGAAGAAGACAAGAAACGCAAAGAGGAGGAGGAGGC"
  21.   temp=my_string.split("c.")[1]
  22.   index1=temp.split("_")[0]
  23.   print(index1)
  24.   fullIndex2=temp.split("_")[1]
  25.   if "del" in fullIndex2:
  26.     index2=fullIndex2.split("del")[0]
  27.     print(index2)
  28.     delete(index1, index2)
  29.   elif "ins" in fullIndex2:
  30.     index2=fullIndex2.split("ins")[0]
  31.     print(index2)
  32.     code=fullIndex2.split("ins")[1]
  33.     print(code)
  34.     addition(index1, code)
  35.  
  36. my_string="c.1100_1145insTTAAGGAGGAGGAAGAAGACAAGAAACGCAAAGAGGAGGAGGAGGC"
  37. temp=my_string.split("c.")[1]
  38. index1=temp.split("_")[0]
  39. print(index1)
  40. fullIndex2=temp.split("_")[1]
  41. if "del" in fullIndex2:
  42.   index2=fullIndex2.split("del")[0]
  43.   print(index2)
  44.   delete(index1, index2)
  45. elif "ins" in fullIndex2:
  46.   index2=fullIndex2.split("ins")[0]
  47.   print(index2)
  48.   code=fullIndex2.split("ins")[1]
  49.   print(code)
  50.   addition(index1, code)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement