基因家族分析时,提取cds与pep序列报错

应该如何解决?代码与报错如下,数据来源:ncbi的参考基因组注释文件。

/share/work/biosoft/TransDecoder/latest/util/gff3_file_to_proteins.pl  --gff3  $species.longest_isoform.gff3 --fasta $genome  --seqType CDS >$species.cds.fa
#提取pep序列
/share/work/biosoft/TransDecoder/latest/util/gff3_file_to_proteins.pl --gff3 $species.longest_isoform.gff3 --fasta $genome --seqType prot >$species.pep.fa

报错如下:

  1. No such file or directory at /share/work/biosoft/TransDecoder/TransDecoder-v5.5.0/util/../PerlLib/GFF3_utils2.pm line 40, <$filehandle> line 3167.
    No such file or directory at /share/work/biosoft/TransDecoder/TransDecoder-v5.5.0/util/../PerlLib/GFF3_utils2.pm line 40, <$filehandle> line 3167.
请先 登录 后评论

1 个回答

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

提示文件路径不对,你检查变量对应的文件路径是否正确

请先 登录 后评论
  • 1 关注
  • 0 收藏,952 浏览
  • A.C Ver2 提出于 2023-11-16 19:36

相似问题