@follows( mkdir("plots") )
@merge("/ifs/projects/proj007/web/capseq7/gc/*testes-cap.replicated.gc.export", ("plots/testes_capseq_gc_content.pdf","plots/testes_control_gc_content.pdf") )
def plotFigure1cGCContent( infiles, outfiles):
    capseq_out, control_out = outfiles
    indir = os.path.dirname(infiles[0])
    scriptsdir = PARAMS["scriptsdir"]
    R('''source("%(scriptsdir)s/R/proj007/proj007.R") ''' % locals())
    R('''speciesPlot(dir="%(indir)s", pattern="*testes-cap.replicated.gc.export", main="Testes CAPseq", xlab="GC Content", filename="%(capseq_out)s", plotcol=2, xlimit=c(0,1), ylimit=c(0,15))''' % locals() )
    R('''speciesPlot(dir="%(indir)s", pattern="*testes-cap.replicated.gc.export", main="Testes Control", xlab="GC Content", filename="%(control_out)s", plotcol=3, xlimit=c(0,1), ylimit=c(0,15))''' % locals() )
  
  
############################################################
############################################################
############################################################
## Pipeline organisation
Figure 1c: density plots of GC content
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_species_conservation.py/plotFigure1cGCContent