@transform( annotateGenesetCapseqOverlap, suffix(".genes_capseq_overlap"), ".genes_capseq_overlap.load" )
def loadGenesetCapseqOverlap( infile, outfile ):
    geneset_name = PARAMS["geneset_name"]
    track= P.snip( os.path.basename(infile), ".genes_capseq_overlap").replace(".","_").replace("-","_")
    statement = """cat %(infile)s | python ~/src/csv2db.py 
                         --database=%(database)s
                         --table=%(track)s_%(geneset_name)s_genes_capseq_overlap
                         --index=gene_id
                 > %(outfile)s; """
    P.run()
            
############################################################
############################################################
## Section 1b: Count overlap of CAPseq intervals with gene/transcript TSSs
:type n: int
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_proj007.py/loadGenesetCapseqOverlap--crossover with--dataset/ETHPy150Open kamyu104/LeetCode/Python/integer-break.py/Solution.integerBreak