Generate a blast genome from wormbase genomic sequences
# Download ref genomes/contigs
lftp ftp.wormbase.org -e 'mget /pub/wormbase/species/**/sequence/genomic/*.current.genomic.fa.gz'
function get_genomes() {
for i in `ls -1 *canonical*fa.gz`; do
IFS='.' read -r -a genome <<< "$i";
zcat ${i} | sed "s/>/> ${genome[0]} :: /g"
done;
}
get_genomes() | makeblastdb -in - -out wormbase_blast_genome -title wormbase_blast -dbtype=nucl