Comments (5)
Yes, ABRA will likely have problems with the lack of sequence diversity also. We do have an update in progress that should help with this. If you are feeling adventurous, we could try get you going on an alpha version of the new software.
from abra.
Thanks for your reply, I would be very grateful to try the new functionality. I have a specific case where the variant caller is unable to find a 20bp deletion but it can be seen clearly in igv with soft clipping enabled. How do I get my hands on the new version.
Thanks
Matt
from abra.
An alpha release of ABRA2 is available. See: https://github.com/mozack/abra2/releases
Download and extract abra2_0.1.tgz
Sample command line:
JAR=abra2/release/test/abra2_0.1-SNAPSHOT-jar-with-dependencies.jar
# Path to directory containing the .so files
ABRA_LIBDIR=abra2/release/test
REF=ref/hg19.fa
NORMAL=normal.sort.bam
TUMOR=tumor.sort.bam
BED=uncseq5.bed
java -Xmx4G -Djava.library.path=$ABRA_LIBDIR -jar $JAR --ref $REF --in $NORMAL,$TUMOR --out normal.abra.bam,tumor.abra.bam --working abra_temp_dir --targets $BED --threads 8 --sa > abra2.log 2>&1
from abra.
Thanks Lisle ABRA worked well. See IGV screenshot below. Top pane is ABRA, bottom is original bwa.
Could you explain how the amplicon function works please? Are you still using a de bruijn graph?
Thanks
Matt
from abra.
The --sa param you specified disables the de Bruijn graph assembly. In ABRA2, we've introduced local Smith Waterman alignment for reads with high base quality soft clipping. Optimal use cases for one method over the other (or both combined) has not yet been determined (although cases such as yours will clearly be better served by Smith Waterman).
from abra.
Related Issues (20)
- Excessive run time HOT 2
- [fwrite] Remote I/O error HOT 5
- polyA/polyT/lowcomplexity region realignments for single ended reads HOT 3
- Exception in thread Error HOT 2
- Error in Cadabra HOT 4
- [Request] Replace System.out.print* statements with System.err.print* for stream compatibility HOT 3
- xlC compatibility - C++ STL referenced in .c files HOT 4
- hidden dependency on bwa mem HOT 2
- Abra error HOT 2
- Change header in bam HOT 1
- Is YA tag deterministic? HOT 4
- Abra will realign duplicates? HOT 5
- Testing practices HOT 2
- Abra processes unique reads? HOT 2
- Is it possible to get the assembled contig? HOT 2
- No contigs assembled: no space left on device HOT 1
- abra2 error java.lang.NumberFormatException HOT 3
- build issue
- java.lang.IllegalStateException
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from abra.