BlastN2.2.1Search DBs for Nucleotide Sequence similarityAltschul, Madden, Schaeffer, Zhang, Miller, LipmanAltschul, Stephen F., Thomas L. Madden, Alejandro A. Schaeffer,Jinghui Zhang, Zheng Zhang, Webb Miller, and David J. Lipman (1997), Gapped BLAST and PSI-BLAST: a new generation of protein database search programs, Nucleic Acids Res. 25:3389-3402.Nucleic Acid Sequenceblastnblast_initperl " " -10scheduler_conf3scheduler.confperl "node_exclusive=1\\n" .
"nodes=1\\n" .
"threads_per_process=16\\n" querySequence Fileperl"blast_wrapper -query query.fasta" query.fasta3blast2perl"-out blast_output.txt"9blastn compares a nucleotide query sequence against a nucleotide sequence databasenum_threadsperl"-num_threads 16"6runtime1scheduler.confMaximum Hours to Run (up to 168 hours)0.5The maximum hours to run must be less than 168perl$runtime > 168.0The maximum hours to run must be greater than 0.05perl$runtime < 0.05perl"runhours=$value\\n"The job will run on 16 processors as configured. If it runs for the entire configured time, it will consume 16 x $runtime cpu hoursperldefined $nucleotid_dbEstimate the maximum time your job will need to run. We recommend testimg initially with a < 0.5hr test run because Jobs set for 0.5 h or less depedendably run immediately in the "debug" queue.
Once you are sure the configuration is correct, you then increase the time. The reason is that jobs > 0.5 h are submitted to the "normal" queue, where jobs configured for 1 or a few hours times may
run sooner than jobs configured for the full 168 hours.
use_epaThis search is for Evolutionary placementperl""nucleotid_dbnucleotide dbperl "-db /projects/ps-ngbt/databases/blast/blast_ref_dbs/$value.expt" 45_8Sitsnuclsunucssumitssurpb1rpb2UNITE_public_tbasunaligned_itsunaligned_lsuunaligned_its_lsuunaligned_ssuunaligned_mitssuunaligned_rpb1unaligned_rpb2Sorry, unaligned databases are not valid for EPAperl$use_epa && ($value eq "unaligned_its" || $value eq "unaligned_lsu" || $value eq "unaligned_ssu" || $value eq "unaligned_mitssu" || $value eq "unaligned_rpb1" || $value eq "unaligend_rpb2") Please select a databaseperl!defined $nucleotid_dbchoose a nucleotide db for blastn, tblastn or tblastxgap_open_penaltyCost to open a gap (-gapopen; default=5)perl(defined $value)? " -gapopen $value":""5for blastn default=54gap_extend_penaltyCost to extend a gap (-gapextend; default=2)perl(defined $value)? " -gapextend $value":""25for blastn default=2; for blastp, blastx and tblastn default=1Limited values for gap existence and extension are supported for these three programs. Some supported and suggested values are:Existence Extension10 110 211 18 29 2(source: NCBI Blast page)all_results*