makeblastdb check for correct formatting
1
0
Entering edit mode
5 months ago

hello, I made a blast database with

makeblastdb -in GRCh38.fa -dbtype nucl -out GRCh38.fa

the execution is fine but is there a way to check if the formatting is correct? Thank you

software error blastdb sanity check • 244 views
ADD COMMENT
0
Entering edit mode

out of curiosity : why do you want to check that? do you suspect something not being OK?

ADD REPLY
0
Entering edit mode

Yes, I am using a program that calls blast internally and I am getting errors such as [main_samview] random alignment retrieval only works for indexed BAM or CRAM files. and FileNotFoundError: [Errno 2] No such file or directory: '/home/gigiux/Documents/Transposons/Results/A0N-vs-Human endogenous retrovirus.bam.9'. I would like to exclude basic mistakes before debugging further...

ADD REPLY
0
Entering edit mode

the samview error/warning on the retrieval you can resolve by indexing your bam files (this has nothing to do with blast itself though)

ADD REPLY
1
Entering edit mode
5 months ago
GenoMax 102k

From blast+ package:

   $  blastdbcheck -h
    USAGE
      blastdbcheck [-h] [-help] [-db DbName] [-dbtype molecule_type]
        [-dir DirName] [-recursive] [-verbosity DefaultKey] [-full]
        [-stride StrideLength] [-random NumSequences] [-ends NumSequences]
        [-no_isam] [-legacy] [-must_have_taxids] [-cdd_delta] [-version]

    DESCRIPTION
       BLAST database integrity and validity checking application
ADD COMMENT
0
Entering edit mode

I got this:

$ blastdbcheck -db ~/Documents/Transposons/GRCh38.fa -dbtype nucl -full
Writing messages to <stdout> at verbosity (Summary)
ISAM testing is ENABLED.
Legacy testing is DISABLED.
TaxID testing is DISABLED.
Using `full' mode: every OID will be tested.

Testing 1 volume(s).
  /home/gigiux/Documents/Transposons/GRCh38.fa / MetaData:   [ERROR] caught exception.
 Result=FAILURE. 1 errors reported in 1 volume(s).
Testing 0 alias(es).
 Result=SUCCESS. No errors reported for 0 alias(es).

Total errors: 1
ADD REPLY
0
Entering edit mode

what is the exact error you're getting from running this?

side note: it's not because this check fails there is for sure something wrong with your blast DBs. (it could for instance complain on missing some of the taxid files, which is not a vital thing for running blasts)

ADD REPLY
1
Entering edit mode

MetaData: [ERROR] caught exception.

See: Blastdbcheck Error and http://seqanswers.com/forums/archive/index.php/t-34826.html

It is related to absence of taxonomy DB, which may be fine in this case.

ADD REPLY

Login before adding your answer.

Traffic: 1616 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6