1.将fastq 文件转换成fasta 文件 seqtk seq -A input.fastq > output.fasta 2.得到反向互补序列 seqtk seq -Ar input.fastq > output.fasta 3.seqtk comp: 得到fastq/fasta 文件的碱基组成 (输出格式:序列id 序列长度 A C G T ) seqtk comp in.fa > out.fa 4.subseq 根据name.list(不带>符号)提取...
seqtk seq -Ar input.fastq > output.fasta 3.seqtk comp: 得到fastq/fasta 文件的碱基组成(输出格式:序列id 序列长度 A C G T )seqtk comp in.fa > out.fa4.subseq 根据name.list(不带>符号)提取子序列 -l可设定输出的每行长度seqtk subseq -l 40 B1_NR100nl.fasta name.list > out.fa...
seqtk seq -A input.fastq > output.fasta input.fastq的内容: @NB001 ATGCACAAAACCCC+/// output.fasta 的内容: >NB001 ATGCACAAAACCCC 2)得到反向互补序列 seqtk seq -Ar input.fastq > output.fasta output.fasta的内容为: >NB001 GGGGTTTTGTGCAT seqtk comp: 得到fastq/fasta 文件的碱基组成 seqtk c...
3.2、seqtk comp $seqtk comp Usage: seqtk comp [-u] [-r in.bed] <in.fa> # 获取fa的碱基的组成信息,用-r参数可以输出bed中的给定区间的序列 Output format: chr, length, #A, #C, #G, #T, #2, #3, #4, #CpG, #tv, #ts, #CpG-ts # 输出格式:序列号 序列长度 A C G T ...
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 {...
Options: -t TAB delimited output# 输出以tab分割 -l INT sequence line length [0]# 输出序列以长度INT换行 Note: Use 'samtools faidx' if only a few regions are intended.#注意:如果只有少数几个区域,请使用'samtools faidx' $seqtk mergepe ...