AnnotateVariants.RmdPlotting the number of variants by annotation category. One needs to first annotate carients using SnpEff. This may take a while.
library(pacbiowdlR)
variantsfile <- "location/small_variants.ann.txt.gz"
plot_snpEff(variantsfile, cores = 16)