|
MUMmer是一个快速比对两个基因组序列的软件。在一个普通的linux机器上,比对两个3~4Mb的基因组所需要的时间只要30秒不到。其余还有很多功能将在下面详细说明。
根据不同的情况,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 <prefix>.out -p <plot prefix>"
gnuplot <plot prefix>.gp" 需要图示结果的话,需要安装gnuplot
上一篇:序列比对软件及blast 下一篇:第八章 多序列比对的实际应用
|