[Bioperl-l] Truncate sequence with features
heikki at sanbi.ac.za
Fri Jun 9 03:35:12 EDT 2006
The definitive document describing the locations is the feature table
but you probably know that already.
Two questions come to mind:
1. Can you parse your joint location using bioperl without errors?
2. Is there a practical advantage in including a location which has no
relevance to the sequence in hand?
I notice that the /partial qualifier is deprecated and the docs suggest using
</> signs to indicate that the sequence is partial, so I guess what you are
doing is correct.
On Thursday 08 June 2006 18:31, Roy Chaudhuri wrote:
> Hi all.
> I've been playing around with a subroutine to truncate a sequence and
> adjust the coordinates of any features that overlap the specified
> region- something that according to the comments in
> Bio::Location::Simple has been abortively worked on in the past.
> I've submitted the subroutine as an enhancement in Bugzilla. It's a bit
> hacky but works for what I needed it for. However I'm a bit unsure on
> the best way to deal with split locations where one of the sublocations
> is entirely outside the truncated region. My current method results in
> locations like:
> join(1..500, >1000..>1000)
> which is quite ugly and possibly invalid, but kind of makes sense. Does
> anyone know what would be the correct behaviour for this situation?
> Dr. Roy Chaudhuri
> Bioinformatics Research Fellow
> Division of Immunity and Infection
> University of Birmingham, U.K.
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
______ _/ _/_____________________________________________________
_/ _/ _/ Heikki Lehvaslaiho heikki at_sanbi _ac _za
_/_/_/_/_/ Associate Professor skype: heikki_lehvaslaiho
_/ _/ _/ SANBI, South African National Bioinformatics Institute
_/ _/ _/ University of Western Cape, South Africa
_/ Phone: +27 21 959 2096 FAX: +27 21 959 2512
More information about the Bioperl-l