I have created a Blast database using a reference genome. Then, I have performed a local blast search in command line using a gene of interest. I have obtained some hits with the usual Blasting information. Now, I want to extract the exact matching sequence from the blast search with their corresponding start and stop position. I have a
.gff file for the genes as well for the reference genome.
Can you give me any suggestions on how to do that?
Thanks a lot.