Correctly compute nb of recombinations_ in DnaReplicationReport

src/libaevol/io/DnaReplicationReport.cpp

@@ -303,6 +303,14 @@ void DnaReplicationReport::compute_stats()
            mut->mut_type() == S_DEL);
     nb_mut_[mut->mut_type()]++;
   }
+
+#ifdef __EUKARYOTE
+  for (const auto& recomb : recombinations_) {
+    assert(recomb->mut_type() == RECOMB);
+    nb_mut_[recomb->mut_type()]++;
+  }
+#endif
+
 }
 
 void DnaReplicationReport::write_to_tree_file(gzFile tree_file) const {