[Bioperl-l] Taxa Id from blast report

Jason Stajich jason.stajich at gmail.com
Fri Apr 19 19:38:00 UTC 2013


Did you provide -parse_seqids in the header?  

Peter dealt with related things here:
http://blastedbio.blogspot.com/2012/10/my-ids-not-good-enough-for-ncbi-blast.html

Jason
On Apr 19, 2013, at 1:05 PM, shalu sharma <sharmashalu.bio at gmail.com> wrote:

> Hi,
>    Thanks everyone for you inputs.
> @Peter:
> I got really excited when i saw that you can even get super kingdom, but
> when i tried to test it i just got taxa ids but not the super kingdom. Do
> you have any idea whats going wrong?
> my command:
> blastx  -query test.fas -db /db/ncbiblast/refseq/latest/refseq_protein
> -max_target_seqs 1 -outfmt "6 staxids sskingdoms"
> 
> output:
> 246200    N/A
> 246200    N/A
> 
> Thanks
> Shalu
> 
> 
> On Thu, Apr 18, 2013 at 3:52 PM, Smithies, Russell <
> Russell.Smithies at agresearch.co.nz> wrote:
> 
>> I agree they have finally listened and added features requested by users
>> but I've been suggesting they have a compressed output format available
>> from eutils or genbank for years but have made no headway ;- (
>> What's so hard about gzip'ping the output? I'm sure it would go a long way
>> toward solving all the problems we get with truncated replies from queries!!
>> 
>> --Russell
>> 
>> -----Original Message-----
>> From: bioperl-l-bounces at lists.open-bio.org [mailto:
>> bioperl-l-bounces at lists.open-bio.org] On Behalf Of Fields, Christopher J
>> Sent: Friday, 19 April 2013 6:26 a.m.
>> To: Peter Cock
>> Cc: bioperl-l at lists.open-bio.org; shalu sharma; shalabh sharma
>> Subject: Re: [Bioperl-l] Taxa Id from blast report
>> 
>> On Apr 18, 2013, at 11:48 AM, Peter Cock <p.j.a.cock at googlemail.com>
>> wrote:
>> 
>>> On Thu, Apr 18, 2013 at 5:32 PM, shalabh sharma
>>> <shalabh.sharma7 at gmail.com> wrote:
>>>> Hey Peter,
>>>>      Thanks a lot, I really appreciate it. I wanted these things
>>>> implemented in blast from long time.
>>>> 
>>>> Thanks
>>>> Shalabh
>>> 
>>> Me too. You can get the descriptions from the plain text BLAST or XML
>>> output already of course, but they're not so nice to work with.
>>> 
>>> Peter
>> 
>> NCBI has been much more receptive of user input over the last several
>> years, much more so than in the past.  I understand the reasoning for
>> dropping BLAST support (though there were definitely needless bumps in that
>> process).
>> 
>> chris
>> 
>> 
>> _______________________________________________
>> 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.
>> =======================================================================
>> 
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l

Jason Stajich
jason.stajich at gmail.com
jason at bioperl.org





More information about the Bioperl-l mailing list