楼主 | 收藏 | 举报 2018-03-06 00:00   浏览:104   回复:0

MUMmer使用介绍

MUMmer是一个快速比对两个基因组序列的软件。在一个普通的linux机器上,比对两个3~4Mb的基因组所需要的时间只要30秒不到。其余还有很多功能将在下面详细说明。

MUMmer的主页为http://www.tigr.org/software/mummer/,只要填写好一个licenes,就可以免费下载该软件。安装是将下载的MUMmer2.13.tar.gz解压缩,进入解压缩目录后,键入make install,就完成安装,如果有报错信息,则根据报错信息对Makefile进行修改。

根据不同的情况,MUMmer有不同的程序进行比对,下面将一一说明

run-mummer1.csh

用法 run-mummer1.csh <seq 1> <seq 2> <tag>

seq1 是一条fasta格式的序列

seq2 是另一条fasta格式的序列

tag 是输出的文件名,输出有tag.out tag.align tag.gaps

该程序是最适合用在两个高度相似的基因组,并且没有很多的重排的片断

run-mummer2.csh

用法 run-mummer2.csh <reference> <query> <prefix>

refercence 是参考序列

query 是比对序列

prefix 是输出名字

输出文件有

  • <prefix>.out 列出比对上的位置和长度
  • <prefix>.gaps 列出了比对结果中的gap分布
  • <prefix>.errorsgaps 列出了上面gap分布结果中可能错误的gap
  • <prefix>.align 和.gaps文件相似,多了gap前后的位置和gap长度的比率

该程序是在mummer1的基础上改进了算法,增加了比对上的序列聚类的功能,适合有一定重排情况的序列之间的比对

nucmer

用法: nucmer [options] <reference> <query>

reference 参考序列

query 比对序列

options 参数太多,无法在这里一一说明,可以参考MUMmer目录下的doc文档

该程序最大的特点是支持多序列对多序列比对,适合差异较大的序列之间的比对

promer

用法: promer [options] <reference> <query>

reference 参考序列

query 比对序列

options 参数太多,无法在这里一一说明,可以参考MUMmer目录下的doc文档

该程序是多序列对多序列进行6个相位比对

运行下面两步,直接运行即可,路径已经设好

mummerplot

mummerplot <prefix>.out -p <plot prefix>"

gnuplot <plot prefix>.gp"

需要图示结果的话,需要安装gnuplot

原文来自:http://lwzhanghz.blog.163.com/blog/static/1368263692013128104648/

打赏
网站首页 | 关于我们 | 联系方式 | 使用协议 | 版权隐私 | 网站地图  |  排名推广  |  广告服务  |  积分换礼  |  网站留言  |  RSS订阅  |  违规举报
 
免责声明:本站有部分内容来自互联网,如无意中侵犯了某个媒体 、公司 、企业或个人等的知识产权,请来电或致函告之,本网站将在规定时间内给予删除等相关处理。