Question: Error: NCBI C++ Exception: Error pre-fetching sequence data
0
gravatar for 1291016966
7 months ago by
12910169660
12910169660 wrote:

Hi all, Recently, I run blastn on ubuntu. The command just as below:

nohup /analysis/Software/ncbi-blast-2.6.0+/bin/blastn -query data/S08-C-1-K.1.fasta -db /database/blast_db/h_abfpv_genome -out blastn_m6.out -outfmt 6 -evalue 1e-5 -num_threads 10 -max_target_seqs 1 &

But there is an error.

Error: NCBI C++ Exception:
    T0 "/home/coremake/release_build/build/PrepareRelease_Linux64-Centos_JSID_01_350334_130.14.22.10_9008__PrepareRelease_Linux64-Centos_1481139955/c++/compilers/unix/../../src/corelib/ncbiobj.cpp", line 977: Critical: ncbi::CObject::ThrowNullPointerException() - Attempt to access NULL pointer.
    T0 "/home/coremake/release_build/build/PrepareRelease_Linux64-Centos_JSID_01_350334_130.14.22.10_9008__PrepareRelease_Linux64-Centos_1481139955/c++/compilers/unix/../../src/app/blast/blast_app_util.cpp", line 766: Critical: BLAST::ncbi::BlastFormatter_PreFetchSequenceData() - Error pre-fetching sequence data 
     Stack trace:
      /analysis/Software/ncbi-blast-2.6.0+/bin/blastn ???:0 ncbi::CStackTrace::CStackTrace(std::string const&) offset=0x8B addr=0x1a9587b
      /analysis/Software/ncbi-blast-2.6.0+/bin/blastn ???:0 ncbi::CException::SetSeverity(ncbi::EDiagSev) offset=0xEB addr=0x1bc0d8b
      /analysis/Software/ncbi-blast-2.6.0+/bin/blastn ???:0 ncbi::CObject::ThrowNullPointerException() offset=0x2EA addr=0x1c2db5a
      /analysis/Software/ncbi-blast-2.6.0+/bin/blastn ???:0 ncbi::objects::CBioseq_Handle::IsSetInst_Fuzz() const offset=0x0 addr=0x143e9a0
      /analysis/Software/ncbi-blast-2.6.0+/bin/blastn ???:0 ncbi::blast::LoadSequencesToScope(std::vector<ncbi::objects::CSeq_id_Handle, std::allocator<ncbi::objects::CSeq_id_Handle> >&, std::vector<ncbi::CRange<unsigned int>, std::allocator<ncbi::CRange<unsigned int> > >&, ncbi::CRef<ncbi::objects::CScope, ncbi::CObjectCounterLocker>&) offset=0x3BF addr=0xcc93ff
      /analysis/Software/ncbi-blast-2.6.0+/bin/blastn ???:0 ncbi::BlastFormatter_PreFetchSequenceData(ncbi::blast::CSearchResultSet const&, ncbi::CRef<ncbi::objects::CScope, ncbi::CObjectCounterLocker>, ncbi::blast::CFormattingArgs::EOutputFormat) offset=0x82 addr=0x990d52
      /analysis/Software/ncbi-blast-2.6.0+/bin/blastn ???:0 CBlastnApp::Run() offset=0xF51 addr=0x987141
      /analysis/Software/ncbi-blast-2.6.0+/bin/blastn ???:0 ncbi::CNcbiApplication::x_TryMain(ncbi::EAppDiagStream, char const*, int*, bool*) offset=0x92A addr=0x1a9ddfa
      /analysis/Software/ncbi-blast-2.6.0+/bin/blastn ???:0 ncbi::CNcbiApplication::AppMain(int, char const* const*, char const* const*, ncbi::EAppDiagStream, char const*, std::string const&) offset=0x773 addr=0x1a9c653
      /analysis/Software/ncbi-blast-2.6.0+/bin/blastn ???:0 main offset=0x2CB addr=0x985f3b
      /lib/x86_64-linux-gnu/libc.so.6 ???:0 __libc_start_main offset=0xF0 addr=0x7faa19d5e830
      /analysis/Software/ncbi-blast-2.6.0+/bin/blastn ???:0 /analysis/Software/ncbi-blast-2.6.0+/bin/blastn() [0x985b29] offset=0x0 addr=0x985b29

Does anyone know how to solve it?

Looking for the reply! Thank you very much!

Myshu

ADD COMMENTlink modified 5 months ago by nemo0 • written 7 months ago by 12910169660

I think the problem is with your query fasta. Did the error occur right away, or the BLAST run for a while before crashing? You can try to troubleshoot your input fasta by sppliting it in smaller chunks and testing then separately.

edit: see pointers at Question: ncbi error C++ exception.

ADD REPLYlink modified 7 months ago • written 7 months ago by h.mon20k

The error occur right away. The result is null.

Yes, I tried to split the input fasta into small chunks...But even I split into 3000 chunks, and each of file only has about 600 sequences, the error still occurred....

Do you know what the largest number of input sequence is?

ADD REPLYlink written 7 months ago by 12910169660

Select one (preferably short) sequence that causes the error and post it here.

ADD REPLYlink written 7 months ago by h.mon20k

Thanks for your reply. And I tried to select several sequences to do the test.

And I found that the sequence below causes the error...

>ST-E00291:181:HFVCTALXX:1:1101:8420:1872/1
GTTTGTATTTCTGTGGGATCAGTGGTGATATCCTCTTTAGCATTTTTTATTGCATCTAATTGATTCTTCTCTATTATTCTTGCTAGTGGTCTACCAATTTTATCTTTTAAAAAAACCAGCTCCTCGATTGATCTTTTTGAAGGGTT

But why? I can't find any different....

ADD REPLYlink modified 7 months ago • written 7 months ago by 12910169660

Hi 1291016966,

I tried with your sequence using the same command except for the database with blastn v2.6.0+. It worked well.

ADD REPLYlink written 5 months ago by Sishuo Wang140

Hi 1291016966,

Did you solve the above problem? I got similar error when I run tblastn. In my case, the sequence below causes the error.

gi|1072238524|gb|OEU23251.1| MKQSVAAIVGTNNNADNLAVSASHLYSYESAVNALLSPGLHQSITKEDILKSSLRRTKTV ADMRYYWNKILKYNRHTNNVNTNDDDQKKNKKKPLLIHITGTKGKGSTACMCENILRSNG YKTGLFTSPHLINIRERIRYNGQPINEKLFSDVYWKIRQAFEVDESDETSDEEVDPPPKL PGYFRMLTLMGMYTFLYELSDDVDVIILEVGMGGRYDATNFLDTTNTTTTTTTAATTNSN NYCSYFQRVVYGVTLLDLDHTRILGTTLKQIAWEKGGIFSMNKLNPNGISSSTSSTTSTS KLLSEPKSNSKSTTTEKEETDDKDEPSSSTSSTSSSSSNKKYYILDRHDFRMVPPSVTLN ALSSAKWPGRCQTLINKKKNVRFYLDGAHTPQSLDATVEWFRSKSNSSSSSNGTGIFNCS HERNPVELLELLLKHQSGFSRVYFAQSDTSRPSPIAKASAESLLKERGITIREELLADEG KKTDDDSEGEIEGKMNTKTTTWQETLAIIWKHLLIQSTADTKSASTSASALVSLSQTRDQ DSKVFVTGSLYLVGSFLTALGWTEESSPTSPSPSSSSSSTLASTTIQ--

ADD REPLYlink written 5 months ago by nemo0

You might try to run the blast with a diff number of threads parameter (eg. 5 ) we've noticed that sometimes certain values of the num_threads are not working well, no idea why though.

Otherwise, perhaps activating/disabling low-complexity filter might help?

ADD REPLYlink written 5 months ago by lieven.sterck2.6k

Thanks for your suggestion. I tried to run the blast with 5 threads and 10 threads. but I found a same error. I will try your other suggestion.

ADD REPLYlink written 5 months ago by nemo0
1

Hi nemo,

I got the following message: CFastaReader: Hyphens are invalid and will be ignored around line 2

Not sure but could it be because of the presence of the hyphens?

ADD REPLYlink modified 5 months ago • written 5 months ago by Sishuo Wang140
1

Excellent observation! Did not notice those

@ nemo : where do the hyphens come from? They should not be present in a fasta file.

ADD REPLYlink modified 5 months ago • written 5 months ago by lieven.sterck2.6k

Hi Thank you for your reply. hyphens come from our in-house developed program. yes, I noticed the messeage " CFastaReader: Hyphens are invalid and will be ignored around line 2". but this error is not the problem. I still have the same error (Error: NCBI C++ Exception:...) as 1291016966 mentioned above.

ADD REPLYlink written 5 months ago by nemo0
1

OK, so you did try when you removed the hyphens, right?

Did the formatting of your blastDB went ok? is there anything particular in the makeblastdb log file?

Are you able to blast anything with your 2.6.0 version?

ADD REPLYlink written 5 months ago by lieven.sterck2.6k

Hi lieven,

yes, it seems there was no problem with makeblatdb. thank you for your suggestion. When I used old ver.(2.2.29), blast search was completed without any error!

Thank you,

ADD REPLYlink written 4 months ago by nemo0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1751 users visited in the last hour