-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathFrequencyMap Function
More file actions
30 lines (29 loc) · 819 Bytes
/
FrequencyMap Function
File metadata and controls
30 lines (29 loc) · 819 Bytes
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
def FrequentWords( str, int ):
numb = {}
for i in range(len(str)-int+1):
if str[i:i+int] not in numb:
numb[str[i:i+int]] = 0
numb[str[i:i+int]] += 1
maxim = max(numb.values())
res = []
for kmers in numb:
if numb[kmers] == maxim:
res.append(kmers)
return res
str = 'ACGTTACGTTACGTCGACTGCACTACG'
int = 3
print(FrequentWords( str, int ))
# Copy your FrequencyMap() function here.
def FrequencyMap(DNA, int):
freq = {}
n = len(DNA)
for i in range(n-int+1):
Pattern = DNA[i:i+int]
freq[Pattern] = 0
for i in range(n-int+1):
if (Pattern == DNA[i:i+int]):
freq[Pattern]=freq[Pattern]+1
return freq
DNA= 'ACGATTTACGAACGTTACGTAGCTAGCTACGTAC'
int = 3
print (FrequencyMap(DNA, int))