SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
What tools can convert sequence file from tabular format to fasta format? yangjianhunt Bioinformatics 5 03-26-2014 02:48 PM
Merging Blast results - xml or tabular format mstagliamonte Bioinformatics 3 08-01-2013 03:02 AM
Weird Tabular BLAST Occurrence tboothby Bioinformatics 5 03-06-2012 04:20 PM
Convert SOAPsplice v1.4 output to SAM cdb3 RNA Sequencing 0 08-25-2011 11:21 AM
convert blast output rururara Bioinformatics 1 04-08-2011 01:48 AM

Reply
 
Thread Tools
Old 08-27-2013, 06:52 AM   #1
Rsola
Junior Member
 
Location: Europe

Join Date: Aug 2013
Posts: 2
Default Convert Blast output to tabular

Hi,

Does anybody know how to convert Blast Output as below to a tabular (.csv) output or how to extract top hits? In the online tool that I use I can't select for tabular output.

Thanks!

Query= OsSNP0002
(142 letters)

Database: ../uniprot-oryza-sativa-indica.fasta
42,240 sequences; 15,872,719 total letters

Searching..................................................done



Score E
Sequences producing significant alignments: (bits) Value

tr|A2WJN5|A2WJN5_ORYSI Putative uncharacterized protein OS=Oryza... 55 2e-18
tr|B8AKP7|B8AKP7_ORYSI Putative uncharacterized protein OS=Oryza... 53 6e-17
tr|B8AEL6|B8AEL6_ORYSI Putative uncharacterized protein OS=Oryza... 28 1.5
tr|B8B114|B8B114_ORYSI Putative uncharacterized protein OS=Oryza... 28 2.0
tr|B8AQS5|B8AQS5_ORYSI Putative uncharacterized protein OS=Oryza... 27 2.6
tr|B8ADP7|B8ADP7_ORYSI Putative uncharacterized protein OS=Oryza... 27 4.4
tr|A2ZKU4|A2ZKU4_ORYSI Putative uncharacterized protein OS=Oryza... 26 7.5
tr|A2Y774|A2Y774_ORYSI Putative uncharacterized protein OS=Oryza... 26 7.5
tr|A2WXP8|A2WXP8_ORYSI Putative uncharacterized protein OS=Oryza... 25 9.7 [/SIZE][/SIZE]
Rsola is offline   Reply With Quote
Old 08-27-2013, 06:59 AM   #2
lran2008
Member
 
Location: quebec

Join Date: Apr 2013
Posts: 35
Default

You could choose the Blast output format by setting -outfmt. Try -outfmt 6 .
lran2008 is offline   Reply With Quote
Old 08-27-2013, 07:05 AM   #3
Rsola
Junior Member
 
Location: Europe

Join Date: Aug 2013
Posts: 2
Default

Thank you for your reply, but like I indicated, I can't set the output format in the online tool that I used (no command lines in the tool). This is the output I have to deal with.
Rsola is offline   Reply With Quote
Old 08-27-2013, 07:22 AM   #4
lran2008
Member
 
Location: quebec

Join Date: Apr 2013
Posts: 35
Default

Quote:
Originally Posted by Rsola View Post
Thank you for your reply, but like I indicated, I can't set the output format in the online tool that I used (no command lines in the tool). This is the output I have to deal with.
Sorry, I didn't notice it. I don't know what's the output format by online blast. But if it can be opened by excel or be imported in to R, you could get what you want then.
lran2008 is offline   Reply With Quote
Old 08-27-2013, 07:29 AM   #5
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,087
Default

Are you able to use the command line/have access to python? Following may work in that case

https://bitbucket.org/ctb/ngs-course....py?at=default

Extended in this post:

http://haotu.wordpress.com/2012/03/0...-to-csv-excel/
GenoMax is offline   Reply With Quote
Old 08-27-2013, 07:33 AM   #6
rhinoceros
Senior Member
 
Location: sub-surface moon base

Join Date: Apr 2013
Posts: 372
Default

I'm not really following how your output is formatted, but I'm sure clever use of e.g. sed could solve your problem. That or using blast_formatter (if your output is asn.1)..
rhinoceros is offline   Reply With Quote
Old 08-27-2013, 07:35 AM   #7
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,087
Default

If you used NCBI to do the search then you can download the results as "hit table.csv" file under the "download" option on the results page.
Attached Images
File Type: png ncbi_download.PNG (26.8 KB, 25 views)
GenoMax is offline   Reply With Quote
Old 08-01-2014, 09:57 AM   #8
SDPA_Pet
Senior Member
 
Location: US

Join Date: Apr 2013
Posts: 222
Default

Hi GenoMax,

Actually, it is more complicated. The blast output was extracted from MEGAN. It was originally blasted using -m 0 format. I extracted some blast output of some pathways in order to look at the details. I would like to report the highest hit score species in my paper. So, it would be easy for me to convert it to tabular format.
SDPA_Pet is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 05:14 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO