diff --git a/build_pop.py b/build_pop.py
index 845660c47cbd7dca3740b9c0d97c5228631c3379..1528dda4d57f5f8739f8ae92ac143170dd3f57c2 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))