@files( buildGeneAnnotations, "annotations_genes.counts" )
def makeGeneCounts( infile, outfile ):
    
    statement = """
    cat < %(infile)s |\
    python %(scriptsdir)s/gtf2table.py \
        --genome-file=genome \
        --counter=length \
        --log=%(outfile)s.log \
    > %(outfile)s
    """ 
    P.run()

###################################################################
###################################################################
###################################################################
coun gene exon statistics.
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_snps.py/makeGeneCounts