diff --git a/src/bin/kmerclust_plot.R b/src/bin/kmerclust_plot.R
new file mode 100644
index 0000000000000000000000000000000000000000..57d6408ce00ba8b603fb19faf58fedf28f6fc497
--- /dev/null
+++ b/src/bin/kmerclust_plot.R
@@ -0,0 +1,22 @@
+library(kmerclust)
+library(ggplot2)
+
+args <- commandArgs(trailingOnly = TRUE)
+print(args)
+load(file = paste0(args[3], ".Rdata"))
+load(file = paste0(args[3], "_boot.Rdata"))
+
+res %>%
+  ggplot(aes(x = name, y = BIC)) +
+  geom_violin() %>%
+  theme_bw()
+ggsave(paste0(args[3], "_BIC.pdf"))
+res %>%
+  ggplot(aes(x = name, y = WSS_f / BSS)) +
+  geom_violin() %>%
+  theme_bw()
+ggsave(paste0(args[3], "_WSS_BSS.pdf"))
+count %>%
+  ggplot(aes(x = count_m, y = count_f)) %>%
+  theme_bw()
+ggsave(paste0(args[3], "_scatter.pdf"))