Error When Using Makeblastdb On A Large File

Question

[vijay]

Dear All,

I have created a stand alone blast on my system, to work with NGS data and it is working fine. However, I have got a sequence file which is about 2400MB, which is two large to create a local database and I am facing a error when i use the makeblastdb command.

Please can someone advise me on how to create a local database out of this file.

-Vijay

[Istvan Albert]

please indicate the command that you issue and the error that you are getting

Answer 1

[Sujai Kumar]

I agree that the OP should specify the full command and error (thanks Istvan) - but just in case this helps:

makeblastdb gives an error if there are any tabs or other weird characters in the fasta file. Check your sequence headers, and replace all weird spaces with ordinary spaces.

I use the following:

perl -i -plne 's/\s+/ /g' sequence.fasta

Answer 2

[neal.platt]

That is not an overly large file to make a blastdb from. Most vertebrate genomes are that size...granted you will have millions of reads compared to dozens of chromosomes.

I tested making a blastdb on a short read file of similar size and it was fine. Nothing special..

formatdb -pF -i short.read.fa

What type of error are you getting specifically?