Difference between revisions of "Parsing Blast Results from Your Favorite Database"
From GcatWiki
Line 1: | Line 1: | ||
− | This tutorial assumes the user has blast version 2.2.24 installed and has already made their local blast database (on their computer). It is also written for Macintosh users; however, all scripts and tools are Windows compatible or have similar programs for Windows. | + | This tutorial demonstrates how to isolate data from large blast hits. |
+ | |||
+ | NOTE: This tutorial assumes the user has blast version 2.2.24 installed and has already made their local blast database (on their computer). It is also written for Macintosh users; however, all scripts and tools are Windows compatible or have similar programs for Windows. | ||
+ | |||
#Run your search using the command<pre>/usr/local/ncbi/blast/bin/blastn -query querySequence.fasta -db dataBase.fasta -outfmt "7 qacc sacc evalue qstart qend sstart send" -out blast_output.txt</pre> | #Run your search using the command<pre>/usr/local/ncbi/blast/bin/blastn -query querySequence.fasta -db dataBase.fasta -outfmt "7 qacc sacc evalue qstart qend sstart send" -out blast_output.txt</pre> |
Revision as of 20:34, 23 February 2011
This tutorial demonstrates how to isolate data from large blast hits.
NOTE: This tutorial assumes the user has blast version 2.2.24 installed and has already made their local blast database (on their computer). It is also written for Macintosh users; however, all scripts and tools are Windows compatible or have similar programs for Windows.
- Run your search using the command
/usr/local/ncbi/blast/bin/blastn -query querySequence.fasta -db dataBase.fasta -outfmt "7 qacc sacc evalue qstart qend sstart send" -out blast_output.txt
- What does this command do? "-" indicates a command. The text that follows is the actual command parameters.
- -query = query file
- -db = database to search
- -outfmt = output format DO NOT CHANGE THIS UNLESS YOU KNOW HOW TO EDIT blastParse.py
- -out = output file for blast results
- What does this command do? "-" indicates a command. The text that follows is the actual command parameters.
- Download blastParse.py
- Place blastParse.py into the same folder as your blast results file
- Open terminal and navigate into the folder containing blastParse.py using following Unix commands
cd OR ls
- Run blastParse.py using the command
python blastParse.py
- Follow prompts of the BLASTPARSE program
- Results will be saved as tab delimited data in text files. If you would like to visualize the data, open files in excel to make graphs (right click > open with > excel).