应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 生物医药区 » 分子生物 » 生物信息 » 【请教大家】blast序列比对问题
51/1返回列表
查看: 1880  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-08-30 17:44:23
已阅   回复此楼   关注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的回帖
查看全部 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客户端
赞一下(1人) 回复此楼
科研是将金钱转换为知识的过程，而创新则是将知识转换为金钱的过程。
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序列找的不对，应该有完全匹配的

另外，你的回复很有启发性
赞一下 回复此楼
5楼2015-09-01 09:55:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 309求调剂 +4 gajsj 2026-03-25 5/250 2026-03-26 00:27 by Dyhoer
[考研] 086000生物与医药292求调剂 +4 小小陈小小 2026-03-22 7/350 2026-03-25 19:07 by 星空星月
[考研] 284求调剂 +15 Zhao anqi 2026-03-22 15/750 2026-03-25 12:51 by wht0531
[考研] 食品专硕 一志愿双一流 328 +3 xiaom99 2026-03-21 4/200 2026-03-24 21:20 by lailaisimei
[考研] 307求调剂 +5 超级伊昂大王 2026-03-24 5/250 2026-03-24 15:46 by 星空星月
[考研] 277分求调剂，跨调材料 +3 考研调剂lxh 2026-03-24 3/150 2026-03-24 13:52 by JourneyLucky
[考研] 【双一流院校新能源、环境材料，材料加工与模拟招收大量调剂】 +4 Higraduate 2026-03-22 7/350 2026-03-24 11:23 by 种大树
[考研] 341求调剂(一志愿湖南大学070300) +5 番茄头--- 2026-03-22 6/300 2026-03-23 23:45 by Txy@872106
[考研] 一志愿武理材料工程348求调剂 +6 ￣^￣゜汗 2026-03-19 9/450 2026-03-23 19:53 by pswait
[考研] 350求调剂 +6 weudhdk 2026-03-19 6/300 2026-03-23 15:47 by tangyuan0840221
[考研] 328求调剂 +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
[考研] 石河子大学（211、双一流）硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 材料与化工085600，总分304，本科有两篇sci参与，求调剂 +4 幸运的酱酱 2026-03-22 5/250 2026-03-22 20:15 by edmund7
[考研] 一志愿中南化学（0703）总分337求调剂 +9 niko- 2026-03-19 10/500 2026-03-22 16:08 by ColorlessPI
[考研] 275求调剂 +6 shansx 2026-03-22 8/400 2026-03-22 15:27 by barlinike
[考研] 303求调剂 +5 安忆灵 2026-03-22 6/300 2026-03-22 12:46 by 素颜倾城1988
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 261求B区调剂，科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
信息提示
关闭
请填处理意见
关闭