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

登录

该问题已被锁定！

1: 关注

1171: 浏览

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 人关注

相关问题

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

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

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

可否用sanger测序测某一基因的上下游序列: 705 浏览 2 关注 1 回答 0 评论

如何提取比对上某一个基因或者某一段序列的所有reads: 560 浏览 2 关注 1 回答 0 评论

有哪些常用的有参序列拼接工具？: 683 浏览 3 关注 1 回答 0 评论

如何知道自己测序的adapter序列？: 966 浏览 3 关注 1 回答 0 评论

知道一批基因的具体位置，如何批量的从基因组中取出基因序列？: 870 浏览 3 关注 2 回答 0 评论

blastp后，我想提取score的序列，请问使用python应该怎么写这个脚本: 698 浏览 2 关注 1 回答 0 评论

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

推荐内容

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

问如何按一个列表对基因型文件进行过滤，剔除不需要的样本？: 841 浏览 2 关注 1 回答 0 评论

问 Bowtie2使用中遇到的内存不足问题: 732 浏览 1 关注 2 回答 0 评论

问 linux下非root用户设置所运行任务的CPU占用率和线程数: 1232 浏览 2 关注 3 回答 0 评论

问 GAPIT包导出的GWAS结果如何添加新的阈值线？以及GAPIT的结果文件中的nobs、H&B.P.Value、Effect分别是什么意思？: 1101 浏览 2 关注 1 回答 0 评论

问用GWAS筛选受选择基因，样本数量不够，请问可以用选择消除分析吗: 677 浏览 2 关注 1 回答 0 评论

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

问关于RSEM和RPKM: 666 浏览 2 关注 1 回答 0 评论

问使用deeptools生成ChIP-seq信号热图与谱图: 2085 浏览 3 关注 2 回答 0 评论

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

All Rights Reserved Powered BY WeCenter V4.1.0 © 2024