diff --git a/Snakefile b/Snakefile
index 8a1cf44..9296446 100644
--- a/Snakefile
+++ b/Snakefile
@@ -309,6 +309,7 @@ rule all:
 
 
 if (config["bam_rmdup_realign_indels"]
+    or config["merge_vcfs_per_dataset"]
     or config["mlRho"]
     or config["pca"]
     or config["ROH"]
diff --git a/workflow/report/missingness_filtered_vcf_multiqc.rst b/workflow/report/missingness_filtered_vcf_multiqc.rst
new file mode 100644
index 0000000..08f02d3
--- /dev/null
+++ b/workflow/report/missingness_filtered_vcf_multiqc.rst
@@ -0,0 +1 @@
+MultiQC report, summarizing statistics for merged and filtered VCF files including {{ snakemake.wildcards.dataset }} samples.
\ No newline at end of file
diff --git a/workflow/rules/9_merge_vcfs.smk b/workflow/rules/9_merge_vcfs.smk
index 7a998a3..3766baf 100644
--- a/workflow/rules/9_merge_vcfs.smk
+++ b/workflow/rules/9_merge_vcfs.smk
@@ -612,7 +612,10 @@ rule missingness_filtered_vcf_multiqc:
     input:
         missingness_filtered_vcf_multiqc_inputs,
     output:
-        stats="results/{dataset}/vcf/" + REF_NAME + "/stats/vcf_merged_missing/multiqc/multiqc_report.html",
+        stats=report(
+            "results/{dataset}/vcf/" + REF_NAME + "/stats/vcf_merged_missing/multiqc/multiqc_report.html",
+            caption="../report/missingness_filtered_vcf_multiqc.rst",
+            category="VCF file processing",),
     params:
         indir="results/{dataset}/vcf/" + REF_NAME + "/stats/vcf_merged_missing/",
         outdir="results/{dataset}/vcf/" + REF_NAME + "/stats/vcf_merged_missing/multiqc",