Advertisement
Guest User

Untitled

a guest
Nov 22nd, 2017
78
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 0.44 KB | None | 0 0
  1. # Reverse Complement
  2.  
  3.  
  4. def rc(read):
  5.     seq = read[1]
  6.     c = ""  # this will be the reverse complement of the input seq
  7.     r = seq[::-1]  # string is reversed here
  8.     for b in r:
  9.         if b == "A":
  10.             c += "T"
  11.         elif b == "T":
  12.             c += "A"
  13.         elif b == "C":
  14.             c += "G"
  15.         elif b == "G":
  16.             c += "C"
  17.     print(seq)
  18.     print(r)
  19.     print(c)
  20.     return c
  21.  
  22.  
  23. rc("ATCGCCCCCC")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement