文献笔记七十一: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文件。可以参考。
欢迎大家关注我的公众号
小明的数据分析笔记本
- 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 数组属性和方法
- JavaWeb——JSON语法讲解与Jackson解析器完成JSON数据与Java对象的转换(应用Ajax与JSON实现校验用户名是否在的功能)
- JavaWeb——JQuery之高级案例实战(打开网页自动弹出广告效果、抽奖效果实现)
- JavaWeb——JQuery之高级操作应用及实践案例总结(动画、遍历、事件绑定)
- Java——简单Java类深入(数据表与简单Java类、一对多映射、双向一对多映射、多对多映射)
- Java——代码块(普通块、构造块、静态块)
- JavaWeb——JQuery之五种选择器的应用及实践案例总结(基本选择器、层级选择器、属性选择器、过滤选择器、表单过滤选择器)
- Java——引用传递实例分析(进阶分析、对象比较、类与类的关联实现)
- JavaWeb——JQuery之基本概述与快速入门实践总结(JQuery各版本区别、JQuery对象和JS对象的区别与转换)
- Java——动态代理设计模型概述(实现步骤、增强方式)与代理商采购电脑模拟程序实战
- Java——this关键字(调用本类属性、调用本类方法、表示当前对象)
- 5分钟玩转Lighthouse|零基础也能拥有WordPress个人博客
- Java——String类常用方法总结,看这一篇就够啦(比较、查找、截取、替换、拆分、其他)
- Java——深入分析类与对象(封装性、构造方法与匿名对象、简单Java类开发原则)
- Java——类与对象(基本概念、定义、内存分析、引用传递)
- Java——Stream数据流