[Bioperl-l] Re: Sim4

Ewan Birney birney@ebi.ac.uk
Fri, 28 Jul 2000 19:46:01 +0100 (BST)


On Fri, 28 Jul 2000 hilmar.lapp@pharma.novartis.com wrote:

> 
> 
> 
> 
> 
> It is also using the outdated Bio::SeqFeature::Homol. I'll change that to
> the newer Bio::SeqFeature::FeaturePair.
> 
> 
>      Fine. This was my idea, too.
> 
>      There's a little semantic problem left with the notion of Exon for
>      Sim4 alignments and the notion of feature and hit in FeaturePair. For
>      Sim4, you put the shorter sequence first, so that some people might be
>      referring to this as the query. In the FeaturePair object, query will
>      however refer to the sequence having the Exons, that is, the second
>      (genomic) sequence, because you obviously expect
>           $exon_start = $exon->start()
>      return the right result. $exon->hstart() would then return the
>      corresponding coordinate on the 'other' sequence.
> 

Indeed.

>      I'll check that Sim4/Result.pm correctly initializes the Exon object,
>      correctly meaning in the above sense.
> 
>      Then there's still ExonSet which is missing, but I haven't been able
>      to find any special functionality this is supposed to provide. Do you
>      remember the rationale behind this module? Could it be replaced by a
>      2-dimensional array of Exon objects? Or, easier, and maybe more
>      elegant, to handle, an array of Bio::Seq objects (without actual
>      sequences) each having an array of Exon objects as Feature.
> 


I very much dislike the Bio::Seq idea.


>      Another possibility, which I tend to favor: process the Sim4 result in
>      a similar way to multiple sequence files, that is, the user is
>      requested to call repeatedly something like parse_next_alignment(),
>      with each call returning an array of Exon objects.
> 
>      If there are no objections I'll go for that.


sounds good to me.

> 
>           Hilmar
> 
> 
> 
> 
> 
> 
> _______________________________________________
> 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>. 
-----------------------------------------------------------------