Skip to content
Snippets Groups Projects
Commit ac36dfde authored by Cecile Triay's avatar Cecile Triay
Browse files

re-correct output file names

parent 0224a23c
No related branches found
No related tags found
No related merge requests found
......@@ -209,7 +209,7 @@ def generate_individuals (nb_individuals, chromosome_size, marker_density, mean_
#markers_positions = sorted(random.sample(range(chromosome_size),size))
markers_positions = np.linspace(1, chromosome_size, markers_nb, dtype="int")
with open('Breakpoints_3x.csv', 'w') as breakpointFile:
with open('Breakpoints.csv', 'w') as breakpointFile:
breakpointFile.write(",".join(["sample","average_bkp_position", "bkp_start_position", "bkp_stop_position", "transitionType"]) + "\n");
for i in range(nb_individuals):
print("individual " + str(i+1))
......@@ -276,7 +276,7 @@ header = [
## Write the VCF files (with and without errors)
with open('test_3x.vcf', 'w') as file:
with open('test.vcf', 'w') as file:
# Write the header of the file
file.writelines(header)
......@@ -292,7 +292,7 @@ with open('test_3x.vcf', 'w') as file:
row.append("1/1:"+ str(round(mean_depth)) + ":0," + str(round(mean_depth)) + ":.:.:.:.")
file.write("\t".join(row) + "\n");
with open('test_error_3x.vcf', 'w') as file:
with open('test_error.vcf', 'w') as file:
# Write the header of the file
file.writelines(header)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment