diff --git a/src/bin/kmerclust_merge.R b/src/bin/kmerclust_merge.R
index cadf100beea16e7fcb4f034dea5e4efa18129c0f..675b3da2e7ee20067c7329b4b69d8218c640be02 100644
--- a/src/bin/kmerclust_merge.R
+++ b/src/bin/kmerclust_merge.R
@@ -3,10 +3,11 @@ library(kmerclust)
 args <- commandArgs(trailingOnly = TRUE)
 print(args)
 res_files <- list.files(pattern = "Rdata$")
+print(res_files)
 res_total <- NULL
 for (res_file in res_files) {
   load(res_file)
   res_total <- rbind(res_total, res)
 }
 res <- res_total
-save(res, file = paste0(args[1], "_boot", ".Rdata"))
+save(res, file = paste0(args[1], "_boot.Rdata"))
diff --git a/src/bin/kmerclust_plot.R b/src/bin/kmerclust_plot.R
index 42d1aab5720639d0b9d557aea20b267cdef1b6ad..9f7bca8e8fbd19f9f2a02ddd2fcb8fb7d2378f64 100644
--- a/src/bin/kmerclust_plot.R
+++ b/src/bin/kmerclust_plot.R
@@ -6,17 +6,18 @@ print(args)
 load(file = paste0(args[1], ".Rdata"))
 load(file = paste0(args[1], "_boot.Rdata"))
 
+print(res)
 res %>%
   ggplot(aes(x = name, y = BIC)) +
-  geom_violin() %>%
+  geom_violin() +
   theme_bw()
 ggsave(paste0(args[1], "_BIC.pdf"))
 res %>%
   ggplot(aes(x = name, y = WSS_f / BSS)) +
-  geom_violin() %>%
+  geom_violin() +
   theme_bw()
 ggsave(paste0(args[1], "_WSS_BSS.pdf"))
 count %>%
-  ggplot(aes(x = count_m, y = count_f)) %>%
+  ggplot(aes(x = count_m, y = count_f)) +
   theme_bw()
 ggsave(paste0(args[1], "_scatter.pdf"))
diff --git a/src/main.nf b/src/main.nf
index 423c4d3382d6c430b0560741348725f3246fa0e5..174ca7592279aede954ed4d5528e0fc4b65f1860 100644
--- a/src/main.nf
+++ b/src/main.nf
@@ -33,6 +33,6 @@ workflow {
   COLLATEKMER(MERGEKMER.out.csv.map{it -> [it[0].specie, it[1]] }.groupTuple())
   KMERCLUST_LOAD(COLLATEKMER.out.csv, params_csv.collect())
   KMERCLUST(Channel.of(1..params.bootstrap).combine(KMERCLUST_LOAD.out.rdata))
-  KMERCLUST_MERGE(KMERCLUST.out.rdata.map{it -> [it[1], it[2]]}.groupTuple())
-  KMERCLUST_PLOT(KMERCLUST_MERGE.out.rdata)
+  KMERCLUST_MERGE(KMERCLUST.out.rdata.groupTuple())
+  KMERCLUST_PLOT(KMERCLUST_MERGE.out.rdata.mix(KMERCLUST_LOAD.out.rdata).groupTuple())
 }