如序列长度 GC $ seqkit fx2tab hairpin.fa.gz -l -g -n -i -H | head -n 4 | csvtk -t -C '&' pretty#name seq qual length GCcel-let-7 99 43.43cel-lin-4 94 54.26cel-mir-1 96 40.62#两种形式转换zcat hairpin.fa.gz | seqkit fx2tab | seqkit tab2fx#按照长度排列序列seqkit sort ...
如序列长度 GC $ seqkit fx2tab hairpin.fa.gz-l-g-n-i-H|head-n4|csvtk-t-C'&'pretty#name seq qual length GCcel-let-79943.43cel-lin-49454.26cel-mir-19640.62#两种形式转换zcat hairpin.fa.gz|seqkit fx2tab|seqkit tab2fx#按照长度排列序列seqkit sort-l hairpin.fa.gz#得到前1000条readsseqkit fx...
统计每个序列的GC情况seqkit fx2tab --name --only-id --gc *.gz ==根据id号提取序列== ### 随机创造个id列表 seqkit sample-p0.001duplicated-reads.fq.gz|seqkit seq--name--only-id>id.txt ###根据id列表提取序列 seqkit grep--pattern-file id.txt duplicated-reads.fq.gz 或者-f参数 ### 提取...
如序列长度 GC $ seqkit fx2tab hairpin.fa.gz -l -g -n -i -H | head -n 4 | csvtk -t -C '&' pretty#name seq qual length GCcel-let-7 99 43.43cel-lin-4 94 54.26cel-mir-1 96 40.62#两种形式转换zcat hairpin.fa.gz | seqkit fx2tab | seqkit tab2fx#按照长度排列序列seqkit sort ...