Instead of using an ncbi database, I ran the makeblastdb command on a fasta file I downloaded from ncbi.

The following script will produce...