老师您好!我在使用ParaAT是报错,提示Generating homologous group files for alignments: 0,请问这个如何解决?

具体的错误提示如下:

****************************************************************

Checking input parameters: [OK]

Checking Epal2nal.pl: [OK]

Checking multiple sequence aligner: muscle[OK]

Checking output format: [OK]

Checking output folder: [OK]


Outputing input parameters:

Homologous groups = homo_pairs.txt1

Amino acids sequences = Input_pep.fa

Nucleotide sequences = Input_cds.fa

Process file = proc.txt

Output folder = align_out

Multiple sequence aligner = muscle

Output format = axt

Genetic code = 1-The Standard Code 

KaKs_Calculator used = FALSE

Remove gap = FALSE

Remove mismatched codons = FALSE

Verbose output = FALSE


Reading homologous groups from homo_pairs.txt1: 100 groups

Reading nucleotide sequences from Input_cds.fa: 118339 nucleotide sequences

Reading amino acid sequences from Input_pep.fa: 118339 amino acid sequences

Generating homologous group files for alignments: 0 from 100 groups

Aligning amino acid sequences for 0 homologous groups:


Translating amino acid alignments into codon alignments:

Cleaning temporary files:

Mission Accomplished (Time used: 0s).


请先 登录 后评论

2 个回答

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

Reading homologous groups from homo_pairs.txt1: 100 groups

Reading nucleotide sequences from Input_cds.fa: 118339 nucleotide sequences

Reading amino acid sequences from Input_pep.fa: 118339 amino acid sequences 

看看这些文件里面的序列ID是否一致

100 行基因的ID是否在序列里面

请先 登录 后评论

老师您好,是有的,我拿homo_pairs.txt1里的序列ID尝试了下,可以搜索到。序列的ID也确认了一致,但是仍然不行。

grep  'GWHTABJY009817'  Input_cds.fa 

>GWHTABJY009817



请先 登录 后评论

相似问题