-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathFind a Median String.py
More file actions
43 lines (33 loc) · 1.08 KB
/
Find a Median String.py
File metadata and controls
43 lines (33 loc) · 1.08 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
import sys
def hammingDistance(s, t):
return len([a for (a,b) in zip(s,t) if a!=b])
def distanceBetweenPatternAndStrings (pattern,dna):
def hamming(pattern,genome):
return min([hammingDistance(pattern, genome[i:i + len(pattern)])
for i in range(len(genome)-len(pattern)+1)])
return sum([hamming(pattern,motif) for motif in dna])
def k_mers(k):
if k<=0:
return ['']
else:
result=[]
for ks in k_mers(k-1):
for l in ['T','G','C','A']:
result.append(ks+l)
return result
def medianString(k,dna):
def findClosest(d):
distance=sys.float_info.max
closest=None
for k_mer in k_mers(k):
if distance>d(k_mer,dna):
distance=d(k_mer,dna)
closest=k_mer
return closest
return findClosest(distanceBetweenPatternAndStrings)
k=6
f = open("C:/Users/iMan/PycharmProjects/RosalindExercises/rosalind_ba2b.txt","r")
text = f.read()
dna =[]
dna = text.split("\n")
print(medianString(k,dna))