seqtk comp in.fa > out.fa4.subseq 根据name.list(不带>符号)提取子序列 -l可设定输出的每行长度seqtk subseq -l 40 B1_NR100nl.fasta name.list > out.fa5.随机抽取序列,按照比例或者数量 -s设定随机种子,便于重复seqtk sample -s 10 test.fq 0.4 #比例 seqtk sample -s 10 test.fq 100 #数量6...
可以观察到第一个参数是源文件,第二个参数是对应键名文件,我们根据name.list去提取文件. seqtk subseq genome.fa name.list | less -N 我们可以改变name.list的文件内容,让subseq提取不同位置的碱基.代码保持不变,获得的碱基不同了. 这里有一个点需要科普.我们在文件里列的是1~10,但最后展示出...
Usage: seqtk subseq [options] <in.fa> <in.bed>|<name.list> #提取name.list中指定名称的fa序列, Options: -t TAB delimited output# 输出以tab分割 -l INT sequence line length [0]# 输出序列以长度INT换行 Note: Use 'samtools faidx' if only a few regions are intended.#注意:如果只有少数几个...
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 {...