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.

Why GATK 3.7 VariantAnnotator - Processed 0 loci

I was trying to use GATK 3.7 VariantAnnotator to annotate a vcf.gz file. Here is my code:

/u/local/apps/java/jdk1.8.0_111/bin/java -jar /u/local/apps/gatk/3.7/GenomeAnalysisTK.jar \
-T VariantsToBinaryPed \
-V ${NAME}_dbsnp-GRCh37.vcf \
-m ${HOME}/raw/0.2_SCZ.fam \
-bed ${NAME}.bed \
-bim ${NAME}.bim \
-fam ${NAME}.fam \
-mgq 0 \

But the output file was extremely small, and then in the error log I found the following:

Note that "INFO 15:07:00,074 VariantAnnotator - Processed 0 loci." at the end of all the records of ProgressMeter.

How shall I fix this?


  • Sorry, I was using the following script:

    /u/local/apps/java/jdk1.8.0_111/bin/java -jar /u/local/apps/gatk/3.7/GenomeAnalysisTK.jar \
    -R $FASTA \
    -T VariantAnnotator \
    -V $VCF \
    -o ${NAME}_dbsnp-GRCh37.vcf \
    --alwaysAppendDbsnpId \
    --dbsnp $DBSNP \

  • SheilaSheila Broad InstituteMember, Broadie, Moderator admin

    Hi Jenny,

    Can you check the output VCF and make sure the dbsnp rsid field is populated?


Sign In or Register to comment.