How To Extract Max Score Blast Hit Among A Huge Data Set?

Hi everyone, I want to extract max score blast sequence among 50 genomic data set : my approach like this:

1. I did all the denovo assembly of all 50 genome.

2. Next i predicted all the ORF using Prodigal: (A Microbial Gene Prediction Software)

3. Now i want to set up a local blast with a query gene sequence vs all 50 genomic ORF.
4. Finally extract the best blast score sequence among all 50 genome. Can anyone suggest me how to proceed step 3 and step 4

Thank you advance.