文献笔记七十一:REDO根据vcf文件检测植物细胞器基因组RNA编辑位点
时间:2022-07-28
本文章向大家介绍文献笔记七十一:REDO根据vcf文件检测植物细胞器基因组RNA编辑位点,主要内容包括其使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
REDO: RNA editing detection in plant organelles based on variant calling results
期刊
Journal of computational biology
2018
四区, 影响因子1点多
完成单位
Chinese Academy of sciences
Beijing Institute of Genomics
主要功能
利用转录组数据比对到细胞器参考基因组得到vcf文件,比对工具使用GSNAP或者BWA,检测变异使用GATK或者SAMtools.
然后REDO这个软件 利用得到的vcf文件检测可能的RNA编辑位点。要求的输入文件是 参考序列fasta格式,注释文件 tbl格式,vcf格式的变异文件,还需要指定输出文件的后缀名,然后就是很多过滤参数,有默认设置,也可以自己指定。但这些参数的意思我还得在仔细看看。
软件是perl脚本,论文中写道 鉴定 注释和统计RNA编辑位点使用perl,画图调用的是R语言。
软件的下载链接
https://sourceforge.net/projects/redo/
直接解压出来就可以使用
使用的时候可能会遇到报错
Can't locate Text/NSP/Measures/2D/Fisher/left.pm in @INC (you may need to install the Text::NSP::Measures::2D::Fisher::left module) (@INC contains: /etc/perl /usr/local/lib/x86_64-linux-gnu/perl/5.26.1 /usr/local/share/perl/5.26.1 /usr/lib/x86_64-linux-gnu/perl5/5.26 /usr/share/perl5 /usr/lib/x86_64-linux-gnu/perl/5.26 /usr/share/perl/5.26 /usr/local/lib/site_perl /usr/lib/x86_64-linux-gnu/perl-base) at REDO.pl line 11.
BEGIN failed--compilation aborted at REDO.pl line 11.
是因为缺少模块 Text::NSP::Measures::2D::Fisher::left
直接使用命令
cpan install Text::NSP::Measures::2D::Fisher::left
安装就可以了
解压出来的文件带了测试数据,试了一下
perl REDO.pl -g example/cp/data/cp.fsa -v example/cp/data/SRR1063407_cp.vcf -t example/cp/data/cp.tbl -o AA
当前目录就会多出来7个文件,7个文件怎们看我还得再研究一下。文件里有个图
这个图应该怎么理解,也得花时间想一想
软件里还有一个readme.txt文件,还介绍了如何使用bwa+samtools得到vcf文件。可以参考。
欢迎大家关注我的公众号
小明的数据分析笔记本
- Jetty入门
- java序列化反序列化深入探究
- 前后端分离开发模式下后端质量的保证 —— 单元测试
- java如何获取一个对象的大小
- JDK1.7新特性(3):java语言动态性之脚本语言API
- JDK1.7新特性(4):java语言动态性之反射API
- ASP.NET Core集成现有系统认证
- Redis(2):常用命令详解
- C#移动跨平台开发(2)Xamarin移动跨平台解决方案是如何工作的?
- Ruby(3):基本语法中
- Python(3):文件读写与异常
- 向ASP.NET Core迁移
- Gitlab CI 自动部署 asp.net core web api 到Docker容器
- 从XMLHttpRequest请求响应里getResponseHeader(header)报错:Refused to get unsafe header "**" 问题解决
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- 数据源管理 | 搜索引擎框架,ElasticSearch集群模式
- PHP的CLI命令行运行模式浅析
- 基于Pytorch构建三值化网络TWN
- 从零学Paddle系列-1 Paddle框架CNN相关API详解
- 智能搜索模型预估框架的建设与实践
- 1,Jupyter NoteBook 常用魔法命令
- 60行代码徒手实现深度神经网络
- 30行代码徒手实现logistic回归
- 8,模型的训练
- 在腾讯云上部署科学计算软件Amber
- 手把手教你搭建一个灰度发布环境
- Kibana: 如何使用 Search Bar
- 「PHP」以nginx、php-cgi为例,把nginx、php-cgi安装为Windows系统服务
- 聊聊dubbo-go的GenericFilter
- 知新 | koa框架入门到熟练第二章