As I continue to add steps to my SNP/Indel discovery workflow, the latest recommendation is local realignment around Indels using GATK following the initial alignment step. I have just commenced the step which generates the target intervals for realignment (RealignerTargetCreator) and it looks like it will take an hour to complete, with the realignment still required after that. My test data set is a single sample of approx 5 million paired end 100bp reads.
For an upcoming project, my plan is to run 150 similarly sized samples. Therefore the addition of such time-consuming steps will have a major impact on timelines. Can anyone with experience in this area comment on the time required for Indel realignment vs the benefits received? Is it worth it?
1 answer
Read realignment certainly appears to be an essential step. Do you have access to a computer cluster? Working through the samples in parallel batches would certainly be less painful.
As an aside, I have had a good experience using Nils Homer's SRMA package. Perhaps you could do a time trial. SRMA just works over the whole batch of reads, so could be faster?
Log in to answer this question.
-
DNAseq pipeline steps
written by shuksi1984 •
-
GATK; Variant calling RNA-Seq best practices: bam processing steps
written by user230613
-
IndelRealigner and RealignerTargetCreator for amplicon sequencing
written by kulvait
-
Coverage depth and breadth after recaliberation
written by nikkihathi •
-
Making a Consensus Sequence of a Gene
written by User000
-
Gatk Indel Realignment Taking Forever! Help
written by Ds89 •
-
Error When Running Gatk
written by huboqiang •
-
How Do You Do Local Realignment Around Indels For Exome Sequencing?
written by Travis
-
Filtering Ngs Genomic Alignments
written by Travis
-
Snp Indel Discovery - Critique A Proposed Workflow
written by Travis
