knownGene = GenomeIntervalTree.from_table(url=knownGene_localurl, decompress=True) # Py3 downloads .gz files to local files with names not ending with .gz
assert len(knownGene) == 82960
- result = knownGene[b'chr1'].search(100000, 138529)
+ result = knownGene[b'chr1'].overlap(100000, 138529)