From 3790e7da275900ad333d83c7248111dd4d884eaf Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Flor=C3=A9al=20Cabanettes?= <floreal.cabanettes@inra.fr> Date: Wed, 19 Apr 2017 16:22:38 +0200 Subject: [PATCH] Fix build pop script --- build_pop.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/build_pop.py b/build_pop.py index 845660c..1528dda 100755 --- a/build_pop.py +++ b/build_pop.py @@ -226,11 +226,11 @@ print("GENERATE RANDOM READS FOR EACH INDIVIDUAL FROM GENOME...\n") # Generate reads for all individuals: if not haploid: cmd = str("{4}pirs/pirs simulate -z -x {3} -d -B {4}pirs/Profiles/Base-Calling_Profiles/humNew.PE100.matrix.gz " - "-I {4}pirs/Profiles/InDel_Profiles/phixv2.InDel.matrix -l {5} -m {6} -v {7}" + "-I {4}pirs/Profiles/InDel_Profiles/phixv2.InDel.matrix -l {5} -m {6} -v {7} " "-G {4}pirs/Profiles/GC-depth_Profiles/humNew.gcdep_100.dat -o {0} {1} {2}") else: cmd = str("{4}pirs/pirs simulate -z -x {3} -B {4}pirs/Profiles/Base-Calling_Profiles/humNew.PE100.matrix.gz " - "-I {4}pirs/Profiles/InDel_Profiles/phixv2.InDel.matrix -l {5} -m {6} -v {7}" + "-I {4}pirs/Profiles/InDel_Profiles/phixv2.InDel.matrix -l {5} -m {6} -v {7} " "-G {4}pirs/Profiles/GC-depth_Profiles/humNew.gcdep_100.dat -o {0} {1}") for i in range(1, nb_inds+1): prefix = os.path.join(output_dir, "INDIV_" + str(i)) -- GitLab