Entering edit mode
7.3 years ago
ylaguka
•
0
Hello, I am trying to use blast in R. It requires downloading the database to a local file. I'm having trouble finding the database for human genome grch37 (I will settle for grch38 too). I've tried refseq_genomic.242 and human_genomic_transcript, none of it returned the right result.
I've read the Manual for the database but couldn't find more detailed explanations on the exact database I'm looking for. The source where I was looking is ftp://ftp.ncbi.nlm.nih.gov/blast/db
Been stuck on it for a while, any help is appreciated, thanks!
You could get the hg19 genome sequence from UCSC here. Uncompress and then
cat
individual chromosome files into a single large file. Create blast indexes usingmakeblastdb
and then use those indexes for your program.Thank you for the suggestion. I was thinking about doing it that way. As it seems difficult to locate the the exact grch37 database the blast website is using, I will use the makeblastdb. Thanks again.