[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