Average kc over all segments in chain

marbl · Nov 30, 2023 · 83ab90f · 83ab90f
1 parent be6dbb6
commit 83ab90f
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/src/map/include/computeMap.hpp b/src/map/include/computeMap.hpp
@@ -1638,6 +1638,12 @@ namespace skch
                   ) / it->n_merged; // this would scale directly by the number of mappings in the chain
 
 
+              //Mean kmer complexity of all mappings in the chain
+              it->kmerComplexity = ( std::accumulate(
+                                      it, it_end, 0.0,
+                                      [](double x, MappingResult &e){ return x + e.kmerComplexity; })
+                  ) / it->n_merged; // this would scale directly by the number of mappings in the chain
+
               //Discard other mappings of this chain
               std::for_each( std::next(it), it_end, [&](MappingResult &e){ e.discard = 1; });