Bandage 序列提取实战：3 种方法从组装图中直接获取目标序列

Bandage 序列提取实战3 种方法从组装图中直接获取目标序列【免费下载链接】Bandagea Bioinformatics Application for Navigating De novo Assembly Graphs Easily项目地址: https://gitcode.com/gh_mirrors/ba/Bandage在基因组组装完成后如何快速从组装图中取出目标序列Bandage 序列提取功能就是为此设计的。Bandage 是一款专为生物信息学设计的开源图形化工具全称是 a Bioinformatics Application for Navigating De novo Assembly Graphs Easily它能将 Velvet、SPAdes、MEGAHIT 等组装软件生成的 de novo 组装图可视化让普通研究者也能直观地浏览 contig 之间的连接关系并直接从图中提取序列数据。本文将用 3 种实战方法教你如何从组装图中快速获取目标序列从新手到命令行高手一次讲清。一、方法一图形界面选中节点一键导出 FASTA 序列 ✅这是最直观、最适合新手的 Bandage 序列提取方式全程鼠标操作不需要写任何命令。操作步骤打开 Bandage加载组装图文件支持.fastg、.gfa、.LastGraph、.csv等多种格式。在图上用鼠标框选或点选你感兴趣的节点每个节点代表一条 contig。点击菜单Copy selected node sequences to clipboard序列即复制到剪贴板或点击Save selected node sequences to FASTA弹出保存对话框直接导出为标准的 FASTA 文件默认文件名为selected_sequences.fasta。这套逻辑在源码中实现得非常简洁清晰核心导出函数是saveSelectedSequencesToFile()它遍历所有选中的节点逐个调用getFasta()写入文件mainwindow.cpp。小技巧如果只想要某一条序列选中单个节点即可想要多条就按住 Ctrl 多选或框选一片区域一次导出全部。二、方法二Specify exact path 精确指定路径提取完整目标序列 当你的目标序列横跨多个节点、需要按特定顺序拼接时光靠选中就不够了。Bandage 提供了Specify exact path工具让你手动输入节点路径字符串从而精准定义序列提取范围。最快操作方法在菜单中找到Specify exact path for copy/save打开路径指定对话框按方向节点ID的格式输入路径例如5 -7 9表示先正向经过节点 5再反向经过节点 7最后正向经过节点 9输入的同时图形窗口会实时高亮显示这条路径方便你确认无误对应实现见 pathspecifydialog.cpp确认路径有效后即可Copy path sequence to clipboard或Save path sequence to FASTA得到这条完整路径拼接后的序列。这个功能特别适合从环状质粒、重复区等复杂结构中提取完整序列。路径序列的生成与保存逻辑位于 mainwindow.cpp 中路径对象会按顺序把每个节点的序列拼接起来输出为一条完整 FASTA 记录。三、方法三命令行 querypaths批量 BLAST 提取目标序列 当查询序列有成百上千条、或者你想把搜索提取自动化时图形界面就力不从心了。此时可以使用 Bandage 内置的命令行子命令querypaths它先用 BLAST 在组装图中搜索每条查询序列再自动提取命中的路径序列非常适合批量处理。最快配置方法Bandage querypaths graph文件 查询FASTA 输出前缀 --pathfasta --hitsfasta命令完成后会生成 3 类结果文件前缀.tsv查询、路径、长度、覆盖度、平均一致性等详细统计表格前缀_paths.fasta每条查询对应的完整路径序列多 FASTA 格式前缀_hits.fasta路径中所有 BLAST 命中的节点序列。其中--pathfasta与--hitsfasta两个参数控制 FASTA 文件的输出完整的参数说明与实现位于 querypaths.cpp。相关命令还包含info、image、reduce等子命令全部集中在 command_line 目录下方便脚本化调用。实用提示querypaths 的输出文件不允许覆盖已存在的文件重复运行时需先清理旧结果避免报错。总结三种方法怎么选场景推荐方法上手难度看图随手取几条序列方法一选中节点导出⭐ 极简提取跨节点的完整目标序列方法二Specify exact path⭐⭐ 简单批量搜索并提取大量序列方法三querypaths 命令行⭐⭐⭐ 进阶总而言之Bandage 把看图和取序列打通了图形界面适合交互式探索与单条序列提取命令行则适合大规模批处理与流程集成。无论你是刚接触组装的生信新手还是需要批量处理数据的资深研究者掌握这 3 种 Bandage 序列提取方法都能显著提升从组装图中获取目标序列的效率。快去试试吧【免费下载链接】Bandagea Bioinformatics Application for Navigating De novo Assembly Graphs Easily项目地址: https://gitcode.com/gh_mirrors/ba/Bandage创作声明：本文部分内容由AI辅助生成（AIGC），仅供参考