I keep getting a segmentation fault message whenever trying to run a tblastn command:
tblastn -db genome.fna -query protein.fas -out protein.out -num_threads 14 -outfmt 7.
I have done the following so far: I downloaded the latest linux version of BLAST+ from here https://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/LATEST/ , compiled and installed the program. Afterwards I created a database from a somewhat large genome (over 15 gbp). The database consisted of two files as it was larger than 4GB.
makeblastdb -in genome.fna -parse_seqids -blastdb_version 5 -title "genome" -dbtype nucl -max_file_sz 4GB
I then tried to use tblastn with the following command
tblastn -db genome.fna -query protein.fas -out protein.out -num_threads 14 -outfmt 7
, which resulted in the aforementioned segmentation fault. (blastn runs normally) tblastn works fine online, when searching for the same protein and same genome used in the above command line.
The computer on which the commands are running on has a 128GB RAM and 14 core processor, so I doubt the hardware is to blame.
I wonder what the cause of this error could be.
Thank you in advance.