RNAduplex sequence input format

Question

[ete]

Hello!

I am using RNAduplex from Vienna package and trying to find out in which format the two sequences has to be passed to the cmd. The manual only say "reads two RNA sequences from stdin or <filename>". I used a filename but this might lead to IOErrors if my query sequence is very long. Does anybody know how to give the two sequences?

Thanks in advance!

[ete]

Hello Asaf!

I want to start RNAduplex from a cmd in a Python script by entering the sequence directly (not a file).

I tried already but without success:

'RNAduplex ATG GTC'
'RNAduplex ATG\nGTC'

Thanks

[Asaf]

echo "AGCT\nACGT\n" | RNAduplex

or I can send you a script that open a process, write sequences to the process and read the results, works with threading as well if you have a lot of sequences

[ete]

Sorry, Asaf I am still having problems:

In Python I do something like this:

cmd = 'RNAduplex < duplex.seq'
p = subprocess.Popen(cmd,shell=True,stdout=subprocess.PIPE)

But I want to give the raw sequences as parameters. I tried something like:

cmd = "AGCT\nACGT\n" RNAduplex

but it does not work. Any hints?

[Asaf]

cmd = 'echo "AGCT\nACGT\n" | RNAduplex'

you can either have stdin=subprocess.PIPE and have:

p.stdin.write("ACGACT\nACGAGCTG\n")
print p.stdout.next()

Answer 1

[Asaf]

You should give two sequences with or without the headers (>seq1...)

e.g.

$RNAduplex

Input two sequences (one line each); @ to quit
....,....1....,....2....,....3....,....4....,....5....,....6....,....7....,....8
AGCGAGCTGTGTGACTGTGACTGATGCTG
ACGAGCTGATGCTGATCTATGCGTATGCTGATGC
lengths = 29,34
((((.((.(((((((.(((.(....(((.&.))).).))).).)).)))))).)))).   1,29  :   3,30  (-10.00)

If your sequences are in a file you should execute:

$RNAduplex < input_file

each pair of RNA sequences will be evaluated for their binding.