Now computing add_indel seed from global dataset seed.

lib/dataset.ml

@@ -30,5 +30,7 @@ let add_indels_to_dataset d ~seed =
   let model_prefix = sprintf "%s_0.33_i" d.model_prefix in
   let tree_prefix  = d.tree_prefix in
   let is_real      = d.is_real in
-    let dataset      = Ready_dataset.add_indels_to_ready_dataset ~p ~seed d.dataset in
-    {model_prefix; tree_prefix; is_real; dataset; seed = Hashtbl.hash d.seed}
+  let new_seed     = Hashtbl.hash seed in
+  let indel_seed   = Hashtbl.hash new_seed in
+  let dataset      = Ready_dataset.add_indels_to_ready_dataset ~p ~seed:indel_seed d.dataset in
+  {model_prefix; tree_prefix; is_real; dataset; seed = new_seed}

lib/ready_dataset.ml

@@ -86,7 +86,7 @@ let add_indels_to_fna ~(p:float) ~(seed:int) (fna:nucleotide_fasta workflow) : n
 let add_indels_to_ready_dataset ~p ~seed d =
   let r_d = to_raw d in
   (* VL: one seed per fna. FIXME:should be computed from global seed and task info instead *)
-  let fna = add_indels_to_fna ~p ~seed:(Random.int Int.max_value) r_d.fna in
+  let fna = add_indels_to_fna ~p ~seed r_d.fna in
   let fna_infos = r_d.fna_infos in
   let ready_dataset = of_raw {Raw_dataset.input_tree=r_d.input_tree ; fna; fna_infos} in
   ready_dataset