可以做蛋白的比对,在转换成cds的比对:https://www.omicsclass.com/article/1939
老师我在做KAKS分析时,cds的fa序列是三的倍数,但生成的.axt文件对比序列不是三的倍数,我按操作手册的方法,perl ../script/KAKS_SHAIXUAN.pl -in1 WRKY_cds_confirmed.fa.fai -in2 WRKY_cds_confirmed_blast.out -out duplication_gene.out,中提高的筛选条件,甚至我调高到90,输出的结果中还是包括不是三的倍数的比对序列,这种情况应该怎么办呢,
还有,同样的问题,由于基因组注释的问题,我的CDS序列没有起始密码子,提高筛选条件后并没有筛掉这些序列,这种情况怎么分析呢?还是说直接把这一对丢掉不要了?谢谢老师