Question: Ngs: Indel Analysis Tool Order
0
gravatar for David K Crossman
7.3 years ago by
United States
David K Crossman130 wrote:
Hello! I have a question about the NGS: Indel Analysis toolset in Galaxy. I have aligned my samples from Illumina's HiSeq2000 to the reference genome using BWA. I've called SNPs using SAMTools and now need to call indels. Under the NGS: Indel Analysis toolset, I see two options: "Filter Indels for SAM" and "Extract indels from SAM." Since the descriptions are a little vague, I would presume that I start with the "Filter Indels for SAM" first on the BWA SAM file and then run "Extract indels from SAM" on the "Filter" output SAM file (i.e. BWA SAM -> Filter Indels for SAM -> Extract indels from SAM). Is this the correct order? Or would I skip the "Filter Indels for SAM" step (since the BWA SAM file technically already contains indels so there would be no need to filter) and just go straight to "Extract indels from SAM" (i.e. BWA SAM -> Extract indels from SAM)? I've tried both ways and get different results. For example: 1. BWA SAM -> Filter Indels for SAM -> Extract indels from SAM -> gives me 26,417 regions of interest 2. BWA SAM -> Extract indels from SAM -> gives me 93,974 regions of interest Which way is correct? Any help/info would be greatly appreciated! Thanks, David
bwa alignment samtools bam • 1.4k views
ADD COMMENTlink modified 7.3 years ago by Jennifer Hillman Jackson25k • written 7.3 years ago by David K Crossman130
0
gravatar for Jennifer Hillman Jackson
7.3 years ago by
United States
Jennifer Hillman Jackson25k wrote:
Hello, The "Filter Indels" tool is intended to act as a quality screen before "Extract Indels" is used. Whether you want to use this step is something you will need to determine, as are the exact criteria for filtering. The best advice is to examine a sample of the differences between the two methods and tune or exclude the Filter step. Hopefully this helps, Best, Jen Galaxy team -- Jennifer Jackson http://usegalaxy.org http://galaxyproject.org/Support
ADD COMMENTlink written 7.3 years ago by Jennifer Hillman Jackson25k
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 16.09
Traffic: 175 users visited in the last hour