Question: How Do You Do Local Realignment Around Indels For Exome Sequencing?
3
gravatar for Travis
8.0 years ago by
Travis2.8k
USA
Travis2.8k wrote:

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!

ADD COMMENTlink modified 8.0 years ago by Docroberson280 • written 8.0 years ago by Travis2.8k

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.

ADD REPLYlink modified 6.7 years ago • written 6.7 years ago by Biomonika (Noolean)3.1k
1
gravatar for Drio
8.0 years ago by
Drio910
United States
Drio910 wrote:

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.

ADD COMMENTlink written 8.0 years ago by Drio910

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

ADD REPLYlink written 8.0 years ago by Travis2.8k
0
gravatar for Docroberson
8.0 years ago by
Docroberson280
the lab
Docroberson280 wrote:

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.

ADD COMMENTlink written 8.0 years ago by Docroberson280
1

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

ADD REPLYlink written 8.0 years ago by Travis2.8k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1631 users visited in the last hour