Running BLAST optimally

From IBERS Bioinformatics and HPC Wiki
Revision as of 13:48, 7 August 2013 by Mjv08 (talk | contribs) (Created page with "This is a brief description of a couple of experiments that I (mjv08) have tried out using blast in order to determine the best way in which to run it. Okay, so you have a fa...")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

This is a brief description of a couple of experiments that I (mjv08) have tried out using blast in order to determine the best way in which to run it.

Okay, so you have a fasta file you wish to blast. The file I had was supplied by Russ (rom) and contained 18788 sequences in a fasta file. He wanted to run the default settings for blast using blast2go. These settings using blast+ on the command line are; 

blastx -db nr -query input.fasta -out output.xml -evalue 1e-3 -outfmt 5 -show_gis -num_threads 1

(I'm not sure the number of threads the blast2go use when accessing the NCBI www-blast)

So,

Retrieved from "https://bioinformatics.ibers.aber.ac.uk/wiki/index.php?title=Running_BLAST_optimally&oldid=124"