@files( ( os.path.join( PARAMS["geneset_dir"],PARAMS["geneset_gene_tss"] ), PARAMS["bed_ucsc_cgi"]), "cgi/cgi.gene_tss.overlap.count" )
def getCGIGeneTSSOverlapCount( infiles, outfile ):
    tss, cgi = infiles
    to_cluster = True
    statement = """echo "Predicted CGIs overlapping 1 or more TSS" > %(outfile)s; intersectBed -a %(cgi)s -b %(tss)s -u | wc -l >> %(outfile)s; 
                   echo "Predicted CGIs not overlapping any TSS" >> %(outfile)s; intersectBed -a %(cgi)s -b %(tss)s -v | wc -l >> %(outfile)s; 
                   echo "TSS overlapped by 1 or more CGI" >> %(outfile)s; intersectBed -a %(tss)s -b %(cgi)s -u | wc -l >> %(outfile)s; 
                   echo "TSS not overlapped by any predicted CGI" >> %(outfile)s; intersectBed -a %(tss)s -b %(cgi)s -v | wc -l >> %(outfile)s; 
                   sed -i '{N;s/\\n/\\t/g}' %(outfile)s; """
    P.run()

############################################################
Change the stored password
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_proj007.py/getCGIGeneTSSOverlapCount--crossover with--dataset/ETHPy150Open MacSysadmin/pymacadmin/lib/PyMacAdmin/Security/Keychain.py/GenericPassword.update_password