用mega软件进行进化分析,
最大似然法,boot值设定以前,使用泊松法。
使用iqtree进行进化分析
使用perl脚本提取蛋白质序列过程中,遇到Can't locate Bio/SeqIO.pm in
尝试重新安装bioperl无果
尝试cpan安装bio/seqIO模块无果
无意中重启ubuntu解决问题。
安装iqtree
conda install iqtree
model finder寻找最优模型
iqtree -s 16s_muscle_trim.fas -m MF -nt 2
出现错误,序列差异太大,使用muscle对齐序列
muscle in srf_pep_rmw.fa -out srf_pep_rmw_muscle.fa
再次寻找最佳模型得到最佳模型 VT+F+R10
建树
iqtree -s 16s_muscle_trim.fas -m TIM3+F+R4 -bb 1000 -nt AUTO
- Figtree可视化树文件
ps:在进行构建进化树之前可以用pfam数据库鉴定一下提取的蛋白序列文件是否全部是目标蛋白质结构。