-
Notifications
You must be signed in to change notification settings - Fork 2
Description
Hello,
I am trying to run taxor and taxor profile gives a core dump. Created a custom hixf with refseq archaea
taxor build --input-file refseq-a/2025-11-17_14-13-41/refseq_accessions_taxonomy.csv --input-sequence-dir refseq-a/2025-11-17_14-13-41/files --output-filename refseq-a-k22-s12.hixf --threads 6 --kmer-size 22 --syncmer-size 12 --use-syncmer
Steps to recreate:
taxor profile --search-file SAMPLE.search.txt --cami-report-file SAMPLE.report \ --seq-abundance-file SAMPLE.abundance --binning-file SAMPLE.binning --sample-id SAMPLE
Error:
Parsing search results...done
Remove matches to nonunique references...done
Remove low confidence references...done
Filter associated references...done
Start EM algorithm...started
Initialize prior probabilities ...done
Starting EM iteration 0
Calculate Log Likelihoods ...done
Update prior probabilities ...done
final update...
done
done
Starting EM iteration 1
Calculate Log Likelihoods ...done
Update prior probabilities ...done
final update...
done
done
Number of EM steps needed: 1
done
Calculate higher rank sequence abundances..terminate called after throwing an instance of 'std::out_of_range'
what(): basic_string::substr: __pos (which is 3) > this->size() (which is 0)
Aborted (core dumped)
Versions:
taxor version: 0.2.0
SeqAn version: 3.4.0-rc.1
I downloaded version 0.2.1 but the --version command displays the above version
Thanks in advance!
-Prottasha