Heads up:
We’re moving the GATK website, docs and forum to a new platform. Read the full story and breakdown of key changes on this blog.
Notice:
If you happen to see a question you know the answer to, please do chime in and help your fellow community members. We encourage our fourm members to be more involved, jump in and help out your fellow researchers with their questions. GATK forum is a community forum and helping each other with using GATK tools and research is the cornerstone of our success as a genomics research community.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.
Attention:
We will be out of the office for a Broad Institute event from Dec 10th to Dec 11th 2019. We will be back to monitor the GATK forum on Dec 12th 2019. In the meantime we encourage you to help out other community members with their queries.
Thank you for your patience!

Error stack trace Error Message: For input string: "R"

Hi, I've run Samtools1.0 & bcftools1.0 which output vcf4.2 which caused issuse. So I've updated to the latest version of GATK & am now getting this Error which I cannot resolve.

Program Args: --variant SPARK.common.vcf -R Documents/RefSeqIndexAllPrograms/GATK/hg19_1stM_unmask_ran_all.fa -T SelectVariants -o 24999.SPARK.common.vcf -sn 24999 -env -nt 4

ERROR stack trace

java.lang.NumberFormatException: For input string: "R"
ERROR A GATK RUNTIME ERROR has occurred (version 3.2-2-gec30cee):
ERROR MESSAGE: For input string: "R"

Any help would be gratefully received.
Thanks alig

Answers

  • airtimeairtime Member

    Hi,

    it looks like your pathway of R is wrong or your reference file is not in correct order?
    Did you check this?

    Best Air

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie admin

    @alig Can you please post the complete stack trace (meaning the full error message)? This will help us narrow down the source of the error.

  • aligalig Member

    Hi, Here is the complete stack trace, thanks alig

    INFO 17:31:34,943 HelpFormatter - --------------------------------------------------------------------------------
    INFO 17:31:34,946 HelpFormatter - The Genome Analysis Toolkit (GATK) v3.2-2-gec30cee, Compiled 2014/07/17 15:22:03
    INFO 17:31:34,946 HelpFormatter - Copyright (c) 2010 The Broad Institute
    INFO 17:31:34,946 HelpFormatter - For support and documentation go to http://www.broadinstitute.org/gatk
    INFO 17:31:34,950 HelpFormatter - Program Args: --variant SPARKS.common.vcf -R /home/neuro/Documents/RefSeqIndexAllPrograms/GATK/hg19_1stM_unmask_ran_all.fa -T SelectVariants -o 24919.SPARKS.common.vcf -sn 24919 -env -nt 4
    INFO 17:31:34,957 HelpFormatter - Executing as [email protected] on Linux 3.2.0-67-generic amd64; OpenJDK 64-Bit Server VM 1.7.0_65-b32.
    INFO 17:31:34,957 HelpFormatter - Date/Time: 2014/09/01 17:31:34
    INFO 17:31:34,957 HelpFormatter - --------------------------------------------------------------------------------
    INFO 17:31:34,958 HelpFormatter - --------------------------------------------------------------------------------
    INFO 17:31:35,041 GenomeAnalysisEngine - Strictness is SILENT
    INFO 17:31:35,201 GenomeAnalysisEngine - Downsampling Settings: Method: BY_SAMPLE, Target Coverage: 1000
    INFO 17:31:37,130 GATKRunReport - Uploaded run statistics report to AWS S3

    ERROR ------------------------------------------------------------------------------------------
    ERROR stack trace

    java.lang.NumberFormatException: For input string: "R"
    at java.lang.NumberFormatException.forInputString(NumberFormatException.java:65)
    at java.lang.Integer.parseInt(Integer.java:492)
    at java.lang.Integer.valueOf(Integer.java:582)
    at htsjdk.variant.vcf.VCFCompoundHeaderLine.(VCFCompoundHeaderLine.java:171)
    at htsjdk.variant.vcf.VCFInfoHeaderLine.(VCFInfoHeaderLine.java:46)
    at htsjdk.variant.vcf.AbstractVCFCodec.parseHeaderFromLines(AbstractVCFCodec.java:205)
    at htsjdk.variant.vcf.VCFCodec.readActualHeader(VCFCodec.java:111)
    at htsjdk.tribble.AsciiFeatureCodec.readHeader(AsciiFeatureCodec.java:88)
    at htsjdk.tribble.AsciiFeatureCodec.readHeader(AsciiFeatureCodec.java:41)
    at htsjdk.tribble.index.IndexFactory$FeatureIterator.readHeader(IndexFactory.java:413)
    at htsjdk.tribble.index.IndexFactory$FeatureIterator.(IndexFactory.java:401)
    at htsjdk.tribble.index.IndexFactory.createDynamicIndex(IndexFactory.java:312)
    at org.broadinstitute.gatk.engine.refdata.tracks.RMDTrackBuilder.createIndexInMemory(RMDTrackBuilder.java:402)
    at org.broadinstitute.gatk.engine.refdata.tracks.RMDTrackBuilder.loadIndex(RMDTrackBuilder.java:288)
    at org.broadinstitute.gatk.engine.refdata.tracks.RMDTrackBuilder.getFeatureSource(RMDTrackBuilder.java:225)
    at org.broadinstitute.gatk.engine.refdata.tracks.RMDTrackBuilder.createInstanceOfTrack(RMDTrackBuilder.java:148)
    at org.broadinstitute.gatk.engine.datasources.rmd.ReferenceOrderedQueryDataPool.(ReferenceOrderedDataSource.java:208)
    at org.broadinstitute.gatk.engine.datasources.rmd.ReferenceOrderedDataSource.(ReferenceOrderedDataSource.java:88)
    at org.broadinstitute.gatk.engine.GenomeAnalysisEngine.getReferenceOrderedDataSources(GenomeAnalysisEngine.java:990)
    at org.broadinstitute.gatk.engine.GenomeAnalysisEngine.initializeDataSources(GenomeAnalysisEngine.java:772)
    at org.broadinstitute.gatk.engine.GenomeAnalysisEngine.execute(GenomeAnalysisEngine.java:285)
    at org.broadinstitute.gatk.engine.CommandLineExecutable.execute(CommandLineExecutable.java:121)
    at org.broadinstitute.gatk.utils.commandline.CommandLineProgram.start(CommandLineProgram.java:248)
    at org.broadinstitute.gatk.utils.commandline.CommandLineProgram.start(CommandLineProgram.java:155)
    at org.broadinstitute.gatk.engine.CommandLineGATK.main(CommandLineGATK.java:107)

    ERROR ------------------------------------------------------------------------------------------
    ERROR A GATK RUNTIME ERROR has occurred (version 3.2-2-gec30cee):
    ERROR
    ERROR This might be a bug. Please check the documentation guide to see if this is a known problem.
    ERROR If not, please post the error message, with stack trace, to the GATK forum.
    ERROR Visit our website and forum for extensive documentation and answers to
    ERROR commonly asked questions http://www.broadinstitute.org/gatk
    ERROR
    ERROR MESSAGE: For input string: "R"
    ERROR --------------------------------
  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie admin

    OK, so this tells us the error occurs while parsing a VCF header. Since you're only passing in the one VCF file, SPARK.common.vcf, the error must be in that one. It looks like one of the annotation field definitions may be responsible; there is a "R" character where the program expects to see a number. If you can't find it yourself, post the file header here and we'll try to help you identify it.

  • aligalig Member

    Thank you, I have found the R character in one of the annotation fields & replaced with a number, alig

Sign In or Register to comment.