pr2database icon indicating copy to clipboard operation
pr2database copied to clipboard

makeblastdb fails with full database

Open nvpatin opened this issue 2 years ago • 0 comments

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

nvpatin avatar Oct 18 '23 21:10 nvpatin