@jobs_limit(2)
    @transform( "*.pileup.gz", suffix(".pileup.gz"), "_pileup.load")
    def loadPileup( infile, outfile ):

        to_cluster = False
        tablename = outfile[:-len(".load")]
        statement = '''
        gunzip < %(infile)s |
        awk 'BEGIN {printf("contig\\tpos\\treference\\tgenotype\\n")} 
                   {printf("%%s\\t%%s\\t%%s\\t%%s\\n",$1,$2,$3,$4)}' |
           python %(scriptsdir)s/csv2db.py %(csv2db_options)s \
            --index=contig,pos \
            --table=%(tablename)s
        > %(outfile)s
        '''

        P.run()
load pileup information.
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_snps.py/loadPileup