Parsing RNAhybrid output
1
0
Entering edit mode
5.2 years ago
baurumon ▴ 30

Hi, I have used RNAhybride to predict miRna to target desired sequence. Can anyone help me in the way to parse output by column as below.

target: NC_031978.2:30702661-30704116 | miRNA : oni-miR-10613 | mfe: -26.5 kcal/mol | p-value: 0.095468

enter image description here

RNA-Seq • 1.4k views
ADD COMMENT
1
Entering edit mode

Please do not add screenshots of plain text - it is counterproductive. Instead, copy-paste the text and use the Code Formatting option to present the text well.

code_formatting

ADD REPLY
2
Entering edit mode
5.2 years ago
JC 13k

It is easy on Perl (or Python):

#!/usr/bin/perl

use strict;
use warnings;

while (<>) {
    chomp;
    if (/^target/ or /^miRNA/ or /^mfe/) { print "$_ | "; }
    elsif (/^p-value/) { print "$_\n"; }
}

use as: perl parseRNAfold < IN > OUT

ADD COMMENT

Login before adding your answer.

Traffic: 1830 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