diff --git a/noisyplotype.py b/noisyplotype.py index dd6320d3f66ac097ae7c50fa85c75145ca46e0b9..542f40457839627184f650bd042da1e911879e05 100644 --- a/noisyplotype.py +++ b/noisyplotype.py @@ -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)