Entering edit mode
6.9 years ago
lidanqing1990
▴
10
hi,everyone! I use blast+ to make annotation of RNA-seq data against swissprot When I run blast, it always stopped by error. I used blast 2.3.0 and my commands is:
blastx -query /home/aubdyl/purple_nutsedge/purple_nutsedge_trinity/Trinity.fasta -out /home/aubdyl/purple_nutsedge/purple_nutsedge_trinity/purple_nutsedge_swissprot -db /home/aubdyl/data/blast/swissprot.00 -outfmt 6 -evalue 1e-5 -num_descriptions 10 -num_threads 4
and the error is:
Critical: (110.6) CNcbiRegistry: Syntax error in system-wide configuration file: NCBI C++ Exception:
T0 "/home/coremake/release_build/build/PrepareRelease_Linux64-Centos_JSID_01_250088_130.14.22.10_9008__PrepareRelease_Linux64-Centos_1448906370/c++/compilers/unix/../../src/corelib/ncbireg.cpp", line 675: Error: ncbi::IRWRegistry::x_Read() - Invalid registry entry format: '$ cat .ncbirc ' (m_Pos = 1)
Warning: [blastx] The parameter -num_descriptions is ignored for output formats > 4 . Use -max_target_seqs to control output
Error: NCBI C++ Exception:
T0 "/home/coremake/release_build/build/PrepareRelease_Linux64-Centos_JSID_01_250088_130.14.22.10_9008__PrepareRelease_Linux64-Centos_1448906370/c++/compilers/unix/../../src/corelib/ncbireg.cpp", line 675: Error: ncbi::IRWRegistry::x_Read() - Invalid registry entry format: '$ cat .ncbirc ' (m_Pos = 1)
I have changed another edition of blast, but it still doesn't work. Could anyone give me some help ? Thank you so much!
Did you edit the
.ncbirc
file manually and make any changes?Also try changing
-db /home/aubdyl/data/blast/swissprot.00
to-db /home/aubdyl/data/blast/swissprot
.Thanks, I didn't edit .ncbirc file manually and -db /home/aubdyl/data/blast/swissprot doesn't work. Could you tell me where could I find .ncbirc file ?
It should be in your home directory. See if you are able to
cat ~/.ncbirc
?Yes, I saw it. Should I delete this file? Blast+ Error seems similar to my question.
Move the file to a new name e.g.
mv ~/.ncbirc ~/.ncbirc.ORIG
and try.Many thanks, I deleted that file and commands are running well at present.