Found this python script for quick computation of overlapping features in genomic interval databases. The script can be downloaded from the author's FTP site indicated in this paper. It's a small script and very easy to use. It works with standard GFF format.