Hi,
sorry if the question isn't precise enough.
Here's the command:
tools/gatk-4.0.11.0/gatk VariantRecalibrator -R data/humanRefGenome/hs37d5.fa \
--variant data/mongolian/Mongolian_genome_sorted_noChr2.snp.vcf \
-O data/mongolian/mongolian.recal \
-tranches-file data/mongolian/mongolian.tranches \
-rscript-file data/mongolian/mongolian_plots.R \
--resource hapmap,known=false,training=true,truth=true,prior=15.0:data/mongolian/hapmap_3.3.hg19.sites_noChr2.vcf.gz \
--resource omni,known=false,training=true,truth=true,prior=12.0:data/mongolian/1000G_omni2.5.hg19.sites_noChr2.vcf.gz \
--resource 1000G,known=false,training=true,truth=false,prior=10.0:data/mongolian/1000G_phase1.snps.high_confidence.hg19.sites_noChr2.vcf.gz \
--resource dbsnp,known=true,training=false,truth=false,prior=2.0:data/mongolian/dbsnp_138.hg19_noChr2.vcf.gz \
-an QD -an MQ -an MQRankSum -an ReadPosRankSum -an FS -an DP \
-tranche 90.0 -tranche 99.0 -tranche 99.9 -tranche 100.0 \
-mode SNP
and that's the error:
21:10:45.445 INFO VariantRecalibrator - Writing out recalibration table...
21:10:45.459 INFO VariantRecalibrator - Shutting down engine
[November 15, 2018 9:10:45 PM CET] org.broadinstitute.hellbender.tools.walkers.vqsr.VariantRecalibrator done. Elapsed time: 33.89 minutes.
Runtime.totalMemory()=2090336256
java.lang.IllegalArgumentException: Can't do comparison because Locatables' contigs not found in sequence dictionary
at org.broadinstitute.hellbender.utils.IntervalUtils.compareContigs(IntervalUtils.java:149)
at org.broadinstitute.hellbender.utils.IntervalUtils.compareLocatables(IntervalUtils.java:85)
at org.broadinstitute.hellbender.tools.walkers.vqsr.VariantDatum.lambda$getComparator$2(VariantDatum.java:49)
at java.util.TimSort.binarySort(TimSort.java:296)
at java.util.TimSort.sort(TimSort.java:239)
at java.util.Arrays.sort(Arrays.java:1512)
at java.util.ArrayList.sort(ArrayList.java:1462)
at java.util.Collections.sort(Collections.java:175)
at org.broadinstitute.hellbender.tools.walkers.vqsr.VariantDataManager.writeOutRecalibrationTable(VariantDataManager.java:456)
at org.broadinstitute.hellbender.tools.walkers.vqsr.VariantRecalibrator.onTraversalSuccess(VariantRecalibrator.java:695)
at org.broadinstitute.hellbender.engine.GATKTool.doWork(GATKTool.java:968)
at org.broadinstitute.hellbender.cmdline.CommandLineProgram.runTool(CommandLineProgram.java:139)
at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMainPostParseArgs(CommandLineProgram.java:192)
at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:211)
at org.broadinstitute.hellbender.Main.runCommandLineProgram(Main.java:160)
at org.broadinstitute.hellbender.Main.mainEntry(Main.java:203)
at org.broadinstitute.hellbender.Main.main(Main.java:289)
The output file (mongolian.recal) doesn't have any variant:
...
##contig=<ID=6,length=171115067,assembly=hs37d5.fa>
##contig=<ID=7,length=159138663,assembly=hs37d5.fa>
##contig=<ID=8,length=146364022,assembly=hs37d5.fa>
##contig=<ID=9,length=141213431,assembly=hs37d5.fa>
##contig=<ID=M,length=16571,assembly=hs37d5.fa>
##contig=<ID=X,length=155270560,assembly=hs37d5.fa>
##contig=<ID=Y,length=59373566,assembly=hs37d5.fa>
##reference=hs37d5
##source=VariantRecalibrator
#CHROM POS ID REF ALT QUAL FILTER INFO
END OF FILE (there isn't other rows).
What am I doing wrong?