[Bioperl-l] bug
    Stefan Kirov 
    skirov at utk.edu
       
    Thu Oct  6 19:33:20 EDT 2005
    
    
  
Matt,
try this:
#!/usr/bin/perl
use warnings;
use strict;
use Bio::SeqIO;
#use Benchmark;
my $file = shift;
#my $t0 = new Benchmark;
my $s15;
# get sequence
{
my $fasta_in = Bio::SeqIO->new(  -file => "<$file", -format => 'fasta' );
my $seq_obj = $fasta_in->next_seq();
 $s15 = $seq_obj->seq();
}
use lib "bioperl-branch-1.4";
my $fasta_in = Bio::SeqIO->new(  -file => "<$file", -format => 'fasta' );
my $seq_obj = $fasta_in->next_seq();
my $s14 = $seq_obj->seq();
print "same\n" if ($s15 eq $s14);
exit();
    
    
More information about the Bioperl-l
mailing list