@transform( exportLiverTestesTSSTranscriptList, suffix( ".transcript.tss.distance.export"), ".transcript.tss.distance.export.load" )
def loadLiverTestesTSSTranscriptList( infile, outfile ):
    geneset_name = PARAMS["geneset_name"]
    statement = """cat %(infile)s | python %(scriptsdir)s/csv2db.py 
                         --database=%(database)s
                         --table=liver_testes_merged_%(geneset_name)s_interval_transcript_mapping
                         --index=gene_id
                         --index=interval_id
                 > %(outfile)s; """
    P.run()
        
############################################################
Load CAPseq interval annotations: distance to transcript transcription start sites
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_proj007.py/loadLiverTestesTSSTranscriptList