Hi, I am trying to blast to the shared genbank database but am getting the following error. Perhaps the database was not indexed properly?
Error:
Inactive Modules:
- py-cython
Due to MODULEPATH changes, the following have been reloaded:
- bzip2/1.0.8
The following have been reloaded with a version change:
- gcc/11.3.0 => gcc/9.2.0 4) perl/5.36.0 => perl/5.30.0
- htslib/1.17 => htslib/1.10.2 5) python/3.11.3 => python/3.7.6
- openblas/0.3.21 => openblas/0.3.7 6) zlib/1.2.12 => zlib/18.0.4
BLAST Database error: No alias or index file found for nucleotide database [/project/biodb/genbank/2023-03-01/FASTA/nt] in search path [/scratch1/mruggeri/Aagar_denovo/Ref/subset10_cdh_alltags::]
Job script:
#!/bin/bash
#SBATCH --ntasks=1
#SBATCH --cpus-per-task=16
#SBATCH --mem-per-cpu=2G
#SBATCH --time=06:00:00
#SBATCH --error=%x_%j.err
#SBATCH --output=%x_%j.out
#SBATCH --mail-type=END
#SBATCH --mail-user=mruggeri@usc.edu
module load gcc/9.2.0
module load blast-plus/2.11.0
blastn -query subset10_cdh_alltags.fasta -db /project/biodb/genbank/2023-03-01/FASTA/nt -evalue 0.0001 -num_threads 16 -num_descriptions 5 -num_alignments 5 -out subset10_cdh_alltags.br