using GATK IndelRealigner

guillaumers · 07-20-2016, 01:52 AM

Hi,

I have some difficulties using IndelRealigner on my data. After comparing bam files before and after realignement process I find some regions where a correction should have occured. (an igv capture of an ambiguous region that isnt resolved after realignement : https://postimg.org/image/7yuy0hz81/ ). Is it normal that all my indels "problems" are not solved? And are there parameters to adjust to solve them? Have I made mistakes running the commands?

(here are my command lines ) java -Xmx12000m -jar ~/Soft/GATK/GenomeAnalysisTK.jar -T RealignerTargetCreator -nt 15 -R ~/igv/genomes/genome.fa -I sample.bam -o temp.indel.intervals java -Xmx12000m -jar ~/Soft/GATK/GenomeAnalysisTK.jar -T IndelRealigner -R ~/genome.fa -I sample.bam -targetIntervals temp.indel.intervals -o sample_realigned.bam

I'm willing to take any ideas!

tristan dubos · 07-20-2016, 02:45 AM

Hi Guillaume,

The correction appends if there is a better alignment possible, how do you know you should have a different alignment in this position ?

Tristan

wdecoster · 07-23-2016, 03:28 PM

Not an answer to your question, but interesting read considering IndelRealigner: https://software.broadinstitute.org/gatk/blog?id=7847

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
GATK IndelRealigner error	shocker8786	Bioinformatics	1	12-21-2014 07:56 AM
GATK IndelRealigner error	kga1978	Bioinformatics	7	12-16-2011 09:58 AM
GATK IndelRealigner error	Heisman	Bioinformatics	2	09-14-2011 04:25 PM
GATK IndelRealigner error	carljo	Bioinformatics	1	04-06-2011 12:25 PM
GATK IndelRealigner Error:Code Exception	HongxiangZheng	Bioinformatics	1	12-14-2010 08:21 PM

07-20-2016, 01:52 AM	#1
guillaumers Junior Member Location: Colmar (France) Join Date: Oct 2014 Posts: 1	using GATK IndelRealigner Hi, I have some difficulties using IndelRealigner on my data. After comparing bam files before and after realignement process I find some regions where a correction should have occured. (an igv capture of an ambiguous region that isnt resolved after realignement : https://postimg.org/image/7yuy0hz81/ ). Is it normal that all my indels "problems" are not solved? And are there parameters to adjust to solve them? Have I made mistakes running the commands? (here are my command lines ) java -Xmx12000m -jar ~/Soft/GATK/GenomeAnalysisTK.jar -T RealignerTargetCreator -nt 15 -R ~/igv/genomes/genome.fa -I sample.bam -o temp.indel.intervals java -Xmx12000m -jar ~/Soft/GATK/GenomeAnalysisTK.jar -T IndelRealigner -R ~/genome.fa -I sample.bam -targetIntervals temp.indel.intervals -o sample_realigned.bam I'm willing to take any ideas!

07-20-2016, 02:45 AM	#2
tristan dubos Member Location: France Join Date: Dec 2015 Posts: 39	Hi Guillaume, The correction appends if there is a better alignment possible, how do you know you should have a different alignment in this position ? Tristan

07-23-2016, 03:28 PM	#3
wdecoster Member Location: Antwerp, Belgium Join Date: Oct 2015 Posts: 97	Not an answer to your question, but interesting read considering IndelRealigner: https://software.broadinstitute.org/gatk/blog?id=7847