编者按

欢迎关注微信公众号:基因学苑,更多精彩内容等你发掘!

前面写了专题《手把手教你生物信息分析平台搭建》,然后又介绍了很多《生物神奇网站》资源,也介绍了《生物信息之独孤九剑》Linux操作。那么万事俱备,就开始学习生物信息吧。所以,我们开始新的篇章——《生物信息百jia软件》。百Jia是什么意思呢?可以是百佳,也可以是百家,还可以是百加。从100家中选择100款优秀软件,掌握这些软件,就可以扩展出更多内容,这就是百Jia。

一、功能分类:
原核生物基因预测

二、软件官网:
compbio.ornl.gov/prodigal/

三、软件介绍:
prodigal的全称是Prokaryotic Dynamic Programming Genefinding Algorithm,原核的动态编程基因查找算法,prodigal主要应用于细菌和古生菌的基因预测,不能用于真核生物,如果要对meta样品做基因预测,prodigal还专门提供了meta的版本。 除此之外,prodigal还支持在线提交序列的方式来预测基因预测。也非常的易于使用。而且相对与glimmer基因预测工具,prodigal更加好用,只需一步即可,而且,软件可以直接输出基因的核酸序列并翻译出的相应的氨基酸序列,这对很多初学者来说是非常方便的。

四、下载安装:
wget https://codeload.github.com/hyattpd/Prodigal/tar.gz/v2.60
tar -zxvf Prodigal-2.6.1.tar.gz
make install
五、软件使用:
软件的使用比较容易,比glimmer要容易很多,只需一步即可,直接敲prodigal命令就会弹出软件帮助信息。

-a 是输出氨基酸文件-c 不允许基因一边断开,也就是要求完整的orf,有起始和终止结构
-d 输出预测基因的序列文件
-f 选择输出文件格式,有gbk,gff,和sco格式可供选择
-g 指定密码子,原核为第11套
-i 输入文件,即需要预测的基因组序列文件
-m 屏蔽基因组中的N碱基
-o 输出文件,默认为屏幕输出
-p 选择方式,是单菌还是meta样品
-q 不输错错误信息到屏幕
-t 指定训练集
-s 输出所有潜在基因以及分值到一个文件中

六、使用案例:
prodigal -a ref.pep -d ref.cds -f gff -g 11 -o ref.gff -p single -s ref.stat -i ref.fna >prodigal.log
七、注意事项:
1、对于原核生物基因预测,我们还需要注意一个问题,就是当样品为支原体时,在密码子选择上要修改一下,因为在支原体中遗传密码中的密码子UGA能够编译成色胺酸,而不是一般的乳白色终止码。

生物信息百Jia软件(六):prodigal相关推荐

  1. 生物信息百Jia软件(七):glimmer

    编者按 欢迎关注微信公众号:基因学苑,更多精彩内容等你发掘! 前面写了专题<手把手教你生物信息分析平台搭建>,然后又介绍了很多<生物神奇网站>资源,也介绍了<生物信息之独 ...

  2. 生物信息百Jia软件(十七):flash

    欢迎订阅WX众号:基因学苑,更多精彩内容等你发掘! 基因学苑Q群:32798724 编者按 前面写了专题<手把手教你生物信息分析平台搭建>,然后又介绍了很多<生物神奇网站>资源 ...

  3. 生物信息百jia软件(25):quast

    欢迎订阅WX众号:基因学苑,更多精彩内容等你发掘! 基因学苑Q群:32798724 去年这时候制定了一个计划,写100篇生物软件的教程推文,取名字"生物信息百jia软件",本来以为 ...

  4. 生物信息百Jia软件(22):featureCounts

    欢迎订阅WX众号:基因学苑,更多精彩内容等你发掘! 基因学苑Q群:32798724 编者按 虽然短序列比对工具subread的使用没有bwa和hisat2流行,但是软件包中的featureCounts ...

  5. 生物信息百Jia软件(二十):lumpy

    欢迎订阅WX众号:基因学苑,更多精彩内容等你发掘! 基因学苑Q群:32798724 通哥点评 基因组结构变异SV的检测一直都是当前基因组研究中非常困难的课题,由于二代测序多长过短,因此无论在优化算法, ...

  6. 生物信息百Jia软件(十四):velvet

    欢迎订阅WX众号:基因学苑,更多精彩内容等你发掘! 基因学苑Q群:32798724 通哥点评 velvet是最早的用于二代测序的序列拼接工具,拼接效果还可以,不过我觉得不如SOAPdenovo效果好. ...

  7. 生物信息百Jia软件(21):hisat2

    欢迎订阅WX众号:基因学苑,更多精彩内容等你发掘! 基因学苑Q群:32798724 编者按 现在的软件更新实在太快了,光是转录组用来比对的软件就有很多,bowtie, tophat,star,现在又有 ...

  8. 生物信息百Jia软件(28):canu

    欢迎订阅WX众号:基因学苑,更多精彩内容等你发掘! 基因学苑Q群:32798724 上一次介绍了三代拼接工具falcon,但是falcon这款无论是从安装还是使用都比较麻烦,对于新手是一个很大的考验, ...

  9. 阿里巴巴开源技术汇总:115个软件(六)

    阿里巴巴开源技术汇总:115个软件(六) 摘要: 今天我们的阿里巴巴开源技术汇总更新到了第六期,也算是这个系列的最后一期了.回顾这一百多个开源项目,它们凝聚了许多阿里人无数的心血和汗水,但我们相信阿里 ...

最新文章

  1. CentOS7.4 安装mongodb
  2. Android开发之Service与IntentService的区别与使用场景(源代码剖析)
  3. HBase在淘宝的应用和优化
  4. 2015年秋季个人阅读计划
  5. CTabCtrl控件标签的相关设置
  6. 轻松学习理解ACL访问控制列表(转)
  7. 计算机光驱故障分析,光驱故障示例解决
  8. js----------------------运算符
  9. 软考时间管理思维导图
  10. 天猫商城在线购物系统
  11. 表单验证工具类ValidationUtils
  12. 设置360浏览器的背景为护眼模式(浅豆绿色)
  13. 微信二维码活码制作管理系统源码+独立版网站
  14. java String工具类/字符串工具类 StringUtil
  15. Pycharm中 No module ‘自定义模块‘ 的问题(mmrotate)
  16. 从事嵌入式行业年薪有多少,你和高薪究竟差了哪些东西?
  17. 电子学习用品商城小程序开发功能需求分析
  18. 梦一样——十一月英语总结
  19. 由浅入深玩转华为WLAN—23 Traffic-filter(ACL)在WLAN无线场景的应用
  20. Java 函数的参数引用,通过传入参数带回函数运行结果

热门文章

  1. 当天期货开户次日就可以交易
  2. 惠普HP Prime可编程计算器之添加坐标表及写入与读取
  3. 如何精确的设置行与行之间的间距
  4. mingw linux c 开发环境,MSYS2和MINGW的开发环境安装
  5. 未来的计算机事业作文600字,模块计算机考试的感想
  6. 通过DataEase行列权限设置实现数据权限管控
  7. 3dmax 2016 2015 2014 2013 vary渲染 视频教程 从入门到精通
  8. 将 Jar 包打成一个 Jar 包
  9. 【文末福利】用Python画了一幅《海上生明月》的画
  10. 【电商】电商后台设计—促销模块(上)