比较基因组学03.Phylo_Tree/aln_out

比较基因组学第三步比对完有的.cds_aln.fasta为0,没有内容


导致#合并成supergene
seqkit concat aln_out/*.cds_aln.fasta.renamed  >  single_copy.cds_msa.fasta  为0



请先 登录 后评论

2 个回答

Song

attachments-2024-10-XIDWDQFx66fe0f6723fb6.png

请先 登录 后评论
Ti Amo

看一下比对的日志文件里面有没有报错说”将氨基酸序列(pep)转换为核酸序列(nuc)时发现了不一致性“,如果有的话,你找到对应基因ID的cds和蛋白序列检查他们之间是不是3倍关系。如果数量不多,可以把这样的基因对删除。

需要paraAT比对的log文件报错信息进一步确认。

请先 登录 后评论
  • 1 关注
  • 0 收藏,282 浏览
  • Song 提出于 2024-10-03 11:27

相似问题