Question: Facing error while trying to build SNP recalibration model using GATK
0
gravatar for adnan.haider.ra
3.2 years ago by
Pakistan
adnan.haider.ra10 wrote:

Hello everyone,

Kindly help me in fixing following error while trying to build SNP recalibration model uing GATK,

$ java -jar GenomeAnalysisTK.jar -T VariantRecalibrator -R Reference/hg19.fa -input raw_variants.vcf -resource:hapmap,known=false,training=true,truth=true,prior=15.0 /Reference/hapmap_3.3.hg19.vcf -resource:omni,known=false,training=true,truth=true,prior=12.0 /Reference/1000G_omni2.5.hg19.vcf -resource:1000G,known=false,training=true,truth=false,prior=10.0 /Reference/1000G_phase1.snps.high_confidence.hg19.vcf -resource:dbsnp,known=true,training=false,truth=false,prior=2.0 /Reference/dbsnp_138.hg19.vcf -an QD -an FS -an SOR -an MQ -an MQRankSum -an ReadPosRankSum -mode SNP -tranche 100.0 -tranche 99.9 -tranche 99.0 -tranche 90.0 -recalFile recalibrate_SNP.recal -tranchesFile recalibrate_SNP.tranches -rscriptFile recalibrate_SNP_plots.R

ERROR ------------------------------------------------------------------------------------------

ERROR A USER ERROR has occurred (version 3.5-0-g36282e4):

ERROR

ERROR This means that one or more arguments or inputs in your command are incorrect.

ERROR The error message below tells you what is the problem.

ERROR

ERROR If the problem is an invalid argument, please check the online documentation guide

ERROR (or rerun your command with --help) to view allowable command-line arguments for this tool.

ERROR

ERROR Visit our website and forum for extensive documentation and answers to

ERROR commonly asked questions http://www.broadinstitute.org/gatk

ERROR

ERROR Please do NOT post this error to the GATK forum unless you have really tried to fix it yourself.

ERROR

ERROR MESSAGE: Invalid argument value '-resource:hapmap known=false training=true truth=true prior=15.0' at position 6.

ERROR Invalid argument value '/Reference/hapmap_3.3.hg19.vcf' at position 7.

ERROR Invalid argument value '-resource:omni known=false training=true truth=true prior=12.0' at position 8.

ERROR Invalid argument value '/Reference/1000G_omni2.5.hg19.vcf' at position 9.

ERROR Invalid argument value '-resource:1000G known=false training=true truth=false prior=10.0' at position 10.

ERROR Invalid argument value '/Reference/1000G_phase1.snps.high_confidence.hg19.vcf' at position 11.

ERROR Invalid argument value '-resource:dbsnp known=true training=false truth=false prior=2.0' at position 12.

ERROR Invalid argument value '/Reference/dbsnp_138.hg19.vcf' at position 13.

ERROR ------------------------------------------------------------------------------------------

Thanks and Best Regards,
Adnan

recalibration gatk • 1.5k views
ADD COMMENTlink modified 3.2 years ago by Pierre Lindenbaum118k • written 3.2 years ago by adnan.haider.ra10
0
gravatar for andrew.j.skelton73
3.2 years ago by
London
andrew.j.skelton735.5k wrote:

Pretty sure that this:

-resource:hapmap known=false training=true truth=true prior=15.0

Should be comma separated to 

-resource:hap map,known=false,training=true,truth=true,prior=15.0
ADD COMMENTlink written 3.2 years ago by andrew.j.skelton735.5k

Yeah pretty sure about it, I am providing these comma seperated arguments as mentioned in GATK guide, 

-resource:hapmap,known=false,training=true,truth=true,prior=15.0

but still facing the same error
ADD REPLYlink written 3.2 years ago by adnan.haider.ra10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2326 users visited in the last hour