[Bioperl-l] CDS across contigs

Mike Muratet muratetm@torchconcepts.com
Mon, 08 Jul 2002 12:57:37 -0500


Greetings

I recall this problem has surfaced before and maybe been answered, but I
can't find where.

One occasionally comes across CDS features defined as:

     CDS
join(AB005248.1:83409..83646,1..236,315..437,525..641, 718..771)

(This one is from A. thaliana.) The location() method will return a
SplitLocationI object, but you can't get the sequence by the subseq()
method, of course. Is the solution to get the id from the sublocation
and fetch the sequence somehow?

Thanks

Mike