基于RNAfold预测RNA的二级结构（命令行版）

RNAfold子程序实质上是封装在ViennaRNA软件包中。

从https://www.tbi.univie.ac.at/RNA/#download网站下载ViennaRNA（version1.8.5）源码包，编译安装。安装过程不再细讲。不会的可参照我之前的博客。
运行以下命令批量预测和绘制RNA二级结构图

/usr/local/bin/RNAfold -p < test.fa > test.res
# /usr/local/bin/RNAplot -t 1 -o svg --pre aaa < test.res # 基于1的结果绘图
/biodata/02.software/ViennaRNA-1.8.5/Utils/relplot.pl test_sequenc_ss.ps test_sequenc_dp.ps > 5S_rss.ps # 可靠性信息以颜色注释的形式呈现给RNA二级结构图（beauty）
gv 5S_rss.ps # 查看二级结构
# 将postscript文件转成pdf图（另还可用ghostscript转化）
magick convert -density 300 5S_rss.ps 5S_rss.pdf # 生成二级结构pdf图
magick convert -density 300 5S_rss.pdf 5S_rss.png

说明：RNAfold会生成两个文件。一个为postscript文件，这种文件可以转换为PDF格式。

生成的二级结构图如下：

另产生的str文件内容如下：
该结果采用dot-bracket表示法标记二级结构，上述用法只给出了最佳的二级结构预测结果和对应的自由能。
关于dot-bracket表示法解释：
该表示方法就称之为dot-bracket notation, 其核心思想是利用配对的括号来表示碱基的互补配对，用连续的点号来表示茎环结构，对于下图所示的二级结构

RNA二级结构表示法Dot-Bracket notation如何理解：

对应的表达式如下

(((..((((...)))).)))

从第一个黑色圆点对应的碱基开始，一开始是3个配对碱基，所以先用3个(表示，接下来是茎环结构中未配对的两个碱基，用2个.表示，然后是4个配对碱基，再然后是未配对的3个碱基，再往后的配对碱基与前面左括号(表示的碱基相配对，所以用右括号)表示。

这种表示方式只采用了两种符号，而且都是计算机可以识别的符号，所以软件也可以识别，在很多软件和数据库中，给出的二级结构都会用这种方式来表示。

基于RNAfold预测RNA的二级结构（命令行版）相关推荐

RNAfold预测RNA的二级结构
欢迎关注"生信修炼手册"! 在mirdeep软件的分析结果中,会提供miRNA前体的二级结构,这个结果实际上是通过调用RNAfold来实现的,该软件是一个经典的预测RNA二级结构的 ...
csvtk：高效命令行版极简dplyr
写在前面什么时候写 csvtk 呀,csvtk 也借鉴了些 datamash 的东西. 之前写 datamash 的使用教程 linux 极简统计分析工具 datamash 必看教程,收到了一位读者 ...
Python 实现有道翻译命令行版
一.个人需求由于一直用Linux系统,对于词典的支持特别不好,对于我这英语渣渣的人来说,当看英文文档就一直卡壳,之前用惯了有道词典,感觉很不错,虽然有网页版的但是对于全站英文的网页来说并不支持.索性 ...
jdk解压版_命令行版的斗地主你玩过没？
相信大家都玩过斗地主游戏,或在现实中斗地主或在电脑和手机上斗地主,但你想过用命令行界面进行斗地主吗? 先来张图体验一下: 是不是觉得挺有意思,下面就带大家一起玩一下吧~部署命令行版斗地主 1 环境准备 ...
python刷网易云_牛逼了！用Python开发的命令行版网易云音乐，Github获8300颗星！...
大家好,我是程序员G哥最近在逛Github发现了一个非常有趣的库musicbox,是用纯Python打造的,收获了8300颗星.Python语言简单易学,好玩有趣,身边越来越多的小伙伴都开始学习Py ...
渗透工具TPLogScan：ThinkPHP3和5日志扫描工具，提供命令行版和BurpSuite插件版，尽可能发掘网站日志信息
Git地址: r3change/TPLogScan: ThinkPHP全日志扫描工具,命令行版和BurpSuite插件版 (github.com) ThinkPHP3和5日志扫描工具,提供命令行版和B ...
java编辑2048小游戏_Java 制作命令行版 2048小游戏
Java 命令行版 2048小游戏(2020年8月14日) 制作背景大二即将开学,从头开始学习了一个多月的java,对二维数组的操作稍微熟悉了一些.于是想做一个简单的2048来试一试. 众所周知,2 ...
康威生命游戏java_Java Python 康威生命游戏 - 命令行版
Java & Python 康威生命游戏 - 命令行版(2020年7月23日) 制作背景高二的时候看霍金的<大设计>最后几页的时候看到里面提到了康威生命游戏,介绍了它的规则,感觉 ...
Visual C C++ studio2019 自制鼠标点击器,窗口版和命令行版 210325
一窗口版 1.h和cpp h: DesktopMouseClick1.h #pragma once#include "resource.h" #include <threa ...

基于RNAfold预测RNA的二级结构（命令行版）

RNAfold子程序实质上是封装在ViennaRNA软件包中。

RNA二级结构表示法Dot-Bracket notation如何理解：

基于RNAfold预测RNA的二级结构（命令行版）相关推荐

最新文章

热门文章