24小时热门版块排行榜    

查看: 4162  |  回复: 6

xiami1943

铁虫 (小有名气)


[交流] 【求助/交流】关于ncbi上比对结果的问题 已有5人参与

PCR了一个16s在ncbi上比对,出了一些参数,不明白什么意思,max score,totle score,query coverage,e value,max ident,是什么意义?为什么扩同样菌种(买来的)的16s,在ncbi上比对以后query coverage值也不是100%呢?有哪位有经验的大虾指点我一下吧,不胜感激(还指望能筛到一个新菌毕业呢,16s比来比去成一样的了,/(ㄒoㄒ)/~~)
回复此楼

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rural2084

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
看天(金币+1):谢谢 交流 2010-07-14 13:00:03
NCBI上有的是16s的全序列,你测序得到的序列没有包含16s全序列两端的序列,所以query coverage不可能是100%,比对序列主要看query coverage,e value相似性,其他的我也不是很清楚
2楼2010-07-14 11:09:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cnlics

木虫 (小有名气)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
lstt09nk(金币+3):感谢交流~· 2010-07-14 12:46:33
只需关注query coverage,e value,max ident即可,你还要好好体味下gapped与ungapped、redundancy的含义。

gapped的意思是允许程序在比较序列的时候将序列割裂产生空缺进行比较,这样就可以找到差异较大的相似序列,高度相似的序列被认为在进化上是同源的;相反,ungapped要求在序列比较过程中保持query完整,不允许程序引入空缺即gap。

“冗余-redundancy”:数据库中有很多序列非常相似,有些序列之间只有个别碱基或氨基酸不同,甚至一模一样,这种状态就叫冗余。这导致blast或fasta查询的结果中,常常有很多冗余性的结果。有些nr库是non-redundant的,即非冗余的,是经过人工处理后,将没有明显意义的冗余序列简化了。实际上,生物医学数据库中,并没有严格意义上的nr库,都有一定的冗余度。

可以简单地理解为,e value值越小,比对的条件越严格,在这样严格的前提下,identity的百分比越高,(默认的从进化角度)那么序列之间同源关系越近或者说亲缘关系越密切。

澄清以上概念,你的结果可以很好理解。

既然你做菌种的进化研究,你应该明白,由于数量巨大,细菌的变异是极其常见的。即使单克隆来源的细菌,经过若干代的培养,基因组DNA上也会出现很多突变。

综合序列比对涉及的因素,加上生物学背景知识,你的16SRNA的结果是非常正常的。
3楼2010-07-14 11:28:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiami1943

铁虫 (小有名气)


引用回帖:
Originally posted by rural2084 at 2010-07-14 11:09:39:
NCBI上有的是16s的全序列,你测序得到的序列没有包含16s全序列两端的序列,所以query coverage不可能是100%,比对序列主要看query coverage,e value相似性,其他的我也不是很清楚

谢谢
4楼2010-07-15 10:06:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiami1943

铁虫 (小有名气)


引用回帖:
Originally posted by cnlics at 2010-07-14 11:28:41:
只需关注query coverage,e value,max ident即可,你还要好好体味下gapped与ungapped、redundancy的含义。

gapped的意思是允许程序在比较序列的时候将序列割裂产生空缺进行比较,这样就可以找到差异较大的相似序 ...

看着这么多字有点晕,咕~~(╯﹏╰)b,谢谢你回了这么详细的解释
5楼2010-07-15 10:08:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fusiminhz

金虫 (小有名气)

不错。祝福!也粘粘喜气。

不错。祝福!也粘粘喜气。
6楼2010-07-15 10:15:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

longlong166

新虫 (小有名气)

送红花一朵
学习了
7楼2013-05-06 01:15:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiami1943 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:39 by CgyNCDVNhGVg
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-15 16:26 by CgyNCDVNhGVg
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-15 15:35 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-15 15:27 by CgyNCDVNhGVg
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-15 15:14 by CgyNCDVNhGVg
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:34 by VMDqgnlKimtZ
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 14:25 by VMDqgnlKimtZ
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:13 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 09:38 by VMDqgnlKimtZ
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:34 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-15 06:14 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 21:37 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见