Abra for indel realignment
Hi all,
Is anyone using Abra for realignment of exome/wgs data? Looks like a nice tool, and paper looks good too. I'd be interested in hearing peoples experience of it.
Repo: https://github.com/mozack/abra
cheers
realignment
abra
gatk
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Recommended Tool for Somatic SV Discovery in Tumor-Only WES Data
written by George •Hi all, I'm currently working on an analysis using tumor-only whole exome sequencing (WES) data and am looking for a reliable tool for structural variant …
-
how to run the compare_genomes for comparative analysis
written by rj.rezwan •Hi, I am interested to compare the 10 genomes for comparative analysis. I have genome, cds, protein and gff files for this analysis. I want …
-
Exon phase visualisation tools
written by oakhamwolf •Hi all Hope someone can help. I'm looking for a way to visualise exon phase like [here][1] Ultimately I'd like to be able to ingest …
-
autosomal ancestry inference
written by J.F.JiangHi all, I am very interested in ancestry inference. 1. autosomal ancestry I've followed the tutorial from https://github.com/pcgoddard/Burchardlab_Tutorials/wiki/ADMIXTURE#popfile. However, decision on cluster information K is …
-
Forum: Interesting/useful/beautiful bio sequence viz tools
written by JoeBelow is the previous post content, but I've just reformatted all this to a list now for easier future reference. > Looking for suggestions for …
-
Processing Whole Genome Sequencing data (shallow) and call Copy Number Variations
written by Alejandro Jimenez SanchezDear all, We are proposing a project where we want to perform shallow Whole Genome Sequencing (<1x) on DNA obtained from tumour FFPE human samples. …
-
Whole Genome Sequences somatic MAF files at TCGA
written by always_learningHi All, I was trying to get WGS MAF files from TCGA. I was going through this wonderful post here at https://www.biostars.org/p/69222/ and it looks …
-
BamUtil for Base Qual Recalibration
written by DanielkHi all, Does anyone have experience with or has benchmarked BamUtil for Base Quality Score Recalibration of NGS data? Looks like an interesting alternative to …
-
Alternative to Sage Synapse
written by DanielkHi all, [Sage Synapse][1] is a very nice tool for data sharing within project as well as making data public, like in the TCGA pan …
-
Srma (Short Read Micro-Aligner) Optimising For Speed
written by Travis<p>Hi all, </p> <p>I am wondering if anyone has experience using this package for local realignment of reads? I am wondering how best to tweak …
Just a thought, but... if you align with BBMap, you won't need to realign indels, because they'll be correctly aligned already.
Brian Bushnell This comment is not constructive in its current form. Why would BBMap help? What does BBMap have to do with the original question?
As I understand it, Abra is a tool that takes mapped reads and realigns them around indels, because the reads originating over the indel were not aligned correctly. Locally realigning over suspected indels can help improve alignments, if the mapping (approximate location) is correct. However, it will not place reads on indels if they were mapped somewhere else in the first place, because the mapper could not tolerate long indels and thus did not consider that a valid mapping location. So mapping with a program that cannot tolerate long indels, and subsequently realigning with a tool that can, will generally yield inferior results to mapping and aligning with a single tool that does tolerate long indels... because in the second case, there will be a lower rate of reads being mismapped or left unmapped.
If the goal of the thread is solely to discuss Abra, then my comment is irrelevant. But if the goal is to find tools that do a good job of aligning reads around indels, then it is relevant.
We've had good results using Abra. Of course, we developed it, so my opinion might be biased :)
We use it extensively for cancer gene panels and exomes. We've not used it heavily for WGS. In our WGS test case, it was necessary to split processing up by chromosome.
There is also now a somatic indel caller that specifically works with Abra realignments. Will post details on this on the Abra website soon.