public
Created

length memoization

  • Download Gist
length.pl
Perl
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23
#!/usr/bin/env perl
 
use lib '/home/ben/l/git/bioperl-live';
 
use 5.006;
use strict;
use warnings;
use Bio::Seq;
use Bio::SeqFeature::SubSeq;
 
 
my $reference_sequence = Bio::Seq->new( -id => 'blah', -seq => 'ACGT'x1000000 );
 
#$reference->length( $reference_sequence->length );
 
for my $i (1..4000){
my $newseq = Bio::SeqFeature::SubSeq->new(
-template => $reference_sequence,
-start => 200000+$i,
-end => 200000 + 500+$i,
);
my $shiny_seq = $newseq->seq;
}

Please sign in to comment on this gist.

Something went wrong with that request. Please try again.