帮助文档
https://bioinf.shenwei.me/seqkit/
安装我直接使用conda安装了
conda install seqkit
统计fastq文件的信息
seqkit stats SRR6236885.sra.fastq -a
选取部分fastq
按比例
seqkit sample SRR6236885.sra.fastq -p 0.01 -o sample.fastq
按数量
seqkit sample SRR6236885.sra.fastq -n 100 -o sample-1.fastq
但是这里不是精确的输入你自己指定的数字，这个帮助文档也提到了
不去看具体原因了先
fastq 转换fasta
seqkit fq2fa sample-1.fastq -o sample-1.fasta
根据id 提取序列 fasta
fastq好像也可以
seqkit grep sample-1.fasta -f id.list
id只是一部分好像也可以
比如我这里fasta文件 的完整id是 SRR6236885.sra.9047 3a19e708-9d65-4f29-a332-d2e9a2b39234_Basecall_Alignment length=2375
id.list文件里的id是SRR6236885.sra.9047
这样也可以的
欢迎大家关注我的公众号
小明的数据分析笔记本
小明的数据分析笔记本 公众号 主要分享：1、R语言和python做数据分析和数据可视化的简单小例子；2、园艺植物相关转录组学、基因组学、群体遗传学文献阅读笔记；3、生物信息学入门学习资料及自己的学习笔记！