@jobs_limit(1, "R")
@transform(runMetavelvet
           , suffix(".contigs.fa")
           , ".stats.pdf")
def plotCoverageHistogram(infile, outfile):
    inf = P.snip(infile, ".contigs.fa") + ".stats.txt"
    outf = P.snip(inf, ".txt") + ".pdf"
    R('''library(plotrix)''')
    R('''data = read.table("%s", header=TRUE)''' % inf)
    R('''pdf("%s", height = 7, width = 7 )''' % outf)
    R('''weighted.hist(data$short1_cov, data$lgth, breaks=seq(0, 200, by=1))''')
    R["dev.off"]()

###################################################################                                                                                                                                                                          
###################################################################                                                                                                                                                                          
###################################################################
splithost('//host[:port]/path') --> 'host[:port]', '/path'.
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_genomeassembly.py/plotCoverageHistogram--crossover with--dataset/ETHPy150Open ofermend/medicare-demo/socialite/jython/Lib/urllib.py/splithost