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

GATK working, but -ERC command not found

Answered
0

3 comments

  • Avatar
    Genevieve Brandt (she/her)

    Hi Ann W,

    Thanks for trying out GATK! Thankfully, I think this error can be easily solved. It looks like it's just a small typo in your command line so that the -ERC line is read as a different command than part of the GATK command line. 

    Look for a small issue in your script, potentially a newline character or something else after the -O $HaplotypeCaller_dir/$tm.g.vcf.gz argument and before the -ERC GVCF argument. 

    Let me know if you have any other questions.

    Best,

    Genevieve

    1
    Comment actions Permalink
  • Avatar
    Ann W

    That was the issue! Thanks so much, Gen! For reference, here is the new bash script that worked. 

    #!/bin/bash
    lee=/pine/scr/a/n/anndo1/Lee_NatureGenetics_2017_Germline
    align_dir=$lee/bulk_align
    HaplotypeCaller_dir=$lee/HaplotypeCaller

    tm=blah_case_blah    

    module load gatk
    gatk --java-options "-Xmx4g" HaplotypeCaller \
        -R /proj/yuchaojlab/lib/ucsc.hg19.fasta \
        -I $align_dir/*$tm.sorted.rg.dedup.realigned.recal.bam \
        -O $HaplotypeCaller_dir/$tm.g.vcf.gz \
        -ERC GVCF

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

    Thank you for the update, glad it worked!

    0
    Comment actions Permalink

Please sign in to leave a comment.

Powered by Zendesk