[Bioperl-l] SeqIO/scf.pm bug (with fix)

Heikki Lehvaslaiho heikki@ebi.ac.uk
29 Oct 2002 17:08:53 +0000


I've commmited the fix into the main trunk. 

Thanks, Tony.

	-Heikki

On Tue, 2002-10-29 at 15:57, Tony Cox wrote:
> 
> There is a bug in the SeqIO SCF module that means that SCF files with low
> dynamic range data that have their trace data values encoded as 8-bit unsigned
> values (rather than the normal 16-bit) are incorrectly parsed. This caused the
> unpack function to overrun the data and extract garbage.
> 
> I'm hoping that somebody can add the patch below to the main trunck SeqIO/scf.pm
> module. Here is the diff to fix it:
> 
> 
> 147,151c158
> <             my $byte = "n";
> <             if ($self->{sample_size} == 1){
> <                 $byte = "c";
> <             }
> <               @read = unpack "${byte}${length}",$buffer;
> ---
> >                       @read = unpack "n$length",$buffer;
> 
> 
> regards
> 
> Tony
> 
> 
> ******************************************************
> Tony Cox			Email:avc@sanger.ac.uk
> Sanger Institute		WWW:www.sanger.ac.uk
> Wellcome Trust Genome Campus	Webmaster
> Hinxton				Tel: +44 1223 834244
> Cambs. CB10 1SA			Fax: +44 1223 494919
> ******************************************************
> 
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l@bioperl.org
> http://bioperl.org/mailman/listinfo/bioperl-l
-- 
______ _/      _/_____________________________________________________
      _/      _/                      http://www.ebi.ac.uk/mutations/
     _/  _/  _/  Heikki Lehvaslaiho          heikki@ebi.ac.uk
    _/_/_/_/_/  EMBL Outstation, European Bioinformatics Institute
   _/  _/  _/  Wellcome Trust Genome Campus, Hinxton
  _/  _/  _/  Cambs. CB10 1SD, United Kingdom
     _/      Phone: +44 (0)1223 494 644   FAX: +44 (0)1223 494 468
___ _/_/_/_/_/________________________________________________________