Rosalind Challenge 9: Finding a Motif in DNA # http://rosalind.info/problems/subs/ # #------------------------------------------------# print("Rosalind Challenge 9: Finding a Motif in DNA\n") s = Seq("ACAAGCTGAGTTGCTGAGTCGCTGAGTCCGCTGAGTATCGGCTGAGTACTTGCTGAGTCAGCTGAGTGCTGAGTCGCTGAGTTAATATGGCTGAGTAGCTGAGTTACTAGCTGAGTCGCTGAGTGAGCCGCTGAGTCACCGCATCGTGCATACGCTGAGTACGCGCTGAGTATGCTGAGTGCTGAGTCGGGAGGAAATAAAGCTGAGTTGCTGAGTGCTGAGTTTCTGCTGAGTGCTGAGTGCTGAGTGGCTGAGTAAGAAGCTGAGTAGTGCTGAGTTTGTCGCTGAGTGCTGAGTTCGGCTGAGTGCTGAGTGCTGAGTGCTGAGTTACAGCTGAGTAGAAAGGGGCTGAGTGCTGAGTAGCTGAGTCGTGCTGAGTCGCTGAGTGCTGAGTAGCTGAGTAGGCTGAGTAGCTGAGTGCTGAGTGCTGAGTAAGCTGAGTGTATCGGTGGCTGAGTAAGCTGAGTGCTGAGTTGCGCCGCTGAGTCGCTGAGTTGTTGCTGAGTCGCTGAGTGCTGAGTTGGCTGAGTTGAGCTGAGTCTTGCTGAGTTTCGCTGAGTGCTGAGTGCTGAGTAGAGGGGCTGAGTACTGCTGAGTGCTGAGTTGTATGCTGAGTGCTGAGTGGGCCGCTGAGTGCTGAGTGAGGGGCTGAGTGCTGAGTTGCTGAGTAGCTGAGTGCTGAGTGCTGAGTGCTGAGTGCTGAGTGCTGAGTGCTGGCTGAGTAGGCTGAGTTGCGCTGAGTGCTGAGTTGGTTTGCTGAGTAGCTGAGTGGCTGAGTATGCTGAGTTCGCTGAGTGATCCTCGGCTGAGTGCTGAGTCCGCTGAGTGGCTGAGTGCTGAGTATACTGCTGAGTGCTGAGTTGCTGAGTTTAGCTGAGTGCTGAGTGGCTGAGTAGCCTGCTGAGTGGACGGTGCTGAGTCCCAGCTGAGTAGCTGAGTGCTGAGTAGCTGAGTGGCTGAGTATTCAAGCTGAGTGCTGAGTGCTGAGTGCTGAGT") t = Seq("GCTGAGTGC") locations = [] count = 0 loc = 0 while (count < len(s)): loc = s.find(t, start=count) if (loc > 0): count = loc + 1 locations.append(loc + 1) # adding plus 1 since the result needs to be in count from 1 format else: count = count + 1 print(' '.join(map(str,locations))) print("\n---\n")
Preview:
downloadDownload PNG
downloadDownload JPEG
downloadDownload SVG
Tip: You can change the style, width & colours of the snippet with the inspect tool before clicking Download!
Click to optimize width for Twitter