To work on query related to overlapping genomic features: use IntervalTree python pkg
http://informatics.malariagen.net/2011/07/07/using-interval-trees-to-query-genome-annotations-by-position/
http://hackmap.blogspot.de/2008/11/python-interval-tree.html
http://en.wikipedia.org/wiki/Interval_tree#Construction