initial commit

LindoNkambule · Aug 10, 2019 · 41a1b45 · 41a1b45
1 parent a022abf
commit 41a1b45
Showing 1 changed file with 15 additions and 0 deletions.
diff --git a/VCFVisualize.R b/VCFVisualize.R
@@ -0,0 +1,15 @@
+library(data.table)
+library(packHV) # required for plotting the histogram and boxplot together
+
+file <- fread('gatk.vcf')
+print (file)
+quality <- file$QUAL
+position <- file$POS
+
+quality_histogram <- function(x){
+  hist_boxplot(x,
+       main = "Quality Distribution",
+       xlab = "QUAL",
+       col = 'cyan')
+}
+quality_histogram(quality)