@transform( exportCapseqNoncodingTSSGeneList, suffix( ".noncoding.tss.distance.export"), ".noncoding.tss.distance.export.load" )
def loadCapseqNoncodingTSSGeneList( infile, outfile ):
    track = P.snip( os.path.basename( infile ), ".noncoding.tss.distance.export" ).replace("-","_").replace(".","_")
    geneset_name = PARAMS["geneset_name"]
    statement = """cat %(infile)s | python %(scriptsdir)s/csv2db.py 
                         --database=%(database)s
                         --table=%(track)s_%(geneset_name)s_interval_noncoding_mapping
                         --index=gene_id
                         --index=interval_id
                   > %(outfile)s; """
    P.run()

############################################################    
############################################################
## External linCRNA datasets
(internal) Got an answer to a previously sent message.
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_proj007.py/loadCapseqNoncodingTSSGeneList--crossover with--dataset/ETHPy150Open nlloyd/SubliminalCollaborator/libs/twisted/spread/pb.py/Broker.proto_answer