Strange BLAST output
0
0
Entering edit mode
7.5 years ago
Prasad ★ 1.6k

hi,

I have used standalone blast to search few of transcript sequences against NCBI NR database. For few transcripts (~88) i have found strange output as

c93530_g1_i1    gi|281206185    0.00    0       0       0       0       0       0       0       4e-19   94.4
c93530_g1_i1    gi|281206185    0.00    0       0       0       0       0       0       0       1e-17   90.1
c93530_g1_i1    gi|281206185    0.00    0       0       0       0       0       0       0       4e-17   88.6

how come everything else is 0 except evalue and score. command used,

blastx -query transcript.fasta -db /nr_db_oct_16/nr -evalue 0.0000000001 -outfmt 6 -out blast_out -num_threads 20 -max_target_seqs 1

I did online blast for the same transcript against NR getting proper result. Any idea what might be the problem?

Thanks

blast similarity sequence • 1.8k views
ADD COMMENT
0
Entering edit mode

Which version of blast are you using? Do you get the same results for this query if you use another output format?

ADD REPLY
0
Entering edit mode

ncbi-blast-2.2.30. I haven't tried in other format. let me try that

ADD REPLY
0
Entering edit mode

It was released over 2 years ago. You should see if you can reproduce the problem with the latest release.

ADD REPLY
0
Entering edit mode

i did the blast again with these transcripts, now getting proper results. now just wondering what could be the reason. was this because of a glitch in server?

ADD REPLY

Login before adding your answer.

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