生信初学者使用GeneMark_ES脚本经验

最新推荐文章于 2024-01-04 00:01:49 发布
最新推荐文章于 2024-01-04 00:01:49 发布
阅读量1.4k 收藏 2
点赞数 1
文章标签： 大数据 perl
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/sunwb83/article/details/125794845
版权

GeneMark-ES perl脚本使用

  1. ES需要perl module可以在install文件中查询。

  2. 在安装module时分自动和手动安装。自动安装指令为: 按装cpan：$sudo perl -MCPAN -e install Spiffy; 安装perl模块：$sudo cpan install DBI; 验证是否安装成功: $perl -e 'use DBI'  （没报错即成功）。手动安装，先从module网站（http://matecpan.org或者http://search.cpan.org/ ）下载module包。加压后在module目录下运行perl Makefile.PL指令检测mudole包，先后用make指令和sudo make install指令安装module。包安装好后在gm_key文件夹下运行指令cp gm_key ~/.gm_key指令拷贝key文件（要将key文件名称改为gm_key） 。genemark_ES解压目录gmes_linux文件夹下运行perl check_install.bash脚本检查ES运行的所以包和环境是否完整。

  3. 在GeneMark-E-test文件下新建ES文件夹；在ES文件夹内新建input文件夹存放要分析的数据；在新建test文件夹，在test文件夹目录下执行$ ../../../gmes_petap.pl --sequence ../input/genomic.fna --ES --fungus指令就可以成功预测目标基因组内基因。（../表示上一级文件夹，有多少上一级文件就有多少../）

  4. 注意：直接在gmes_petap.pl脚本所在文件夹目录下运行此脚本会报错code not excute in installation folder。所以必须在脚本所在文件下新建文件夹内才可以正常运行。

sunwb83
关注
  • 1
    点赞
  • 踩
  • 2
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
GeneMark-ES:真核生物编码基因预测软件
庐州月光的博客
08-07 4661
欢迎关注"生信修炼手册"!GeneMark-ES软件用于预测真核生物中的蛋白编码基因，和其他预测基因结构的软件不同，它采用的是非监督算法，可以不依赖训练集进行预测。官网如...
GeneMark进行基因组无监督从头预测
m0_51499191的博客
12-03 475
参考基因结构注释（1）：从头注释
使用MAKER进行基因注释(高级篇之GeneMark-ET模型训练）
xuzhougeng blog
11-04 5085
GeneMarkGeorgia Institute of Technology开发的一系列基因预测工具。真核生物基因组预测主要会用到GeneMark-ES/ET, 其中GeneMark-ES可用于无监督自训练，也就是只要提供一个基因组序列即可，而GeneMark-ET则是在GeneMark-ES的基础上整合了高通量的RNA-Seq转录本数据，工作流程如下 工作流程 如果...
GeneMark-ES/ET/EP的安装
Sunday1016的博客
11-06 1293
首先看一下官网的GeneMark的说明他有好几个专用版本，例如GeneMark-ES，GeneMarkS，MetaGeneMark，GeneMark.hmm with Heuristic Models，GeneMark.hmm Prokaryotic，GeneMark.hmm Eukaryotic，GeneMark。需要哪个下载那个我的样品是真核的刚开始看文章用了GeneMark.hmm，下载后发现库很少，后来发现GeneMark-ES里面也有GeneMark.hmm。
java生成海报实例源码-BRAKER:BRAKER是使用GeneMark-ES/ET和AUGUSTUS在新型真核基因组中全自动预测蛋白质编码
06-15
java生成海报实例源码刹车用户指南 BRAKER 的 Github 存储库的联系方式： Katharina J. Hoff，德国格赖夫斯瓦尔德大学，+49 3834 420 4624 Tomas Bruna，乔治亚理工学院，美国， BRAKER 的作者 凯瑟琳娜·J·霍夫, 西蒙·兰格, 亚历山大·洛姆萨泽, 托马斯布鲁纳, 马克·博罗多夫斯基 , 马里奥·斯坦克 [a]格赖夫斯瓦尔德大学，数学与计算机科学研究所，Walther-Rathenau-Str。 47, 17489 德国格赖夫斯瓦尔德 [b]格赖夫斯瓦尔德大学，微生物功能基因组学中心，Felix-Hausdorff-Str。 8, 17489 德国格赖夫斯瓦尔德 [c]佐治亚理工学院和埃默里大学 Wallace H Coulter 联合生物医学工程系，美国亚特兰大 30332 [d]计算科学与工程学院，30332 亚特兰大，美国 [e]莫斯科物理技术学院，莫斯科地区 141701，多尔戈普鲁德内，俄罗斯 图 1：当前 BRAKER 的作者，从左到右：Mario Stanke、Alexandre Lomsadze、Katha
funannotate安装
weixin_44573946的博客
05-05 2554
funannotate安装
生信分析服务器配置脚本
08-13
生信分析linux环境配置，安装分析软件环境，需要root权限
生信知识点-串讲-富集分析_R生信_生物信息_富集分析_R富集分析_生信分析_
10-02
使用`clusterProfiler`，我们首先需要获取差异表达基因列表，然后执行富集分析，最后绘制富集图和条形图以可视化结果。 以下是一般步骤： 1. **数据准备**：导入基因表达数据和相应的差异表达基因列表。这可能来自...
生信分析脚本
04-16
RNA数据格式转换脚本 在linux系统上使用python进行RNA数据格式转换
生信分析论文套路R语言代码
12-05
TCGA GEO数据处理，基因注释，差异分析，GO，KEGG，GSEA富集分析，肿瘤突变负荷，免疫浸润，LASSO回归，随机森林，SVM-RFE，COX回归，WGCNA网络，共识聚类分析，药物敏感性分析，干性指数，免疫浸润指数，预后模型等...
Perl-for-Bioinformatics：尝试帮助有兴趣将Perl用于生物信息学的任何人
02-05
用Perl编写的生物信息学工具 这里的大多数脚本是在我从事不同项目时编写的，我认为这对其他人将很有用，并且可以根据需要进行扩展/修改。 IO ::常规 脚本使用自定义 Perl模块。 请通过浏览目录查看安装说明。 如果要安装lncRNApipe Pipeline，则会自动安装IO::Routine模块。 需要Bio::SeqIO模块已安装且可用。 nc lncRNA管道 从头开始提取推定的新型lncRNA的管道，其中提供了从深度测序数据（例如：RNA-Seq）和注释数据组装而成的GTF格式的转录本列表。 转到目录以获取脚本列表。 安装lncRNApipe及其所有依赖项（Mac和
GeneMarkerHID_V1.90_Demo_Setup.EXE
08-11
用于基因分型检测结果分析 微卫星STR荧光电泳结果分析
生信常见数据格式.pdf
03-28
生信常见数据格式.pdf生信常见数据格式.pdf生信常见数据格式.pdf生信常见数据格式.pdf生信常见数据格式.pdf生信常见数据格式.pdf
生物信息学技术：如何推动生物研究和医疗创新
最新发布
禅与计算机程序设计艺术
01-04 484
1.背景介绍 生物信息学技术(Bioinformatics)是一门结合生物学、计算机科学和数学的跨学科学科，其主要目标是研究生物数据的存储、管理、分析和应用。随着生物科学的发展，生物信息学技术在生物研究和医疗创新中发挥着越来越重要的作用。 生物信息学技术的发展受到了生物科学、计算机科学和数学等多个领域的支持。生物科学为生物信息学提供了研究对象和背景知识，计算机科学为生物信息学提供了数据处理和分...
BGI-College生信入门——6、Perl小骆驼学习笔记（一）
bgicollege的博客
07-19 340
Perl 入门
2020.12.19丨根据差异基因ID匹配注释文本脚本
12-19 586
上周遇到一个比较麻烦的项目。物种是一种酵母菌，参考基因组是组装的，并没有像样的gtf文件，使用genemark-ES生成注释文件后，需要对差异基因进行注释。本来之前有一个脚本能够很流畅地处理这个步骤。然而，由于genemark-ES自动生成的geneID在perl脚本中存在bug，个人对perl又不是很熟悉，因此重新写了一个脚本，用于差异基因与注释文本的匹配注释，下面直接上脚本。 import csv genome_file = open('C:/Users/bbplayer/Downloads/ge
2018-6-25转录组学习3 参考基因组和基因注释
L_yivs的博客
06-25 4252
1. 参考基因组下载 目前常用的国际通用的人类参考基因组在NCBI，UCSC和ENSEMBL的版本和对应情况如下： NCBI UCSC ENSEMBL GRCH36 hg18 ENSEMBL release_52 GRCH37 hg19 ENSEMBL release_59/61/64/68/69/75 GRCH38 hg38 ...
各种组织的marker gene
zengwanqin的博客
04-21 9901
最近做单细胞分群，总结一下我自己查文献查到的各种组织有哪些细胞，这些细胞的marker gene 是什么。 一、bladder 膀胱 bladder genes_EC = c("Pecam1", "CD31") #endothelial cells genes_IC = c("Cd14")#immune cells genes_BEC = c("Epcam","Upk3a","Upk1b")# bladder epithelial cells genes_UIC = c("Grhl3")#umbrella
一文读懂参考基因组和基因组注释+最全下载方法
热门推荐
Baimoc
06-02 5万+
文章目录一、什么是参考基因组和基因组注释？二、参考基因组版本命名1、常用人参考基因组对应表2、常用小鼠参考基因组对应表三、下载1、NCBI2、Ensemble3、GENCODE4、UCSC5、iGenomes四、其他参考基因组信息 一、什么是参考基因组和基因组注释？ 先来理一理参考基因组，基因组注释文件，测试数据间的关系。 自从 1990 启动的家喻户晓的人类基因组计划开始，全世界的科学家竭尽全力破译了第一个完整的人类基因组，从那时开始人类拿到了一本只有 ATCG 四个碱基书写的天书。后续人们逐步完善了基因
abricate生信分析软件使用方式以及输出代码怎么写
06-09
Fabricate是一个用于生物信息学数据分析的Python软件包。它主要用于基因组学和转录组学数据分析。下面是Fabricate的使用方式和代码输出方式： 使用方式： 1. 安装Fabricate 您可以使用pip安装Fabricate。在终端中输入以下命令即可： ``` pip install fabricate ``` 2. 编写分析脚本 Fabricate的使用方式与常规的Python脚本类似。您需要编写Python脚本来完成您的生物信息学数据分析。 例如，以下是一个简单的Fabricate分析脚本： ```python from fabricate import * def build(): # 运行您的数据分析代码 run('python analysis.py') def clean(): # 清理分析结果 run('rm -rf results') def all(): # 执行所有任务 clean() build() if __name__ == '__main__': main() ``` 在上面的脚本中，`build()`函数运行您的数据分析代码，`clean()`函数清理分析结果，`all()`函数则是执行所有任务。`main()`函数将根据命令行参数来执行不同的任务。 3. 运行分析脚本 在终端中运行以下命令来执行您的分析脚本： ``` python your_script.py all ``` 其中`your_script.py`是您的Fabricate分析脚本的文件名，`all`是要执行的任务。 输出代码： Fabricate的输出代码方式与常规的Python脚本类似。您可以使用print语句将结果输出到终端或将结果写入文件。 例如，以下是一个将分析结果写入文件的示例代码： ```python from fabricate import * def build(): # 运行您的数据分析代码 run('python analysis.py > results.txt') def clean(): # 清理分析结果 run('rm -rf results') def all(): # 执行所有任务 clean() build() if __name__ == '__main__': main() ``` 在上面的脚本中，`build()`函数将分析结果写入`results.txt`文件。在您的数据分析代码中，您可以将结果输出到标准输出，例如： ```python print('分析结果：', result) ``` 在运行分析脚本时，您可以将结果重定向到文件中，例如： ``` python your_script.py all > results.txt ``` 其中`your_script.py`是您的Fabricate分析脚本的文件名，`all`是要执行的任务，`results.txt`是要将结果写入的文件名。
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值