序列比对软件 MUMmer 简单上手(一)

MUMmer 4 是一个快速比对 DNA 或蛋白序列的软件,最常用的 program 就是 nucmerpromer,分别是比对核苷酸和氨基酸序列的,接下来以 nucmer 为例详细介绍。

安装

wget  https://github.com/mummer4/mummer/releases/download/v4.0.0rc1/mummer-4.0.0rc1.tar.gz

tar  -zxvf  mummer-4.0.0rc1.tar.gz

cd  mummer-4.0.0rc1

./configure  --prefix=/path/to/installation

make

make install

简单上手

取一个参考基因组 ref.fa,和另一个要与之比对的参考基因组 qry.fa 就可以上手进行比对了:

./nucmer -p <prefix> ref.fa qry.fa

上述命令,会生成如下文件:

<prefix>.delta

如果要查看比对的坐标,可以使用如下命令:

show-coords <prefix>.delta > <prefix>.coords

如果想要可视化,在安装了 gnuplot 的前提下,可以直接运行 perl 脚本 mummerplot,输入文件为第一步比对的结果文件:

./mummerplot -l <prefix>.delta


如果想要进一步了解该软件的使用,可以进一步阅读我写的其它分享;

序列比对软件 MUMmer 高级使用(二)

序列比对软件 MUMmer 结果文件解读(三)

序列比对软件 MUMmer 结果可读化处理(四)

要是觉得有用可以登录一下账号,点个赞,以表支持!

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容