###can used
mafft --auto --inputorder /fa/$a > mafft/${a%.*}.faa
iqtree -s ../$a -bb 1000 -nt AUTO
########################################################################################
mafft --auto --inputorder GmARF8a_8b_Medtr5g076270.1.txt.fa > out2
mafft --auto --clustalout --inputorder GmARF8a_8b_Medtr5g076270.1.txt.fa > out2
###多序列比对&&建树#############################################################################################################################
#muscle 生成phy格式
#muscle -in mouse_J.pro -physout seqs.phy
#ClustalW2 生成phy格式###不好用啊
#echo "1\nyour.fa\n2\n9\n1\n4\n\n1\nyour.phy\nyour.dnd\nX\n\nX\n" |clustalw2
#MAFFT 生成phy格式
mafft --auto --phylipout --inputorder final.faa > final.phy
###muscle结果可能用到使用BMGE对序列剪齐
#java -jar BMGE.jar -i 16s_muscle.fas -t DNA -of 16s_muscle_trim.fas -m DNAPAM250:4 –g 0.5
iqtree2 -s 16s_muscle_trim.fas -T AUTO -B 1000
iqtree -s 16s_muscle_trim.fas -bb 1000 -nt AUTO
#使用ModelFinder寻找最优模型###得到最佳模型TIM3+F+R4
iqtree -s 16s_muscle_trim.fas -m MF -nt 2
#使用UFBoot2建树
iqtree -s 16s_muscle_trim.fas -m TIM3+F+R4 -bb 1000 -nt AUTO
################################################################################################################################################
多序列比对&&建树(mafft && iqtree2)
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
