INNER CODE UNIT · C

b

lh3/seqtk · seqtk.c:1021

			int b = seq->seq.s[i];
			if (bitcnt_table[seq_nt16_table[b]] == 2)
				printf("%s\t%ld\t%c\n", seq->name.s, i+1, b);
		}
	}
	kseq_destroy(seq);
	gzclose(fp);
	return 0;
}

int stk_split(int argc, char *argv[])
{
	gzFile fp;
	kseq_t *seq;
	int64_t l;
	int c, i, n = 10, len = 0;
	char *prefix, *fn;
	FILE **out;

View source record →