Skip to content
Snippets Groups Projects
Commit cfea4c00 authored by Rohmer Coralie's avatar Rohmer Coralie
Browse files

Sorting in the total output

parent 47a82eb5
Branches
No related tags found
No related merge requests found
......@@ -104,8 +104,8 @@ for threshold in files :
tab_line=line.split(",")
if not (tab_line[0]=="MSA") :
msa=tab_line[0]
base_pair=tab_line[1]
nb_read=tab_line[2]
base_pair=int(tab_line[1])
nb_read=int(tab_line[2])
if (base_pair not in data[exp_name][dir]):
data[exp_name][dir][base_pair] = {}
......@@ -178,10 +178,14 @@ for threshold in files :
output_all.write("\n")
# Writing data
bases=sorted(bases)
MSA=sorted(MSA)
reads=sorted(reads)
print(reads,bases,MSA)
for base in bases :
for read in reads :
output_mean.write(base + "," + read + ",")
output_all.write(base + "," + read + ",")
output_mean.write(str(base) + "," + str(read) + ",")
output_all.write(str(base) + "," + str(read) + ",")
for exp_name in data["order"]:
for dir in data[exp_name]["order"]:
if dir == "":
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment