Parsing Blast Results from Your Favorite Database
From GcatWiki
This tutorial assumes the user has blast version __________ installed and has already made their local blast database (on their computer). It is also written for Macintosh users, however all scripts and tools mentioned here are windows compatible with Windows.
- Run your search using the command
/usr/local/ncbi/blast/bin/blastn -query scf1453.fasta -db bb_latest_assembly.fasta -outfmt "7 qacc sacc evalue qstart qend sstart send"
- What does this command do? "-" indicates a command. The text that follows is the actual command.
- -query = query file
- -db = database to search
- -outfmt = output format DO NOT TOUCH THIS UNLESS YOU KNOW HOW TO EDIT blastParse.py
- What does this command do? "-" indicates a command. The text that follows is the actual command.
- 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.