24小时热门版块排行榜    

查看: 3705  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 311求调剂 +5 冬十三 2026-03-15 5/250 2026-03-15 18:38 by 无际的草原
[考研] 0703化学调剂 ,六级已过,有科研经历 +4 曦熙兮 2026-03-15 4/200 2026-03-15 18:01 by JourneyLucky
[考研] 290求调剂 +3 孔志浩 2026-03-12 8/400 2026-03-15 15:30 by 孔志浩
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 267一志愿南京工业大学0817化工求调剂 +5 SUICHILD 2026-03-12 5/250 2026-03-14 14:53 by jean5056
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 306求调剂 +4 唐薏薏 2026-03-09 4/200 2026-03-14 01:19 by JourneyLucky
[考研] 环境调剂 +6 晓看天暮看云 2026-03-09 6/300 2026-03-14 01:16 by JourneyLucky
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 一志愿中科院,化学方向,295求调剂 +4 一氧二氮 2026-03-11 4/200 2026-03-13 22:35 by JourneyLucky
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 26考研求调剂 +5 丶宏Sir 2026-03-13 5/250 2026-03-13 13:05 by JourneyLucky
[考研] 304求调剂(085602一志愿985) +12 化工人999 2026-03-09 12/600 2026-03-13 12:02 by JourneyLucky
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考研] 279求调剂 +3 莫xiao 2026-03-10 4/200 2026-03-11 08:06 by 斩魂滴兔子!
[考研] 家人们 调剂不迷路 看这里 +8 likeihood 2026-03-09 13/650 2026-03-10 08:09 by likeihood
信息提示
请填处理意见