-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathSUBS-Finding a Motif in DNA.py
20 lines (17 loc) · 1.2 KB
/
SUBS-Finding a Motif in DNA.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
'''
Rosalind Problems: [SUBS] Finding a Motif in DNA
'''
def locate_motif(dnaseq, sub):
sub_len = len(sub)
dna_len = len(dnaseq)
frame = dna_len - sub_len +1
location = []
for i in range(frame):
if dnaseq[i:i+sub_len] == sub:
location.append(i+1)
return location
dna = 'CTGTCAGCCTGTCAGTCGCCTTTGACTGTCAGCCTGTCAGGCTGTCAGACTGTCAGTCAGCTGTCAGGACCTGTCAGATCTGTCAGCATTTTCTGTCAGGGGAATCCCCTGTCAGGTGGGTTCTGTCAGACCAACCTGTCAGGTGCTGTCAGTACTGTCAGCTGTCAGCTCTGTCAGTATCATGCGCTGTCAGCTGTCAGCTGTCAGCTGTCAGTCTGTCAGCACTGTCAGACTGTCAGACTCAGCTGTCAGGAATTACTACTGTCAGCGGTCGTGCCTGTCAGCCAGACTGTCAGGCTGTCAGCTGTCAGCTGTCAGCAACTGTCAGTGGGGCCCCGCTGTCAGCGCGCTGTCAGGATCTGTCAGGGTTGCTGTCAGCTGTCAGCTGTCAGTCCTTTTTAGACTGTCAGACCTGTCAGCCTGTCAGCGCCCAGGTCTGTCAGGCTGTCAGCAGGGTGTGAACTGTCAGCCTGTCAGACCTGTCAGACCCTGTCAGCTGTCAGCTGTCAGCCGCCTGTCAGCGTCCCTGGCCTGTCAGGGCCCTCTGTCAGCTGTCAGGCCTGTCAGACTGTCAGGGATCGCCTGTCAGCTGTCAGACTGTCAGCTGTCAGCCTGTCAGGTTTGACTGTCAGCTGTCAGAAACTGTCAGACTGTCAGTCTGTCAGCTCTGTCAGCTGTCAGGCTGTCAGCAACTGTCAGGCTACTGTCAGCTGTCAGGCGGCACTGTCAGGCTGTCAGCTGTCAGCTGTCAGGTCACTGTCAGTCTGTCAGGACTGTCAGCTGTCAGATGGCTGTCAGCATACTGTCAGACCTGTCAG'
substring = 'CTGTCAGCT'
location_list = locate_motif(dna,substring)
for i in location_list:
print(i, end=' ')