Return the k-mer coverage histogram. */
| 10 | |
| 11 | /** Return the k-mer coverage histogram. */ |
| 12 | static inline |
| 13 | Histogram coverageHistogram(const SequenceCollectionHash& c) |
| 14 | { |
| 15 | typedef SequenceCollectionHash Graph; |
| 16 | |
| 17 | Histogram h; |
| 18 | for (Graph::const_iterator it = c.begin(); |
| 19 | it != c.end(); ++it) { |
| 20 | if (it->second.deleted()) |
| 21 | continue; |
| 22 | h.insert(it->second.getMultiplicity()); |
| 23 | } |
| 24 | return h; |
| 25 | } |
| 26 | |
| 27 | /** Calculate a k-mer coverage threshold from the given k-mer coverage |
| 28 | * histogram. */ |
no test coverage detected