Notice:
If you happen to see a question you know the answer to, please do chime in and help your fellow community members. We appreciate your help!

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.

BaseRecalibrator only doing recalibrations on substitutions, not indels

Hello, I am trying out BaseRecalibrator. My command is:

gatk BaseRecalibrator \
      --reference reference.fa \
      --input input.bam \
      --known-sites good_filtered_snps.vcf \
      --known-sites good_filtered_indels.vcf \
      --output recal_table

The output recalibration table only include information for substitutions (event type "M"), none for insertions or deletions, which other people's recalibration tables seem to have. I cannot find anything in the BaseRecalibrator help or in its output that explains why I'm only getting a third of the recalibration table.

Thoughts?

Best Answer

Answers

Sign In or Register to comment.