BackBLAST.py

Description:

A Biopython program that takes a list of query proteins and uses local BLASTp to search for highly similer proteins within a local blast database (usally a local db of a target proteome). The program then BLASTps backwards from the found subject proteins to the query proteome to confirm gene orthology.

Requirements:

• This program requires the Biopython module: http://biopython.org/wiki/Download
• This script requires BLAST+ 2.2.9 or later.
• All operations are done with protien sequences.
• All query proteins must be from sequenced genomes in order to facilitate backwards BLAST.
• MakeBlastDB must be used to create BLASTp databases for both query and subject proteomes.
• BLAST databases require the FASTA file they were made from be in the same directory.

Usage: BackBLAST.py <queryGeneList.faa> <queryBLASTDB.faa> <subjectBLASTDB.faa>

Example: BackBLAST.py queryGeneList.faa AL123456.3.faa AUUJ00000000.faa