[Bioperl-l] Bio::Tools::SeqStat

Heikki Lehvaslaiho heikki at ebi.ac.uk
Fri Feb 28 09:50:34 EST 2003


Peter,

I noticed that SeqStat module choked on some Bio::RichSeq entries coming
from EMBL/Genbank. The reason is that SeqStat imposes strict alphabets
which for RNA sequences is (A,U,G,C) but the incoming sequence contains
Ts.

I fixed that in the HEAD so that  _is_alphabet_strict(),
get_mol_wt() and count_monomers() internally convert s/T/U/g.
However, I am not 100% sure that this is what users expect. Your call.

I ran into this while cleaning scripts. Try:

bala ~/src/bioperl-live/scripts/seqstats> perl gccalc.PLS -f embl
../../t/data/ar.embl


	-Heikki


-- 
______ _/      _/_____________________________________________________
      _/      _/                      http://www.ebi.ac.uk/mutations/
     _/  _/  _/  Heikki Lehvaslaiho          heikki at 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
___ _/_/_/_/_/________________________________________________________



More information about the Bioperl-l mailing list