ASR祖先序列重建，最后一步使用PAML时出现一些问题 - 做专业的生信问答网站

登录

该问题已被锁定！

1: 关注

4244: 浏览

ASR祖先序列重建，最后一步使用PAML时出现一些问题

这是我目前的ASR流程

1.在refseq搜集相关蛋白质序列，并cd-hit以0.95相似度去冗余剩下大概四五百条
2.mafft进行多序列比对，参数--auto，trimal剪辑对齐比对结果
3.iqtree建系统发育树
4.paml下面codeml ASR分析

关键卡在了第四步linux版本的codeml需要输入codeml.ctl配置文件，aaRatefile=氨基酸替代模型iqtree预测的是LG+F+R6，但是paml自带的.dat氨基酸替换模型没有这个组合，请问我该如何生成.dat文件。下面是我的配置文件比对序列时氨基酸序列，希望各位大佬检查一下有没有上面问题。

seqfile = supergene.phy
treefile = ff.fa.treefile
outfile = mlc
getSE = 0
noisy = 9
verbose = 1
seqtype = 2
runmode = 0
CodonFreq = 2
clock = 0
aaDist = 0
model = 2
aaRatefile =
icode = 0
Mgene = 0
fix_kappa = 0
kappa = 2
fix_omega = 0
omega = .4
fix_alpha = 1
alpha = 0
Malpha = 0
ncatG = 10
RateAncestor = 1
Small_Diff = .5e-6
method = 0

好问题 0 评论收藏举报

查看全部 3 个回答

BYXU123 注册会员用户来自于: 湖北省武汉市
2023-05-31 20:13

ProTest3可以生成氨基酸替换模型吗？不是只能预测吗？

我没有找到相关参数

赞同 0 0评论

关于作者

: BYXU123 注册会员
这家伙很懒，还没有设置简介

1: 回答

0: 文章

5: 问题

问题动态

发布时间: 2023-05-31 16:19

更新时间: 2023-05-31 20:49

关注人数: 1 人关注

相关问题

从多个.fa文件中提取以“poptri.”开头的蛋白序列: 1718 浏览 2 关注 1 回答 0 评论

请问一下，如何根据基因名批量下载CDS区的序列？: 2358 浏览 2 关注 1 回答 0 评论

如何从NCBI上分别下载所有的RNA病毒和DNA病毒的序列: 3396 浏览 2 关注 2 回答 0 评论

常用的对测序序列进行聚类分析的软件有什么？及其使用的了什么算法？: 2933 浏览 2 关注 2 回答 0 评论

翻译序列: 1613 浏览 2 关注 1 回答 0 评论

使用bedtools提取序列时的deyond length问题；bed负值: 3251 浏览 2 关注 2 回答 0 评论

根据已知序列查找其第三代基因组编号: 1860 浏览 2 关注 1 回答 0 评论

序列的二级结构呈现一下: 2112 浏览 2 关注 1 回答 0 评论

如果在基因组中鉴定到病毒序列，如何从生信角度判断它是EVEs还是HGTs: 1859 浏览 2 关注 1 回答 0 评论

植物基因组组装过程中如何去除质体序列: 2314 浏览 2 关注 1 回答 0 评论

推荐内容

问请问在不考虑测序质量的情况下，如何根据bam文件获得比对到指定position的碱基？: 4010 浏览 2 关注 2 回答 0 评论

问使用tophat2和bowtie1寻找环形RNA时报错: 2942 浏览 2 关注 2 回答 0 评论

问 RNA-seq比对region：exonic/intronic/intergenic 比例异常: 2951 浏览 2 关注 2 回答 0 评论

问使用Tracking Tumor ImmunoPhenotype(TIP)网站分析TCGA的BLCA_tpm数据: 2466 浏览 2 关注 2 回答 0 评论

问生存分析KM-plot交叉问题: 2831 浏览 2 关注 2 回答 0 评论

问请问如何查询某基因的详细的Biological Function？: 1767 浏览 2 关注 1 回答 0 评论

问 cox: 1962 浏览 2 关注 1 回答 0 评论

问 linux下使用convert出现报错，可能是什么原因？如何解决？: 2686 浏览 2 关注 1 回答 0 评论

问 chipseq分析，利用deeptools的computeMatrix reference-point画样本的peak center信号图: 3813 浏览 3 关注 2 回答 0 评论

问动态库存在但调用报错问题: 2099 浏览 2 关注 1 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2025