Malisa malisas
malisas
/ reciprocal_best_hits_for_blastp.py
Created
April 10, 2016 07:56
A script for finding pairs of orthologous genes between two species, based on their blastp alignment results.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| #!/usr/bin/python | |
| #Author: Malisa Smith 08-12-2015 | |
| import optparse | |
| from operator import itemgetter | |
| """ | |
| This program takes as input the results of 2 blastp runs to identify pairs of orthologous genes | |
| between two species. File 1 is the result of performing blastp with your query organism against | |
| the subject database. File 2 is the result of performing blastp with your subject organism against | |
| query database. |