Bug in Bio::EnsEMBL::DBSQL::SequenceAdaptor
Jonathan Barber
jon at compbio.dundee.ac.uk
Fri Jul 4 12:12:12 UTC 2003
The method _reverse_comp is buggy as it assigns to $_, which messes up
$_ (which is global is scope).
Diff is as follows:
Index: SequenceAdaptor.pm
===================================================================
RCS file: /cvsroot/CVSmaster/ensembl/modules/Bio/EnsEMBL/DBSQL/SequenceAdaptor.pm,v
retrieving revision 1.17
diff -r1.17 SequenceAdaptor.pm
299,304c299,302
< my $self = shift;
< my $seq = shift;
<
< $_ = reverse( $seq );
< tr/CGTAcgta/GCATgcat/;
< return $_;
---
> my ($self, $seq) = @_;
> $seq = reverse( $seq );
> $seq =~ tr/CGTAcgta/GCATgcat/;
> return $seq;
--
Jon
More information about the EMBOSS
mailing list