@transform( annotateLiverTestesMergedTranscriptTSSDistance, suffix( ".transcript.tss.distance"), ".transcript.tss.distance.load" )
def loadLiverTestesMergedTranscriptTSSDistance( infile, outfile ):
    track= P.snip( os.path.basename(infile), ".transcript.tss.distance").replace(".","_").replace("-","_")
    geneset_name = PARAMS["geneset_name"]
    statement = """cat %(infile)s | python ~/src/csv2db.py 
                         --database=%(database)s
                         --table=liver_testes_merged_%(geneset_name)s_transcript_tss_distance
                         --index=gene_id
                         --index=closest_id
                         --index=id5
                         --index=id3
                 > %(outfile)s; """
    P.run()

############################################################
Load CAPseq interval annotations: distance to transcript transcription start sites
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_proj007.py/loadLiverTestesMergedTranscriptTSSDistance