Attention:
The frontline support team will be unavailable to answer questions on April 15th and 17th 2019. We will be back soon after. Thank you for your patience and we apologize for any inconvenience!

Indel Variant recalibration is showing an error

NeethuNeethu IndiaMember

I have done joint genotyping of 200 samples and I have a multi sample vcf file. SNP recalibration worked but when I did Indel quality score recalibration it shows an error "MESSAGE: Bad input: Values for QD annotation not detected for ANY training variant in the input callset. VariantAnnotator may be used to add these annotations"

I have used the below command
java -Xmx8g -XX:ParallelGCThreads=20 -jar /mnt/exome/Softwares/GenomeAnalysisTK.jar -T VariantRecalibrator -R /mnt/exome/ReferenceFiles/human_g1k_v37.fasta -input recalibrated_snps_raw_indels.vcf -resource:mills,known=true,training=true,truth=true,prior=12.0 /mnt/exome/ReferenceFiles/Resource/Mills_and_1000G_gold_standard.indels.b37.vcf -an QD -an DP -an FS -an SOR -an ReadPosRankSum -an MQRankSum -mode INDEL -tranche 100.0 -tranche 99.9 -tranche 99.0 -tranche 90.0 --maxGaussians 4 -recalFile recalibrate_INDEL.recal -tranchesFile recalibrate_INDEL.tranches -rscriptFile recalibrate_INDEL_plots.R

It would be really great if some one can help me to trouble shoot this.

Answers

  • SheilaSheila Broad InstituteMember, Broadie, Moderator admin
    edited April 2017

    @Neethu
    Hi,

    Can you try using dbSNP as well? It looks like you set the Mills resource file as known=true,training=true,truth=true which may be producing the error. Have a look at this article for information on what we recommend for indels in VQSR.

    -Sheila

Sign In or Register to comment.