diff --git a/deep_sequencing_unit/source/bash/CreatePerCycleReadQuality.R b/deep_sequencing_unit/source/bash/CreatePerCycleReadQuality.R deleted file mode 100644 index 3baeb2ba76eb1af0e2007b39f6d43a994a4af286..0000000000000000000000000000000000000000 --- a/deep_sequencing_unit/source/bash/CreatePerCycleReadQuality.R +++ /dev/null @@ -1,11 +0,0 @@ -library(ShortRead) -args <- commandArgs(TRUE) -reads <- readFastq(args[1]) -if (length(reads) > 1000000) { - reads <- sample(reads,1000000) # if more than a million reads, sample randomly -} -qual <- FastqQuality(quality(quality(reads))) # get quality scores -readM <- as(qual, "matrix") # convert scores to matrix -pdf(file=paste(args[1],"boxplot.pdf", sep="_")) # Save box plot as boxplot.pdf in current folder -boxplot(readM, outline = FALSE, main="Per Cycle Read Quality", sub=args[1], xlab="Cycle", ylab="Phred Quality", col="grey") -dev.off() \ No newline at end of file