Bioinformatics/Subsequence: Difference between revisions
julia example
(julia example) |
|||
Line 72:
312..316
</pre>
=={{header|Julia}}==
<lang julia>DNArand(n, bases=['A', 'T', 'C', 'G']) = String(rand(bases, n))
DNAsearch(needle, haystack, lap=true) = findall(needle, haystack, overlap=lap)
const rand_string = DNArand(200)
const subseq = DNArand(4)
println("Search sequence:\n$rand_string\nfor substring $subseq. Found at positions: ")
foreach(p -> print(rpad(p[2], 8), p[1] % 10 == 0 ? "\n" : ""), DNAsearch(subseq, rand_string))
</lang>{{out}}
<pre>
Search sequence:
AAGAGTACGTCGCCAGGGAAGCTTCGGAACGTGCCCGGTGCCAAGCGTCACTACGTGGCAGAGTATATTCATGCTGCAAGGAAATTATTAATCGGGTACTGTGCGCAGCTTGTGCCGCGATATTTTGTACTCTCTCGGAATAGGCAACGCTGGATCATGCGTAGACTGTTGATCGGACGAGTTTGTTACGGATAATAAGG
for substring TGTT. Found at positions:
168 185
</pre>
=={{header|Phix}}==
| |||