-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathfinding_motif_in_dna.py
6 lines (5 loc) · 1.08 KB
/
finding_motif_in_dna.py
1
2
3
4
5
6
seq = 'GCGAACAGTCGTGACGAAGCGAACAGCGAACACGCGAACAGGCAGTGAGATCCACAGGGCGAACAGAGCGAACAGGGCGAACAAGTATGGGCGAACAGCTTTGCGAACACGCGAACAGTGCGAACAGCGAACAGCGAACAATGCTCATCCGGGCTCAGCGAACAAGCGAACAGCGAACAGCGAACAATCGGCGAACAGTTGCGAACATCCCAGCGAACAGAGGCGAACACCGCGAACACCCACGCGAACACCGGGCGAACAGCGAACAGCGAACATTAAGCGAACACGTTAGCGAACAGAGCGAACACGCGAACAACGCGAACATGGCGAACATCCACGTTAGCGAACAGGTCCGGATGGCGAACAGCTGATCGGGACAATGCGAACACAGGCGAACAGCTTTTGCGAACAAGCGAACAGCGAACACTAGCGAACACGCGAACAGTTTCGGCGAACACTGCGAACAGCGAACACGCGAACAATGCGAACATACCGCAACGCGAACAAGCGAACAAAAATGCGAACAGCGAACAATAAGACCTAGCGAACAGCCATCCGGGCGAACAGGCCGAGCGAACATTATTGTGCGCGAACAGTAGCGAACACGGCGAACAGCGAACAGCGAACAGGCGAACATAGCGAACAGCGGCGAACACAAATAACTGCTGCGAACAGCGAACAGGGCGAACACGCGAACAGCGAACAGCGAACAGCGAACAGCGAACATAGAGCGAACAACAAGCGAACAGTATCTGCGAACACTACGCGAACAGCGAACATACGCGAACAACGTCGCCGCTGCGAACATGGGGCGAACAATAGCGCGAACAGCGAACAGTCAGCGAACAGCGAACATTTCATAGCGAACAGATCTCATAGGCGAACATCGGATCTAAGGGCGAACAAGCGAACAGGAGTAGCCGCGAACAAGCGCCTTTGCGAACAGCGAACAGCGCGAACAGCGAACATGCGTAGCGAACACAGCC'
motif = 'GCGAACAGC'
for i in range(len(seq)):
if motif in seq[i:i+len(motif)]:
print(i+1)