Bioinformatics/Sequence mutation: Difference between revisions

Content added Content deleted
Line 475: Line 475:
=={{header|Julia}}==
=={{header|Julia}}==
<lang julia>dnabases = ['A', 'C', 'G', 'T']
<lang julia>dnabases = ['A', 'C', 'G', 'T']
randpos(seq) = rand(1:length(seq) # 1
randpos(seq) = rand(1:length(seq)) # 1
mutateat(pos, seq) = (s = seq[:]; s[pos] = rand(dnabases); s) # 2-1
mutateat(pos, seq) = (s = seq[:]; s[pos] = rand(dnabases); s) # 2-1
deleteat(pos, seq) = [seq[1:pos-1]; seq[pos+1:end]] # 2-2
deleteat(pos, seq) = [seq[1:pos-1]; seq[pos+1:end]] # 2-2