INNER CODE UNIT · C

c

lh3/seqtk · seqtk.c:654

			c = i < len? seq_nt6_table[(uint8_t)seq->seq.s[i]] : 5;
			if (i == len || (c >= 1 && c <= 4)) {
				if (l > 0 && l >= min_size)
					printf("%s\t%ld\t%ld\n", seq->name.s, i - l, i);
				l = 0;
			} else ++l;
		}
	}
	kseq_destroy(seq);
	gzclose(fp);
	return 0;
}

/* subseq */

int stk_subseq(int argc, char *argv[])
{
	khash_t(reg) *h = kh_init(reg);

View source record →