Gracefully exit when contig taxonomy doesn't find any hits in contig

[genomewalker]

Hi

Expected Behavior

When using taxonomy to annotate a contig, mmseqs2 should gracefully exit if there are not hits when doing the search

Current Behavior

At the moment it exits with an error when Query database size: 0
Full log here

Query database size: 0 type: Aminoacid
Estimated memory consumption: 151G
Target database size: 43032414 type: Aminoacid
Index table k-mer threshold: 163 at k-mer size 7 
Index table: counting k-mers
[=================================================================] 43.03M 4m 11s 993ms
Index table: Masked residues: 327874438
Index table: fill
[=================================================================] 43.03M 11m 31s 69ms
Index statistics
Entries:          10015666182
DB size:          67075 MB
Avg k-mer size:   7.824739
Top 10 k-mers
    DFEQLPH	65289
    NVPGWSP	65148
    FRYAFPS	64963
    RYYVLGW	64817
    VDRLNTQ	63812
    WRLDFLN	62949
    TVDGDFS	62609
    NKTDFVQ	61859
    QDWVQIP	61316
    LDGAYVP	59471
Time for index table init: 0h 16m 9s 488ms
Time for processing: 0h 16m 22s 538ms
lcaalign /dev/shm/tmp/contig-taxonomy/e77760e07c3167ef94f937365175d68b.tax.UniRef50/12595762918832005520/taxonomy_tmp/13584421311535647719/orfs_filter DB/uniref50-tax/UniRef50 /dev/shm/tmp/contig-taxonomy/e77760e07c3167ef94f937365175d68b.tax.UniRef50/12595762918832005520/taxonomy_tmp/13584421311535647719/tmp_taxonomy/4936078414881223269/tmp_hsp1/16551820070899017371/pref_0 /dev/shm/tmp/contig-taxonomy/e77760e07c3167ef94f937365175d68b.tax.UniRef50/12595762918832005520/taxonomy_tmp/13584421311535647719/tmp_taxonomy/4936078414881223269/first --sub-mat nucl:nucleotide.out,aa:blosum62.out -a 0 --alignment-mode 1 --alignment-output-mode 0 --wrapped-scoring 0 -e 1 --min-seq-id 0 --min-aln-len 0 --seq-id-mode 0 --alt-ali 0 -c 0 --cov-mode 0 --max-seq-len 65535 --comp-bias-corr 1 --max-rejected 5 --max-accept 30 --add-self-matches 0 --db-load-mode 0 --pca 1 --pcb 1.5 --score-bias 0 --realign 0 --realign-score-bias -0.2 --realign-max-seqs 2147483647 --gap-open nucl:5,aa:11 --gap-extend nucl:2,aa:1 --zdrop 40 --threads 32 --compressed 0 -v 3 

Input /dev/shm/tmp/contig-taxonomy/e77760e07c3167ef94f937365175d68b.tax.UniRef50/12595762918832005520/taxonomy_tmp/13584421311535647719/tmp_taxonomy/4936078414881223269/tmp_hsp1/16551820070899017371/pref_0 does not exist
Error: Alignment died
Error: First search died
Error: taxonomy died
Error: Search died

Many thanks
Antonio

Your Environment

Include as many relevant details about the environment you experienced the bug in.

• Git commit used (The string after "MMseqs Version:" when you execute MMseqs without any parameters): 5daca42
• Which MMseqs version was used (Statically-compiled, self-compiled, Homebrew, etc.): Self-compiled
• For self-compiled and Homebrew: Compiler and Cmake versions used and their invocation: GCC 7.5.0 cmake 3.10.2
• Server specifications (especially CPU support for AVX2/SSE and amount of system memory): AVX2/SSE supported, 256G
• Operating system and version: Ubuntu 18.04

[milot-mirdita]

Thanks, I think the issues should be fixed now.