从植物sRNA-seq数据中de novo预测miRNA的工具--miRDeep-P2
写在前面
- 植物miRNA与动物miRNA在多种特征上存在区别,所以直接使用miRDeep或miRDeep-2去分析植物miRNA数据是非常不负责任的。
- miRDeep-P2是de novo预测,不是通过同源序列匹配
- 得到的结果需要依据已有的miRNA家族命名,获得的新的miRNA家族也需要命名,这部分分析代码并不包含在miRDeep-P2中
- 关于miRDeep-P2及miRNA的问题欢迎交流qq:734256624
软件下载地址:https://sourceforge.net/projects/mirdp2/files/,推荐下载最新版
参数
NECESSARY:
-g/–genome reference genome file in fasta format
-x/–index /prefix bowtie index file corresponding to reference genome(是bowtie的index,不是bowtie2)
-f/–fasta formatted sRNA-seq file in fasta format
-q/–fastq unformatted sRNA-seq file in fastq
-i/–input input file(s), multiple files are seperated by comma (,)
-b/–batch list file of input files
-o/–output path to output folder (输出文件夹路径最后不要加斜线"/",输出文件夹需要自己提前构建)
OPTIONAL:
-L/–locate THRESHOLD of different locations one read can map to, default is 15 (个人认为,如果是4倍体,这里的参数可以设置为30试一试)
-N/–length THRESHOLD of length of precursor candidates, default is 300 (根据目前的植物miRNA标准,这里不建议修改)
-M/–mismatch THRESHOLD of allowed mismatches for bowtie mapping, default is 0
-R/–rpm THRESHOLD of reads rpm when filtering, default is 10 (这里可以适当放低,RPM为5也是可以接受的)
-p/–thread number of thread used for RNAfold, default is 1
–large-index use this option when using large bowtie index, e.g. when using *.ebtwl files (如果基因组很大,bowtie构建的index后缀为 ebtwl,则需要加上 --large-index参数)
-h/–help print this help
示例:
- 构建bowtie index(不是bowtie2)
bowtie genome.fasta genome.index
- 运行miRDeep-P2
i. 如果基因组不大,则可以看到bowtie构建的index后缀为ebtw,不需要加–large-index参数;如果index后缀为ebtwl,则需要加上–large-index参数。
ii. 输出文件夹需要提前构建
iii. -o后加的输出文件夹最后不要加斜线"/"mkdir miRDP2_out miRDP2-v1.1.4_pipeline.bash -g genome.fasta -x genome.index -q -i sRNA-seq_trimmed.fq -o miRDP2_out -R 5
输出文件夹中后缀为P_prediction的文件为miRNA的信息文件
从植物sRNA-seq数据中de novo预测miRNA的工具--miRDeep-P2相关推荐
- Nanopore测序在基因组 de novo中的应用
Nanopore测序在基因组 de novo中的应用 自1977年第一代sanger测序问世来,经过几十年的发展,测序技术得到了极大的发展. 从第一代测序到第二代测序再到第三代测序,测序技术的每一次变 ...
- The power of single molecule real-time sequencing technology in the de novo assembly of a eukaryotic
The power of single molecule real-time sequencing technology in the de novo assembly of a eukaryotic ...
- Efficient Hybrid De Novo Error Correction and Assembly for Long Reads
Efficient Hybrid De Novo Error Correction and Assembly for Long Reads 长read的高效的混合从头纠错和装配 Abstr ...
- Trinity简介(1)--用于无参考基因组的转录组de novo组装
一. Trinity简介 Trinity,是由 the Broad Institute 开发的转录组de novo组装软件,由三个独立的软件模块组成: Inchworm,Chrysalis和Butte ...
- NGS项目四:高通量测序在植物生物胁迫研究中的应用
在植物生物胁迫研究方面,本文通过作者所精密相关课题组的研究,比较传统研究方法和高通量测序方法的优劣. Sun X, Tan Q, Nie Z,et al. Differential Expressio ...
- 三代测序之微生物基因组 de novo 测序
三代测序之微生物基因组 de novo 测序 2016-05-20 编辑:诺禾致源 微生物基因组 de novo 测序 产品优势 • 单分子,长序列,完美跨越细菌全部重复单元 • 真菌 20Kb ...
- Tools and Strategies for Long-Read Sequencing and De Novo Assembly of Plant Genomes
Tools and Strategies for Long-Read Sequencing and De Novo Assembly of Plant Genomes 用于植物基因组长读测序和从头组装 ...
- Evaluation of hybrid and non-hybrid methods for de novo assembly of nanopore reads
Evaluation of hybrid and non-hybrid methods for de novo assembly of nanopore reads 混合和非混合方法对nanopore ...
- Oxford Nanopore sequencing, hybrid error correction, and de novo assembly of a eukaryotic genome
Oxford Nanopore sequencing, hybrid error correction, and de novo assembly of a eukaryotic genome 牛津纳 ...
最新文章
- LINQ根据某字段GroupBy
- 交流继电器rc吸收电路_使用固态继电器的注意事项
- Linux 命令之 declare -- 声明或显示 shell 变量
- 使用final修饰局部变量???
- java基础知识点整理一
- boost::bind with ros topic,ros中subscribe用boost::bind绑定多个参数
- Java基础篇:什么是FileWriter
- Go database/sql 教程
- adminLTE的介绍
- 微信提示在客户端提交验证_微信提示非常用设备要求输入短信验证码解决方法...
- 带你逐步深入了解SSM框架——淘淘商城项目之redis缓存
- 什么叫死区时间_关于pwm死区时间的介绍
- 网络原理4 数据链路层
- 【vue3.0学习】父传子 子传父
- google官方提供的颜色表。
- 大卫艾克~雄狮醒来2:实相的本质(中字)
- 【软件应用】word数学公式插件mathtype安装
- MonkeyPatchWarning: Monkey-patching ssl after ssl has already been imported may lead to errors
- 【Linux认证考试分数线多少】
- 使用cmd命令行或运行框进行关机重启操作