[BioPython] FeatureLocation

Peter peter at maubp.freeserve.co.uk
Tue Apr 15 11:17:14 UTC 2008


On Tue, Apr 15, 2008 at 11:30 AM, Timothy Wu <2huggie at gmail.com> wrote:
> Hi,
>  I'm trying to access locations in a Genbank file like this:
>
>  mRNA            join(5172..5189,6254..6277,6744..6851,7266..7321,
>                      8515..8749,10663..10826,11626..11807,12784..12879,
>                      16701..17002,18265..18330,20541..20678,21188..21334,
>                      23789..23873,24587..24653,26071..26197,31165..31341,
>                      32922..33089,35124..35304,39185..39305,39545..39680,
>                      40530..40749,40901..41045,43199..43355,44727..44926,
>                      46662..46788,48405..48506,49040..49183,50001..50159,
>                      51315..51481,54580..54774,56455..56746)
>
>  I intend to get all the positions. However as it stands, it looks to me like
>  only the starting and ending position is supported in the FeatureLocation
>  class. Is there a way to obtain the coordinates?

They should be there as sub-features (a list property of the SeqFeature).

Peter



More information about the Biopython mailing list