cblaster icon indicating copy to clipboard operation
cblaster copied to clipboard

Diamond database not made when making local database

Open renVer01 opened this issue 10 months ago • 5 comments

Hi,

I was trying to make a local database with the cblaster makedb command from genbank files. Everything seems to run fine, except no diamond database is made. I saw in another thread that an incorrect installation of the diamond package might be the problem, but after reinstalling I still have the same issue.

Do you have any idea what else could be the problem?

Thanks!

renVer01 avatar Apr 04 '25 12:04 renVer01

I am having the same issue.

Kiratrea avatar Apr 08 '25 22:04 Kiratrea

Ok I think I solved it. I had to check my fasta file for gaps. Once I got rid of the gaps I was able to make the diamond database directly with diamond and use it with cblaster going forward.

Kiratrea avatar Apr 08 '25 22:04 Kiratrea

Having the same results also. It looks like that the fasta file output from the command contained some empty sequences and diamond could not create the database. Cblaster did not handle the error from diamond correctly and therefore there was no database created.

fjell-dev avatar Apr 09 '25 06:04 fjell-dev

Yes, removing empty sequences fixed the issue for me as well. Thanks!

renVer01 avatar Apr 09 '25 09:04 renVer01

I used this command to remove the empty sequences in my large FASTA file sed -i 'N;/\n$/d' large_file.fasta (use at your own risk) and use the diamond makedb command successfully. The cblaster search command worked with the "manually" created diamond database, apparently.

fjell-dev avatar Apr 09 '25 09:04 fjell-dev

I have encountered the same problem. cblaster makedb outputs the faster file and sqlite file but not the dmnd file. I did as fjell-dev suggested but using the command sed 'N;/\n$/d' [large_file.fasta] > [large.fasta].

[large_file.fasta] is the fasta output of cblaster [large.fasta] is the input for diamond after cleaning the empty lines, as follows: diamond makedb --in large.fasta --db myDb.dmnd

I tested cblaster search with myDb.dmnd and all worked well

mvalvano avatar Sep 29 '25 13:09 mvalvano