-h Give a short summary of command line options.
-bhit_number
Maximal number of hits to show. hit_number hits with increasing minimum free energy (reminder:
larger energies are worse) are shown, unless the -e option is used and the energy cut-off has been
exceeded (see -e option below) or there are no more hits. Hits may only overlap at dangling bases
(5' or 3' unpaired end of target).
-c Produce compact output. For each target/query pair one line of output is generated. Each line is a
colon (:) separated list of the following fields: target name, query name, minimum free energy,
position in target, alignment line 1, line 2, line 3, line 4. If a target or a query is given on
command line (ie. no -t or -q respectively), its name in the output will be "command line".
-dxi,theta
xi and theta are the position and shape parameters, respectively, of an extreme value distribution
(evd). p-values of duplex energies are assumed to be distributed according to such an evd. For a
length normalised energy en, we have P[X <= en] = 1 - exp(-exp(-(-en-xi)/theta)), where en =
e/log(m*n) with m and n being the lengths of the target and the query, respectively. If the -d
option is omitted, xi and theta are estimated from the maximal duplex energy of the query,
assuming a linear dependence. The parameters of this linear dependence are coded into the program,
but the option -s has to be given to choose from the appropriate set. Note that the evd is
mirrored, since good mfes are large negative values.
-sset_name
Used for quick estimate of extreme value distribution parameters (see -d option above). Tells
RNAhybrid which target dataset to assume. Valid parameters are 3utr_fly, 3utr_worm and 3utr_human.
-eenergy_cutoff
Hits with increasing minimum free energy (reminder: larger energies are worse) less than or equal
to energy_cutoff are shown, unless the -b option is used and the number of already reported hits
has reached the maximal hit_number (see -b option above). Hits may only overlap at dangling bases
(5' or 3' unpaired end of target).
-pp-value_cutoff
Only hits with p-values not larger than p-value_cutoff are reported. See also options -d and -s.
-ffrom,to
Forces all structures to have a helix from position from to position to with respect to the query.
The first base has position 1.
-mmax_target_length
The maximum allowed length of a target sequence. The default value is 2000. This option only has
an effect if a target file is given with the -t option (see below).
-nmax_query_length
The maximum allowed length of a query sequence. The default value is 30. This option only has an
effect if a query file is given with the -q option (see below).
-uiloop_upper_limit
The maximally allowed number of unpaired nucleotides in either side of an internal loop.
-vbloop_upper_limit
The maximally allowed number of unpaired nucleotides in a bulge loop.
-g(ps|png|jpg|all)
Produce a plot of the hybridisation, either in ps, png or jpg format, or for all formats together.
The plots are saved in files whose names are created from the target and query names
("command_line" if given on the command line). This option only works, if the appropriate graphics
libraries are present.
-ttarget_file
Each of the target sequences in target_file is subject to hybridisation with each of the queries
(which either are from the query_file or is the one query given on command line; see -q below).
The sequences in the target_file have to be in FASTA format, ie. one line starting with a > and
directly followed by a name, then one or more following lines with the sequence itself. Each
individual sequence line must not have more than 1000 characters. If no -t is given, either the
last argument (if a -q is given) or the second last argument (if no -q is given) to RNAhybrid is
taken as a target.
-qquery_file
See -t option above.
Install Now
PNG Icons
