Attention:
The frontline support team will be unavailable to answer questions on April 15th and 17th 2019. We will be back soon after. Thank you for your patience and we apologize for any inconvenience!

HaplotypeCaller haploid GVCF format error?

Hi all,

I prepared a clean Bam file following GATK Best Practice and used GATK4 HaplotypeCaller to create a gvcf with ploidy1 option:

gatk-4.0.2.1/gatk HaplotypeCaller --native-pair-hmm-threads 24 -I KU_filtered_sorted_mdup.bam -O HC.KU.raw.snps.indels.g.vcf -R ref.fasta -ploidy 1 --emit-ref-confidence GVCF

When I validated the gvcf, ValidateVariants threw errors at the end:

<br />11:27:55.681 INFO  ProgressMeter - Traversal complete. Processed 124689522 total variants in 3.8 minutes.
11:27:55.681 INFO  ValidateVariants - Shutting down engine
[April 10, 2018 11:27:55 AM JST] org.broadinstitute.hellbender.tools.walkers.variantutils.ValidateVariants done. Elapsed time: 3.82 minutes.
Runtime.totalMemory()=4682940416
java.lang.IllegalArgumentException: Illegal character in path at index 15:HC.KU.raw.snps.indels.g.vcf
    at java.net.URI.create(URI.java:852)
    at org.broadinstitute.hellbender.engine.FeatureInput.makeIntoAbsolutePath(FeatureInput.java:242)
    at org.broadinstitute.hellbender.engine.FeatureInput.toString(FeatureInput.java:314)
    at java.util.Formatter$FormatSpecifier.printString(Formatter.java:2886)
    at java.util.Formatter$FormatSpecifier.print(Formatter.java:2763)
    at java.util.Formatter.format(Formatter.java:2520)
    at java.util.Formatter.format(Formatter.java:2455)
    at java.lang.String.format(String.java:2940)
    at org.broadinstitute.hellbender.engine.FeatureDataSource.close(FeatureDataSource.java:589)
    at org.broadinstitute.hellbender.engine.FeatureManager.lambda$close$9(FeatureManager.java:505)
    at java.util.LinkedHashMap$LinkedValues.forEach(LinkedHashMap.java:608)
    at org.broadinstitute.hellbender.engine.FeatureManager.close(FeatureManager.java:505)
    at org.broadinstitute.hellbender.engine.GATKTool.onShutdown(GATKTool.java:857)
    at org.broadinstitute.hellbender.engine.VariantWalker.onShutdown(VariantWalker.java:95)
    at org.broadinstitute.hellbender.cmdline.CommandLineProgram.runTool(CommandLineProgram.java:138)
    at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMainPostParseArgs(CommandLineProgram.java:180)
    at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:199)
    at org.broadinstitute.hellbender.Main.runCommandLineProgram(Main.java:159)
    at org.broadinstitute.hellbender.Main.mainEntry(Main.java:202)
    at org.broadinstitute.hellbender.Main.main(Main.java:288)
Caused by: java.net.URISyntaxException: Illegal character in path at index 15: /media/yoshi/My Book/Aet_v4.0_ChrSeqSplit/HC.KU-2103.raw.snps.indels.g.vcf
    at java.net.URI$Parser.fail(URI.java:2848)
    at java.net.URI$Parser.checkChars(URI.java:3021)
    at java.net.URI$Parser.parseHierarchical(URI.java:3105)
    at java.net.URI$Parser.parse(URI.java:3063)
    at java.net.URI.<init>(URI.java:588)
    at java.net.URI.create(URI.java:850)
    ... 19 more

When I ran GenotypeGVCFs using the gvcf, it ran to completion, but threw the same "java.lang.IllegalArgumentException" errors at the end.

My java is Java runtime: OpenJDK 64-Bit Server VM v1.8.0_162-8u162-b12-0ubuntu0.16.04.2-b12.

Could you please help me with this problem?
Many thanks.

Best Answer

Answers

Sign In or Register to comment.