Commit 7d70df1d authored by Mikaël Salson's avatar Mikaël Salson
Browse files

KmerAffect: We should not have same affectation with different lengths

At least for now this should not happen since we index all a germline with the same seed
parent 22489f9f
......@@ -22,6 +22,7 @@
*/
#include "kmeraffect.h"
#include <cassert>
int affect_strand(const affect_t &affect) {
return (affect.c & (1 << 7)) ? 1 : -1;
......@@ -120,8 +121,10 @@ KmerAffect &KmerAffect::operator+=(const KmerAffect &kmer) {
// strand is ambiguous but not the label, but we don't have enough space
// in 1 byte…
*this = AFFECT_AMBIGUOUS;
else
else {
assert(affect.c != kmer.affect.c || getLength() == kmer.getLength());
*this = AFFECT_AMBIGUOUS;
}
}
return *this;
}
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment