Benjamin Buchfink
Benjamin Buchfink
Hi Andreas, unfortunately there is no error message from Diamond in these logs, only this: ``` + diamond blastp --threads 1 --db test/C.faa.diamond --query test/E.faa --out /tmp/tmp.dUhbhJeGjW/E.faa.vs.C.faa.diamond -e 1e-05 --outfmt...
Hi @tillea , I was able to track down the issue using the travis ci. Will post a new release in the next days.
Hi @tillea my new attempt: https://github.com/bbuchfink/diamond/releases/tag/v0.9.35
Last time the big endian architectures failed, which are working now, so its not related to the endianess any more. I'll see if I can do something about the mips64el...
Thanks for posting the backtrace, but unfortunately it's not clear to me why this crash happens. It does not occur on other platforms, including when compiled with an address sanitizer....
The taxonomic format does not output alignments. It merges all alignments for one query into one classification using the LCA algorithm.
It's possible that the targets they align against have no taxonomy assignment, you can check for that.
No, the NCBI files do not map Uniprot accessions. I'll have to look into providing taxonomy mapping for uniprot.
The nr database would be the most obvious choice for an extensive database, and I think most sequences should have a taxonomic assignment. AnnoTree could also be an option, see...
Diamond was not designed for such short sequences. With some tuning of the parameters that may work but also be pretty slow if you use such a big database as...