[Bioperl-l] blast question

Smithies, Russell Russell.Smithies at agresearch.co.nz
Sun Nov 4 21:00:37 UTC 2012


What version of blast are you using?
There have been quite a few bug fixes and I suspect any responses from NCBI will suggest upgrading to the current version of blast+


--Russell

-----Original Message-----
From: bioperl-l-bounces at lists.open-bio.org [mailto:bioperl-l-bounces at lists.open-bio.org] On Behalf Of shalabh sharma
Sent: Saturday, 3 November 2012 3:50 a.m.
To: Fields, Christopher J
Cc: bioperl-l; blast-help at ncbi.nlm.nih.gov
Subject: Re: [Bioperl-l] blast question

I know, i am really worried about my past analysis now.
Thanks a lot for cc'ing this mail Chris.

-Shalabh

On Thu, Nov 1, 2012 at 5:41 PM, Fields, Christopher J <cjfields at illinois.edu
> wrote:

> That's a scary error, but the best place to submit this would be the 
> BLAST help list at NCBI (cc'd)
>
> chris
>
> On Nov 1, 2012, at 2:36 PM, shalabh sharma <shalabh.sharma7 at gmail.com>
> wrote:
>
> > Hi All,
> >          First of all i am really very sorry for posting blast 
> > question
> in
> > this forum, I am not sure if this is the right place.
> > I will really appreciate if anyone can guide me to the right direction.
> >
> > I am using blastall to get a top hit from a database so i am using 
> > -v 1
> -b
> > 1 (i hope this is right).
> > But the strange part is that i am getting wrong results.
> >
> > for example: if i use -v 1 -b 1 then for one of the hit i am getting
> this:
> >
> >
> > Sequences producing significant alignments:                      (bits)
> > Value
> >
> > fig|6666666.11092.peg.487 Thiol:disulfide oxidoreductase TlpA          38
> > 4e-04
> >
> >
> > If i use -v 3 -b 3 then i am getting this for the same query:
> >
> > Sequences producing significant alignments:                      (bits)
> > Value
> >
> > fig|6666666.11092.peg.1134 COG3118: Thioredoxin domain-containin...   570
> > e-167
> > fig|6666666.11092.peg.487 Thiol:disulfide oxidoreductase TlpA          38
> > 9e-07
> > fig|6666666.11092.peg.1133 Exodeoxyribonuclease III (EC 3.1.11.2...    18
> > 1.0
> >
> > As you can see the top hit in the first case is totally wrong.
> >
> > I would really appreciate if someone can help me out, or direct to 
> > in the right direction.
> >
> > Thanks
> > Shalabh
> >
> >
> >
> > --
> > Shalabh Sharma
> > Scientific Computing Professional Associate (Bioinformatics 
> > Specialist) Department of Marine Sciences University of Georgia 
> > Athens, GA 30602-3636 
> > _______________________________________________
> > Bioperl-l mailing list
> > Bioperl-l at lists.open-bio.org
> > http://lists.open-bio.org/mailman/listinfo/bioperl-l
>
>


--
Shalabh Sharma
Scientific Computing Professional Associate (Bioinformatics Specialist) Department of Marine Sciences University of Georgia Athens, GA 30602-3636 _______________________________________________
Bioperl-l mailing list
Bioperl-l at lists.open-bio.org
http://lists.open-bio.org/mailman/listinfo/bioperl-l

=======================================================================
Attention: The information contained in this message and/or attachments
from AgResearch Limited is intended only for the persons or entities
to which it is addressed and may contain confidential and/or privileged
material. Any review, retransmission, dissemination or other use of, or
taking of any action in reliance upon, this information by persons or
entities other than the intended recipients is prohibited by AgResearch
Limited. If you have received this message in error, please notify the
sender immediately.
=======================================================================




More information about the Bioperl-l mailing list