[Bioperl-l] Bio::SeqFeature::Gene::Exon throws exception when encountering split location (Bio::Location::Split)
cjfields at uiuc.edu
Tue Jun 12 12:42:29 EDT 2007
On Jun 12, 2007, at 11:10 AM, Mark Johnson wrote:
> On 6/12/07, Torsten Seemann
> <torsten.seemann at infotech.monash.edu.au> wrote:
>> Can you use the ->spliced_seq() method to do this?
>> --Torsten Seemann
>> --Victorian Bioinformatics Consortium, Monash University
>> --Tel +61 3 9905 9010
> Actually, I'd forgotten about spliced_seq(). That seems like it
> will Do The Right Thing. It's just up to the invoker to call
> spliced_seq() instead of seq() as appropriate.
> So, is there any other code that will break if I modify
> Bio::SeqFeature::Gene::Exon::location to not throw an exception when
> encountering Bio::Location::SplitLocationI? I'm wondering if it's
> just a paranoid check or if it's there to guard against something. If
> the latter, I need to know what code to fix. I'll dig and look, but
> if anybody knows or has an idea, save me some time. I suppose I can
> just change it and see what tests start failing. 8)
I'm wondering why you want to use Bio::SeqFeature::Gene::Exon to
describe the 'wrap-around' genes. The SeqFeature::Gene::Exon docs
state that the Exon class is used to specifically describe exons, as
the name implies. Exons are primarily eukaryotic in origin, so you
shouldn't encounter wraparounds, and should not have split locations
by definition (which likely explains the exception).
Wouldn't a SeqFeature::Generic work just as well using a split location?
More information about the Bioperl-l