This is a test version of Biostars. For the public version, visit https://www.biostars.org.
gatk vqsr --recal-file error

In GATK 4.0.6.0 I am getting the below error using the command:

    gatk VariantRecalibrator \
        -R /home/cmccabe/Desktop/fastq/resources/ucsc.hg19.fasta \
        -V NA12878.vcf.gz \
   --resource hapmap,known=false,training=true,truth=true,prior=15.0:/home/cmccabe/Desktop/fastq/resources/hapmap_3.3.hg19.sites.vcf.gz \
   --resource omni,known=false,training=true,truth=false,prior=12.0:/home/cmccabe/Desktop/fastq/resources/1000G_omni2.5.hg19.sites.vcf.gz \
   --resource 1000G,known=false,training=true,truth=false,prior=10.0:/home/cmccabe/Desktop/fastq/resources/1000G_phase1.snps.high_confidence.hg19.sites.vcf.gz \
   --resource dbsnp,known=true,training=false,truth=false,prior=2.0:/home/cmccabe/Desktop/fastq/resources/dbsnp_138.hg19.vcf.gz \
   -an QD -an MQ -an MQRankSum -an ReadPosRankSum -an FS -an SOR \
   -mode SNP \
   --recal-file /home/cmccabe/Desktop/fastq/NA12878.recal \
   --tranches-file /home/cmccabe/Desktop/fastq/NA12878.tranches \
   --rscript-file /home/cmccabe/Desktop/fastq/NA12878.plots.R


A USER ERROR has occurred: recal-file is not a recognized option

Is that option not supported, I am following the best practices docs and think its correct or is there something else. Thank you :).

gatk vqsr
• 2,908 views • link
updated by Santosh Anand Santosh Anand • written by bioguy24 bioguy24

1 answer

Probably the argument has been deleted as I could not find anything on the argument lists

https://software.broadinstitute.org/gatk/documentation/tooldocs/4.0.3.0/org_broadinstitute_hellbender_tools_walkers_vqsr_VariantRecalibrator.php

Also, there is a thread on GATK reg. this and you may also like to follow there

https://gatkforums.broadinstitute.org/gatk/discussion/12449/gatk-4-variantrecalibrator-example-using-unavailable-option-recal-file

• 1 views • link
written by Santosh Anand Santosh Anand

Thank you very much :).

• 1 views • link
written by bioguy24 bioguy24

Log in to answer this question.

More posts like this
  • Help wanted for a struggling bioinfmatician! GATK Variantrecalibration
    written by chrisgr • chrisgr
    Hello, any answer on any the following questions will be much appreciated! I'm playing around with gatk's VariantRecalibration tool. I have a few questions that …
  • Can not generate tranches plot after VQSR
    written by Gerlumen • Gerlumen
    Hello! I have encountered with a strange issue when trying to analyze VariantRecalibrator output. I can not generate or find an R-script file of *.tranche …
  • gatk VariantRecalibrator positional argument error
    written by michael.flower.14 • michael.flower.14
    I'm trying to use recalibrate my vcf using gatk VariantRecalibrator, but keep getting an error "Illegal argument value: Positional arguments were provided". But I don't …
  • GATK : VariantRecalibrator no positional argument is defined for this tool.
    written by Arslan Tariq • Arslan Tariq
    Hi, I am trying to run the following command: gatk VariantRecalibrator \ -R genome.fa \ -V all.Sample.SNP.vcf.gz \ --trust-all-polymorphic -tranche 100.0 -tranche 99.95 -tranche 99.9 …
  • VQSR no r file being created
    written by geneart$$ • geneart$$
    Hi , I have r module loaded in the path I am running my scripts , yet the `--rscript all_jointcalls.AS.plots.R` is not even created ! …
  • Tranche Plot after VSQR
    written by BAGeno BAGeno
    Hi I am doing exome sequencing of 15 patients. After variant calling I run `VariantRecalibrator` using this command gatk VariantRecalibrator -R new_hg38.fa \ -V allsample.vcf …
  • GATK4 VariantRecalibrator resource error
    written by samreen.jasvi • samreen.jasvi
    My following command java -Xmx10g -Djava.io.tmpdir=Varcalling/cohort.vcf_VQSR -jar gatk-package-4.0.1.2-local.jar VariantRecalibrator -R /xxx/REFERENCES/hg19//bwa_index/hg19.fa -V Varcalling/cohort.vcf -resource hapmap,VCF,known=false,training=true,truth=true,prior=15.0 /xxx/REFERENCES/hg19/hapmap_3.3.hg19.sites.reorder.vcf -resource omni,known=false,training=true,truth=false,prior=12.0 /xxx/REFERENCES/hg19//1000G_omni2.5.hg19.sites.reorder.vcf -resource dbsnp,known=true,training=false,truth=false,prior=6.0 /xxx//REFERENCES/hg19//dbsnp_138.hg19.reorder.vcf -resource 1000G,known=false,training=true,truth=false,prior=10.0 /xxx//REFERENCES/hg19//1000G_phase1.snps.high_confidence.hg19.sites.reorder.vcf …
  • error using bqsr and incompatible contigs
    written by bioguy24 bioguy24
    Using the below command in GATK 4.0.6.0, I am getting: A USER ERROR has occurred: Input files reference and features have incompatible contigs: No overlapping …
  • multithread support for gatk
    written by bioguy24 bioguy24
    The below command is meant to use bqsr with multithreading: gatk BaseRecalibrator --num_cpu_threads_per_data_thread=8 \ -I /home/cmccabe/Desktop/fastq/NA12878_marked_duplicates.bam \ -R /home/cmccabe/Desktop/fastq/resources/ucsc.hg19.fasta \ --known-sites /home/cmccabe/Desktop/fastq/fastq/resources/dbsnp_138.b37.vcf \ --known-sites /home/cmccabe/Desktop/fastq/resources/Mills_and_1000G_gold_standard.indels.hg19.vcf …
  • Facing error while trying to build SNP recalibration model using GATK
    written by adnan.haider.ra • adnan.haider.ra
    <p>Hello everyone,</p> <p>Kindly help me in fixing following error while trying to build SNP recalibration model uing GATK,</p> <pre> $ java -jar GenomeAnalysisTK.jar -T VariantRecalibrator …