参考基因组建立索引报错

今天合并gVCF文件的时候,出现下面的错误

A USER ERROR has occurred: Fasta dict file file:///public/jychu/refs/ncbi/GCF_000002315.5_GRCg6a_genomic.dict for reference file:///public/jychu/refs/ncbi/GCF_000002315.5_GRCg6a_genomic.fa does not exist. Please see http://gatkforums.broadinstitute.org/discussion/1601/how-can-i-prepare-a-fasta-file-to-use-as-reference for help creating it.
我试着用samtools faidxsamtools dict去建立索引,还是出现一样的错误

后来我用picard去建立索引,如下

picard CreateSequenceDictionary R=GCF_000002315.5_GRCg6a_genomic.fa O=GCF_000002315.5_GRCg6a_genomic.dict

结果运行成果!


image.png
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容