ParaAT 多序列比对 报错 2.0

老师好,我写了sh后台运行ParaAT 多序列比对
查看运行完成的日志文件时,一共有148个这样的报错情况:

在将氨基酸序列(pep)转换为核酸序列(nuc)时发现了不一致性

我查找了其中三个ID对应的cds和pep序列,碱基/氨基酸个数为:cds=pep*3+3,是3倍关系

——————

并且继续用这个结果进行下一步的合并分析时  #合并成supergene

seqkit concat ./aln_out/*.pep_aln.renamed  >  single_copy.pep_msa.fasta

输出的结果是空的……

用蛋白质或cds比对结果构建进化树都出现这个情况

——————

前面都是按照课程的脚本跑的

请问该怎么处理这个报错呢?

attachments-2024-04-p5m4onFh661ce4145c4ba.png


请先 登录 后评论

1 个回答

omicsgene - 生物信息
擅长:重测序,遗传进化,转录组,GWAS

这个需要测试查看你的输入文件是不有问题;光从你的描述我也不好推断原因;

这些异常的基因都去看看删除一下试试,还有就是文件不要在windows下用记事本修改导致文件兼容问题;编辑用专业的工具并设置:https://www.omicsclass.com/article/395

请先 登录 后评论