[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