Bioinformatics/Sequence mutation: Difference between revisions

{{trans|Python}}

F nonrandom(n)

:seed = 1664525 * :seed + 1013904223

print(‘ #10 @#.’.format(kind, index))

print()

{{out}}

=={{header|Ada}}==

with Ada.Numerics.Discrete_Random;

with Ada.Text_Io;

Pretty_Print (Sequence);

{{out}}

<pre>Initial sequence:

=={{header|Arturo}}==

dna: map 1..200 => [sample bases]

print "------------------------------"

prettyPrint dna

{{out}}

=={{header|C}}==

Adenine ( A ) is always swapped for Thymine ( T ) and vice versa. Similarly with Cytosine ( C ) and Guanine ( G ).

#include<stdlib.h>

#include<stdio.h>

return 0;

}

Sample run :

<pre>

=={{header|C++}}==

#include <iomanip>

#include <iostream>

sequence_generator::print_sequence(std::cout, sequence);

return 0;

{{out}}

:: :genome <i><Genome Sequence></i>)

<b>All keys are optional. <i><Genome length></i> is discarded when :genome is set.</b>

(defun random_base ()

(random 4))

(t (delete_base genome)))))

(output_genome_info genome "MUTATED"))

{{out}}

<pre>

=={{header|Factor}}==

macros math math.statistics namespaces prettyprint quotations

random sequences sorting ;

[ mutate ] curry times nl "MUTATED " write show-dna ;

{{out}}

<pre>

=={{header|Go}}==

import (

fmt.Println()

prettyPrint(dna, 50)

{{out}}

</pre>

=={{header|Haskell}}==

import Data.List.Split (chunksOf)

import System.Random (Random, randomR, random, newStdGen, randoms, getStdRandom)

showSequence = mapM_ (uncurry (printf "%3d: %s\n")) . chunkedDNASequence

showBaseCounts = mapM_ (uncurry (printf "%s: %3d\n")) . baseCounts

{{out}}

<pre>Initial Sequence:

=={{header|J}}==

MUTS=: ;: 'del ins mut'

)

{{out}}

=={{header|Java}}==

import java.util.Random;

private static final int OP_COUNT = 3;

private static final char[] BASES = {'A', 'C', 'G', 'T'};

{{out}}

=={{header|JavaScript}}==

const numBases = 250

const numMutations = 30

console.log('\nMUTATED BASE COUNTS:')

printBases(mut);

{{out}}

=={{header|Julia}}==

randpos(seq) = rand(1:length(seq)) # 1

mutateat(pos, seq) = (s = seq[:]; s[pos] = rand(dnabases); s) # 2-1

testbioseq()

<pre>

500nt DNA sequence:

=={{header|Lua}}==

Using the <code>prettyprint()</code> function from [[Bioinformatics/base_count#Lua]] (not replicated here)

bases = {"A","C","T","G"}

function randbase() return bases[math.random(#bases)] end

print("MUTATED:")

{{out}}

<pre>ORIGINAL:

=={{header|Mathematica}} / {{header|Wolfram Language}}==

BioSequence is a fundamental data type in Mathematica:

seq = BioSequence["DNA", "ATAAACGTACGTTTTTAGGCT"];

randompos = RandomInteger[seq["SequenceLength"]];

ends = Rest[Accumulate[Prepend[StringLength /@ parts, 0]]];

StringRiffle[MapThread[ToString[#1] <> "-" <> ToString[#2] <> ": " <> #3 &, {begins, ends, parts}], "\n"]

{{out}}

<pre>1-50: TAGCAGGGGAATTGTCGACTCCCGGGTTTCAATTGCCAACCAAGCATATT

=={{header|Nim}}==

import strformat

import strutils

echo "\nMutated sequence"

echo "————————————————\n"

{{out}}

=={{header|Perl}}==

{{trans|Raku}}

use warnings;

use feature 'say';

say "Total bases: ". length $mutate;

say "$_: $cnt{$_}" for @bases;

{{out}}

<pre>Original DNA strand:

=={{header|Phix}}==

<span style="color: #004080;">string</span> <span style="color: #000000;">dna</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">repeat</span><span style="color: #0000FF;">(</span><span style="color: #008000;">' '</span><span style="color: #0000FF;">,</span><span style="color: #000000;">200</span><span style="color: #0000FF;">+</span><span style="color: #7060A8;">rand</span><span style="color: #0000FF;">(</span><span style="color: #000000;">300</span><span style="color: #0000FF;">))</span>

<span style="color: #008080;">for</span> <span style="color: #000000;">i</span><span style="color: #0000FF;">=</span><span style="color: #000000;">1</span> <span style="color: #008080;">to</span> <span style="color: #7060A8;">length</span><span style="color: #0000FF;">(</span><span style="color: #000000;">dna</span><span style="color: #0000FF;">)</span> <span style="color: #008080;">do</span> <span style="color: #000000;">dna</span><span style="color: #0000FF;">[</span><span style="color: #000000;">i</span><span style="color: #0000FF;">]</span> <span style="color: #0000FF;">=</span> <span style="color: #008000;">"ACGT"</span><span style="color: #0000FF;">[</span><span style="color: #7060A8;">rand</span><span style="color: #0000FF;">(</span><span style="color: #000000;">4</span><span style="color: #0000FF;">)]</span> <span style="color: #008080;">end</span> <span style="color: #008080;">for</span>

<span style="color: #000000;">mutate</span><span style="color: #0000FF;">()</span>

<span style="color: #000000;">show</span><span style="color: #0000FF;">()</span>

{{out}}

<pre>

=={{header|PureBasic}}==

#SEQLEN=200

#PROTOCOL=#True

PrintN("After 10 mutations:")

pprint()

{{out}}

<pre>Initial sequence:

Similarly parameter choice is chosen from to give the base for substitution or insertion - the more any base appears, the more likely it is to be chosen in any insertion/substitution.

from collections import Counter

print(f" {kind:>10} @{index}")

print()

{{out}}

<code>prettyprint</code> and <code>tallybases</code> are defined at [[Bioinformatics/base count#Quackery]].

[ $ "" swap times

cr cr say "Mutating..." cr

10 times mutate

{{out}}

=={{header|Racket}}==

(define current-S-weight (make-parameter 1))

(define s+d (parameterize ((current-D-weight 5)) (for/fold ((s initial-sequence)) ((_ 10)) (mutate s))))

(newline)

{{out}}

# The DNA strand

sub diff ($orig, $repl) {

($orig.comb Z $repl.comb).map( -> ($o, $r) { $o eq $r ?? $o !! $r.lc }).join

{{out}}

<pre>ORIGINAL DNA STRAND:

=={{header|Ring}}==

row = 0

dnaList = []

total = dnaBase["A"] + dnaBase["T"] + dnaBase["C"] + dnaBase["G"]

see "Total: " + total+ nl + nl

{{out}}

<pre>

=={{header|Ruby}}==

attr_accessor :seq

test.delete

puts test

{{out}}

<pre> 0 TAAGGTGAGGAGTGTGATGGAGTTCGGTGGCTAGCCACAAATACAACACA

=={{header|Swift}}==

enum Action: CaseIterable {

}

{{out}}

=={{header|Vlang}}==

{{trans|Go}}

import rand.seed

println('')

pretty_print(dna, 50)

{{out}}

{{libheader|Wren-sort}}

{{libheader|Wren-fmt}}

import "/fmt" for Fmt

import "/sort" for Sort

for (i in 0...muts) dna = mutate.call(dna, w)

System.print()

{{out}}

=={{header|Yabasic}}==

{{trans|Phix}}

// by Galileo, 07/2022

show()

mutate()

{{out}}

<pre>1: TCCATCGTGG GATCGCTCTA GCGGTATGCT ATCATTCCTA TAGCAATTCT

=={{header|zkl}}==

dna:=(190).pump(Data().howza(3),(0).random.fp(0,4),bases.get); // bucket of bytes

[0..*,50].zipWith(fcn(n,bases){ println("%6d: %s".fmt(n,bases.concat())) },

dna.walker().walk.fp(50)).pump(Void); // .pump forces the iterator

{{out}}

<pre>

=={{header|Rust}}==

use rand::prelude::*;

use std::collections::HashMap;

println!("\nMutated sequence:\n{}", seq);

}

{{out}}

<pre>