Merge pull request #214 from bkille/fix-kmer-complexity

Average kmer complexities over chain in mergeMappingsInRange

src/map/include/computeMap.hpp

@@ -1655,6 +1655,11 @@ namespace skch
  [](double x, MappingResult &e){ return x + e.nucIdentity; })
  ) / it->n_merged; // this would scale directly by the number of mappings in the chain
 
+ //Mean identity of all kmer complexities 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; });