INNER CODE UNIT · C
beg
lh3/seqtk · seqtk.c:401
beg = left; end = seq->seq.l - right;
if (beg >= end) beg = end = 0;
if (fixed_len > 0 && end - beg > fixed_len) end = beg + fixed_len;
} else if (seq->qual.l > min_len) {
for (i = 0, beg = tmp = 0, end = seq->qual.l, s = max = 0.; i < seq->qual.l; ++i) {
int q = seq->qual.s[i];
if (q < 36) q = 36;
if (q > 127) q = 127;
s += param - q_int2real[q];
if (s > max) max = s, beg = tmp, end = i + 1;
if (s < 0) s = 0, tmp = i + 1;
}
/* max never set; all low qual, just give first min_len bp */
if (max == 0.) beg = 0, end = min_len;
if (end - beg < min_len) { // window-based
int is, imax;