@follows( mkdir("plots"), getReplicatedTranscriptTSSProfileCapseq )
@transform("tss-profile/*.replicated.transcript.tss-profile.capseq.counts.tsv.gz", regex(r"tss-profile/(\S+).replicated.transcript.tss-profile.capseq.counts.tsv.gz"), r"plots/\1.combined.tss-profile.pdf")
def plotFigure2b( infile, outfile):
    capseq = infile
    scriptsdir = PARAMS["scriptsdir"]
    nocapseq = capseq.replace("capseq", "nocapseq")
    R('''source("%(scriptsdir)s/R/proj007/proj007.R") ''' % locals() )
    R('''combinedTSSPlot(capseqfile="%(capseq)s", nocapseqfile="%(nocapseq)s", outfile="%(outfile)s", ylimit=c(0,10), scale=1)''' % locals() )
    
############################################################
Test that with replace_existing=True, a new job replaces another with the same id.
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_proj007.py/plotFigure2b--crossover with--dataset/ETHPy150Open agronholm/apscheduler/tests/test_schedulers.py/TestBaseScheduler.test_add_job_replace