Gleeson-Lab / wxs_pipeline

Starting with BAMs and FASTQs, follow GATK 4.0 Best Practices up to generating a joint-genotyped VCF
1 stars 1 forks source link

Adjust VQSR if Given a Small amount of Data #11

Open brcopeland opened 2 years ago

brcopeland commented 2 years ago

For example as per https://gatk.broadinstitute.org/hc/en-us/articles/360035531112--How-to-Filter-variants-either-with-VQSR-or-by-hard-filtering, it can be necessary to reduce the --max-gaussians parameter. Even doing that, small numbers of variants can cause the step to fail entirely. One possible solution is that VQSR could be run only if the total size of the inputs is sufficiently large.