Question: GATK4 VariantRecalibrator resource error
0
gravatar for samreen.jasvi
3 months ago by
samreen.jasvi0 wrote:

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 -an SOR -an QD -an MQRankSum -an ReadPosRankSum -an MQ  -an FS -an DP -mode SNP -O Varcalling/cohort.vcf.SNP.recal --tranches-file Varcalling/cohort.vcf.SNP.tranches --rscript-file Varcalling/cohort.vcf.SNP.plots.R

gives an error for the resource

 `***********************************************************************

A USER ERROR has occurred: Invalid argument '/sc/orga/projects/PBG/REFERENCES/hg19/hapmap_3.3.hg19.sites.reorder.vcf'.

***********************************************************************

I've tried changing it all the following but keep getting errors

 --resource hapmap,VCF,known=false,training=true,truth=true,prior=15.0 /xxx/REFERENCES/hg19/hapmap_3.3.hg19.sites.reorder.vcf

same error for above.

OR

 --resource:hapmap,VCF,known=false,training=true,truth=true,prior=15.0 /xxx/REFERENCES/hg19/hapmap_3.3.hg19.sites.reorder.vcf

the error here is

A USER ERROR has occurred: The argument: "resource/resource" does not accept tags: "hapmap,VCF,known=false,training=true,truth=true,prior=15.0"

Anyone know what the mistake is that I am making here?

wes ngs gatk4 • 171 views
ADD COMMENTlink written 3 months ago by samreen.jasvi0

Have a look at this thread

Basically you need a : between resource parameters and the resource path:

--resource hapmap,VCF,known=false,training=true,truth=true,prior=15.0:/xxx/REFERENCES/hg19/hapmap_3.3.hg19.sites.reorder.vcf

Although this also depends on your version of GATK...

ADD REPLYlink modified 3 months ago • written 3 months ago by bruce.moran620

It worked for me with a a singly hyphen before resource and no 'VCF' in the tags, like so, -resource hapmap,known=false,training=true,truth=true,prior=15.0:$HapMap_Sites

I have a similar syntax question, if its acceptable to ask here as a follow-up. HaplotypeCaller is supposed to be able to take a VCF as a comp track, I've been unable to do so without error. Only my -dbSNP argument works. from the documentation -> e.g. -comp:FOO will have 'FOO' in the INFO field). Records that are filtered in the comp track will be ignored. Note that 'dbSNP' has been special-cased (see the --dbsnp argument).

$GATK HaplotypeCaller \
        -I $INP \
        -O $OUTPUT \
        -R $REF \
        -L $Targets \
        -ERC GVCF \
        --dbsnp $DBSNP138 \
       -comp:hapmap3.3 $HapMapV3VCF \
        $infofields

Here too I've tried various combinations of spaces and colons in between -comp and the variable. Any thoughts. Is this even an important resource to use during Haplotype calling?

ADD REPLYlink written 3 months ago by samreen.jasvi0

Don't do much germline mutation calling, but I think ExAC might have superseded HapMap? No idea WRT -comp. I highly recommend the GATK website, it's a bit of a mess because of versioning but the users and mods are very active, and because it's GATK your error will have occurred for others and almost certainly will have a Q/A on there.

ADD REPLYlink written 3 months ago by bruce.moran620

Thanks, I'll go through their forum again, I've looked through before and couldn't find my query there. I can't register to post a question on GATK website, since I don't have a uni/org email (am a contractor).

ADD REPLYlink written 3 months ago by samreen.jasvi0

Hello samreen.jasvi!

We believe that this post does not fit the main topic of this site.

Question solved.

For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.

If you disagree please tell us why in a reply below, we'll be happy to talk about it.

Cheers!

ADD REPLYlink written 3 months ago by samreen.jasvi0
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: 1615 users visited in the last hour