[Bioperl-l] split seq feature and fuzzy feature proposal
Fri, 19 Jan 2001 13:14:53 +0100
Hilmar Lapp wrote:
> Yeah, that's the really hairy case. We probably should define
> first what we would like to be able to do with compound locations.
> This is a strong call for feedback: what do people out there using
> the package intend to do with compound locations? E.g. if you draw
> annotations, would you just draw the part referring to the
> attached seq? Ensembl people, any experience/wishlists for this?
I hope do not mind me giving some comments on this issue.
I am writing some programs to automatically extract genes and intergenic
regions from DNA sequences. So, I am mostly interested in the type of a
feature and also its start and end position in the sequence. The main problem
I am facing is that sometimes a feature is not extracted from the sequence
because it has a fuzzy location.
eg. if the location of a CDS is described as "join(AL101010.1:1..201,123..245)
this CDS is not add to the list of feature and it is impossible to do anything
usefull with this sequence for me.
In my opinion, I think it is important that a feature is created even if the
location is fuzzy. When there is a problem, it should be possible to access
the description of the location.
+ Gert Thijs
+ email: firstname.lastname@example.org
+ homepage: http://www.esat.kuleuven.ac.be/~thijs
+ Kasteelpark Arenberg 10
+ B-3001 Leuven-Heverlee
+ Tel : +32 16 32 18 84
+ Fax : +32 16 32 19 70