[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