多序列比对

最经典和广为熟知的多序列比对软件是 clustalw 。 但是现有的多序列比对软件较多，有文献报道：比对速度（Muscle>MAFFT>ClustalW>T-Coffee），比对准确性（MAFFT>Muscle>T-Coffee>ClustalW）。因此，推荐使用 MAFFT 软件进行多序列比对。【来自陈连福博客】

MAFFT

一、安装

wget https://mafft.cbrc.jp/alignment/software/mafft-7.453-with-extensions-src.tgz
tar zxvf mafft-7.453-with-extensions-src.tgz

vi修改extensions和core文件夹中MakeFile文件的第一行prefix路径到个人文件夹

cd mafft-7.453-with-extensions/extensions
make clean
make
make install

cd ../core
make clean
make
make install
#mafft程序在scripts里

二、运行

1. 精确比对

最准确的方法。适合于 <200 条序列，且序列长度 <~2000 aa/nt 的比对

mafft --maxiterate 1000 --localpair input.fa > output.fa 或者 linsi input.fa > output.fa

适合于序列长度相似的多序列比对。序列条数 <200, 序列长度 <~2000 aa/nt

mafft --maxiterate 1000 --genafpair  input.fa > output.fa 或者 einsi input.fa > output.fa

适合序列中包含较大的非匹配区域。序列条数 <200, 序列长度 <~2000 aa/nt

mafft --maxiterate 1000 --globalpair input.fa > output.fa 或者 ginsi input.fa > output.fa

2. 节约时间

减少迭代次数，最大迭代次数减为 2

mafft --retree 2 --maxiterate 2 input.fa > output.fa 或者 fftnsi input.fa > output.fa

最大迭代次数减为 0

mafft --retree 2 --maxiterate 0 input.fa > output.fa 或者 fftns input.fa > output.fa

此方法非常快速，适合 >2000 条序列的多序列比对

mafft --retree 1 --maxiterate 0 input.fa > output.fa

迭代过程中不进行 FFT aproximation

mafft --retree 2 --maxiterate 2 --nofft input.fa > output.fa 或者 nwnsi input.fa > output.fa
mafft --retree 2 --maxiterate 0 --nofft input.fa > output.fa 或者 nwns input.fa > output.fa

3 个参数都设置为最不消耗时间的类型，适合于 ~10,000 到 ~50,000 条序列的比对

mafft --retree 1 --maxiterate 0 --nofft --parttree input.fa > output.fa