REFOLD
NAME
refold.pl − refold using consensus structure as constraint
SYNOPSIS
refold.pl [−t thresh] file.aln alidot.ps refold.pl [−−turn length] file.aln file.alifold
DESCRIPTION
refold.pl reads an alignment in CLUSTAL format, and a consensus secondary structure, either in the form of the standard output of "RNAalifold", or in the form of a dot plot as produced by "RNAalifold −p". For each sequence in the alignment it writes the name, sequence, and constraint structure to stdout in a format suitable for piping into "RNAfold −C".
The constraint string is produced by removing from the consensus structure all gaps, as well as all pairs not compatible with the particular sequence. If the structure is read from a dot plot file, only pairs with probability > then some threshold (default 0.9) are used.
OPTIONS
−t thresh
use only pairs with p>thresh as constraint. Only applicable when reading consensus structure from a dot plot.
−−turn length
use only pairs closing a hairpin with looplength >length. This option allows for setting minimal loop length of hairpin loops (default: 3)
EXAMPLE
RNAalifold test.aln > test.alifold refold.pl test.aln test.alifold > test.cfold RNAfold −C < test.cfold