How To Get The Result Which Mismatch 1 To The Query Sequence By Nucleotide-Nucleotide Blast 2.2.27+
3
0
Entering edit mode
11.2 years ago
2011101101 ▴ 110

the program is blastn -task blastn-short -query rfam_trimed_short_uniq.fa -db /root/bio-analysis-soft/database/gss-est-db/tage -ungapped -penalty -1 -reward 1 -outfmt 6 -out 2 -num_threads 12

I get an error.Error: NCBI C++ Exception: "/root/bio-analysis-soft/msw/ncbi-blast-2.2.27+-src/c++/src/corelib/ncbiobj.cpp", line 689: Critical: ncbi::CObject::ThrowNullPointerException() - Attempt to access NULL pointer.

I think the error caused by -num_threads 12.can you help me ?

blast • 5.1k views
ADD COMMENT
2
Entering edit mode
10.3 years ago

I guest you compiled your BLAST+ from source tarball.

To use num_threads, you have to add a "--with-mt" argument when run configre to support multithreads.

You can follow this manual: BLAST Command Line Applications User Manual

cd c++
./configure --without-debug --with-strip --with-mt --with-build-root=ReleaseMT
cd ReleaseMT/build
make all_r
ADD COMMENT
0
Entering edit mode

Installing BLAST+ executables is the better choice. Compilation from source costs long time and errors may occur.

ADD REPLY
0
Entering edit mode

Not when they have got the "wrong" compile options. Also if that you claim is true then the Gentoo distribution is a bad idea is general (spoiler: it is not).

ADD REPLY
0
Entering edit mode
11.2 years ago

According to NCBI's BLAST Manual:

num_threads    integer    1    Number of threads (CPUs) to use in blast search.

The command-line option is mostly referred to using CPUs (as opposed to simple threads) in the manual. I suspect you don't have 12, but did you try it with the number of CPUs you actually have got?

ADD COMMENT
0
Entering edit mode

My ,cups 4 ,but 16 core.when num_threads 2,there is still an error

ADD REPLY
0
Entering edit mode
11.2 years ago
Sujai Kumar ▴ 270

I think whenever I've seen that error, the query fasta file has had problems: can you check for

  • blank sequences (i.e. no sequence after a line beginning with ">"); or
  • tabs in the fasta header
  • non-standard line endings (Win/DOS endings on Linux/OSX etc)

On a side note - I seem to remember seeing a tool that checks and cleans Fasta files... might be worth running once on your query file

ADD COMMENT
0
Entering edit mode

when num_threads 1,there is no error

ADD REPLY

Login before adding your answer.

Traffic: 1337 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