2020-07-15记录下今天做进化分析遇到的问题

用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

  1. Figtree可视化树文件

ps:在进行构建进化树之前可以用pfam数据库鉴定一下提取的蛋白序列文件是否全部是目标蛋白质结构。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。