BLAST Reciprocal Best Hits (Linux)
0
0
Entering edit mode
8.0 years ago

I am an undergraduate and relatively new to bioinformatics and am currently trying to perform a reciprocal best hits analysis in linux. I blasted my transcript sequences against both of my genomes of interest (in my case, Megachile rotundata and Apis mellifera). I have the outputs of those blasts in csv tables. I now want to compare the two outputs to determine corresponding genes between the two species but am unsure of how to do this. Everything I have been reading has been written summary explaining this step but I haven't been able to see any example scripts. I was wondering if anyone knew the necessary script for this step or if there is an active example/tutorial that someone could point me towards? Thanks!

blast • 4.2k views
ADD COMMENT
0
Entering edit mode

An example script: rbbh.py. Another: Galaxy RBH.

ADD REPLY

Login before adding your answer.

Traffic: 3485 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6