Hello,
I have around 2000 DNA nucleotide sequences (60 bases long) stored in each row in an excel sheet. I want to run BLAST over each one of them individually and extract the "Description" of the first hit.
Like for Example: Suppose on NCBI BLAST website I insert one query sequence and hit BLAST with Homo sapiens as organisms. I extract the description of the first hit and paste it aside the sequence in the excel sheet (same row different column). Rather than doing this one by one, is there any package in R which helps me capture the description of individual sequences of the 100% query cover hit. By description I mean gene name or gene symbol
Thank you.
Does it need to be a R package ? Otherwise, Blast+ (blast from the command line) is the tool you need to run thousands of blast jobs locally and streamline the output.
I guess you're probably looking for
metablastr
.