22 lines
508 B
Python
22 lines
508 B
Python
def MinimumSkew(Genome):
|
|
positions = []
|
|
skew = SkewArray(Genome)
|
|
minimum = min(skew)
|
|
for i in range(0, len(Genome)):
|
|
if skew[i] == minimum:
|
|
positions.append(i)
|
|
return positions
|
|
|
|
|
|
def SkewArray(Genome):
|
|
Skew = []
|
|
Skew.append(0)
|
|
for i in range(0, len(Genome)):
|
|
if Genome[i] == "G":
|
|
Skew.append(Skew[i] + 1)
|
|
elif Genome[i] == "C":
|
|
Skew.append(Skew[i] - 1)
|
|
else:
|
|
Skew.append(Skew[i])
|
|
return Skew
|