How to get location of sequence in GenomicAnnotations

Specific Domains Biology, Health, and Medicine
1

Hello
I am newcomer in BioJulia and getting interested in it. Now I am starting to use GenomicAnnotations package. For example, I open a genbank file of E. coli K12 as

using GenomicAnnotations
K12 = readgbk(“K12.gbff”)
K12v = [ ]
for record in K12
K12v = [K12v; record]
end
println(K12v[1].genes[2])

and obtain
gene 190…255
/db_xref=“ASAP:ABE-0000006”
/db_xref=“ECOCYC:EG11277”
/db_xref=“GeneID:944742”
/gene=“thrL”
/locus_tag=“b0001”
/gene_synonym=“ECK0001”

I can see that the location is from 190 to 255.
But, how can I get this location in the code as a vector, for example?

Version of Julia : 1.8.5
Version of GenomicAnnotations : 0.3.3

Thank you.

1 Like
2

This is a little bit late, but take a look at this issue:

2 Likes
3

That is just what I have wanted!
Thank you very much.

1 Like