INNER CODE UNIT · C
min_l
lh3/seqtk · seqtk.c:814
min_l = seq[0]->seq.l < seq[1]->seq.l? seq[0]->seq.l : seq[1]->seq.l;
printf(">%s", seq[0]->name.s);
for (l = 0; l < min_l; ++l) {
c[0] = seq[0]->seq.s[l]; c[1] = seq[1]->seq.s[l];
if (seq[0]->qual.l && seq[0]->qual.s[l] - 33 < qual) c[0] = tolower(c[0]);
if (seq[1]->qual.l && seq[1]->qual.s[l] - 33 < qual) c[1] = tolower(c[1]);
if (is_intersect) is_upper = (isupper(c[0]) || isupper(c[1]))? 1 : 0;
else if (is_mask) is_upper = (isupper(c[0]) || isupper(c[1]))? 1 : 0;
else is_upper = (isupper(c[0]) && isupper(c[1]))? 1 : 0;
c[0] = seq_nt16_table[c[0]]; c[1] = seq_nt16_table[c[1]];
if (c[0] == 0) c[0] = 15;
if (c[1] == 0) c[1] = 15;
b[0] = bitcnt_table[c[0]];
b[1] = bitcnt_table[c[1]];
if (is_upper) {
if (b[0] == 1 && b[1] == 1) {
if (c[0] == c[1]) ++cnt[0];
else ++cnt[1];