Hi everyone, I want to extract max score blast sequence among 50 genomic data set : my approach like this:
I did all the denovo assembly of all 50 genome.
Next i predicted all the ORF using Prodigal: (A Microbial Gene Prediction Software)
- Now i want to set up a local blast with a query gene sequence vs all 50 genomic ORF.
- 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.