Change output file name and format of aevol_post_change_size_neutral_mut

src/post_treatments/aevol_post_change_size_neutral_mut.cpp

@@ -127,7 +127,9 @@ int main(int argc, char ** argv) {
   std::cout << " - output file = " << outputFile << std::endl;
   IOJson inputJson(inputFile);
 
-  out::init("result.txt", "mutation.txt");
+  std::string result = "result_seed_" + std::to_string(seed_prng) + ".txt";
+  std::string mutation = "mutation_seed_" + std::to_string(seed_prng) + ".csv";
+  out::init(result.c_str(), mutation.c_str());
 
   if (wanted_size == 0) {
     wanted_size = inputJson.getIndividuals()[0]->amount_of_dna() + delta;

src/post_treatments/libaevol_post/neutral_mutation_output.cpp

@@ -56,8 +56,8 @@ void out::init(const char *name_result,
   mutation_file << "################################################################################\n";
   mutation_file << "#\n";
   mutation_file << "# Header for R\n";
-  mutation_file << "gener gen_unit mut_type pos_0 pos_1 pos_2 pos_3 invert align_score align_score_2 seg_len"
-                   " repl_seg_len GU_len impact nbgenesatbreak nbgenesinseg nbgenesinreplseg"
+  mutation_file << "gener,gen_unit,mut_type,pos_0,pos_1,pos_2,pos_3,invert,align_score,align_score_2,seg_len"
+                   ",repl_seg_len,GU_len,impact,nbgenesatbreak,nbgenesinseg,nbgenesinreplseg"
                 << std::endl;
 }
 
@@ -124,21 +124,21 @@ int32_t mutation_event_length_mutation(MutationEvent &mutationEvent,
 void out::record_mutation_event(MutationEvent &mutationEvent,
                                 unsigned int number_generation,
                                 unsigned int size) {
-  mutation_file << number_generation << " ";
-  mutation_file << 0 << " ";  // GeneticUnit
-  mutation_file << mutationEvent.type() << " ";
-  mutation_file << mutationEvent.pos_1() << " ";
-  mutation_file << mutation_event_pos_2(mutationEvent) << " ";
-  mutation_file << mutation_event_pos_3(mutationEvent) << " ";
-  mutation_file << mutation_event_pos_4(mutationEvent) << " ";
-  mutation_file << mutation_event_invert(mutationEvent) << " ";
-  mutation_file << 0 << " "; // align_score
-  mutation_file << 0 << " "; // align_score2
-  mutation_file << mutation_event_length_mutation(mutationEvent, size) << " ";
-  mutation_file << -1 << " "; // rep_seg_length
-  mutation_file << size << " ";
-  mutation_file << 0 << " "; // Impact of the mutation on the metabolic error
-  mutation_file << 0 << " "; // Nb of coding RNAs possibly disrupted by the breakpoints
-  mutation_file << 0 << " "; // Nb of coding RNAs included in the segment
+  mutation_file << number_generation << ",";
+  mutation_file << 0 << ",";  // GeneticUnit
+  mutation_file << mutationEvent.type() << ",";
+  mutation_file << mutationEvent.pos_1() << ",";
+  mutation_file << mutation_event_pos_2(mutationEvent) << ",";
+  mutation_file << mutation_event_pos_3(mutationEvent) << ",";
+  mutation_file << mutation_event_pos_4(mutationEvent) << ",";
+  mutation_file << mutation_event_invert(mutationEvent) << ",";
+  mutation_file << 0 << ","; // align_score
+  mutation_file << 0 << ","; // align_score2
+  mutation_file << mutation_event_length_mutation(mutationEvent, size) << ",";
+  mutation_file << -1 << ","; // rep_seg_length
+  mutation_file << size << ",";
+  mutation_file << 0 << ","; // Impact of the mutation on the metabolic error
+  mutation_file << 0 << ","; // Nb of coding RNAs possibly disrupted by the breakpoints
+  mutation_file << 0 << ","; // Nb of coding RNAs included in the segment
   mutation_file << 0 << std::endl; // Nb of coding RNAs that were completely included in the replaced segment
 }