基因丰度统计

基因丰度计算

使用soap将clean data比对到非冗余的核酸序列NonRundant.gene.nucl.fasta。根据soap比对结果,获得比对到某基因的read数。Soap比对工具的使用在之前的微信推送稿中介绍过,这里就不做说明。然后,基于read数获得基因的丰度值,计算方法如下:

S1

S2

S3

G1

.

.

.

G2

.

.

.

G3

.

.

.

样本i中基因j表示为SiGj, 样本i中比对到基因j的read数表示为SiGj_read,基因j长度表示为Gj_len。

相对丰度:SiGj丰度值=(SiGj_read/Gj_len)/(SiG1_read/G1_len+SiG2_read/G2_len+SiG3_read/G3_len ...)

参考:宏基因组分析——基因预测篇

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值