[Bioperl-l] Bio::PrimarySeq and Selenocysteine

[Heikki Lehvaslaiho]

Lincoln,

How did you manage to get that error?
This works for me:

#----------------------------------------------
use Bio::PrimarySeq;
$seq = Bio::PrimarySeq->new
    (-seq              => 'UABCDEFGHIJKLMN',
     -display_id       => 'new-id',
     -alphabet         => 'protein',
     -verbose          => 2);
$seq->seq('ABCDEFGHIJKLMNU');
#----------------------------------------------

	-Heikki

On Mon, 2003-07-21 at 20:51, Lincoln Stein wrote:
> Hi,
> 
> I guess I'm missing something.  I'm trying to load a protein that contains a 
> selenocysteine (code "U") and am getting an exception about the protein 
> alphabet not being right.  Is there an alternative alphabet that I should be 
> using?
> 
> Lincoln
-- 
______ _/      _/_____________________________________________________
      _/      _/                      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
___ _/_/_/_/_/________________________________________________________