Hello,
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.
Use
blastdbcmd
: https://www.ncbi.nlm.nih.gov/books/NBK569853/thank you for the useful link