24小时热门版块排行榜    

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

weinerchen

新虫 (初入文坛)

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

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

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注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的回帖
查看全部 10 个回答

daqwu

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
scelab(金币+5): 高手啊,谢谢交流 2011-11-04 15:47:32
普通的blast就可以,把你要比的基因组建一个database,以你要比的scaffold作为query序列比对就行,比对结果,建议加一个参数-m 8,这样出的结果是excel格式,方便看一些。
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,也可以
4楼2011-11-04 20:37:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见