This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Strange BLAST output

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
• 2,362 views • link
written by Prasad Prasad

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

• 0 views • link
written by 5heikki 5heikki

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

• 0 views • link
written by Prasad Prasad

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

• 0 views • link
written by 5heikki 5heikki

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?

• 0 views • link
written by Prasad Prasad

0 answers

No answers yet.

Log in to answer this question.

More posts like this