What Is The Input Format For Rnaplot In The Vienna Package?
2
1
Entering edit mode
11.6 years ago

Hi all,

I am going to use RNAplot (in Vienna RNA package) to draw a figure of my RNA secondary structure. I however, I do not know what should be the input. Could anyone please tell me? Thanks!

• 7.4k views
ADD COMMENT
3
Entering edit mode
11.6 years ago

Here is an example of RNAplot's input format:

>seq1
ACTGAGCTTAGCGATTAGGC
(((..((((...)))).)))

This format is known as "dot-bracket notation", and it is used to represent a RNA secundary structure by many other packages, and is also used by many Vienna scripts. I can not find much documentation online, but here there is some.

Note that you can calculate this structure by using the RNAfold tool:

$: echo ">seq1" > myseq.txt
$: echo "ACTGAGCTTAGCGATTAGGC" >> myseq.txt
$: cat myseq.txt
>seq1
ACTGAGCTTAGCGATTAGGC

$: cat myseq.txt | RNAfold > myfold.txt
$: cat myfold.txt
>seq1
ACUGAGCUUAGCGAUUAGGC
.....(((((.....))))) ( -2.70)

# Note: to use RNAplot, you have to remove the score first:
$: sed -i 's/ \(.*\)$//' myfold.txt
$: cat myfold.txt
>seq1
ACUGAGCUUAGCGAUUAGGC
.....(((((.....)))))

Once you have the input file, you can produce the plot by simply calling RNAplot:

cat myfold.txt | RNAplot -o svg
ADD COMMENT
0
Entering edit mode

Thanks very much!!!!!!

ADD REPLY
0
Entering edit mode
11.6 years ago

http://www.tbi.univie.ac.at/~ronny/RNA/RNAplot.html

RNAplot 2.1.1

Draw RNA Secondary Structures

The program reads RNA sequences and structures from stdin in the format as produced by RNAfold and produces drawings of the secondary structure graph.

ADD COMMENT
0
Entering edit mode

Thanks !a lot for your help

ADD REPLY

Login before adding your answer.

Traffic: 766 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6