Seqkit说明文档
参见：https://bioinf.shenwei.me/seqkit/
常用参数记录
Read and print
# Only print seq (global flag -w defines the output line width, 0 for no wrap)
$ seqkit seq hairpin.fa.gz -s -w 0
Stats
$ seqkit stats *.f{a,q}.gz
file format type num_seqs sum_len min_len avg_len max_len
hairpin.fa.gz FASTA RNA 28,645 2,949,871 39 103 2,354
mature.fa.gz FASTA RNA 35,828 781,222 15 21.8 34
reads_1.fq.gz FASTQ DNA 2,500 567,516 226 227 229
reads_2.fq.gz FASTQ DNA 2,500 560,002 223 224 225
fx2tab
# Print sequence length, GC content, and only print names (no sequences), we could also print title line by flag -H
$ seqkit fx2tab hairpin.fa.gz -l -g -n -i -H | head -n 4 | csvtk -t -C '&' pretty
# name seq qual length GC
cel-let-7 99 43.43
cel-lin-4 94 54.26
cel-mir-1 96 40.62
