Skip to content

Instantly share code, notes, and snippets.

Created December 23, 2015 01:00
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save anonymous/5138c61c41fb985c7837 to your computer and use it in GitHub Desktop.
Save anonymous/5138c61c41fb985c7837 to your computer and use it in GitHub Desktop.
sequence = 'TAAAGACTGCCGAGAGGCCAACACGAGTGCTAGAACGAGGGGCGTAAACGCGGGTCCGAT'
def skew(sequence):
c = 0
g = 0
min_skew = 0
skew_list = []
index = 0
for i in sequence:
index += 1
if i == 'C':
c += 1
if i == 'G':
g += 1
skew = g-c
if skew < min_skew:
skew_list = [index]
min_skew = skew
if skew == min_skew and index not in skew_list:
skew_list.append(index)
print(skew_list)
#with open('data.txt', 'r') as in_file:
#sequence = in_file.readline()
#skew(sequence)
skew(sequence)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment