今天合并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 faidx 和samtools dict去建立索引,还是出现一样的错误
后来我用picard去建立索引,如下
picard CreateSequenceDictionary R=GCF_000002315.5_GRCg6a_genomic.fa O=GCF_000002315.5_GRCg6a_genomic.dict
结果运行成果!

image.png