[Bioperl-l] Shuffling sequences

martin s9904982 at sms.ed.ac.uk
Wed May 26 03:59:47 EDT 2004


shuffleseq from emboss would allow you to specify the number of
replicates you want on the command line:

shuffleseq in.fasta out.fasta -shuffle = 1000

out.fasta then contains the 1000 shuffled replicates.

Martin

On Wed, 2004-05-26 at 05:30, Hilmar Lapp wrote:
> Isn't there a tool in EMBOSS that exactly that and does it pretty fast?
> 
> On Tuesday, May 25, 2004, at 12:16  PM, Jason Stajich wrote:
> 
> > my $i = 1;
> > for my $s ( @randomseqs ) {
> >   my $newseq = new Bio::PrimarySeq(-display_id => "rand.$i",
> > 	                           -seq        => $s);
> >   $out->write_seq($s);
> 
> $i++;
> 
> 	-hilmar
> 
> > }
> >
-- 
Martin Jones
The Nematode Genomics Lab
Institute of Cell, Animal and Population Biology
University of Edinburgh
King's Buildings
West Mains Road
Edinburgh, EH9 3JT




More information about the Bioperl-l mailing list