Create a gist now

Instantly share code, notes, and snippets.

What would you like to do?
length memoization
#!/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;
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment