view mapping_quality_stats.r @ 0:f00479673d47 draft

planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/mapping_quality_stats commit e4b37874b820a2ac48732667128a08e5755b7c4b

## Setup R error handling to go to stderr
options(show.error.messages = FALSE,
        error = function() {
            cat(geterrmessage(), file = stderr())
            q("no", 1, FALSE)
        }
)


warnings()
library(optparse)
library(ggplot2)

option_list <- list(
    make_option(c("-i", "--input"), type = "character", help = "Path to tabular file"),
    make_option(c("-o", "--output"), type = "character", help = "path to the pdf plot")
    )

parser <- OptionParser(usage = "%prog [options] file", option_list = option_list)
args <- parse_args(parser)

# data frame implementation
table <- read.delim(args$input, header = TRUE)
colnames(table) <- c("MAPQ", "Counts")


# Barplot
pdf(file = args$output)
ggplot(table, aes(x = MAPQ, y = Counts)) +
    geom_bar(stat = "identity")
devname <- dev.off()