Skip to content

Instantly share code, notes, and snippets.

@pcantalupo
Created November 5, 2014 18:55
Show Gist options
  • Save pcantalupo/f31c6ef8482c2663253d to your computer and use it in GitHub Desktop.
Save pcantalupo/f31c6ef8482c2663253d to your computer and use it in GitHub Desktop.
Get NCBI Viral RefSeq and make BLAST databases
#!/bin/bash
# get annotations
wget ftp://ftp.ncbi.nlm.nih.gov/refseq/release/viral/viral.1.genomic.gbff.gz
gunzip viral.1.genomic.gbff.gz
wget ftp://ftp.ncbi.nlm.nih.gov/refseq/release/viral/viral.1.protein.gpff.gz
gunzip viral.1.protein.gpff.gz
# download and create BLAST databases for viral refseq genomic
wget ftp://ftp.ncbi.nlm.nih.gov/refseq/release/viral/viral.1.1.genomic.fna.gz
gunzip viral.1.1.genomic.fna.gz
makeblastdb -dbtype nucl \
-parse_seqids \
-in viral.1.1.genomic.fna \
-out viral.1.1.genomic \
-title 'Virus Refseq Genomic'
# download and create BLAST databases for viral refseq protein
wget ftp://ftp.ncbi.nlm.nih.gov/refseq/release/viral/viral.1.protein.faa.gz
gunzip viral.1.protein.faa.gz
makeblastdb -dbtype prot \
-parse_seqids \
-in viral.1.protein.faa \
-out viral.1.protein \
-title 'Virus Refseq Protein'
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment