What it does
For each fragment of input sequence we use RNAshapes to create a set of structures.
The default parameters for example consider for each input fragment again a
window of size 40nt and 150nt with a window shift of 30%. This allows to
consider local structures as well as global structures for a fragment.
From each such RNAshape window we take the top 5 shreps (suboptimal
structures for the top 5 shapes) within 20% of the mfe energy of that window
and convert them into graphs.
As shape level (abstraction level) we use 3 for short sequences and 5 for
sequences >= 80nt. Please see also RNAshapes documentation for all these
terms.
Parameters
- input : Sequences in Fasta format.
- stack : This adds an additional vertex (type P) for each pair of stacked base-pairs and four edges
(type p) from each of the involved bases to the new vertex.
- abstr : Add abstract structure graphs to the single shrep graph instances.
- cue : Crop unpaired ends. If you give this flag, then the unpaired ends of each single structure are ignored. E.g. the structure
...(((...))).. becomes just (((...)))
- M : Max number of shreps that should be taken per window.
- c : Relative energy range, i.e. percentage (%) of MFE energy (RNAshapes) Use only one of -e and -c!
- wins : A list of window sizes to use. If none are given (empty string ''), then the entire sequence is taken with no windows. Each window > 1 required!
- shift : The shift of the window, relative to the window size given in percent. So you give which percent of the window size shall be
used for the shift. Of course the shift is rounded down to the
nearest whole number.
Example 20 % of a window 150 would result in a step size of 30 nt.
It is a relative parameter, as you can give different window sizes.
If you do not give this parameter there is a default shift of 1 nt.
- u : Ignore unstable structures (RNAshapes). This option filters out closed structures with positive free energy.
- seq-graph-t : Add for each 't #' a graph which contains no structure
- group : Combine/group that number of input seqs into 1 gspan file output name is then '<INT>.group.gspan.bz2'