workLog: seqkit 安装与使用

本文介绍了seqkit的安装过程,通过conda进行便捷安装,并展示了使用seqkit对NGS序列进行统计和翻译的操作。在尝试翻译DNA到氨基酸时,注意到需要调整阅读框,并且由于原始文件包含正反两种序列,需要进一步处理。
摘要由CSDN通过智能技术生成

处理
https://bioinf.shenwei.me/seqkit/

安装

$ conda install -c bioconda seqkit

conda 还真是方便啊……

统计了一下上次NGS拼接后截出来的300个seq

$ seqkit stats testM300.fa
file         format  type  num_seqs  sum_len  min_len  avg_len  max_len
testM300.fa  FASTA   DNA        300  127,344      366    424.5      490

试一下NGS的整个文件,也很快出结果

$ seqkit stats outM.fa
file     format  type    num_seqs        sum_len  min_len  avg_len  max_len
outM.fa  FASTA   DNA   10,760,889  4,571,253,828      250    424.8      4
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值