关于基因结构注释,我的一些折腾和思考

我在2018年,写过一篇《如何对基因组序列进行注释》,简书上有4万多人阅读,CSDN上有8万多人阅读,说明确实有不少人有相关的需求。我自己也琢磨了一段时间,刚好过年回家,简单的梳理一下。

最初的时候,我的思路就是EvidenceModeler(后面简称EVM)的思路,训练不同从头预测软件,例如AUGUSTUS,组装RNA-seq数据,比对同源蛋白,最后交给EVM整合。这个思路挺好的,唯一的问题就是要要自己的动手的地方比较多,那个时候我自己写了一套流程做了串联,但是代码质量很差,所以我在找寻其他工具。

我找到一些替代工具,比如说MAKER, BRAKER, Mikado。这些工具各有不同,MAKER可以整合不同来源的证据, BRAKER会先利用已有同源蛋白或者RNA-seq训练AUGUSTUS模型然后预测结果,Mikado则是直接使用转录组作为唯一证据,最后我选择了MAKER。因为BRAKER,是一个perl串联的脚本,没有完善重跑工具,也就是失败了,就可能得从头开始。mikado只用转录本,那么就可能有一些基因会因为表达量低别遗漏,或者结构不完整。

使用MAKER也有问题,就是官方提供SNAP的训练流程,但是AUGUSTUS的模型训练方法则没有,这很让我头疼。一开始,我想的是找现有的解决方案,发现Transdecoder有一个工具,似乎可以直接把预测的ORF结果导出为AUGUSTUS的输入,但是AUGUSTUS的教程说并不需要那么多,大概1000左右就饱和了。由于我水平有限,一下子卡壳了,不得不去找其他的方案。于是我就想着能不能用BRAKER的输出的模型作为AUGUSTUS的输

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/677307.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

微服务学习 | Spring Cloud 中使用 Sentinel 实现服务限流

前些天发现了一个巨牛的人工智能学习网站,通俗易懂,风趣幽默,忍不住分享一下给大家。点击跳转到网站https://www.captainbed.cn/kitie。 目录 前言 通过代码实现限流 定义资源 通过代码定义资源 通过注解方式定义资源 定义限流规则 通过…

nginx 配置404 页面

自己定义404 页面,通过 nginx 配置一下,重新定向 1. http 里面添加配置 proxy_intercept_errors on; fastcgi_intercept_errors on; 2. server 中添加配置 error_page 500 502 503 504 /404.html;error_page 404 /404.html;location /404.html…

cximage在vs2013下使用方法

1.下载源码 Cximage源码官网 CxImage download | SourceForge.net 下载最新版本 702版本 Download cximage702_full.7z (CxImage) 2.编译 vs2013打开CxImageFull_vc10.sln 这个源码版本是vc10的版本,所以vs2013会自动更新项目 因为cximage需要在后面的项目中使…

Nginx实战:2-日志配置

目录 前言 一、访问日志 1.字段配置 2.日志配置 3.默认配置 二、错误日志 前言 Nginx是一款高性能的HTTP和反向代理服务器,同时也是一个IMAP/POP3/SMTP服务器。在Nginx的日常使用中,日志记录是非常重要的一部分,它可以帮助我们监控服务…

2月7日《CS2》终于放大招,玩家激情再次被点燃

2024.2.7号,也就是昨天,V社终于放了大招,对CS2做了高达5个多G的大更新,这次更新内容还是比较多的,说几个比较有意思的点吧。 1、新武器箱:千瓦武器箱! 全新的武器箱千瓦箱,能开出全…

java 回答问题

1. How do you create a variable with the numeric value 5? int x 5; 2. The value of a string variable can be surrounded by single quotes. False 3. Which method can be used to return a string in upper case letters? toUpperCase()

AutoSAR(基础入门篇)5.4-Autosar_汽车ECU项目的开发流程

目录 前言 一、团队介绍 二、AutoSAR开发流程 前言 项目的基础入门篇终于是要讲完了(后面可能还有一些补充内容)。大家耐着性子能看到这里也是真心不容易,AutoSAR很多内容都是理论上的东 西,还需要多加实践。所以如果你能认真看完前面的内容,并加以理解,那么我相信再看…

Netty应用(一) 之 NIO概念 基本编程

目录 第一章 概念引入 1.分布式概念引入 第二章 Netty基础 - NIO 1.引言 1.1 什么是Netty? 1.2 为什么要学习Netty? 2.NIO编程 2.1 传统网络通信中开发方式及问题(BIO) 2.1.1 多线程版网络编程 2.1.2 线程池版的网络编程…

搭建yum仓库服务器

安装 1.安装linux 1.1安装依赖 yum -y install gcc zlib zlib-devel pcre-devel openssl openssl-devel 1.2下载 cd /opt/nginx wget http://nginx.org/download/nginx-1.25.3.tar.gz 1.3解压 tar -xvf nginx-1.25.3.tar.gz 1.4配置 cd nginx-1.25.3 ./configure --pre…

FTP 文件传送协议

目录 1 文件传送协议 FTP 1.1 FTP 的基本工作原理 FTP 特点 主进程的工作步骤 两个连接 两个不同的端口号 NFS 采用另一种思路 1.2 简单文件传送协议 TFTP TFTP 的主要特点 TFTP 的工作(很像停止等待协议) 1 文件传送协议 FTP 文件传送协议 …

Redisson分布式锁 原理 + 运用 记录

Redisson 分布式锁 简单入门 pom <dependency><groupId>org.redisson</groupId><artifactId>redisson</artifactId><version>3.13.6</version></dependency>配置类 package com.hmdp.config;import org.redisson.Redisson;…

NLP_引入注意力机制

文章目录 点积注意力创建两个张量x1和x2计算张量点积&#xff0c; 得到原始权重对原始权重进行归一化求出注意力分布的加权和 缩放点积注意力编码器-解码器注意力定义Attention类重构Decoder类重构Seq2Seq类可视化注意力权重 注意力机制中的 Q、K、V自注意力多头自注意力注意力…

CSS 2D转换 3D动画 3D转换

目录 2D转换(transform): 移动translate: 旋转rotate: 缩放scale&#xff1a; CSS3动画&#xff08;transform&#xff09;&#xff1a; 动画常用的属性&#xff1a; 将长图片利用盒子实现动画的效果&#xff1a; 3D转换&#xff1a; 透视perspective&#xff1a; 旋转r…

动态内存经典笔试题分析

1.代码1 void GetMemory(char *p) { p (char *)malloc(100); } void Test(void) { char *str NULL; GetMemory(str); strcpy(str, "hello world"); printf(str); } int main&#xff08;&#xff09; { Test&#xff08;&#xff09;&#xff1b; return 0&#x…

WWW 2024 | 时间序列(Time Series)和时空数据(Spatial-Temporal)论文总结

WWW 2024已经放榜&#xff0c;本次会议共提交了2008篇文章&#xff0c;research tracks共录用约400多篇论文&#xff0c;录用率为20.2%。本次会议将于2024年5月13日-17日在新加坡举办。 本文总结了WWW 2024有关时间序列&#xff08;Time Series&#xff09;和时空数据&#xf…

使用rem单位制,实现页面适应性窗口元素

目录 REM单位&#xff1a; 媒体查询&#xff1a; 引入资源&#xff1a; REM适配方案&#xff1a; 适配方案一&#xff1a; rem媒体查询less技术 适配方案二&#xff1a; REM单位&#xff1a; rem(rootem)是一个相对单位&#xff0c;类似于em,em是父元素字体大小。 不同的…

elasticsearch下载及可视化工具下载使用

elasticsearch下载及配置、启动 一、下载 Download Elasticsearch | Elastic 二、启动 双击bat即可。 出现如下说明启动成功&#xff1a; 访问测试&#xff1a; 三、注意 &#xff08;1&#xff09;因为es启动默认端口是&#xff1a;9200,所以需要检查此端口是否被占用。…

Linux系统中HTTP代理的常见问题及解决方案

亲爱的Linux用户们&#xff0c;是不是有时候觉得HTTP代理就像是一个魔法盒子&#xff0c;让你在数字世界中自由穿梭&#xff1f;但是&#xff0c;就像所有的魔法物品一样&#xff0c;它也会偶尔出点小状况。今天&#xff0c;我们就来一起探讨一下Linux系统中HTTP代理的常见问题…

DC-9靶机渗透详细流程

信息收集&#xff1a; 1.存活扫描&#xff1a; arp-scan -I eth0 -l 发现靶机ip&#xff1a;192.168.10.132 └─# arp-scan -I eth0 -l 192.168.10.1 00:50:56:c0:00:08 (Unknown) 192.168.10.2 00:50:56:e5:b1:08 (Unknown) 192.168.10.132 //靶机 00:0c…

【漏洞复现】狮子鱼CMS某SQL注入漏洞

Nx01 产品简介 狮子鱼CMS&#xff08;Content Management System&#xff09;是一种网站管理系统&#xff0c;它旨在帮助用户更轻松地创建和管理网站。该系统拥有用户友好的界面和丰富的功能&#xff0c;包括页面管理、博客、新闻、产品展示等。通过简单直观的管理界面&#xf…