diff --git a/Snake.config.json b/Snake.config.json
index 98afdad1c9c6e534e4eaec1d4fb32891e4aad4b7..53cd157fff0b3e35074034d689688465637dbbf5 100644
--- a/Snake.config.json
+++ b/Snake.config.json
@@ -16,6 +16,7 @@
         "50000"     : "utils/variable_bins.GRCh38.50kb.bed",
         "100000"    : "utils/variable_bins.GRCh38.100kb.bed"
     },
+    "R_reference"   : "BSgenome.Hsapiens.UCSC.hg38",
 
     "bp_density"    : {
         "few"       : 0.1,
diff --git a/Snakefile b/Snakefile
index f4688c2fe0ac3a1b05731f314a3694305ac76423..b29a56abc83d11b2af1461c148ad44b43504f8fb 100644
--- a/Snakefile
+++ b/Snakefile
@@ -11,7 +11,6 @@ BAM, = glob_wildcards("bam/{bam}.bam")
 
 rule all:
     input:
-        "log/StrandPhaseR.config",
         expand("plots/" + config["sample"] + ".{window}_fixed.pdf", window = [50000, 100000, 200000, 500000]),
         expand("plots/" + config["sample"] + ".{window}_variable.pdf", window = [50000, 100000]),
         expand("segmentation2/" + config["sample"] + ".{window}_fixed.{bpdens}.txt",
@@ -287,7 +286,7 @@ rule prepare_strandphaser_config_per_chrom:
             print("compareSingleCells = TRUE",     file = f)
             print("callBreaks       = FALSE",    file = f)
             print("exportVCF        = '", config["sample"], ".txt'", sep = "", file = f)
-            print("bsGenome         = 'BSgenome.Hsapiens.UCSC.hg19'", file = f)
+            print("bsGenome         = '", config["R_reference"], "'", file = f)