-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathpattern clusters.py
More file actions
18 lines (16 loc) · 9.44 KB
/
Copy pathpattern clusters.py
File metadata and controls
18 lines (16 loc) · 9.44 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
genome = 'AGGATTGGCAGTACTGCTGCGTGAACACGCGATGCGTTTCGGGTAACGACTTCCACCCGCTTAAAGATAACAAAAATGTCTAATAAACAATCGGTTGTATGAAATTATGCCAGTCTAGCCCAAGAGCATTACATCTGGTTCCATAAATACGGTCGGCTTGTCGCCAGCGGCTTAGTCGCGGATTAGCGTAGTCCGAGGTTACCCTTTCCCTCCCACTCGAACGGACGATTGTGCCTACATCTTGGTCTGCTTACCAATTCCATACGCAGCCATGCCTAATTCAGCTCTGTAGGCTATGGTCACCACGTAGGAACTATTCGACTGTACCCACAAGCGCGTGATTCGGATATTAAACTCAGGAGGTCACTTGGCGTTCATGGATACCACACAGAGCTATTTATGGATATTAAAACCCAGAATACGATACAGGATGTAATTGACCCATACTGCTTATGAAAAATGTGCCTCTGAGCTCATCGACGATAGTTTATGTACTCGGTTTCTGGGTTATGCTGCGGTCCGTAATCCATCTCCGCCGCGGGAGCCTAACTCCCTCCACTAGGAAGTGACCATACGAGTGTGCGGGTTGGGTCACCATAAGGCGCTGGTCAACTGATGAGTCAGACCTGGACGCGGTACCTCTCTAGTTCTCCTTTGCGTGCCCGGTTGAAGATCACGATGTGTTCGTCCCAGGAAACGGCGGCTGGCTGGCTCGGGTAACTAACCCTCAGAGCTCACAGAGGGGTGCTACGCACGGAAGATATAACCTGTGACTTATACGTCAACTCTACTGTATGCTTCGTCCGGCCTCACTGTATTGACCAGCACGTTGCATGGGAAACATGAGGACACAGCGTGACAGACTTTCATAGGCGCCGCCTCTGAGGGGCGCGCCCGGTAGTACCACTAAGCGCCCTCGAAAGCTGTTCTGGACGCCTGGTTCGGGGCGTATTAGGGATTGAATAGGGTAATGCCTGCAAGACTGGCATAGGCTCGAGTTTTTCTTCAGATTAACATCCTTAGACAAAGTAAGGTAAATGTTACGCGTTAGGCTATGCAACGTCCCACAATTAAATGTCCACGTCGAATTGTGGCTCGTGATGTGTACTTCAATCAGGGTAGATAGTCGTAGCCAAGTTGTAGAGCATCACCAACTGTTACCACCAACTGTTACCATAACAACGCCAACTGTTACCGTCGATCCAACTGTTACCGGATGTTCTTAGGACCCGTCGCCAACTGTTACCGTTCCAACTGTTACCTTACCACTAATTGGTTAGCTCCAACTGTTACCTATGGACCAACTGTTACCGCCAACTGTTACCCGCCGCCAACTGTTACCGGAGAAACACCTAACTTCTGGAGAAACACCTCCCCAACTCCAACTGTTACCCTGGTGGAGAAACACCCCGGGTCTTGGAGAAACACCAACTGTTACCTCCAACTGTTACCACGAGGGAACAACCTCCAGCAACCTTGACTGCCACTTACCAACTGTTACCTAATTGTTCCAACCCAACTGTTACCAACCCAACTGTTACCGCTGAATGCCCAACTGTTACCCACTGCCAACCCCAACTGTTACCTAACCCCCAACTGTTACCTTACCGAGAAACACCCTTTCATGGAGAAACACCCACGTAAGGGTGCTGCATCTGCCAACCCAACTGTTACCAACTGTTACCCAACTGTTACCAACACCCACCCCGCCCAACTGTTACCTGCATCTGGAGAAACACCATGTACTGTGGTGTTCTGGAGAAACACCTGATGGAGAAACACTGGTGGAGAAACACCCCCTGGAGAAACACCGGCGGTTGTGTGGAGAAACACCCCTGGAGAAACACCTGCTGCATCTGTGCTGCATCTGAATTGCTGCATCTGTGGAGAAACACCGCATCTGACTGTGCTGCATCTGAAAAAGCTCTGTAGCCCCACTAATTGTCGATGCTGCATCTGAGATTTGTGAGCGCGAGGCGTGCTGCATCTGAGGTGTTGTTCGCTATCTAGTTTGAGCTGATTCGCTGCTGCATCTTGCTGCATCTGATGGATAATGCTGCTGCATCTGATCTTGCTGCATCTGATTGGTTCCTGATGATGCTGCATGCTGCATCTGACTGAGACTTATAGGTGGTCCCCAGTACCATTTAGAATAAAATCGGAACTCCCGCTGCGGCACTGTTCGTCTGGTTTCACAGTATCCGATAGCGCCTCACATAGATCCTATAATAGTATCCTACCCCTTTGGAGCTCTTTGAGGCTATGACCACGTTAATTCCTGAGCATACATCAATCAACGATGACGTACTGATCCGGGACAGGCATATTCTCGGTGATCCGGGACAGTATGTGGCGATCCGGGACGATCCGGGACAGTCCCATGCAGCTGGTAATCCAGTCGTGCTGACCGATCCGGGACAGCCGGGACAGGCAGATTGTCGGATAGGATCCGGGACAGCGGGACAGGCTATGGGCTTCCCGATCCGGGACAGAATGGTGCGCGGATCCGGGACAGTCGAGACTCCGCCGGGGTAGATCCGGGACAGCAGATGTAGGGCCGTAGGAGCTCTGGTACAAGATCCGGGACAGCGATCCGGGACAGGTACACCACCGCTGAGTTGGACCGGATGCCTCGATCCGGATCCGGGACAGTCGGCGGATGCTCACTGTCCTTAGCTTGTTCATGCAGATCCGGGACAGGGGATCCGGGACAGTTAAAGCAACGTTCCCTAATTGGTGGCAAGCAACGTTCCAAAGCAACGTTCCGACAGCGAAGCAACGTTCCCGGGACAGCAAGCAACGTTCCACAGGACAGCGAAGCAAGCAACGTTCCAACGTTCCAAGCAACGAAGCAACGTTCCAAGACGCCGATCCGGAAGCAACGTTCCCTTAAGTGAAGCAACGTTCCATATCTGCCGCGCGCCGCTCAATAATAAGCAACGTTCCTAGGACGTTGCATAGGTACGAGTAGCCTCGAGCTCAGGAATAAGTGAGCAACGCAAAAAGCAACGTTCCTGTAAACTATAGAAAGCAACGTTCCTTCCTATCTCTCTTGAGGAAGCAACGTTCCAGCAACGTTCCCGTTCCATATGGTTTTCCCCTTGTAAAGCAACGTTCCTCCTTCCGAAAAACACGACGATCTTATGAAGCAACGTTCCCGCTAAATACGTAAGCAACGTTCCAAGCAACGTTCCGGAAGCAACGTACTTGGCAGGACACGTTCCTTCCAGACTTGGCAGGACAACGTTCCGAAAAGCAACGTTCCTGGCCCGTATCCGCGCCAAGTAACTTGGCAGGACTCTAGTACATCCGCGCCAAGTCGCGCCAAGTCCAATCCGCGCACACTTGGCAGGACTTGGCAGGACGCCAAGTTAGGAAAACTTCCTGGTACTTGGCAGGACATTCCCTACCGTATACTTGGCAGGACACGACCAAGTAATATAACTTGGCAGGACACTGGCAGGACCTTGGCAGGACGTCGCCAAGTACTTGGCAGGACAGGACCCCCTCTAGGACGAGCTTGTCTCCAACTTGGCAACTTGGCAGGACAGTAGTCCCTCCCCTCTCCGCACTTGGCAGGACCGCGCCAAGTCCATTTCCGCGCCAAGTACTCCGCGCCAAGACTTGGCAGGACCAGGACCCTTGGCAGGACAGACTTGGCAGGACTTGGCAGGACTTGGCAGGACGCGTTGCCTTAGCTGCGATCTCCGCGCCAAGTGCTAGGACACTTGGCAGGACAAGTTATTGAACTGCAACTTGGCAACTTGGCAGGACTTGGCAGGACTTGGCAGGACCGCGCCAAGTTCAAGTAGGTTCCGTCCGCGTCCGCGCCAAGTAAGTGCTGACCACATCACACTGTTTTCCCTCTAGGACGTTTATGGCCACCCTCTAGGACGGCAGCCCTCTAGGACGGCGTTAGTTTATATTCCGCACGTTAGCGGTTACCCCTCTAGGACGAGGACCCCCTCTAGGACGGGACCCCTCTAGGACGGGACGATCCCTCTAGGACGATCCCCTCCCCCGTAGTAGGCGGTAGTAGGCGTCCGTAGTAGGCGTAGCTCCGTACCGTAGTAGGCGGTAGGCGGTATCTGCCCGTAGTAGGCGCGTACCGTAGTAGGCGGTCTTGTCTGGCTTACTCGGGTGCCCGTTCCCGTAGTAGGCGAGGCGACCACTTATATTGCCGTAGTAGGCGACCGTAGTAGGCGACTGAATATGGCCGACTTCGCCTGTTCTGTTCCATGCTCCGTAGTAGGCGCCGATCACGCGCGAGGACGTTCGACTTGTCACCGTAGTAGGCGTCCAACGACCGTAGTCCGTAGTAGGCGCGTACCAAAGCCGTAGTAGCCGTAGTAGGCGCCCCCTAATTGCTAACTGCTACTATCGATTGCCCGTAGTAGGCGCCATTGGACCTGCAAAAATTGTGATACGCCACCGTCCGTAGTAGGCGTAGGCGCGAGACTGTGGGAACCGTAGTAGGCGGTAGTAGGCGGGCGCAACGTAAACCTCGCCTTGACGAGAGTAAAGCTGCAACCGTAGTAGGCGAGGCGGCAGCCGAGCCGTCCGTACGTTAGACTTGCAGCTAGCTTGGTACTAACAGTGTTGATAATCTCCGCCAAGCAACGTTTCACATAGGGGCTGGGGTCAAAAGTTAAAACCCGCGCGAATTTCTCGCACCATGTTTGGCGCGGGTACTGTCCTAAGGTAGTCGGTGGGCTCACGAACTCGGTCCCGTCCTGCTTCGAGGGCTTGTATCCAACATTCCGGAAGCTTGACTCTAGGACTGCCATCCGATAGCGGACTCTAACATAACTGTTTTAAACTTGCGACTCGACCTAGCACGAGACGTCGGCCAGTGATGCTCACACCTACAGTTGCATCCCCCGATCAGAACGGGCGTTTTGTCATGGTCTCGAATATGGCAATAATACCGCGACTACGTGGTCCCGCTAAGTATTCAATTACTAGGTCAGTCATGCGGTATCCCGATAATTCTGACTACGCTTTTAAGGCGCGTGCGAATTAGCCTCTAGGGCTCCTATTCATCTAGACAAAATTAGTCTCGTGCCGCGTCGAAGTACACTCCAGCCGGGTTTTCCACACGACGGATATCCGAGATCCTTAGAAGCCGGGTTCTGTCGATTGTTAAATACTTCGCAAGGGAAGGTTGTCGCGTCGAATACCTCTATCAGGATTGTTGCTCACTATAGTCCGGAAAAGGGATAGTCCGGAAACAAAACTTAGTTAGCAAGACCAATAGTCCGGAAACTTCAAACAGGTAAGATAGTCCGGAAAGAATAGTCCGGAAAGTCCGGAAAGTAGTATTGTCATTTAGGATAGTCCGGAAAAAAAAATGGGACATAGTCCGGAAAATAGTCCGGAAAGGAAAGACAGTACCTCAGCCCACCCAGCAACCGGATTGGTCAGGAGATAGTCCGATAGTCCGGAATAGTCCGGAAAACTCTGTATCATGGGATTTCTCATAGTCCGGAAAGGGGTATCGGAATAGTCCGGAAAGGACATAGTCCGGAAAAAACATTATGCCGTATAGTCCGGAAAGTGACCGGCATAGTCCGGAAATAATAGTCCGGAAAAGGTGTGGATATGTTTTTCACATAGTCCGGAAAATAGTCCGGAAAAATCAAACTCAAACTCGCTCTGTATAGTCCGGAAACACCCAGAATAGTCCGGAAAAATCTGTATAGTCCGGAAACTCTCTAGCAACTATGTCGTATAGTCCGGAAAAAAATTTCCGGTCCCATCCATAGTCCGGAAACTATAGTCCGGAAAGTCCGGAAAGGAATCGAATTCAGTGTGACGCAGACCGATGTTGTTTCCAGGGAGTGTTTCCAGGGAATTGTTCTAATGGTGGATTCAGAGAACTTAGAGTTCTTACCCCGCGTGAACAATGTGCCAGTCTGTTTCCAGGGATCCAGGGAGAACTAAACCTGTTGTTTCCAGGGAGGCTGATTCATAATCAATCGCGTGTTTCCATGTTTCCAGGGATATGTCGGGACTCATAAGATCATCCGGGTAACGTTTTTGTTTCCAGGGAGTTTCTGTTTCCAGGGATTGTTTCCAGGGAGTATATGTTTCCAGGGATCCATGTTTCCAGGGACACACCCGACTAAGTGTTTCTGTTTTGTTTCCAGGGATCTTGTCGTCCCTCTTTCCAGGGATGTCGTCCCTCTTACTGTGAATGTTTGTCTGTCGTCCCTCTTTATGTTTCCAGGGATTGCGACTGTCGTCCCTCTTGTCGTCCCTCTCCTCTCCCTCTCTTCCTCTTCTTGTTTCCAGGGATTCGTGTCGTCCCTCTACGCAGATTACTGTCGTGTCGTCCCTCTCTGTCGTGTTGTCGTCCCTCTAGGTTGTCGTCCCTCTGTCGTCCCTCTTCCCTCTCTCTTCTTAGGTTGCAGCGAGAGTCGATGTCGTTGCCTGACAACTCGGGCTACGCATTCCCCTTGGTCTGTCGTGTCGTCCCTCTCAATAAACAAGCTGTCGTCCCTCTTGTCGTCCCTCTATACCTATTTTCGGTTGGTGACACAAATGGCATTTTGTCGTCCCTCTTGTCGTCCCTCTTCCCTCTTCTGTCGTCCCTCTTTGGATGGTAATCCATTAGCACCTAATCCATTAGCACCATTAATCCATTAGCACTGATGGGCTGTTATCCATTAGCAATCCATTAGATCCATTAGCACGTCCCTATCCATTAGCACCGTACTGTCGTCCCTCTCCATGATGTGATGGCCAATCCATTAGCACATCCATTAGCACATGATGCTGAACCCCTTTATTGCTCTCAGGAAGCCAATCAATCCATTAGCACTAGCACCATCATCCATTAGCACACGGCCATGATGGATATCCATTAGCAATCCATTAGCACCATCCATTAGCACCCCATGATAATCCATTAGCACATGTTTCTACGCCATGATCCATTATCCATTAGCACATGCATGATCCATTAGCACCACGATGGCCATGATGGCCATGATGTGATGCGGAATGGCCATGATGCATATGTGCCGACGCGCTACGTGACCACTCCCATTGTAGCAATCCATTAGCACCATCCATTAGCACCACCCGGCCCCCAACTGTCCTTGCGAGACTATCCATTAGCACATGAGTATCCATTAATCCATTAGCACATGATGCTCTCTATCCATTAGCACGATGGCCATGATGGGCAGTCAATTCTAGAAGCTCTCTCTAGACATACGATCGCTTGCTATAATTATCCTGGGCGAGTAAAGCTCATTAAGAGATACATTGATAGGACATAAGCTAGTCGCTTAGGTATCTTCCGCCACGACCAAAGGTGCATCACGATACTCTTTTAAGAAAGGAGAAAAGCTTGCGCAAAGTCCCTAAGCGTCCGAGAGGCAAGGTTTCTGTCTAACAATATCAAGCCGTCCCCCATGGCAGAGGGAGATCGCCACAGGAGCGCAGGTGATGGTAAGCATCCCCCGGGAACGTGAGACCGCGAACGCCGCACTTACTGGAGGGATTTATCCTTTATCCCGTCCAACAACCATGGAGACAAGAGGAATAGGAATCATCTACCAGCGGTTACTGATCCTGCGATCCCATGTGTTTGGTTTGGGTCTGCCGGATCATCATTCCCGAACCGCCCCATGGAGAGTTTAAAAAGGGAGGTAGTAATACTGGGTCAAGTTGGCCACCTCATAAACCCCTCATAAACTTTTCTGTGTACGGGACGACTGCCCGCTCCCCTCATAAACTTGTACATTACTGATATAATCTAGATGCTCCCCTCATAAACCCTCATAAACTTCTCCTCATAAACCTCATAAACTTGGCCCGCACAGCATAAGCCTCATAAACTTATTAGCGGGAATTTTTTTATTCCCTCATAAACTTTGTAGAAAGGACTTCTTGTTTGCCTTCCTCATAAACTTCCTCATAAACTTCGCAGACACCCTCATAAACTTATAAACTTCATAAAGTCAGCCTCATAAACTTTCATAAACTTTTGCAAACCTCATAAACTTAGCGCTGCGTAAGAAGACTTAAACCTCATAAACTTAATTAGCGCTGCGCGCTGCGTAAGTCATAAACTTGGCGCTGCGTAAGGGTCATACATACAGTAAGGGGGTTCAAGAATAGGAGTAACCAACCGGGTGCCTCATAGCGCTGCGTAAGTGCGTAAGGGCTGCGTAAGATAACCTCATAAACTTAGACGCGCTGCGTAAGTTCGTCTACCCTCATAAACTCCTCATAAACGCGCTGCGTAAGGCGCTGCGTAGCGCTGCGTAAGGCGGCGCTGCGTAAGCTGCGCGCTGCGCGCTGCGTAAGATATTTGTTAGGCGCTGCGTAAGCGTAAGGCGCGCGCTGCGTAAGGCGTAAGTGCGTAAGTCCAAATCATCCCGTGCATGGTTCTCCGCGCTGCGTGCGCTGCGTAAGTAGCGCTGCGTAAGAACAGGGGATTAGTTGGCCCTTCTCAGGTGACTATGCGCTGCGTAAGCTGCGTAAGTACCTGCCCCGCGCTGCGTAAGGCCAGCTGCGCTGCGTAGCGCTGCGTAAGCTGCGTAAGCTTCTGCGGCCCCCACACTGTTCTGCGGCCCCCTGCGGCCCCCCATCTGGTGTCAGACCGACGGCTCTTCTTTCTGCGGCCCCGGCCTTCTGCGGCCCCGAGGATTCAACCATTTTCTGCGGCCCCCCCCCTCCTCGGTTCTGCGGCCCTTCTGCGGCCCCCTTCTGCGGCCCCTTTCTGCGGCCCCCTGCGGCCCCCGCGCTTCTGCGGCCCCACTATTGTAATGCAGATTCTGCGGCTTCTGCGGCCCCCTGCCGGTCTACGTCACGTGTCAAGGGATGGACATTCTGCGGCCCCTTATTTCTTCTGCGGCCCCAGCGTTCGAATCCTCCCAAATCCCGGCTAGACACGTTCTGCGGCCCCAGTTCTGCGGCCCCGCTTCTGCTTCTGCGGCCCCTAAGGTTCTGCGGCCCCGTTCACACATATAGTTCTGCGGCCCCTCGCGCTTTCTGCGGCCCCCGGCCCCAGTAAAAGATTGTTCCAACTCATATCAACACGCAGAAATTCTGCGGCCCCCGGCGTTACCGCCTAGCCATAAAATAGGCGGCCAAGCCGAGTCTACCGCAGTTCTCGTTACTAATACCTGGACACTACTAATGACTGCAGAATAATACGGCCGCCCGAGACTGTAGAAACAGACTCATCTAGGGGCAATGTTGCGCCCCTAAGGGGAGGTCGACAGTTACCTACCCCACTAGACCTTAGGCGATAGGCGTGGGCGATA'
k, L, t = 12, 492, 17
out=[]
for i in range(len(genome)-L+1):
string = genome[i:i+L]
c={}
for j in range(len(string)-k+1):
kmer = string[j:j+k]
if kmer in c:
c[kmer] += 1
else:
c[kmer] = 1
for key, v in c.items():
if v >= t:
out.append(key)
print(' '.join(set(out)))