gffread: 从基因组提取fa
使用命令gffread -w test.fa -g GRCh38.primary_assembly.genome.fa test.gff3即可实现。
输入文件test.gff3如下所示：
记住：第三类可以写：transcript, protein, cds。不能写其他的字符，不然提取不了。
得到的输出文件test.fa如下所示：
本文来自博客园，作者：橙子牛奶糖（陈文燕），转载请注明原文链接：https://www.cnblogs.com/chenwenyan/p/16442529.html
【推荐】国内首个AI IDE，深度理解中文开发场景，立即下载体验Trae
【推荐】编程新体验，更懂你的AI，立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包，你的智能百科全书，全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell：AI 加持，快人一步
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型，支持深度思考和联网搜索！
· 基于 Docker 搭建 FRP 内网穿透开源项目（很简单哒）
· ollama系列1：轻松3步本地部署deepseek，普通电脑可用
· 按钮权限的设计及实现
· Apache Tomcat RCE漏洞复现（CVE-2025-24813）
2021-07-04 本周最新文献速递20210704
2020-07-04 irlba报错“Error in (function (A, nv = 5, nu = nv, maxit = 1000, work = nv + 7, reorth = TRUE, : BLAS/LAPACK 例行程序'DLASCL'产生了错误-4信息”解决方案
2016-07-04 什么是基因间隔区和内含子？