{"action":"opened","issue":{"url":"https://api.github.com/repos/BioJulia/Bio.jl/issues/348","repository_url":"https://api.github.com/repos/BioJulia/Bio.jl","labels_url":"https://api.github.com/repos/BioJulia/Bio.jl/issues/348/labels{/name}","comments_url":"https://api.github.com/repos/BioJulia/Bio.jl/issues/348/comments","events_url":"https://api.github.com/repos/BioJulia/Bio.jl/issues/348/events","html_url":"https://github.com/BioJulia/Bio.jl/issues/348","id":182451175,"number":348,"title":"inconsistent dna\"\" - DNASequence() behavior","user":{"login":"cjprybol","id":3497642,"avatar_url":"https://avatars.githubusercontent.com/u/3497642?v=3","gravatar_id":"","url":"https://api.github.com/users/cjprybol","html_url":"https://github.com/cjprybol","followers_url":"https://api.github.com/users/cjprybol/followers","following_url":"https://api.github.com/users/cjprybol/following{/other_user}","gists_url":"https://api.github.com/users/cjprybol/gists{/gist_id}","starred_url":"https://api.github.com/users/cjprybol/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/cjprybol/subscriptions","organizations_url":"https://api.github.com/users/cjprybol/orgs","repos_url":"https://api.github.com/users/cjprybol/repos","events_url":"https://api.github.com/users/cjprybol/events{/privacy}","received_events_url":"https://api.github.com/users/cjprybol/received_events","type":"User","site_admin":false},"labels":[],"state":"open","locked":false,"assignee":null,"assignees":[],"milestone":null,"comments":0,"created_at":"2016-10-12T07:08:14Z","updated_at":"2016-10-12T07:08:14Z","closed_at":null,"body":"This is an example of simulating RNA transcripts (in DNA format) from a fake transcriptome and then randomly making sequencing errors. There are two for loops, the first utilizes `dna\"\"` to create a new, empty dna sequence and in the second for loop I use the `DNASequence()` function.\r\n\r\ntracing through the steps with `@show`,  `dna\"\"` doesn't seem to release the object in memory from the prior loop. This causes the `new_seq` variable to be continuously extended during the for loop iterations. `DNASequence()` has the behavior I expected, and assigns an empty sequence to the variable `new_seq` during each for loop.\r\n\r\nThe key lines are\r\n```\r\n@show new_seq = dna\"\"\r\n```\r\nvs.\r\n```\r\n@show new_seq = DNASequence()\r\n```\r\n\r\nIs this intended behavior, or potentially a bug in the parsing of the `dna\"\"` constructor?\r\n\r\n**code**\r\n```\r\nusing Bio.Seq\r\nusing StatsBase\r\n\r\nsrand(1);\r\ntranscriptome = [randdnaseq(rand(8:12)) for i=1:3];\r\nnumber_of_reads = 5;\r\nanswers = sample(transcriptome, number_of_reads);\r\n\r\nerror = .10;\r\nweights = WeightVec([1-error, error * 1/3, error * 1/3, error * 1/3]);\r\n\r\nsequenced_reads_dna = []\r\nfor seq in answers\r\n\t@show new_seq = dna\"\"\r\n\tfor nuc in seq\r\n\t\t@show choice = sample(1:4, weights)\r\n\t\tif choice == 1\r\n\t\t\t@show push!(new_seq, nuc)\r\n\t\telseif choice == 2\r\n\t\t\t# deletion, do nothing\r\n\t\t\t@show new_seq\r\n\t\telseif choice == 3\r\n\t\t\t# mutation\r\n\t\t\t@show new_seq *= randdnaseq(1)\r\n\t\telseif choice == 4\r\n\t\t\t# insertion\r\n\t\t\t@show new_seq *= push!(randdnaseq(1), nuc)\r\n\t\tend\r\n\tend\r\n\tpush!(sequenced_reads_dna, new_seq)\r\nend\r\n\r\nsequenced_reads_function = []\r\nfor seq in answers\r\n\t@show new_seq = DNASequence()\r\n\tfor nuc in seq\r\n\t\t@show choice = sample(1:4, weights)\r\n\t\tif choice == 1\r\n\t\t\t@show push!(new_seq, nuc)\r\n\t\telseif choice == 2\r\n\t\t\t# deletion, do nothing\r\n\t\t\t@show new_seq\r\n\t\telseif choice == 3\r\n\t\t\t# mutation\r\n\t\t\t@show new_seq *= randdnaseq(1)\r\n\t\telseif choice == 4\r\n\t\t\t# insertion\r\n\t\t\t@show new_seq *= push!(randdnaseq(1), nuc)\r\n\t\tend\r\n\tend\r\n\tpush!(sequenced_reads_function, new_seq)\r\nend\r\n\r\nanswers\r\nsequenced_reads_dna\r\nsequenced_reads_function\r\n```\r\n\r\n**REPL**\r\n```\r\n               _\r\n   _       _ _(_)_     |  A fresh approach to technical computing\r\n  (_)     | (_) (_)    |  Documentation: http://docs.julialang.org\r\n   _ _   _| |_  __ _   |  Type \"?help\" for help.\r\n  | | | | | | |/ _` |  |\r\n  | | |_| | | | (_| |  |  Version 0.5.1-pre+2 (2016-09-20 03:34 UTC)\r\n _/ |\\__'_|_|_|\\__'_|  |  Commit f0d40ec* (22 days old release-0.5)\r\n|__/                   |  x86_64-apple-darwin16.1.0\r\n\r\njulia> using Bio.Seq\r\n\r\njulia> using StatsBase\r\n\r\njulia> srand(1);\r\n\r\njulia> transcriptome = [randdnaseq(rand(8:12)) for i=1:3];\r\n\r\njulia> number_of_reads = 5;\r\n\r\njulia> answers = sample(transcriptome, number_of_reads);\r\n\r\njulia> error = .10;\r\n\r\njulia> weights = WeightVec([1-error, error * 1/3, error * 1/3, error * 1/3]);\r\n\r\njulia> sequenced_reads_dna = []\r\n0-element Array{Any,1}\r\n\r\njulia> for seq in answers\r\n               @show new_seq = dna\"\"\r\n               for nuc in seq\r\n                       @show choice = sample(1:4, weights)\r\n                       if choice == 1\r\n                               @show push!(new_seq, nuc)\r\n                       elseif choice == 2\r\n                               # deletion, do nothing\r\n                               @show new_seq\r\n                       elseif choice == 3\r\n                               # mutation\r\n                               @show new_seq *= randdnaseq(1)\r\n                       elseif choice == 4\r\n                               # insertion\r\n                               @show new_seq *= push!(randdnaseq(1), nuc)\r\n                       end\r\n               end\r\n               push!(sequenced_reads_dna, new_seq)\r\n       end\r\nnew_seq = @dna_str(\"\") = 0nt DNA Sequence:\r\n\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 1nt DNA Sequence:\r\nA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 2nt DNA Sequence:\r\nAT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 3nt DNA Sequence:\r\nATA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 4nt DNA Sequence:\r\nATAG\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 5nt DNA Sequence:\r\nATAGC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 6nt DNA Sequence:\r\nATAGCA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 7nt DNA Sequence:\r\nATAGCAA\r\nchoice = sample(1:4,weights) = 3\r\nnew_seq *= randdnaseq(1) = 8nt DNA Sequence:\r\nATAGCAAA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 9nt DNA Sequence:\r\nATAGCAAAA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 10nt DNA Sequence:\r\nATAGCAAAAG\r\nnew_seq = @dna_str(\"\") = 7nt DNA Sequence:\r\nATAGCAA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 8nt DNA Sequence:\r\nATAGCAAA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 9nt DNA Sequence:\r\nATAGCAAAT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 10nt DNA Sequence:\r\nATAGCAAATA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 11nt DNA Sequence:\r\nATAGCAAATAG\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 12nt DNA Sequence:\r\nATAGCAAATAGC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 13nt DNA Sequence:\r\nATAGCAAATAGCA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 14nt DNA Sequence:\r\nATAGCAAATAGCAA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 15nt DNA Sequence:\r\nATAGCAAATAGCAAT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 16nt DNA Sequence:\r\nATAGCAAATAGCAATA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 17nt DNA Sequence:\r\nATAGCAAATAGCAATAG\r\nnew_seq = @dna_str(\"\") = 17nt DNA Sequence:\r\nATAGCAAATAGCAATAG\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 18nt DNA Sequence:\r\nATAGCAAATAGCAATAGT\r\nchoice = sample(1:4,weights) = 3\r\nnew_seq *= randdnaseq(1) = 19nt DNA Sequence:\r\nATAGCAAATAGCAATAGTC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 20nt DNA Sequence:\r\nATAGCAAATAGCAATAGTCA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 21nt DNA Sequence:\r\nATAGCAAATAGCAATAGTCAC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 22nt DNA Sequence:\r\nATAGCAAATAGCAATAGTCACA\r\nchoice = sample(1:4,weights) = 4\r\nnew_seq *= push!(randdnaseq(1),nuc) = 24nt DNA Sequence:\r\nATAGCAAATAGCAATAGTCACAGC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 25nt DNA Sequence:\r\nATAGCAAATAGCAATAGTCACAGCT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 26nt DNA Sequence:\r\nATAGCAAATAGCAATAGTCACAGCTA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 27nt DNA Sequence:\r\nATAGCAAATAGCAATAGTCACAGCTAG\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 28nt DNA Sequence:\r\nATAGCAAATAGCAATAGTCACAGCTAGT\r\nchoice = sample(1:4,weights) = 4\r\nnew_seq *= push!(randdnaseq(1),nuc) = 30nt DNA Sequence:\r\nATAGCAAATAGCAATAGTCACAGCTAGTGC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 31nt DNA Sequence:\r\nATAGCAAATAGCAATAGTCACAGCTAGTGCT\r\nnew_seq = @dna_str(\"\") = 18nt DNA Sequence:\r\nATAGCAAATAGCAATAGT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 19nt DNA Sequence:\r\nATAGCAAATAGCAATAGTT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 20nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTC\r\nchoice = sample(1:4,weights) = 3\r\nnew_seq *= randdnaseq(1) = 21nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 22nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCCC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 23nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCCCA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 24nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCCCAC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 25nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCCCACT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 26nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCCCACTA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 27nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCCCACTAG\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 28nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCCCACTAGT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 29nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCCCACTAGTC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 30nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCCCACTAGTCT\r\nnew_seq = @dna_str(\"\") = 20nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 21nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 22nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCTC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 23nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCTCA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 24nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCTCAC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 25nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCTCACA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 26nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCTCACAC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 27nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCTCACACT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 28nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCTCACACTA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 29nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCTCACACTAG\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 30nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCTCACACTAGT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 31nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCTCACACTAGTC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 32nt DNA Sequence:\r\nATAGCAAATAGCAATAGTTCTCACACTAGTCT\r\n\r\njulia> sequenced_reads_function = []\r\n0-element Array{Any,1}\r\n\r\njulia> for seq in answers\r\n               @show new_seq = DNASequence()\r\n               for nuc in seq\r\n                       @show choice = sample(1:4, weights)\r\n                       if choice == 1\r\n                               @show push!(new_seq, nuc)\r\n                       elseif choice == 2\r\n                               # deletion, do nothing\r\n                               @show new_seq\r\n                       elseif choice == 3\r\n                               # mutation\r\n                               @show new_seq *= randdnaseq(1)\r\n                       elseif choice == 4\r\n                               # insertion\r\n                               @show new_seq *= push!(randdnaseq(1), nuc)\r\n                       end\r\n               end\r\n               push!(sequenced_reads_function, new_seq)\r\n       end\r\nnew_seq = DNASequence() = 0nt DNA Sequence:\r\n\r\nchoice = sample(1:4,weights) = 2\r\nnew_seq = 0nt DNA Sequence:\r\n\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 1nt DNA Sequence:\r\nT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 2nt DNA Sequence:\r\nTA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 3nt DNA Sequence:\r\nTAG\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 4nt DNA Sequence:\r\nTAGC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 5nt DNA Sequence:\r\nTAGCA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 6nt DNA Sequence:\r\nTAGCAA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 7nt DNA Sequence:\r\nTAGCAAT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 8nt DNA Sequence:\r\nTAGCAATA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 9nt DNA Sequence:\r\nTAGCAATAG\r\nnew_seq = DNASequence() = 0nt DNA Sequence:\r\n\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 1nt DNA Sequence:\r\nA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 2nt DNA Sequence:\r\nAT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 3nt DNA Sequence:\r\nATA\r\nchoice = sample(1:4,weights) = 2\r\nnew_seq = 3nt DNA Sequence:\r\nATA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 4nt DNA Sequence:\r\nATAC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 5nt DNA Sequence:\r\nATACA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 6nt DNA Sequence:\r\nATACAA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 7nt DNA Sequence:\r\nATACAAT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 8nt DNA Sequence:\r\nATACAATA\r\nchoice = sample(1:4,weights) = 2\r\nnew_seq = 8nt DNA Sequence:\r\nATACAATA\r\nnew_seq = DNASequence() = 0nt DNA Sequence:\r\n\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 1nt DNA Sequence:\r\nT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 2nt DNA Sequence:\r\nTC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 3nt DNA Sequence:\r\nTCA\r\nchoice = sample(1:4,weights) = 4\r\nnew_seq *= push!(randdnaseq(1),nuc) = 5nt DNA Sequence:\r\nTCATC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 6nt DNA Sequence:\r\nTCATCA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 7nt DNA Sequence:\r\nTCATCAC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 8nt DNA Sequence:\r\nTCATCACT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 9nt DNA Sequence:\r\nTCATCACTA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 10nt DNA Sequence:\r\nTCATCACTAG\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 11nt DNA Sequence:\r\nTCATCACTAGT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 12nt DNA Sequence:\r\nTCATCACTAGTC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 13nt DNA Sequence:\r\nTCATCACTAGTCT\r\nnew_seq = DNASequence() = 0nt DNA Sequence:\r\n\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 1nt DNA Sequence:\r\nT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 2nt DNA Sequence:\r\nTC\r\nchoice = sample(1:4,weights) = 4\r\nnew_seq *= push!(randdnaseq(1),nuc) = 4nt DNA Sequence:\r\nTCGA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 5nt DNA Sequence:\r\nTCGAC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 6nt DNA Sequence:\r\nTCGACA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 7nt DNA Sequence:\r\nTCGACAC\r\nchoice = sample(1:4,weights) = 2\r\nnew_seq = 7nt DNA Sequence:\r\nTCGACAC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 8nt DNA Sequence:\r\nTCGACACA\r\nchoice = sample(1:4,weights) = 4\r\nnew_seq *= push!(randdnaseq(1),nuc) = 10nt DNA Sequence:\r\nTCGACACAGG\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 11nt DNA Sequence:\r\nTCGACACAGGT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 12nt DNA Sequence:\r\nTCGACACAGGTC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 13nt DNA Sequence:\r\nTCGACACAGGTCT\r\nnew_seq = DNASequence() = 0nt DNA Sequence:\r\n\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 1nt DNA Sequence:\r\nT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 2nt DNA Sequence:\r\nTC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 3nt DNA Sequence:\r\nTCA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 4nt DNA Sequence:\r\nTCAC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 5nt DNA Sequence:\r\nTCACA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 6nt DNA Sequence:\r\nTCACAC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 7nt DNA Sequence:\r\nTCACACT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 8nt DNA Sequence:\r\nTCACACTA\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 9nt DNA Sequence:\r\nTCACACTAG\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 10nt DNA Sequence:\r\nTCACACTAGT\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 11nt DNA Sequence:\r\nTCACACTAGTC\r\nchoice = sample(1:4,weights) = 1\r\npush!(new_seq,nuc) = 12nt DNA Sequence:\r\nTCACACTAGTCT\r\n\r\njulia> answers\r\n5-element Array{Bio.Seq.BioSequence{Bio.Seq.DNAAlphabet{4}},1}:\r\n ATAGCAATAG\r\n ATAGCAATAG\r\n TCACACTAGTCT\r\n TCACACTAGTCT\r\n TCACACTAGTCT\r\n\r\njulia> sequenced_reads_dna\r\n5-element Array{Any,1}:\r\n ATAGCAAAAG\r\n ATAGCAAATAGCAATAGTTCTCACACTAGTCT\r\n ATAGCAAATAGCAATAGTCACAGCTAGTGCT\r\n ATAGCAAATAGCAATAGTTCCCACTAGTCT\r\n ATAGCAAATAGCAATAGTTCTCACACTAGTCT\r\n\r\njulia> sequenced_reads_function\r\n5-element Array{Any,1}:\r\n TAGCAATAG\r\n ATACAATA\r\n TCATCACTAGTCT\r\n TCGACACAGGTCT\r\n TCACACTAGTCT\r\n```"}}
{"actor":{"display_login":"cjprybol"}}