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尋找最優(yōu)模型

iqtree -s 16s_muscle_trim.fas -m MF -nt 2

出現(xiàn)錯誤,序列差異太大,使用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數(shù)據(jù)庫鑒定一下提取的蛋白序列文件是否全部是目標蛋白質結構。

最后編輯于
?著作權歸作者所有,轉載或內容合作請聯(lián)系作者
【社區(qū)內容提示】社區(qū)部分內容疑似由AI輔助生成,瀏覽時請結合常識與多方信息審慎甄別。
平臺聲明:文章內容(如有圖片或視頻亦包括在內)由作者上傳并發(fā)布,文章內容僅代表作者本人觀點,簡書系信息發(fā)布平臺,僅提供信息存儲服務。

友情鏈接更多精彩內容