[Bioperl-l] Tool to mutate DNA sequence

Ryan Golhar golharam at umdnj.edu
Thu Feb 9 21:19:46 UTC 2006

Thanks all.  The responses I got were definitely more than helpful.  FYI
- I did initially look at msbar.  I glanced over the "Number of times to
perform mutation operations", which is what I was looking for.  

I'm looking to statistically test some simply scoring matrices.  I think
msbar will do.


-----Original Message-----
From: bioperl-l-bounces at lists.open-bio.org
[mailto:bioperl-l-bounces at lists.open-bio.org] On Behalf Of Heikki
Sent: Thursday, February 09, 2006 9:55 AM
To: bioperl-l at lists.open-bio.org; golharam at umdnj.edu
Cc: 'The general forum at Bioinformatics.Org'; 'bioperl-l';
emboss at emboss.open-bio.org
Subject: Re: [Bioperl-l] Tool to mutate DNA sequence


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
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
them in a way that simulates evolutionary principles.

My earlier pseudocode example, for example, should allow mutations in
location. Mutations do occur multiple times in same places as sequences
saturated by mutations. Also, you should decide the relative occurrence
transversions versus transitions. Then there are indels; do you want to
those into account?

Also, check the EMBOSS program 'msbar'.

You did not ask this, but... I remember that during the early days of
one of the tools that enabled them to estimate the feasibility of the
genome shotgun sequence assembly, was a very complete program to
in-silico the whole complexity of the human genome. I have no idea of
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.
> Ryan
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org 
> http://lists.open-bio.org/mailman/listinfo/bioperl-l

______ _/      _/_____________________________________________________
      _/      _/
     _/  _/  _/  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
___ _/_/_/_/_/________________________________________________________
Bioperl-l mailing list
Bioperl-l at lists.open-bio.org

More information about the Bioperl-l mailing list