24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3741  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工(0856)304求B区调剂 +8 邱gl 2026-03-30 14/700 2026-04-01 16:41 by Ccclqqq
[考研] 085600 一志愿9 总分351 求调剂学校 +7 czhcz 2026-03-31 8/400 2026-04-01 16:35 by czhcz
[考研] 一志愿085600中科院宁波所276分求调剂 +21 材料学257求调剂 2026-03-28 22/1100 2026-04-01 15:06 by oooqiao
[考研] 288求调剂 一志愿哈工大 材料与化工 +26 洛神哥哥 2026-03-31 26/1300 2026-04-01 12:54 by 山水有情
[考研] 材料专业求调剂 +7 月月鸟木 2026-04-01 7/350 2026-04-01 11:36 by wangjy2002
[考研] 070300一志愿211,312分求调剂院校 +13 小黄鸭宝 2026-03-30 13/650 2026-04-01 10:59 by oooqiao
[考研] 267求调剂 +13 uiybh 2026-03-31 13/650 2026-04-01 10:25 by 探123
[考研] 0856材料与化工调剂,339 +10 10213207 2026-03-31 10/500 2026-04-01 08:28 by Linda Hu
[考研] 求调剂:一志愿:南京大学 专业:0705 总分320 ,本科985,四六级已过 +3 lfy760306 2026-03-31 3/150 2026-04-01 01:57 by Creta
[考研] 318求调剂 +3 笃行致远. 2026-03-31 3/150 2026-03-31 20:27 by 求调剂zz
[考研] 一志愿a区211,085601-307分求调剂 +10 党嘉豪 2026-03-31 23/1150 2026-03-31 18:44 by JourneyLucky
[考研] 299求调剂 +8 嗯嗯嗯嗯2 2026-03-27 8/400 2026-03-31 18:23 by lizhi8172
[考研] 物理学调剂 +4 小羊36 2026-03-30 4/200 2026-03-31 16:16 by lishahe
[考研] 282求调剂 不挑专业 求收留 +4 Yam. 2026-03-30 5/250 2026-03-31 14:41 by 王亮_大连医科大
[考研] 085600材料与化工调剂 +16 kikiki7 2026-03-30 16/800 2026-03-31 10:03 by 氯化亚硝酰
[考研] 抱歉 +4 田洪有 2026-03-30 4/200 2026-03-30 21:26 by mumin1990
[考研] 085602 化学工程专硕 340分求调剂 +4 qianbai11 2026-03-29 4/200 2026-03-30 11:34 by 唐沐儿
[考研] 295求调剂 +5 1428151015 2026-03-27 6/300 2026-03-28 04:04 by fmesaito
[考研] 352分 化工与材料 +5 海纳百川Ly 2026-03-27 5/250 2026-03-28 03:39 by fmesaito
[考研] 调剂 +4 柚柚yoyo 2026-03-26 4/200 2026-03-26 20:43 by fmesaito
信息提示
请填处理意见