Attention:
The frontline support team will be slow on the forum because we are occupied with the GATK Workshop on March 21st and 22nd 2019. We will be back and more available to answer questions on the forum on March 25th 2019.

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 \
-R $FASTA \
-V ${NAME}_dbsnp-GRCh37.vcf \
-m ${HOME}/raw/0.2_SCZ.fam \
-bed ${NAME}.bed \
-bim ${NAME}.bim \
-fam ${NAME}.fam \
-mgq 0 \
-mode SNP_MAJOR

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?
Thanks,
Jenny

Answers

  • 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 \
    --unsafe LENIENT_VCF_PROCESSING

  • SheilaSheila Broad InstituteMember, Broadie, Moderator admin

    @Jennybear
    Hi Jenny,

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

    Thanks,
    Sheila

Sign In or Register to comment.