By downloading this code you agree to the
Source Code Use License (PDF). |
{ version = 1.03; (* of mutseq.p 2004 Jun 29}
(* begin module describe.mutseq *)
(*
name
mutseq: mutate a sequence
synopsis
mutseq(sequseed: in, mutseqp: in, output: out)
files
sequseed: DNA sequence and seed. The first line contains
a sequence of only the characters [acgt]. The second line
contains a random number seed between 0 and 1.
The sequence may be of any length less than constant maxsequ.
mutseqp: parameters to control the program. The file must contain the
following parameters, one per line:
parameterversion: The version number of the program. This allows the
user to be warned if an old parameter file is used.
If the first line contains 'v' then the program will
show the steps it used, verbose mode.
output: messages to the user
description
The program reads in a DNA sequence and a random number seed (a
real number between 0 and 1) from the sequseed file. It makes a
new random number seed and then uses that to select a base to
modify and a base to modify it to. The result is written back out
to the sequseed file.
examples
An example mutseqp parameter file:
1.00 version of mutseq that this parameter file is designed for.
verbose mode
An example sequseed file:
agttatactc
0.82911751563182178870
documentation
see also
Another random sequence generator program: markov.p
author
Thomas Dana Schneider
bugs
technical notes
maxsequ is a constant in the program.
*)
(* end module describe.mutseq *)
{This manual page was created by makman 1.45}
{created by htmlink 1.62}