diff --git a/R/functional_analysis_library.R b/R/functional_analysis_library.R
index f01579d..e5dc0d0 100644
--- a/R/functional_analysis_library.R
+++ b/R/functional_analysis_library.R
@@ -360,6 +360,8 @@ multi_topGOTest <- function(funsys, genes_list, nodeSize = 5, org_db,
                                  p.value = result,
                                  topNodes = length(result@score), 
                                  format.FUN = function(x, dig, eps){x})
+    colnames(res_table)[6] <- "p.value"
+
     if(clean_parentals)
       res_table <- clean_topGO_parentals(res_table, funsys, p_value_cutoff)
     res_table$fdr <- p.adjust(res_table$p.value, method = "BH", n = length(res_table$p.value))