Hello all,
I am using ncbi-blast-2.2.26+ version on Linux machine (64 bit) with Fedora 16 OS. I have tried tblastx with my local nucleotide database
against ensembl fish cDNA sequence. I have followed these steps
1.makeblastdb -in example.fa -dbtype nucl -out example_nucl -logfile example_nuclblast.log
2.tblastx -db example_nucl -query fish.fa -out tblastx_result.tsv -max_target_seqs 20 -evalue 1e-10 -outfmt "7 qseqid sseqid pident length mismatch gaps qstart qend sstart send sqlen slen evalue bitscore" -num_threads 32
I get error
Error: NCBI C++ Exception:
"/home/coremake/release_build/build/PrepareRelease_Linux64-Centos_JSID_01_1540_130.14.18.6_9051_1329940081/rpmbuild/BUILD/ncbi-blast-2.2.26+/c++/src/objtools/alnmgr/alnmap.cpp", line 733: Error: CAlnVec::x_GetSeqLeftSeg(): Invalid Dense-seg: Row 0 contains gaps only.
Can you tell me why tblastx creating error. It outputs file but stops in between.
I tried looking source file at line no 00733
Any help is appreciated
I am using ncbi-blast-2.2.26+ version on Linux machine (64 bit) with Fedora 16 OS. I have tried tblastx with my local nucleotide database
against ensembl fish cDNA sequence. I have followed these steps
1.makeblastdb -in example.fa -dbtype nucl -out example_nucl -logfile example_nuclblast.log
2.tblastx -db example_nucl -query fish.fa -out tblastx_result.tsv -max_target_seqs 20 -evalue 1e-10 -outfmt "7 qseqid sseqid pident length mismatch gaps qstart qend sstart send sqlen slen evalue bitscore" -num_threads 32
I get error
Error: NCBI C++ Exception:
"/home/coremake/release_build/build/PrepareRelease_Linux64-Centos_JSID_01_1540_130.14.18.6_9051_1329940081/rpmbuild/BUILD/ncbi-blast-2.2.26+/c++/src/objtools/alnmgr/alnmap.cpp", line 733: Error: CAlnVec::x_GetSeqLeftSeg(): Invalid Dense-seg: Row 0 contains gaps only.
Can you tell me why tblastx creating error. It outputs file but stops in between.
I tried looking source file at line no 00733
Any help is appreciated
Comment