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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 生物医药区 » 微生物 » 综合其他 » 微生物基因组比对分析求助
101/1返回列表
查看: 3327  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

weinerchen

新虫 (初入文坛)

[求助] 微生物基因组比对分析求助 已有2人参与

我的一个细菌基因组测序得到精细图（scaffolds）,我想将这些scaffolds和该细菌的近缘种进行序列比对，对scaffold进行排序，初步确定基因组的序列顺利。
请问一下，我该用什么软件进行高效比对？即一个scaffold和近缘种基因组序列进行比对或全部scaffold和基因组序列进行比对。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-11-02 20:25:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daqwu

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
scelab(金币+5): 高手啊，谢谢交流 2011-11-04 15:47:32
普通的blast就可以，把你要比的基因组建一个database,以你要比的scaffold作为query序列比对就行，比对结果，建议加一个参数-m 8，这样出的结果是excel格式，方便看一些。
赞一下(2人) 回复此楼
2楼2011-11-04 10:40:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scelab

荣誉版主 (文坛精英)

小木虫之有关部门负责人
引用回帖:
2楼: Originally posted by daqwu at 2011-11-04 10:40:54:
普通的blast就可以，把你要比的基因组建一个database,以你要比的scaffold作为query序列比对就行，比对结果，建议加一个参数-m 8，这样出的结果是excel格式，方便看一些。

请问参数-m8加在哪里？
能否贴个图上来？
赞一下 回复此楼
小木虫之有关部门负责人
3楼2011-11-04 15:52:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daqwu

木虫 (小有名气)

【答案】应助回帖

★ ★ ★
scelab(金币+3): :sweat:不懂程序。。。不过还是谢谢你 2011-11-04 21:05:34
引用回帖:
3楼: Originally posted by scelab at 2011-11-04 15:52:03:
请问参数-m8加在哪里？
能否贴个图上来？

在linux系统里，可以写下面的命令实现
blastall -p blastp -d db/nr/nr -i *.fasta -m 8 -o *.out ;
在本机windows系统里，一些序列处理软件会有参数输入栏，在那输-m 8，也可以
赞一下(1人) 回复此楼
4楼2011-11-04 20:37:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weinerchen

新虫 (初入文坛)
内容已删除
回复此楼
5楼2012-11-03 20:34:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿南hmily

新虫 (小有名气)

【答案】应助回帖

引用回帖:
5楼: Originally posted by weinerchen at 2012-11-03 20:34:00
各位大侠，有用过MAUVE软件吗？目前正在使用MAUVE进行基因组的比对工作。我获得的序列是一株细菌的精细图，含有55个contigs（大片段）。我提交fasta格式的两个基因组序列文件（目的基因组1.fasta和参考基因组2.fast ...

弱弱的问一下，楼主关于mauve的软件问题解决了吗？
我最近在做基因组比对，安装了基因组比对软件mauve，也安装了java，可是还是运行不了，依次出现如下的界面，求指导！！！或者能给我个这个软件的安装包吗？
微生物基因组比对分析求助
打开界面.jpg


微生物基因组比对分析求助-1
运行界面.jpg
赞一下 回复此楼
6楼2014-12-24 12:50:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

张、小杨

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by daqwu at 2011-11-04 10:40:54
普通的blast就可以，把你要比的基因组建一个database,以你要比的scaffold作为query序列比对就行，比对结果，建议加一个参数-m 8，这样出的结果是excel格式，方便看一些。

你好 我是全基因组的新手，现在从0开始分析数据 什么都不懂，，，感觉你是大神啊~能推荐几个软件吗？ 我刚拿到了一个细菌的全基因组的完成图，想找到一对引物 做该菌测定量pcr 万分感谢了~
赞一下 回复此楼
请叫我张有才、
7楼2016-05-10 16:13:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xp198766

铁杆木虫 (著名写手)

小木虫职业打酱油滴~~！

【答案】应助回帖

我觉得看亲缘关系你可以直接选择16S来进行比对，看一下相似度。
另外，如果你想做同源性，可以找出你的物种与近缘物种的直系同源基因，再看直系同源基因的相似度，或者使用这些基因构建联合基因树
赞一下 回复此楼
8楼2016-05-13 12:47:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xp198766

铁杆木虫 (著名写手)

小木虫职业打酱油滴~~！
引用回帖:
3楼: Originally posted by scelab at 2011-11-04 15:52:03
请问参数-m8加在哪里？
能否贴个图上来？...

如果是新的blast软件的话，参数都变了
变成:
blastn -query your.file -db database -evalue 1e-5 -out out.file -outfmt 6

这里面的-outfmt 6就是老版本blast里面的-m 8
赞一下 回复此楼
9楼2016-05-13 12:51:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xp198766

铁杆木虫 (著名写手)

小木虫职业打酱油滴~~！
我好像挖坟了……我错了……
赞一下 回复此楼
10楼2016-05-13 12:52:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 weinerchen 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 课题组会帮你通过考核么？ +3 好运来～9527 2026-10-10 6/300 2026-10-12 14:55 by 简单化xn
[有机交流] OpenAI又攻克数学难题了，离AI做有机合成还有多远？ +4 wangyl_123 2026-10-10 7/350 2026-10-12 14:47 by chuans0923
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷！
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬，路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[基金申请] 入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍 +5 13108017953 2026-10-08 5/250 2026-10-10 11:15 by wdwd123
[有机交流] 我应该选哪条合成路线让实验成功率高一些？ +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
关闭
请填处理意见
关闭