Question: how to include taxonomy information in blstx xml output ?
1
gravatar for Farbod
3.4 years ago by
Farbod3.3k
Toronto
Farbod3.3k wrote:

Dear Friends, Hi ( I'm not native in English so, be ready for some possible language flaws).

I want to use BLASTGrabber (software), but it needs taxonomy information in blast output results and it is only accepting xml output format (outfmt 5).

How can I include or add taxonomy information in my blast script ?

this is my usual script (using NCBI-blast-2.4.0+ ):

blastx -query Trinity.fasta -db nr -out blastx-Trinity.xml -evalue 1e-6 -num_threads 24 -max_target_seqs 1 -outfmt 5

ADD COMMENTlink modified 3.4 years ago by piet1.7k • written 3.4 years ago by Farbod3.3k
2

Correction: You can't include additional information (like taxonomy) with -outfmt 5. You may need to use one of the other options.

Options 6, 7 and 10 can be additionally configured to produce a custom format specified by space delimited format specifiers.

ADD REPLYlink modified 3.4 years ago • written 3.4 years ago by genomax78k
1

Dear genomax2, Hi.

It was the reason that I have asked this question : the blastGrabber need XML output and Taxonomy information at the same time !

What must I do ?

ADD REPLYlink modified 3.4 years ago • written 3.4 years ago by Farbod3.3k

Blast grabber supports XML and text output files (according to the manual you linked above) so you should be able to use 6,7 or 10 (with -outfmt "6 staxids sscinames" etc).

ADD REPLYlink modified 3.4 years ago • written 3.4 years ago by genomax78k
1

Hi,

I have tried it again. It did not accept blastx.outfmt and blastx.txt files but it accepted blastx.xml !

I will send an email to creator of the software in this regard.

Thank you genomax2

ADD REPLYlink modified 3.4 years ago • written 3.4 years ago by Farbod3.3k
1
gravatar for piet
3.4 years ago by
piet1.7k
planet earth
piet1.7k wrote:

Page 8 of the manual you have linked above clearly describes how taxonomy mapping works. Since you are running remote blast against the NCBI "nr" database, you can use the first method described on page 8, which means mapping gi to taxid. You have to download a mapping file from NCBI FTP server. This file is presumably very huge.

ADD COMMENTlink written 3.4 years ago by piet1.7k
1

Dear Piet, Hi and Thank you. I will try that.

ADD REPLYlink written 3.4 years ago by Farbod3.3k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1906 users visited in the last hour