[Bioperl-l] Subseq that transfers features too

Jason Stajich jason.stajich at duke.edu
Wed Nov 17 10:12:04 EST 2004


Off the top of my head, the easiest way to go about this is load data 
into Bio::DB::GFF, then you can grab slices of the sequence and the 
features will have the right coordinates.  Doing within Bioperl 
natively will require doing something with the coordinate mapping, but 
I don't know if the current implementations there try and solve this.

-jason

On Nov 17, 2004, at 4:28 AM, michael watson (IAH-C) wrote:

> Hi
>
> A quick question - I know trunc() will get me a sub sequence object
> without transferring features across, but is there a function that gets
> me a sub sequence that DOES transfer them?
>
> Thanks
>
> Mick
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at portal.open-bio.org
> http://portal.open-bio.org/mailman/listinfo/bioperl-l
>
>
--
Jason Stajich
jason.stajich at duke.edu
http://www.duke.edu/~jes12/



More information about the Bioperl-l mailing list