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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 生物医药区 » 分子生物 » 生物信息 » 【请教大家】blast序列比对问题
51/1返回列表
查看: 1878  |  回复: 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

银虫 (小有名气)
大家难道都没有遇到这种情况吗？
赞一下 回复此楼
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序列找的不对，应该有完全匹配的

这个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的回帖

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的回帖
相关版块跳转 我要订阅楼主 wardwaads 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 【2026考研调剂】制药工程 284分 求相关专业调剂名额 +4 袁奂奂 2026-03-25 8/400 2026-03-25 14:32 by lbsjt
[考研] 334分 一志愿武理-080500 材料求调剂 +3 李李不服输 2026-03-25 3/150 2026-03-25 13:38 by cocolv
[考研] 081700 调剂 267分 +10 迷人的哈哈 2026-03-23 10/500 2026-03-25 12:12 by userper
[考研] 生物技术与工程 +3 1294608413 2026-03-25 3/150 2026-03-25 10:45 by hxsm
[考研] 287求调剂 +10 晨昏线与星海 2026-03-19 11/550 2026-03-25 10:35 by userper
[考研] 一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。 +4 晨昏线与星海 2026-03-20 4/200 2026-03-25 10:16 by allen-yin
[考研] 340求调剂 +5 话梅糖111 2026-03-24 5/250 2026-03-25 06:53 by ilovexiaobin
[考研] 材料考研调剂生 +3 黄粱一梦千年 2026-03-24 3/150 2026-03-24 17:00 by barlinike
[考研] 291求调剂 +3 HanBeiNingZC 2026-03-24 3/150 2026-03-24 16:34 by barlinike
[材料工程] 一志愿C9材料与化工专业总分300求调剂 +4 曼111 2026-03-24 5/250 2026-03-24 15:44 by 星空星月
[考研] 321求调剂 +4 Ymlll 2026-03-24 4/200 2026-03-24 14:44 by sprinining
[考研] 一志愿吉大化学322求调剂 +4 17501029541 2026-03-23 6/300 2026-03-24 10:21 by 戴围脖的小蚊子
[考研] 求调剂 +7 十三加油 2026-03-21 7/350 2026-03-23 23:48 by 热情沙漠
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +5 葵梓卫队 2026-03-18 7/350 2026-03-23 16:26 by lingjue
[考研] 260求调剂 +3 朱芷琳 2026-03-20 4/200 2026-03-22 15:12 by 朱芷琳
[考研] 材料求调剂 +5 @taotao 2026-03-21 5/250 2026-03-21 20:55 by lbsjt
[基金申请] 学校已经提交到NSFC，还能修改吗？ 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 一志愿南昌大学，327分，材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 288求调剂 +16 于海海海海 2026-03-19 16/800 2026-03-20 22:28 by JourneyLucky
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
信息提示
关闭
请填处理意见
关闭