Set RNG seed for diversity sampling

Author: ahmig

workflow/rules/report.smk

@@ -64,6 +64,7 @@ rule diversity_data:
     params:
         bootstrap_reps = config["DIVERSITY_REPS"],
         aln_reference = config["ALIGNMENT_REFERENCE"],
+        seed = 7291,
     input:
         study_fasta = OUTDIR/"nextalign"/f"{OUTPUT_NAME}.aligned.masked.fasta",
         context_fasta = OUTDIR/"context"/"nextalign"/"context_sequences.aligned.masked.fasta",

workflow/scripts/report/diversity_data.R

@@ -5,6 +5,8 @@ log <- file(snakemake@log[[1]], open = "wt")
 sink(log, type = "message")
 sink(log, type = "output")
 
+set.seed(snakemake@params$seed)
+
 library(ape)
 library(pegas)
 library(future.apply)