opal
opal copied to clipboard
SIMD C/C++ library for massive optimal sequence alignment (local/SW, infix, overlap, global)
Hi from a Martin to another! I came to work on Opal since it's used in STAR, and I'm trying to improve the portability of some bioinformatics tools to new...
Hi, Below example that is causing segfault: ``` ./opal_aligner -o 1 -e 1 -x 2 -a NW -p test_query.fasta test_db.fasta -f test.mat ``` test_query.fasta: ``` >test query AABB ``` test_db.fasta:...
I got this error while trying to execute make > g++ -std=c++11 -march=native -Wall -O3 -c opal.cpp > opal.cpp:1169:45: error: no member named 'abs' in namespace 'std'; did you mean...
I should improve alignment speed. One thing I can do for sure is make alignment algorithm banded. I should also investigate if there are any other ways to make algorithm...
Information about E-value: http://blast.ncbi.nlm.nih.gov/Blast.cgi?CMD=Web&PAGE_TYPE=BlastDocs&DOC_TYPE=FAQ#expect . Implementation of E-value in SW#: http://sourceforge.net/p/swsharp/code/ci/master/tree/swsharp/src/evalue.cu .
Could I add SHW mode? Ivan Sovic said it would be useful.
Check this out and fix it, on match_bonus branch. I tried with gap open 0 and extend 1, with gap open 2 and gap extend 3, and I got errors...
Sometimes user may be confident that score will be in some range. In that case, I could allow him to specify for opal that overflow should not be checked. That...
This method is explained in Rognes's article, under "Processing four consecutive cells along the database sequences". They claim that by processing multiple cells they have better usage of SSE registers...
fixing #6 and perform case insentive comparison of matrix names