BLAST Databases-Sapelo2

From Research Computing Center Wiki
Revision as of 15:19, 2 August 2024 by Jordan (talk | contribs)
Jump to navigation Jump to search


Blast databases

Blast databases are pre-formatted to work with blast commands from BLAST and BLAST+. These databases are similar to what you could build using the makeblastdb command with FASTA files. Each database is located at a directory with the name pattern /db/ncbiblast/YYYYMMDD. For example, databases downloaded on February 01, 2022 will be available in /db/ncbiblast/20220201. Starting with the download on August 2nd, 2024 (i.e., /db/ncbiblast/20240802), the following databases are currently available and stored together with databases downloaded at the same time:

nr
nt
taxdb
refseq_rna
refseq_protein
mouse_genome
human_genome
swissprot
cdd_delta
env_nr
env_nt
tsa_nr
tsa_nt
taxdb
ref_prok_rep_genomes
ref_euk_rep_genomes
ref_viroids_rep_genomes
ref_viruses_rep_genomes
16S_ribosomal_RNA
18S_fungal_sequences
28S_fungal_sequences
LSU_eukaryote_rRNA
LSU_prokaryote_rRNA
SSU_eukaryote_rRNA
Betacoronavirus
ITS_RefSeq_Fungi
ITS_eukaryote_sequences
mito

Note that the database refseq_genomic is no longer available to download as a pre-formatted database.


NCBI BLAST Datasets can be loaded just like software modules. This allows users to replicate results by always being able to use a time stamped version of a database. The time stamp corresponds to the download date of the dataset. To search for available NCBI BLAST dataset modules run the command:

ml spider ncbiblastdb

You can then load a database by running the command module load. For example to load the databases timestamped on 08/02/2024:

module load ncbiblastdb/20240802

Loading this module sets the environment variable BLASTDB to /db/ncbiblast/20240802. This allows BLAST+ to search that directory for databases. You can then use the name of the database you would like to use.(nt, nr, etc.) Here is an example of a shell script to run BLAST+ on the batch queue using the database nt:

#!/bin/bash
#SBATCH --job-name=j_BLAST+
#SBATCH --partition=batch
#SBATCH --mail-type=ALL
#SBATCH --mail-user=username@uga.edu
#SBATCH --ntasks=1
#SBATCH --mem=10gb
#SBATCH --time=08:00:00
#SBATCH --output=BLAST+.%j.out
#SBATCH --error=BLAST+.%j.err

cd $SLURM_SUBMIT_DIR
module load BLAST+/2.14.1-gompi-2023a
module load ncbiblastdb/20240802
blastn -query example.fasta -out results.out -db nt

In your actual submission script, use your own discretion for the Slurm header values.

NCBI's nr and nt datasets are the most commonly used at the GACRC. GACRC staff will update these datasets every 3 months along with other potentially helpful NCBI BLAST datasets.


Blast database - old version 5

According to https://ncbiinsights.ncbi.nlm.nih.gov/tag/blastdbv5/, around Feb. 04, 2020 NCBI dropped the _v5 suffix for the v5 Blast databases.

An old version of the v5 database (download in June 2019) is located at a directory with name pattern as /db/ncbiblast.v5/dbname, where dbname is the name of the database, such as nr, nt, taxdb, refseq_rna, swissprot, etc.

These old versions of the nr and nt databases still have the _v5 suffix, i.e. they are called nr_v5 and nt_v5, respectively.

Back to Top