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!
An example script: rbbh.py. Another: Galaxy RBH.