kSNP3寻找SNPs并构建进化树

1. kSNP3 安装

cd ~/tools

wget https://sourceforge.net/projects/ksnp/files/kSNP3.1_Linux_package.zip

unzip kSNP3.1_Linux_package.zip

加入环境变量：

vim ~/.bashrc

i

export PATH=$PATH:$HOME/tools/kSNP3.1_Linux_package/kSNP3

ESC

shift + ：

wq!

source ~/.bashrc

编辑主程序kSNP3的第8行：

将 set kSNP=/usr/local/kSNP3

改为：set kSNP=/home/lhl/tools/kSNP3.1_Linux_package/kSNP3

注：根据自己的路径进行修改

---

2. 基本用法

将所有的基因组文件放于988_ksnp目录之中，在其上一级目录下运行命令创建输入列表：

MakeKSNP3infile 988_ksnp inlist A

运行命令创建输入序列集合：

MakeFasta inlist fastainput_988

计算最佳K值：

Kchooser fastainput_988

寻找SNPs并构建进化树：

nohup kSNP3 -in inlist -outdir SNPs_20181214 -k 23 -ML -NJ -vcf  -CPU 30 -core -min_frac 0.5 |tee Log_988_20181214.txt &