@transform( "*.genelist", regex( r"(\S+).genelist"), r"\1.genelist.load" )
def loadGeneLists( infile, outfile ):

    track = P.snip( os.path.basename( infile), ".genelist" ).replace(".","_").replace("-","_")
    statement = '''cat %(infile)s
                   | python %(scriptsdir)s/csv2db.py
                       --header=gene_id
                       --index=gene_id
                       --table=%(track)s_genelist
                   > %(outfile)s '''
    P.run()

############################################################
generate a covering set of segments with overlap
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_species_conservation.py/loadGeneLists--crossover with--dataset/ETHPy150Open CGATOxford/cgat/obsolete/benchmark_ranges.py/generate_overlapping_segments