24小时热门版块排行榜    

查看: 1879  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wardwaads

银虫 (小有名气)

[求助] 【请教大家】blast序列比对问题 已有1人参与

请教大家一个问题
我做的高通量测序的结果其中的一条序列拿出来blast比对,结果发现query序列从第22个开始,对应着subject序列(16S rRNA gene, complete sequence)第1个碱基,也就是说query序列前面的碱基没有参与比对,但是blast还是认为这两个序列一致性更高,这是否是blast算法的一种参数控制导致的?
用AlignX对比两个序列发现,这种情况是由于query序列在前面比subject序列多21个碱基造成的,之后的序列完全相同。我困惑的是query序列按理应该就是subject序列对应的微生物,但是为什么会比subject的16S全序列前面还多21个碱基?
请大家帮我分析一下,谢谢!

【请教大家】blast序列比对问题
1.blast比对结果.jpg


【请教大家】blast序列比对问题-1
2.alignX比对结果.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wardwaads

银虫 (小有名气)

引用回帖:
3楼: Originally posted by acvhreinlg at 2015-08-31 21:41:27
你用你前面多出来的序列blast看看,估计你所谓的subject序列找的不对,应该有完全匹配的

另外,你的回复很有启发性
5楼2015-09-01 09:55:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

wardwaads

银虫 (小有名气)

大家难道都没有遇到这种情况吗?
2楼2015-08-31 21:17:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acvhreinlg

银虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
西门吹雪170: 金币+1, 鼓励回帖交流 2015-08-31 23:11:17
你用你前面多出来的序列blast看看,估计你所谓的subject序列找的不对,应该有完全匹配的

发自小木虫Android客户端
科研是将金钱转换为知识的过程,而创新则是将知识转换为金钱的过程。
3楼2015-08-31 21:41:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wardwaads

银虫 (小有名气)

引用回帖:
3楼: Originally posted by acvhreinlg at 2015-08-31 21:41:27
你用你前面多出来的序列blast看看,估计你所谓的subject序列找的不对,应该有完全匹配的

这个subject已经是数据库中这种菌目前最长的序列了
总长度是1480,写的是16S rRNA complete sequence
而我们所知的细菌16S rRNA全长为1542nt
是否有可能这里所谓16 rRNA complete sequence并不完整?
如果是这样的话,我用该序列比对该菌的全基因组DNA序列,是否应该就能比对的上?
4楼2015-09-01 09:54:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 274求调剂 +7 顾九笙要谦虚 2026-03-24 7/350 2026-03-25 15:18 by baoball
[考研] 求调剂 +3 李李不服输 2026-03-25 3/150 2026-03-25 13:03 by cmz0325
[考研] 化学调剂 +6 yzysaa 2026-03-21 6/300 2026-03-25 09:27 by aa331100
[考研] 材料学求调剂 +6 Stella_Yao 2026-03-20 6/300 2026-03-25 00:37 by baoball
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +15 吃吃吃才有意义 2026-03-19 16/800 2026-03-24 19:51 by 了了了了。。
[考研] 资源与环境 调剂申请(333分) +7 holy J 2026-03-21 7/350 2026-03-24 17:24 by xiaohai104
[考研] 材料292调剂 +8 橘颂思美人 2026-03-23 8/400 2026-03-24 16:33 by laoshidan
[考研] 321求调剂 +4 Ymlll 2026-03-24 4/200 2026-03-24 14:44 by sprinining
[考研] 305分求调剂(食品工程) +5 Sxy112 2026-03-21 7/350 2026-03-24 12:27 by 544594351
[考研] 335分 | 材料与化工专硕 | GPA 4.07 | 有科研经历 +4 cccchenso 2026-03-23 4/200 2026-03-23 23:00 by 徐ckkk
[考研] 一志愿陕师大生物学071000,298分,求调剂 +3 SYA! 2026-03-23 3/150 2026-03-23 19:09 by macy2011
[考研] 求调剂 +5 Zhangbod 2026-03-21 7/350 2026-03-22 13:13 by Zhangbod
[考研] 303求调剂 +5 安忆灵 2026-03-22 6/300 2026-03-22 12:46 by 素颜倾城1988
[考研] 085600材料与化工306 +4 z1z2z3879 2026-03-21 4/200 2026-03-21 23:44 by ms629
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 265求调剂 +12 梁梁校校 2026-03-19 14/700 2026-03-21 13:38 by lature00
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
信息提示
请填处理意见