24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3749  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 290求调剂085701 +7 1314捧花 2026-04-02 7/350 2026-04-04 23:33 by lqwchd
[考研] 一志愿郑州大学材料与化工085600,求调剂 +24 吃的不少 2026-04-02 24/1200 2026-04-04 23:20 by 永字号
[考研] 081700,311,求调剂 +11 冬十三 2026-04-04 11/550 2026-04-04 22:24 by 啵啵啵0119
[考研] 材料调剂 +15 一样YWY 2026-04-01 15/750 2026-04-04 22:23 by hemengdong
[考研] 0854求调剂 +4 assdll 2026-04-03 4/200 2026-04-04 22:17 by hemengdong
[考研] 278求调剂 +14 范婷娜 2026-04-04 15/750 2026-04-04 22:15 by lqwchd
[考研] 一志愿华北电力大学(北京),材料科学与工程学硕265,求调剂 +11 yelck 2026-04-03 12/600 2026-04-04 19:52 by dongzh2009
[考研] 320求调剂 +5 振—TZ 2026-04-02 5/250 2026-04-03 14:42 by fxue1114
[考研] 建环,能源,土木老师路过看一看!!! +5 嘿嘿uu 2026-04-01 5/250 2026-04-03 11:47 by znian
[考研] 调剂 +3 osbbx 2026-04-02 3/150 2026-04-03 07:47 by cc8418
[考研] 材料340分调剂 +7 夏夜晚风_long 2026-04-02 9/450 2026-04-02 21:20 by dongzh2009
[考研] 318求调剂,计算材料方向 +10 吸喵有害笙命 2026-04-01 11/550 2026-04-02 16:29 by oooqiao
[考研] 26考研调剂 +4 Wnz.20030617 2026-04-01 5/250 2026-04-02 16:11 by 1939136013狗壮
[考研] 0832食品科学与工程学硕282调剂 +4 鱼在水中游a 2026-04-02 7/350 2026-04-02 14:12 by baoball
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +10 哇呼哼呼哼 2026-04-01 11/550 2026-04-02 11:31 by lnilvy
[考研] 求调剂,一志愿 南京航空航天大学 ,080500材料科学与工程学硕,总分289分 +11 @taotao 2026-03-29 11/550 2026-04-02 10:04 by realme321
[考研] 298求调剂 +4 什么是胖头鱼 2026-03-30 6/300 2026-04-01 22:06 by 客尔美德
[考研] 材料与化工(0856)304求B区调剂 +8 邱gl 2026-03-30 16/800 2026-04-01 17:58 by 邱gl
[考研] 一志愿北交材料工程总分358 +5 cs0106 2026-04-01 7/350 2026-04-01 11:45 by wangjy2002
[考研] 269求调剂 +4 我想读研11 2026-03-31 4/200 2026-03-31 10:04 by cal0306
信息提示
请填处理意见