We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
2 parents 1a405d2 + 9afc465 commit cc49b9fCopy full SHA for cc49b9f
src/map/include/computeMap.hpp
@@ -1638,6 +1638,12 @@ namespace skch
1638
) / it->n_merged; // this would scale directly by the number of mappings in the chain
1639
1640
1641
+ //Mean kmer complexity of all mappings in the chain
1642
+ it->kmerComplexity = ( std::accumulate(
1643
+ it, it_end, 0.0,
1644
+ [](double x, MappingResult &e){ return x + e.kmerComplexity; })
1645
+ ) / it->n_merged; // this would scale directly by the number of mappings in the chain
1646
+
1647
//Discard other mappings of this chain
1648
std::for_each( std::next(it), it_end, [&](MappingResult &e){ e.discard = 1; });
1649
0 commit comments