def buildGeneRegions(infile, outfile): '''annotate genomic regions with reference gene set. Only considers protein coding genes. In case of overlapping genes, only take the longest (in genomic coordinates). Genes not on UCSC contigs are removed. ''' PGeneset.buildGeneRegions(infile, outfile)