[Bioperl-l] Tool to mutate DNA sequence
heikki at sanbi.ac.za
Thu Feb 9 14:54:30 UTC 2006
I should have made this very clear in my first reply:
You have to plan very carefully what rules you use when you mutate your
sequence because it will affect directly the resulting sequences. Of course,
all that depends on what you will be using the sequences for. If you are
going to draw evolutionary conclusions from those sequences, you must mutate
them in a way that simulates evolutionary principles.
My earlier pseudocode example, for example, should allow mutations in every
location. Mutations do occur multiple times in same places as sequences get
saturated by mutations. Also, you should decide the relative occurrence of
transversions versus transitions. Then there are indels; do you want to take
those into account?
Also, check the EMBOSS program 'msbar'.
You did not ask this, but... I remember that during the early days of Celera,
one of the tools that enabled them to estimate the feasibility of the whole
genome shotgun sequence assembly, was a very complete program to 'synthesize'
in-silico the whole complexity of the human genome. I have no idea of that
program is generally available now.
On Thursday 09 February 2006 06:46, Ryan Golhar wrote:
> Does anyone know of tool to mutate a DNA sequence by a specified amount?
> For instance, say I have a DNA sequence 1000 bases long, and I want to
> simulate mutations to make it 75% (or 80%, etc) similar to the original.
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
______ _/ _/_____________________________________________________
_/ _/ _/ Heikki Lehvaslaiho heikki at_sanbi _ac _za
_/_/_/_/_/ Associate Professor skype: heikki_lehvaslaiho
_/ _/ _/ SANBI, South African National Bioinformatics Institute
_/ _/ _/ University of Western Cape, South Africa
_/ Phone: +27 21 959 2096 FAX: +27 21 959 2512
More information about the Bioperl-l