利用picard构建参考基因组的dict索引文件

生信学习小达人

已于 2023-07-06 23:29:47 修改

阅读量1.1k

点赞数 1

分类专栏：报错文章标签：服务器

于 2023-07-06 23:00:36 首次发布

本文链接：https://blog.csdn.net/hgz2020/article/details/131586898

版权

报错专栏收录该内容

6 篇文章 0 订阅

订阅专栏

最近，自己在利用gatk中的Haplotypecaller进行calling-snp时，出现报错：

Runtime.totalMemory()=2304245760
***********************************************************************

A USER ERROR has occurred: Fasta dict file file:///data/home/hgzhong/population/workplace/Part1.variant_calling/03.SNP_indel_gatk/./index/genome.dict for reference file:///data/home/hgzhong/population/workplace/Part1.variant_calling/03.SNP_indel_gatk/./index/genome.fasta 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.

***********************************************************************
Set the system property GATK_STACKTRACE_ON_USER_EXCEPTION (--java-options '-DGATK_STACKTRACE_ON_USER_EXCEPTION=true') to print the stack trace.

报错显示 reference.genome.dict 文件不存在

故利用picard构建参考基因组的索引reference.genome.dict文件

java -jar /home/software/picard/build/libs/picard.jar CreateSequenceDictionary R=Arabidopsis_thaliana.TAIR10.dna.toplevel.fa O=Arabidopsis_thaliana.TAIR10.dna.toplevel.fa.dict

仍继续报错，原因是参数O=reference.dict,而不是O=reference.fasta.dict