We've moved!
This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!

BaseRecalibrator

ilariailaria ItalyMember

Hi,
I have a problem with the depth of coverage obtained after GATK workflow, which included indel realigner, baserecalibrator and unified genotyper. In the last vcf file the depth of coverage is very low (3-10), but if I check on IGV the bam file produced after alignment the reads count ranges between 70 and more than 100. Is it possible that the Baserecalibrator step reduce the read counts?
How can I manage this issue?
Thanks in advance

Answers

  • SheilaSheila Broad InstituteMember, Broadie ✭✭✭✭✭

    @ilaria
    Hi,

    Can you please post the vcf record and the IGV screenshot for the site?

    Also, we do not recommend using Unified Genotyper for variant calling. You should try Haplotype Caller.

    Thanks,
    Sheila

Sign In or Register to comment.