Tools‎ > ‎BLAST‎ > ‎

Generate_database

Create a blastn species database

For faster mapping of many sequences against the same reference database

Example: create E. coli genome DB

# merge all E. coli reference genomes into one ecoli.fasta file
cat genomes/ecoli_strains/*.fna > ecoli.fasta

# create blast database
makeblastdb -in ecoli.fasta -parse_seqids -dbtype nucl

  -parse_seqids is required to identify gene/contig IDs in mapping hits
  -dbtype   specifies the type of sequences: protein 'prot' or nucleotide 'nucl'


# Blast using species DB  (option:  -db)
blastn -query genes.fasta -db ecoli.fasta  -outfmt 6 -evalue 1e-30

http://www.ncbi.nlm.nih.gov/books/NBK279688/



Subpages (1): Public sequences