Is it possible to use VariantRecalibrator with only one sample? (90,000 MuTect2 calls)
1
0
Entering edit mode
9 months ago
wiki97 • 0

Hi, I want to filter variants (called with MuTect2) with VariantRecalibrator and ApplyVQSR GATK tools.

However, I have only one sample, with 90,000 variants called. Would this be enough for GATK VariantRecalibrator and ApplyVQSR tools to apply a good model and properly filter my variants?

I will be grateful for any suggestions :)

MuTect2 VQSR GATK • 376 views
ADD COMMENT
0
Entering edit mode

Hi geneart$$, thank you very much for your answer!

ADD REPLY
0
Entering edit mode
9 months ago
geneart$$ ▴ 40

Hi wiki, VQSR is not recommended to be used when you have less number of samples and definitely not when you have 1 sample. Please see :https://gatk.broadinstitute.org/hc/en-us/articles/360037594511 However , you can pad the samples using samples from good database source like 1000G. GATK explains how to do that. One thing though you may want to dig more into their documentation as they have it all there but sort of needle in the haystack.

Hope this helps !

ADD COMMENT

Login before adding your answer.

Traffic: 2549 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6