Test-drive the GATK tools and Best Practices pipelines on Terra
Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.
BQSR with intervals & cut-off base coverage or interval average coverage
if I'm correct:
when I process in the BQSR an interval with very low coverage (0x - 3x base coverage), then I have the error "RecalibrationReport - Missing read group(s) ..." if I go forward to the HaplotypeCaller I have a new error "USER ERROR ... Interval/contig... loci" ...
When the coverage is high everything is ok.
If the above "relation" between errors and steps is correct... what happen when I have two target regions inside an interval, one with a high coverage and the other one with a very low coverage? I lose both regions data or only of the not well covered?
What is the cut-off of reads/coverage/base coverage of BQSR, or average coverage of the inteval?