INNER CODE UNIT · C++
veclen
soedinglab/plass · lib/flash/combine_reads.cpp:229
int veclen = *len_p & ~0xf;
const char *_seq_1 = seq_1 - veclen;
const char *_seq_2 = seq_2 - veclen;
const char *_qual_1 = qual_1 - veclen;
const char *_qual_2 = qual_2 - veclen;
unsigned _num_mismatches = 0;
unsigned _mismatch_qual_total = 0;
for (int i = 0; i < veclen; i++) {
if (_seq_1[i] != _seq_2[i]) {
_num_mismatches++;
_mismatch_qual_total += min(_qual_1[i], _qual_2[i]);
}
}
assert(num_mismatches == _num_mismatches);
assert(mismatch_qual_total == _mismatch_qual_total);
}
#endif