This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How Do You Do Local Realignment Around Indels For Exome Sequencing?

Hi all,

I am wondering how people generally go about local realignment around indels in their workflow?

The most popular programs seem to be SRMA and GATK - does anyone use anything else?

In my initial tests this seems to be a VERY time-consuming step when compared to the initial alignment step. Does anyone employ methods of speeding up the process? For example, only realigning within exonic coordinates OR perhaps performing variant calling first and then using the coordinates from the VCF to restrict the areas where realignment is performed?

I am trying to set up a SNP/Indel workflow and realignment seems to slow it to the the point where it becomes useless for anything but small numbers of samples. I tested on a single chromosome and it took 4.5 hours. Alignment and variant calling without realignment only took about 1.5 hours for the whole genome!

Thanks in advance!

snp indel next-gen sequencing alignment mapping
• 6,289 views • link
updated by Docroberson • Docroberson • written by Travis Travis

Did you find anything new in past months? I am now facing similiar situation as you were, considering which tool (if any) to use. Thanks.

• 0 views • link
written by Biomonika (Noolean) Biomonika (Noolean)

2 answers

Take a look to this question and my answer (and the comments). SRMA allows realignment of specific regions. That will help you reduce the runtime of the process.

• 2 views • link
updated by Ram Ram • written by Drio Drio

I remember the question - I asked it :) My question here is a more general one.

• 0 views • link
written by Travis Travis

I typically use BWA for alignments to get indels at the same time. Pretty fast (on par with bowtie). Alternatively I've used NOVOALIGN to do alignment and indel calculations in one step also. Very good at finding both substitutions and indels, but VERY slow and resource intensive compared to other software.

• 1 views • link
written by Docroberson • Docroberson

I am more curious about the local realignment step that occurs post BWA etc.

• 0 views • link
written by Travis Travis

Log in to answer this question.

More posts like this