If you happen to see a question you know the answer to, please do chime in and help your fellow community members. We appreciate your help!
Test-drive the GATK tools and Best Practices pipelines on Terra
Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.
Syntax problem with VariantRecalibrator
Running VariantRecalibrator on GATK 18.104.22.168. (I know it's an older version - I need this version due to issues in more recent versions with functionality at a different step in my pipeline from this one. Also, I have tested this step on 22.214.171.124 and get the same error.)
My command line:
java -Dsamjdk.use_async_io_read_samtools=false -Dsamjdk.use_async_io_write_samtools=true -Dsamjdk.use_async_io_write_tribble=false -Dsamjdk.compression_level=2 -Xmx10G -jar /seq/vgb/software/gatk/gatk-126.96.36.199/gatk-package-188.8.131.52-local.jar VariantRecalibrator -R /seq/references/Canis_lupus_familiaris_assembly3/v1/Canis_lupus_familiaris_assembly3.fasta -V /seq/vgb/jphekman/singleBamToVcf/31001701034247.vcf --resource dbsnp,known=true,training=false,truth=false,prior=2.0:/seq/vgb/references/dog/population_variation/dog_dbSNP.vcf.gz -an QD -an MQ -an MQRankSum -an ReadPosRankSum -an FS -an SOR -an InbreedingCoeff -mode SNP -O /seq/vgb/jphekman/singleBamToVcf/31001701034247.recal --tranches-file /seq/vgb/jphekman/singleBamToVcf/31001701034247.tranches --rscript-file /seq/vgb/jphekman/singleBamToVcf/31001701034247.plots.R
A USER ERROR has occurred: Couldn't read file file:///home/unix/jphekman/dbsnp,known=true,training=false,truth=false,prior=2.0:/seq/vgb/references/dog/population_variation/dog_dbSNP.vcf.gz. Error was: It doesn't exist.
Note that both /seq/vgb/references/dog/population_variation/dog_dbSNP.vcf.gz and /seq/vgb/references/dog/population_variation/dog_dbSNP.vcf.gz.tbi do exist.
I've seen this error previously when I was using incorrect syntax in the --resources parameter. However I have perused the doc and I can't see where I'm going wrong. (Apologies if I'm missing something obvious.)
Help - what am I missing?