Question: Is it possible to use VariantRecalibrator with only one sample? (90,000 MuTect2 calls)
0
gravatar for wiki97
6 weeks ago by
wiki970
wiki970 wrote:

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 :)

gatk mutect2 vqsr • 144 views
ADD COMMENTlink modified 5 weeks ago • written 6 weeks ago by wiki970

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

ADD REPLYlink written 5 weeks ago by wiki970
0
gravatar for geneart$$
5 weeks ago by
geneart$$40
United States
geneart$$40 wrote:

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 COMMENTlink modified 5 weeks ago • written 5 weeks ago by geneart$$40
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1285 users visited in the last hour