@transform( (getGenesetCapseqOverlapList, getGenesetCapseqOverlapControlList), suffix(".genelist"), ".gtf.gz" )
def getOverlappedGeneGTF( infile, outfile ):
    gene_file = os.path.join( PARAMS["geneset_dir"], PARAMS["geneset_gene_profile"])
    statement = '''zcat %(gene_file)s 
                   | python %(scriptsdir)s/gtf2gtf.py --filter=gene --apply=%(infile)s --log=%(outfile)s.log
                   | python %(scriptsdir)s/gtf2gtf.py --join-exons --log=%(outfile)s.log
                   | sed s/\\\\ttranscript\\\\t/\\\\texon\\\\t/g 
                   | gzip > %(outfile)s; '''
    P.run()

############################################################
## CAPseq profile over overlapped genes
Return whether this plot is in a given row of a GridPlot.
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_proj007.py/getOverlappedGeneGTF--crossover with--dataset/ETHPy150Open bokeh/bokeh/bokeh/models/plots.py/Plot.row