diff --git a/build_pop.py b/build_pop.py index bc4ec727adeaf86f7158893fd303c2e0e1692006..845660c47cbd7dca3740b9c0d97c5228631c3379 100755 --- a/build_pop.py +++ b/build_pop.py @@ -226,7 +226,7 @@ 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 " + "-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 "