[Bioperl-l] Announce: Bio::Seq::Quality

Heikki Lehvaslaiho heikki at ebi.ac.uk
Mon Jul 11 12:02:52 EDT 2005


Bio::Seq::Quality is a new module that allows you to store per-residue quality 
and trace index values using Bio::Seq::MetaI interface. It replaces 
Bio::Seq::SeqWithQuality which is now deprecated.

Solutions to persistence should focus on storing Bio::Seq::Meta and 
Bio::Seq::Meta::Array objects. It should be easy to stringify most real world 
meta values. Then the persistence could be implemented by storing the 
sequence object and N number of meta strings.

All the functional code is in Bio::Seq::Meta::Array, Bio::Seq::Quality merely 
adds a convenient interface.

The POD contains a discussion of differences from Bio::Seq::SeqWithQuality. 
If the following, or anything else,  is a problem let me know as soon as 
possible:

  The greatest difference to Bio::Seq::SeqWithQuality is that in this
  implementation quality for all sequence residues are automatically
  assigned a value of '0' (zero) unless you set it to something
  else. Length of the quality array always equals the length of the
  sequence. Therefore, length() never returns "DIFFERENT".


Enjoy,
 -Heikki


-- 
______ _/      _/_____________________________________________________
      _/      _/                      http://www.ebi.ac.uk/mutations/
     _/  _/  _/  Heikki Lehvaslaiho    heikki at_ebi _ac _uk
    _/_/_/_/_/  EMBL Outstation, European Bioinformatics Institute
   _/  _/  _/  Wellcome Trust Genome Campus, Hinxton
  _/  _/  _/  Cambridge, CB10 1SD, United Kingdom
     _/      Phone: +44 (0)1223 494 644   FAX: +44 (0)1223 494 468
___ _/_/_/_/_/________________________________________________________


More information about the Bioperl-l mailing list