[Bioperl-l] GCG MSF format alignments
Ewan Birney
birney@ebi.ac.uk
Sat, 10 Mar 2001 13:49:44 +0000 (GMT)
On Fri, 9 Mar 2001, Hilmar Lapp wrote:
> Peter Schattner wrote:
> >
> > > Regarding the ~ characters ... my primitive understanding suggests this
> > > throws an exception in Bio::Primaryseq::seq with a message that '[some
> > > sequence ~~~] does not look healthy'.
> >
> > I've had this problem too in other contexts. Bio::Primaryseq is
> > currently rather strict about what it allows in asequence. I would
> > prefer to see it issue a warning rather than die when it comes across a
> > bad character (what do you think about that Ewan? Hilmar? Jason?)
>
> I actually like the exception being thrown there, even though I bumped
> into this several times, too. However, it always indicated a bug.
>
> I'd actually prefer this being solved by creating a specialized class
> that specifically allows for additional characters as they are used in
> e.g. alignments. This should go under Bio::Seq::<AlignedSeq>.pm
> (replace <...> with whatever you like better). Ideally Bio::PrimarySeq
> (or maybe better Bio::PrimarySeqI?) has a method boolean
> validate_seq(), and depending on its return value seq() throws an
> exception. So, in your deriving implementation you'd only have to
> override validate_seq(). (As a matter of fact, there is no
> validate_seq() right now :|
>
> Does this make sense to anyone?
Yup. The refinement here is that Alignments use LocatableSeq so that's the
place to put it (perhaps LocatableSeq should migrate to AlignedSeq...)
>
> Hilmar
> --
> -------------------------------------------------------------
> Hilmar Lapp email: lapp@gnf.org
> GNF, San Diego, Ca. 92121 phone: +1-858-812-1757
> -------------------------------------------------------------
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l@bioperl.org
> http://bioperl.org/mailman/listinfo/bioperl-l
>
-----------------------------------------------------------------
Ewan Birney. Mobile: +44 (0)7970 151230, Work: +44 1223 494420
<birney@ebi.ac.uk>.
-----------------------------------------------------------------