[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