总是记不住blast的参数设置,现在把参数及用法放在此处,方便以后查阅~
Nucleotide-NucleotideBLAST 2.2.28+
usage
blastn [-h] [-help] [-import_search_strategyfilename]
[-export_search_strategy filename] [-tasktask_name] [-db database_name]
[-dbsize num_letters] [-gilist filename][-seqidlist filename]
[-negative_gilist filename] [-entrez_queryentrez_query]
[-db_soft_mask filtering_algorithm][-db_hard_mask filtering_algorithm]
[-subjectsubject_input_file] [-subject_loc range] [-query input_file]
[-outoutput_file] [-evalueevalue] [-word_size int_value]
[-gapopen open_penalty] [-gapextendextend_penalty]
[-perc_identity float_value] [-xdrop_ungapfloat_value]
[-xdrop_gap float_value] [-xdrop_gap_finalfloat_value]
[-searchsp int_value] [-max_hsps_per_subjectint_value] [-penalty penalty]
[-reward reward] [-no_greedy][-min_raw_gapped_score int_value]
[-template_type type] [-template_lengthint_value] [-dust DUST_options]
[-filtering_db filtering_database]
[-window_masker_taxid window_masker_taxid]
[-window_masker_db window_masker_db][-soft_masking soft_masking]
[-ungapped] [-culling_limit int_value][-best_hit_overhangfloat_value]
[-best_hit_score_edge float_value][-window_size int_value]
[-off_diagonal_range int_value] [-use_indexboolean] [-index_name string]
[-lcase_masking] [-query_loc range] [-strandstrand] [-parse_deflines]
[-outfmt format] [-show_gis] [-num_descriptionsint_value]
[-num_alignments int_value] [-html][-max_target_seqs num_sequences]
[-num_threadsint_value] [-remote] [-version]
OPTIONAL ARGUMENTS
-h
Print USAGE and DESCRIPTION; ignore allother parameters
-help
Print USAGE, DESCRIPTION and ARGUMENTS; ignoreall other parameters
-version
Print version number; ignore otherarguments
*** Input query options
-query
Input file name
Default = `-'
-query_loc
Location on the query sequence in 1-based offsets(Format: start-stop)
-strand
Query strand(s) to search againstdatabase/subject
Default = `both'
*** General search options
-task
'megablast' 'rmblastn' >
Task to execute
Default = `megablast'
-db
BLAST database name
* Incompatible with: subject,subject_loc
-out
Output file name
Default = `-'
-evalue
Expectation value (E) threshold for savinghits
Default = `10'
-word_size =4>
Word size for wordfinder algorithm (length ofbest perfect match)
-gapopen
Cost to open a gap
-gapextend
Cost to extend a gap
-penalty
Penalty for a nucleotide mismatch
-reward =0>
Reward for a nucleotide match
-use_index
Use MegaBLAST database index
-index_name
MegaBLAST database index name
*** BLAST-2-Sequences options
-subject
Subject sequence(s) to search
* Incompatible with: db, gilist,seqidlist, negative_gilist,
db_soft_mask, db_hard_mask
-subject_loc
Location on the subject sequence in 1-basedoffsets (Format: start-stop)
* Incompatible with: db, gilist,seqidlist, negative_gilist,
db_soft_mask, db_hard_mask, remote
*** Formatting options
-outfmt
alignment view options:
0 = pairwise,
1 = query-anchored showingidentities,
2 = query-anchored no identities,
3 = flat query-anchored, showidentities,
4 = flat query-anchored, noidentities,
5 = XML Blast output,
6 = tabular,
7 = tabular with comment lines,
8 = Text ASN.1,
9 = Binary ASN.1,
10 = Comma-separated values,
11 = BLAST archive format (ASN.1)
Options 6, 7, and 10 can be additionallyconfigured to produce
a custom format specified by space delimitedformat specifiers.
The supported format specifiers are:
qseqid means Query Seq-id
qgi means Query GI
qacc means Queryaccesion
qaccver means Queryaccesion.version
qlen means Query sequencelength
sseqid means Subject Seq-id
sallseqid means All subject Seq-id(s), separatedby a ';'
sgi means Subject GI
sallgi means All subject GIs
sacc means Subjectaccession
saccver means Subjectaccession.version
sallacc means All subjectaccessions
slen means Subject sequencelength
qstart means Start of alignment inquery
qend means End of alignmentin query
sstart means Start of alignment insubject
send means End of alignmentin subject
qseq means Aligned part ofquery sequence
sseq means Aligned part ofsubject sequence
evalue means Expect value
bitscore means Bit score
score means Raw score
length means Alignment length
pident means Percentage ofidentical matches
nident means Number of identicalmatches
mismatch means Number of mismatches
positive means Number of positive-scoringmatches
gapopen means Number of gapopenings
gaps means Total number ofgaps
ppos means Percentage ofpositive-scoring matches
frames means Query and subjectframes separated by a '/'
qframe means Query frame
sframe means Subject frame
btop means Blast tracebackoperations (BTOP)
staxids means Subject Taxonomy ID(s),separated by a ';'
sscinames means Subject Scientific Name(s),separated by a ';
scomnames means Subject Common Name(s),separated by a ';'
sblastnames means Subject Blast Name(s),separated by a ';'
(in alphabetical order)
sskingdoms means Subject Super Kingdom(s),separated by a ';'
(in alphabetical order)
stitle means Subject Title
salltitles means All Subject Title(s), separatedby a '<>'
sstrand means Subject Strand
qcovs means Query Coverage PerSubject
qcovhsp means Query Coverage PerHSP
When not provided, the default value is:
'qseqid sseqid pident length mismatch gapopenqstart qend sstart send
evalue bitscore', which is equivalent to thekeyword 'std'
Default = `0'
-show_gis
Show NCBI GIs in deflines?
-num_descriptions =0>
Number of database sequences to show one-linedescriptions for
Not applicable for outfmt > 4
Default = `500'
* Incompatible with: max_target_seqs
-num_alignments =0>
Number of database sequences to show alignmentsfor
Default = `250'
* Incompatible with: max_target_seqs
-html
Produce HTML output?
*** Query filtering options
-dust
Filter query sequence with DUST (Format: 'yes','level window linker', or
'no' to disable)
Default = `20 64 1'
-filtering_db
BLAST database containing filtering elements(i.e.: repeats)
-window_masker_taxid
Enable WindowMasker filtering using a TaxonomicID
-window_masker_db
Enable WindowMasker filtering using this repeatsdatabase.
-soft_masking
Apply filtering locations as soft masks
Default = `true'
-lcase_masking
Use lower case filtering in query and subjectsequence(s)?
*** Restrict search or results
-gilist
Restrict search of database to list of GI's
* Incompatible with: negative_gilist,seqidlist, remote, subject,
subject_loc
-seqidlist
Restrict search of database to list ofSeqId's
* Incompatible with: gilist,negative_gilist, remote, subject,
subject_loc
-negative_gilist
Restrict search of database to everything exceptthe listed GIs
* Incompatible with: gilist, seqidlist,remote, subject, subject_loc
-entrez_query
Restrict search with the given Entrez query
* Requires: remote
-db_soft_mask
Filtering algorithm ID to apply to the BLASTdatabase as soft masking
* Incompatible with: db_hard_mask,subject, subject_loc
-db_hard_mask
Filtering algorithm ID to apply to the BLASTdatabase as hard masking
* Incompatible with: db_soft_mask,subject, subject_loc
-perc_identity
Percent identity
-culling_limit =0>
If the query range of a hit is enveloped by thatof at least this many
higher-scoring hits, delete the hit
* Incompatible with: best_hit_overhang,best_hit_score_edge
-best_hit_overhang =0 and =<0.5)>
Best Hit algorithm overhang value (recommendedvalue: 0.1)
* Incompatible with: culling_limit
-best_hit_score_edge =0 and=<0.5)>
Best Hit algorithm score edge value (recommendedvalue: 0.1)
* Incompatible with: culling_limit
-max_target_seqs =1>
Maximum number of aligned sequences tokeep
Not applicable for outfmt <= 4
Default = `500'
* Incompatible with: num_descriptions,num_alignments
*** Discontiguous MegaBLAST options
-template_type
Discontiguous MegaBLAST template type
* Requires: template_length
-template_length
Discontiguous MegaBLAST template length
* Requires: template_type
*** Statistical options
-dbsize
Effective length of the database
-searchsp =0>
Effective length of the search space
-max_hsps_per_subject =0>
Override maximum number of HSPs per subject tosave for ungapped searches
(0 means do not override)
Default = `0'
*** Search strategy options
-import_search_strategy
Search strategy to use
* Incompatible with: export_search_strategy
-export_search_strategy
File name to record the search strategyused
* Incompatible with: import_search_strategy
*** Extension options
-xdrop_ungap
X-dropoff value (in bits) for ungappedextensions
-xdrop_gap
X-dropoff value (in bits) for preliminary gappedextensions
-xdrop_gap_final
X-dropoff value (in bits) for final gappedalignment
-no_greedy
Use non-greedy dynamic programmingextension
-min_raw_gapped_score
Minimum raw gapped score to keep an alignment inthe preliminary gapped and
traceback stages
-ungapped
Perform ungapped alignment only?
-window_size =0>
Multiple hits window size, use 0 to specify 1-hitalgorithm
-off_diagonal_range =0>
Number of off-diagonals to search for the 2ndhit, use 0 to turn off
Default = `0'
*** Miscellaneous options
-parse_deflines
Should the query and subject defline(s) beparsed?
-num_threads=1>
Number of threads (CPUs) to use in the BLASTsearch
Default = `1'
* Incompatible with: remote
-remote
Execute search remotely?
* Incompatible with: gilist, seqidlist,negative_gilist, subject_loc,
num_threads
|