I have a very long list of coordinates that looks like that: `genome_build:chrN:coordinate`

I want to know if the given coordinate hits a gene or whatever else. What is the most efficient way to do so using python?

Example: having hg19:chr8:119331103 as input, it should return [SAMD12]
Hi,

You could use [Ensembl's rest API][1] for that using requests.

But if you have a very long list I would suggest that you try [pyGeno][2], it's typically the kind of applications if was made for:

```
x1 = 119331103
g = Genome(name = "GRCh37.75")
for gene in g.get(Gene, {"start >=": x1, "end <" : x1, "chromosome.number" : "8"}) :
        print gene
```

 [1]: http://http://rest.ensembl.org/documentation/info/overlap_id
 [2]: http://pygeno.iric.ca/