The current GATK version is 3.6-0
Examples: Monday, today, last week, Mar 26, 3/26/04

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Powered by Vanilla. Made with Bootstrap.

Any known problems with IndelRealigner code in gsa-stable?

evolvedmicrobeevolvedmicrobe MGHMember Posts: 16

Does anyone know of any known issues with the indelrealigner? The GATK is calling 1000s of SNPs on one genome I have due to bad realignments. It appears the target finder identifies regions of the genome that are essentially perfectly aligned, but when the realigner gets to these areas it remits the reads as a new alignment that is a train-wreck compared to what the re-aligner started with. Am going to investigate further but thought I would check if this rings any bells. I am using the flag -model USE_READS and have no known indels to work with.


Best Answer


  • Geraldine_VdAuweraGeraldine_VdAuwera Administrator, Dev Posts: 10,739 admin

    Hi there,

    IndelRealigner is a mature tool and as far as we no there are no outstanding issues. Have you checked that your reference is ok and matches your bam files?

    Geraldine Van der Auwera, PhD

  • evolvedmicrobeevolvedmicrobe MGHMember Posts: 16

    Yep, issue doesn't seem to be the reference. Unfortunately I can't reproduce the error too easily as if I only extract the region of the genome that appears to be problematic, the indel realigner then works fine on the smaller dataset. Will keep looking...

Sign In or Register to comment.