GATK showing error of reference file index
0
0
Entering edit mode
16 months ago
rj.rezwan ▴ 10

Hi, i am using the GATK version gatk/4.2.2.0 for HaplotypeCaller. I have been facing the reference.fa indexing issue. I tried to index the file using the following command

samtools faidx ~/path/PitayaGenomic.fa

#three difference formats of reference files 
-rw-r--r-- 1 tariqr 1.3G Feb 21 17:50 PitayaGenomic.fa
-rw-r--r-- 1 tariqr  34K Feb 21 20:01 PitayaGenomic.fa.fai
-rw-r--r-- 1 tariqr 137K Feb 21 21:05 ref.dict

I tried to use all three formats in GATK but facing the following errors. The first error using the .fai format

java.lang.IllegalArgumentException: File is not a supported reference file type: ~/path/genome_gatk/PitayaGenomic.fa.fai

Second error using .fa file

A USER ERROR has occurred: Fasta dict file file:///path/genome_gatk/PitayaGenomic.dict for reference file:///path/genome_gatk/PitayaGenomic.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.

Third error after using the ref.dict format

java.lang.IllegalArgumentException: File is not a supported reference file type: ~/path/genome_gatk/ref.dict

Please suggest to me what should I do to fix the issue of reference. I really acknowledge your kind suggestion. The code is given here:

gatk HaplotypeCaller -R ~/path/genome_gatk/PitayaGenomic.fa -I ~/path/input.bam -O ~/path/output.g.vcf.gz --bam-output ~/path/output.bam -ERC GVCF --native-pair-hmm-threads 8
GATK HaplotypeCaller • 1.7k views
ADD COMMENT
0
Entering edit mode

The command uses: -R ~/path/PitayaGenomic.fa
However, the error says: path/genome_gatk/PitayaGenomic.fa

I see in your command that there is no genome_gatk

ADD REPLY
0
Entering edit mode

I have modified the question. Please have a look and suggest the answer

ADD REPLY
0
Entering edit mode

They're not 3 different formats. .fai is the index of the .fa. .dict is a different format, sure.

Run this command and then try your GATK command:

cd ~/path/
ln -s ref.dict PitayaGenomic.dict
ADD REPLY
0
Entering edit mode

still same issue

java.lang.IllegalArgumentException: File is not a supported reference file type
ADD REPLY
0
Entering edit mode

I used ~/path/ which seems to be the wrong place. Can you show us the output to:

cd ~/path/genome_gatk/
ln -s ref.dict PitayaGenomic.dict
pwd
ls -l
ADD REPLY

Login before adding your answer.

Traffic: 2464 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6