MUMmer使用介绍

楼主  收藏   举报   帖子创建时间:  2018-03-06 00:00 回复:0 关注量:104

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/