24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3484  |  回复: 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 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见