Genome Analysis Toolkit

Variant Discovery in High-Throughput Sequencing Data

GATK process banner

Need Help?

Search our documentation

Community Forum

Hi, How can we help?

Developed in the Data Sciences Platform at the Broad Institute, the toolkit offers a wide variety of tools with a primary focus on variant discovery and genotyping. Its powerful processing engine and high-performance computing features make it capable of taking on projects of any size. Learn more

Can't use GATK or Picard's LiftoverVcf tool without encountering output format exception

0

5 comments

  • Avatar
    Bhanu Gandham

    Hi, 

     

    It looks like you provided a txt file in the REJECT argument, which should be a vcf file. Try

    REJECT=rejectedRecords.vcf
    1
    Comment actions Permalink
  • Avatar
    A S

     

    Bhanu Gandham Great, thanks, that seems to have worked. But now I am getting another error:

    
    java.lang.IllegalStateException: Key LowQual found in VariantContext field FILTER at 1:10014 but this key isn't defined in the VCFHeader. We require all VCFs to have complete VCF headers by default.
    at htsjdk.variant.vcf.VCFEncoder.fieldIsMissingFromHeaderError(VCFEncoder.java:202)
    at htsjdk.variant.vcf.VCFEncoder.getFilterString(VCFEncoder.java:184)
    at htsjdk.variant.vcf.VCFEncoder.write(VCFEncoder.java:135)
    at htsjdk.variant.variantcontext.writer.VCFWriter.add(VCFWriter.java:248)
    at picard.vcf.LiftoverVcf.rejectVariant(LiftoverVcf.java:456)
    at picard.vcf.LiftoverVcf.doWork(LiftoverVcf.java:357)
    at picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:295)
    at org.broadinstitute.hellbender.cmdline.PicardCommandLineProgramExecutor.instanceMain(PicardCommandLineProgramExecutor.java:25)
    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)
    

     

    I think this is because many of the SNPs have a "LowQual" value in one of the fields like so:

    1 6512957 . A . 26.08 LowQual AC=0;AF=0.00;AN=2;DP=21;MQ=28.76;MQ0=4;1000gALT=A;AF1000g=0.01;AFR_AF=0.03;RM;TS=HPGOM;TSseq=A,A,A,A,A;CAnc=A;GAnc=A;OAnc=A;bSC=788;pSC=0.026;GRP=-0.854; Map20=0 GT:DP:GQ:PL:A:C:G:T:IR 0/0:21:1.48:4,0,362:2,18:0,0:0,0:0,1:0

    Now what do I do about this problem? Should I remove the "LowQual" value? The VCF is low quality because that is the nature of the sample; it is from the remains of an Altai Neanderthal.

     

    0
    Comment actions Permalink
  • Avatar
    Bhanu Gandham

    This looks like a malformed vcf. I recommend either recreating this vcf or running some validations tests on it using https://gatk.broadinstitute.org/hc/en-us/articles/360036347892-ValidateVariants

    0
    Comment actions Permalink
  • Avatar
    Brian Wiley

    Anyway the error message can indicate which it is for.  People will see this and think error to -O param when it is the --REJECT param is issue.

    0
    Comment actions Permalink
  • Avatar
    Genevieve Brandt (she/her)

    Hi Brian,

    Thank you for the feedback. In this case, the reject file is an output file type, and so that is why the error message is as written.

    Thanks,

    Genevieve

    0
    Comment actions Permalink

Please sign in to leave a comment.

Powered by Zendesk