[Bioperl-l] Merging exons in all transcripts for a gene

[Elia Stupka]

> It seems non-trivial but not impossible and so before I attempt it I'd
> appreciate if anyone knows of existing code which does this. Doesn't
> appear to be in BioPerl at the moment.

Quite trivial as long as you represent your genes as seqfeatures (or as
Bio::SeqFeature::Gene::GeneStructure) and then make usre of the RangeI
inherited methods such as:

seqfeature->overlaps(seqfeature) #most useful in your case

and all these...

->contains
->equals
->intersection
->union

Elia

********************************
* http://www.fugu-sg.org/~elia *
* tel:    +65 6874 1467        *
* mobile: +65 9030 7613        *
* fax:    +65 6779 1117        *
********************************