【笔记】FASTGA无参基因组比对快速上手

发布时间:2026/9/15 9:35:47

【笔记】FASTGA无参基因组比对快速上手 FastGA是一款开源的超快速双基因组比对工具能在普通笔记本电脑上几分钟内完成10 Gbp级基因组的比对。一、下载与安装下载地址wget https://github.com/thegenemyers/FASTGA/archive/refs/tags/v1.5.zipwget https://github.com/thegenemyers/FASTGA/archive/refs/tags/v1.5.tar.gz该教程编写于2026.7如后续没法下载可以查看其github发布页https://github.com/thegenemyers/FASTGA/releases解压tar -xzf v1.5.tar.gz进入文件夹cd FASTGA-1.5编译整个项目make根据个人偏好可以将编译后的文件放入bin、添加.bashrc或者使用绝对路径使用在此不多赘述。二、使用主要功能FastGA reference.fasta query.fasta [options] # 执行比对ALNtoPAF input.1aln output.paf # 将结果 .1aln 文件转换为通用 .paf格式ALNplot input.1aln output.ps # 对结果链化并生成点图Photoshop格式其它功能ALNshow以多种格式显示选定的比对结果用于查看和调试。FAtoGDB将 FASTA 文件转换为 FastGA 内部使用的基因组数据库GDB格式。GIXmake为给定的 GDB 构建基因组索引GIX。......如果追求快速开始直接使用 FastGAFastGA reference.fasta query.fasta -1:path/to/output.1aln工具会输出 ALN 格式的文件这是一种二进制压缩文件相比 PAF 能压缩十多倍。随后可以使用 ALNplot 绘制点图ALNplot test.1aln test.ps绘制结果为 Photoshop 格式可以在 PS 中打开并编辑。最后可以使用 ALNtoPAF 将 ALN 格式转化为 PAF 格式ALNtoPAF input.1aln output.paf # 不要忘了重定向转化后的文件可以用于绘制图片或者进行分析注该工具包-paf和-psl参数都会将结果直接输出至“标准输出”屏幕需要通过重定向储存文件。注该脚本没有-h和--help选项直接运行对应脚本就能查看帮助。注在LSF集群使用该脚本时出现无法识别相对路径的情况可以将路径改为绝对路径。FastGA help 文档机翻Usage: FastGA [-vkNS] [-L:log:path] [-Tint(8)] [-Pdir($TMPDIR)] [format(-paf)][-fint(10)] [-cint(85) [-sint(1000)] [-lint(100)] [-ifloat(.7)]source1:path[precursor] (#[mask[.1bed]])*[ source2:path[precursor] (#[mask[.1bed]])* ]format -paf[mxsS]* | -psl | -1:align:path[.1aln]precursor .gix | .1gdb | fa_extn | 1_extnfa_extn (.fa|.fna|.fasta)[.gz]1_extn any valid 1-code sequence file type#: 使用此掩码对前面的基因组进行软掩码-v: 详细模式运行时输出统计信息。-L: 将日志输出到指定文件。-k: 保留任何生成的 .1gdb 和 .gix 文件。-N: 不进行软掩码。-S: 使用对称种子不推荐。-T: 使用的线程数。-P: 用于临时文件的目录。-paf: 流式输出 PAF 格式-pafx: 流式输出 PAF 格式CIGAR 字符串中使用 X-pafm: 流式输出 PAF 格式CIGAR 字符串中使用 -pafs: 流式输出 PAF 格式使用简短形式的 CS 字符串-pafS: 流式输出 PAF 格式使用长形式的 CS 字符串-psl: 流式输出 PSL 格式-1: 生成 1-code 输出到指定文件-f: 自适应种子计数截断值-c: 两个基因组中最小种子链覆盖-s: 启动新种子链的阈值-l: 最小比对长度-i: 最小比对一致性-S: 从两个基因组中种子适配体
延伸阅读

更多相关文章

2026/9/8 12:24:48

Ubuntu 26.04 LTS新特性与开发者体验优化

1. Ubuntu 26.04 LTS 深度体验报告作为长期使用Linux系统的开发者,每次Ubuntu LTS版本发布都像等待一场技术盛宴。这次26.04版本代号"Resolute Raccoon"(坚定的浣熊),确实在多个维度展现了令人惊喜的进化。我在实体机和…

2026/9/15 9:32:01

车辆动力学NMPC控制仿真避坑指南——Matlab项目完整复盘

/* MD / 富文本中的 .toc(含博客园搬家等嵌套结构);.toc-box 在侧栏,不受影响 */#content_views .toc,/* 编辑器常在目录前后插入空 p(:empty 仍占 20px),一并去掉避免顶空隙 */#content_views.markdown_views > p:empty:has(+ .toc),#content_views.markdown_views …

2026/9/15 9:32:01

基于SpringBoot的汽车空调设备管控平台设计思路与实战

/* MD / 富文本中的 .toc(含博客园搬家等嵌套结构);.toc-box 在侧栏,不受影响 */#content_views .toc,/* 编辑器常在目录前后插入空 p(:empty 仍占 20px),一并去掉避免顶空隙 */#content_views.markdown_views > p:empty:has(+ .toc),#content_views.markdown_views …

2026/9/15 9:32:01

网络汇聚分流技术演进:从流量复制到智能数据服务平台

一、引言 过去十几年,网络汇聚分流(Network Packet Broker,NPB)一直承担着网络流量“中转站”的角色。其主要职责是将来自不同链路、不同交换机的网络流量进行统一汇聚、过滤、复制,再按需分发给网络安全、业务监控和运…

2026/9/15 9:27:01

LabVIEW实现UDS安全访问SID 27:Seed-Key时序与算法可插拔设计

/* MD / 富文本中的 .toc(含博客园搬家等嵌套结构);.toc-box 在侧栏,不受影响 */#content_views .toc,/* 编辑器常在目录前后插入空 p(:empty 仍占 20px),一并去掉避免顶空隙 */#content_views.markdown_views > p:empty:has(+ .toc),#content_views.markdown_views …

2026/9/15 4:54:30

拯救者Y7000黑屏故障排查与维修实战指南

1. 项目概述:一台黑屏的拯救者Y7000,到底卡在哪一步? 联想拯救者Y7000系列笔记本,从2018年第一代搭载i5-8300H开始,到后来的i7-9750H、i7-10750H、i5-11400H,再到2023年款的R7-7840HS,它始终是学…

2026/9/15 0:01:16

AI英语单词APP开发:自适应学习算法与移动端优化实践

1. 项目概述 作为一名在移动应用开发领域摸爬滚打多年的老手,我最近完成了一个AI英语单词APP的开发项目。这个项目将传统单词记忆方法与现代AI技术相结合,打造了一款能够智能适应不同用户学习习惯的英语学习工具。 市面上大多数单词APP都存在一个通病&a…

2026/9/15 0:01:16

Flutter与OpenHarmony结合开发手语学习APP实战

1. 项目背景与核心价值作为一名同时接触过Flutter和OpenHarmony的开发者,最近我完成了一个基于Flutter for OpenHarmony的手语学习APP实战项目。这个项目最大的特点在于实现了跨平台框架与国产操作系统深度结合的创新实践——用Flutter开发的应用能完美运行在OpenHa…

2026/9/15 0:01:16

六个月成为机器人工程师:从ROS2到SLAM的实战路径

1. 六个月的紧迫感从哪来:先搞清楚你要成为哪种机器人工程师说实话,六个月的期限并不是一个宽松的时间线。市面上任何一本正经的机器人学教材都超过五百页,ROS2的官方文档可以翻到你怀疑人生,再加上ABB、KUKA这些工业机器人厂家动…

2026/9/14 11:59:31

USB Type-C PCB布局分区设计:电源、高速信号与PD协议全攻略

做硬件这行,Type-C接口算是典型的“看着简单,做起来全坑”的东西。光引脚就24个,高低速信号、电源、控制线全部塞在一个小小的连接器里,如果PCB布局不做规划,打样回来基本就是“插上没反应”、“高速掉线”、“静电一打…

2026/9/14 13:53:59

系统编程学习原型如何补齐稳定性边界

系统编程学习原型如何补齐稳定性边界预算有限时&#xff0c;我先优化明显多余的复制&#xff0c;而不是猜测性地换容器。用借用传递只读数据通常就能减少分配&#xff1a; fn parse(line: &str) -> Result<Item, Error> { /* ... */ }用基准确认热点确实在分配&am…

2026/9/14 11:22:57

雨花区哪家财务公司代理记账比较好?

在雨花区&#xff0c;企业处理财税事务常常面临诸多挑战&#xff0c;选择一家靠谱的财务公司至关重要。湖南巨勤财务管理咨询有限公司就是本地正规实体财税服务机构&#xff0c;深耕本地工商财税行业多年&#xff0c;熟悉当地工商局、税务局最新政策与申报流程。主营公司注册、…

还想了解更多?直接咨询顾问

免费诊断 + 免费方案 + 透明报价。

全国咨询热线400-8866-253
免费获取方案
咨询二维码