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