makeblastdb fails with full database
Hello, I downloaded the full PR2 database as a fasta file (" pr2_version_5.0.0_SSU_taxo_long.fasta") which according to the website is suitable for making a local BLAST database. However, when I run makeblastdb using BLAST+ version 2.6.0, it fails to generate a database.
This command
makeblastdb -in pr2_version_5.0.0_SSU_taxo_long.fasta -dbtype nucl -title "PR2 Full Database"
generates only .nin, .nhr, and .nsq files, and shows the following output:
Building a new DB, current time: 10/18/2023 14:05:45 New DB name: /Users/nastassiapatin/DBs/PR2/PR2_full-MagicBlastDB/pr2_version_5.0.0_SSU_taxo_long.fasta New DB title: PR2 Full Database Sequence type: Nucleotide Keep MBits: T Maximum file size: 1000000000B Adding sequences from FASTA; added 221085 sequences in 7.21831 seconds.
While this command (with -parse_seqids)
makeblastdb -in pr2_version_5.0.0_SSU_taxo_long.fasta -dbtype nucl -parse_seqids -title "PR2 Full Database"
also generates only .nin, .nhr, and .nsq files, but these seem to be incomplete as they have ".00" before the .n** suffixes. This command produces the following output:
Building a new DB, current time: 10/18/2023 14:04:02 New DB name: /Users/nastassiapatin/DBs/PR2/PR2_full-MagicBlastDB/pr2_version_5.0.0_SSU_taxo_long.fasta New DB title: PR2 Full Database Sequence type: Nucleotide Keep MBits: T Maximum file size: 1000000000B [1] 6917 abort makeblastdb -in pr2_version_5.0.0_SSU_taxo_long.fasta -dbtype nucl -title
Do you have any ideas why the full PR2 fasta might be failing to build a BLAST database? When I build one from only a select group of taxa (order Dinophyceae) it works fine.
Thanks, Nastassia