@follows(copyCapseqReplicatedBedFiles, mkdir("liver_vs_testes") )
@files( (PARAMS["compare_liver_pattern"]+".replicated.bed", PARAMS["compare_testes_pattern"]+".replicated.bed"), ("liver_vs_testes/liver_nmi.liver.testes.shared.bed", "liver_vs_testes/testes_nmi.liver.testes.shared.bed", "liver_vs_testes/liver_nmi.liver.testes.uniq.bed", "liver_vs_testes/testes_nmi.liver.testes.uniq.bed") )
def liverTestesCompare(infiles, outfile):
    liver, testes = infiles
    to_cluster = False
    
    statement = '''intersectBed -a %(liver)s -b %(testes)s -u > liver_vs_testes/liver_nmi.liver.testes.shared.bed;
                   intersectBed -a %(testes)s -b %(liver)s -u > liver_vs_testes/testes_nmi.liver.testes.shared.bed;
                   intersectBed -a %(liver)s -b %(testes)s -v > liver_vs_testes/liver_nmi.liver.testes.uniq.bed;
                   intersectBed -a %(testes)s -b %(liver)s -v > liver_vs_testes/testes_nmi.liver.testes.uniq.bed; ''' % locals()
    P.run()
        
############################################################
Returns all keywords that match a SQL pattern
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_proj007.py/liverTestesCompare--crossover with--dataset/ETHPy150Open boakley/robotframework-workbench/rwb/lib/keywordtable.py/KeywordTable.get_keywords