We've moved!
This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!

VariantEval error IllegalStateException: The elements of the input Iterators are not sorte

I am able to run VariantEval without a --known argument on my VCF, but it doesn't work when I give the --known argument specifying dbsnp_138.b37.vcf.

Specifically this command:

gatk VariantEval -R Homo_sapiens.GRCh37.dna.primary_assembly.fa -O all.vcf.eval.grp --eval all.vcf.gz --dbsnp dbsnp_138.b37.vcf

(GATK version 4.1.2.0)

As you can see, the reference is GRCh37. The reference fasta, and thus also my all.vcf.gz, contain exactly the same contigs as dbSnp, but in a different order. The analysis terminates after analysing chromosomes: 1, 10, 11, ... 19, 2, 20, 21, 22, 3, 4, .... 9:

15:02:54.016 INFO  ProgressMeter -           9:87555662             47.8              44864000         938981.7
15:03:04.031 INFO  ProgressMeter -           9:97271429             47.9              45106000         940760.2
15:03:14.052 INFO  ProgressMeter -          9:106958387             48.1              45342000         942399.6
15:03:24.055 INFO  ProgressMeter -          9:115867814             48.3              45572000         943909.5
15:03:34.065 INFO  ProgressMeter -          9:126084717             48.4              45820000         945777.7
15:03:44.086 INFO  ProgressMeter -          9:134157685             48.6              46030000         946848.2
15:03:52.808 INFO  VariantEval - Shutting down engine
[November 13, 2019 3:03:52 PM CET] org.broadinstitute.hellbender.tools.walkers.varianteval.VariantEval done. Elapsed time: 48.83 minutes.
Runtime.totalMemory()=12447121408
java.lang.IllegalStateException: The elements of the input Iterators are not sorted according to the comparator htsjdk.variant.variantcontext.VariantContextComparator
    at htsjdk.samtools.util.MergingIterator.next(MergingIterator.java:107)
    at java.util.Iterator.forEachRemaining(Iterator.java:116)
    at java.util.Spliterators$IteratorSpliterator.forEachRemaining(Spliterators.java:1801)
    at java.util.stream.AbstractPipeline.copyInto(AbstractPipeline.java:481)
    at java.util.stream.AbstractPipeline.wrapAndCopyInto(AbstractPipeline.java:471)
    at java.util.stream.ForEachOps$ForEachOp.evaluateSequential(ForEachOps.java:151)
    at java.util.stream.ForEachOps$ForEachOp$OfRef.evaluateSequential(ForEachOps.java:174)
    at java.util.stream.AbstractPipeline.evaluate(AbstractPipeline.java:234)
    at java.util.stream.ReferencePipeline.forEach(ReferencePipeline.java:418)
    at org.broadinstitute.hellbender.engine.MultiVariantWalker.traverse(MultiVariantWalker.java:117)
    at org.broadinstitute.hellbender.engine.GATKTool.doWork(GATKTool.java:1039)
    at org.broadinstitute.hellbender.cmdline.CommandLineProgram.runTool(CommandLineProgram.java:139)
    at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMainPostParseArgs(CommandLineProgram.java:191)
    at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:210)
    at org.broadinstitute.hellbender.Main.runCommandLineProgram(Main.java:162)
    at org.broadinstitute.hellbender.Main.mainEntry(Main.java:205)
    at org.broadinstitute.hellbender.Main.main(Main.java:291)

What could be the reason that I get this result when specifying dbSnp as a source of known variants?

Tagged:

Answers

Sign In or Register to comment.