@transform( getCapseqRNAseqTSSDistance, suffix( ".rnaseq.tss.distance"), ".rnaseq.tss.distance.load" )
def loadCapseqRNAseqTSSDistance( infile, outfile ):
    track= P.snip( os.path.basename(infile), ".rnaseq.tss.distance").replace(".","_").replace("-","_")
    statement = """cat %(infile)s | python ~/src/csv2db.py 
                         --database=%(database)s
                         --table=%(track)s_rnaseq_tss_distance
                         --index=gene_id
                         --index=closest_id
                         --index=id5
                         --index=id3
                 > %(outfile)s; """
    P.run()

############################################################
Load interval annotations: distance to lncRNA transcription start sites
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_proj007.py/loadCapseqRNAseqTSSDistance