This is a test version of Biostars. For the public version, visit https://www.biostars.org.
VariantRecalibrator in GATK

Excuse me:

I used unifiedGenotyper in GATK to call variants and filterd variants with VQSR. There were errors:

java -jar GenomeAnalysisTK.jar \ 
    -T ApplyRecalibration \ 
    -R reference.fa \ 
    -input raw_variants.vcf \ 
    -mode SNP \ 
    --ts_filter_level 99.0 \ 
    -recalFile recalibrate_SNP.recal \ 
    -tranchesFile recalibrate_SNP.tranches \ 
    -o recalibrated_snps_raw_indels.vcf

The log said

##### ERROR MESSAGE: Argument with name '--resource' (-resource) is missing.
##### ERROR Argument with name '--use_annotation' (-an) is missing.

Was the error in UnifiedGenotyper? When I used HaplotypeCaller, it did not print out error.

Any thing that helped me understand the issue would be much appreciated.

alignment sequence next-gen-sequencing
• 4,307 views • link
updated by Ram Ram • written by 897598644 897598644

I have found the error.

Thanks again.

• 0 views • link
updated by Ram Ram • written by 897598644 897598644

This is not an answer to the original question, and it carries no information for other's who come upon this question.

For the benefit of the community, you must tell us what you found, or fixed.

• 0 views • link
updated by Ram Ram • written by karl.stamm karl.stamm

Thanks to Karl's suggestion. It was the path error that I did not find.

• 0 views • link
updated by Ram Ram • written by 897598644 897598644

1 answer

You need to specify the truth and training datasets, as well as the annotations, when you run VariantRecalibrator. You build the model under your set of criteria and then apply it to your dataset.

For more help, check out this GATK guide.

• 0 views • link
updated by Ram Ram • written by Brice Sarver Brice Sarver

Log in to answer this question.

More posts like this