seqtk comp in.fa > out.fa 4.subseq 根据name.list(不带>符号)提取子序列 -l可设定输出的每行长度 seqtk subseq -l 40 B1_NR100nl.fasta name.list > out.fa 5.随机抽取序列,按照比例或者数量 -s设定随机种子,便于重复 seqtk sample -s 10 test.fq 0.4 #比例seqtk sample -s 10 test.fq 100 #...
可以观察到第一个参数是源文件,第二个参数是对应键名文件,我们根据name.list去提取文件. seqtk subseq genome.fa name.list | less -N 我们可以改变name.list的文件内容,让subseq提取不同位置的碱基.代码保持不变,获得的碱基不同了. 这里有一个点需要科普.我们在文件里列的是1~10,但最后展示出...
We read every piece of feedback, and take your input very seriously. Include my email address so I can be contacted Cancel Submit feedback Saved searches Use saved searches to filter your results more quickly Cancel Create saved search Sign in Sign up Reseting focus {...
int stk_subseq(int argc, char *argv[]) { khash_t(reg) *h = kh_init(reg); gzFile fp; kseq_t *seq; int l, i, j, c, is_tab = 0, line = 1024; khint_t k; while ((c = getopt(argc, argv, "tl:")) >= 0) {