bio/Macs2: fix call wrt to temp dir

macs2 has a weird and hidden constraint on having enough space in tmp space, see macs3-project/MACS#122
pveber · May 11, 2022 · 49c7a0e · 49c7a0e
1 parent d363bd2
commit 49c7a0e
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/lib/bio/macs2.ml b/lib/bio/macs2.ml
@@ -54,6 +54,7 @@ let callpeak_gen
     ?broad ?pvalue ?qvalue ?gsize ?call_summits
     ?fix_bimodal ?mfold ?extsize ?nomodel ?bdg ?control ?keep_dup format treatment =
   Workflow.shell ~descr:"macs2.callpeak" ~img [
+    mkdir_p tmp ;
     macs2 "callpeak" [
       option (flag string "--broad") broad ;
       opt "--outdir" Fn.id dest ;
@@ -71,6 +72,7 @@ let callpeak_gen
       option (opt "--keep-dup" keep_dup_expr) keep_dup ;
       option (opt "--control" (list ~sep:" " dep)) control ;
       opt "--treatment" (list ~sep:" " dep) treatment ;
+      opt "--tempdir" Fun.id tmp ;
     ]
   ]