NAME
hybrid-plot, hybrid-plot-ng — plot base pair probabilities
SYNOPSIS
hybrid-plot [OPTION]... PREFIX.br hybrid-plot-ng [OPTION]... PREFIX
DESCRIPTION
hybrid-plot and hybrid-plot-ng display plots of base pair probabilities generated with hybrid or hybrid-ss. hybrid-plot is an interactive plotting program that uses glut and OpenGL, while hybrid-plot-ng is a command-line program that produces output in PNG or Postscript form. In either case, a grid is displayed where the size of the dot at (i, j) indicates the likelihood that base i pairs with base j. (For plots generated by hybrid, i indexes the first sequence and j the second.) PREFIX specifies the sequence(s) that have been hybridized or folded and the .plot files to display.
The length of a side is increased from the length of the longer sequence to a round number, using the "four fifths rule": the length of a side is the smallest multiple of the largest power of 10 so that the length of the longer sequence is at least four fifths of the length of the side. For example, if the longer sequence is 873 characters long the plot will be 1000 by 1000; if the longer sequence is 773 characters long the plot will be 800 by 800.
In the plot and zoom windows of hybrid-plot, the left mouse button selects a pair. The indices and bases of the selected pair are displayed at the bottom of the window along with the probability of the selected pairing.
The middle mouse button is used to zoom. Dragging out a selection causes the zoom window to appear and display the selected area. The zoom area is square and snaps to the nearest pair.
The right mouse button displays a menu that can be used to toggle the zoom window, choose a temperature for which to plot, set a cutoff or end the program. (Setting a cutoff prevents pairs with probabilities below the cutoff from being displayed.) The menu is also used to select color modes.
COLORS
hybrid-plot and hybrid-plot-ng display points in colors corresponding to the probability of the base pair: red indicates the most likely pairs and magenta the least likely. When the color mode is set to "linear", the range [0, 1] is divided linearly into 26 subintervals and each is assigned to a color. When the color mode is set to "logarithmic", the range (epsilon, 1] is divided logarithmically into 25 subintervals and each is assigned to a color. ([0, epsilon] is the 26th interval.) This has the effect of emphasizing the distinctions between low-probability base pairs. Finally, when the color mode is set to "double logarithmic", the intervals [0, epsilon] and [1-epsilon, 1] correspond to pure magenta and pure red, respectively. The other 24 colors are arranged proportionally to log(2p) when p <= 0.5 and to -log(2-2p) when p > 0.5. This has the effect of emphasizing the distinctions among very low-probability base pairs and very high-probability base pairs. For the latter two cases, the value of epsilon can be set by the user.
OPTIONS
With the exception of --grid, the options below apply only to hybrid-plot-ng; hybrid-plot uses a menu-driven interface to change parameters.
- -t, --temperature=REAL
- plot probabilities at REAL °C. PREFIX.REAL.plot must exist. Default is 37.
- -c, --colors=linear|log|double
- set the method of assigning colors to basepairs to linear, logarithmic or double logarithmic. Default is double.
- -e, --epsilon=REAL
- set epsilon (for logarithmic or double logarithmic colors) to REAL. Default is 0.01.
- -g, --grid=NUMBER
- display grid lines every NUMBER bases. A value of 0 turns off the grid. Default is a reasonable value for the size of the plot.
- -d, --dot=NUMBER
- set (maximum) dot size to NUMBER. Default is a reasonable value for the size of the plot.
- -u, --top=NUMBER
- set the upper edge of the plot to base NUMBER. Default is 1.
- -l, --left=NUMBER
- set the left edge of the plot to base NUMBER. Default is 1.
- -s, --size=NUMBER
- set the size of the square to be displayed to NUMBER. By combining --top, --left and --size, it is possible to "zoom in" on any area of the plot.
- -f, --format=ps|png|gif|jpeg
- set the output format to Postscript or a raster image format. Default is ps.
- -i, --i=NUMBER
- in conjunction with --j, choose a basepair to appear selected, as though it were clicked in hybrid-plot.
- -j, --j=NUMBER
- in conjunction with --i, choose a basepair to appear selected, as though it were clicked in hybrid-plot.
- -p, --title=STRING
- set the plot title to STRING. Default is "'prefix'" or "'prefix1' vs. 'prefix2'".
- -r, --filter=on|off
- activate or deactivate the filter. With the filter enabled, a basepair (i, j) is only plotted if none of the basepairs (i - 1, j), (i - 1, j - 1), (i, j - 1), (i + 1, j), (i + 1, j + 1) or (i, j + 1) has a higher probability than (i, j). Default is off.
- -o, --cutoff=REAL
- When a cutoff is set, basepairs with probability less than REAL are not plotted.
OBSCURE OPTIONS
- -m, --machine
REFERENCES
Markham, N. R. and Zuker, M. (2008) UNAFold: software for nucleic acid folding and hybridization. In Keith, J. M., editor, Bioinformatics, Volume II. Structure, Functions and Applications, number 453 in Methods in Molecular Biology, chapter 1, pages 3-31. Humana Press, Totowa, NJ. ISBN 978-1-60327-428-9.
Other references which may be useful may be found at http://www.unafold.org/Dinamelt/dinamelt-references.php
SEE ALSO
AUTHORS
Nick Markham <markham@alum.rpi.edu> and Michael Zuker <zukerm@alum.mit.edu>
COPYRIGHT
Copyright (c) 2006-2021, Rensselaer Polytechnic Institute.
AVAILABILITY
UNAFold is available from http://www.unafold.org/Dinamelt/software/obtaining-unafold.php. Commercial use requires a license; see https://ipo.rpi.edu/invention/unafold-version-40.