BLAST本地检索实践 - SpiralKernel/bioinfo-wiki-migrated GitHub Wiki

BLAST本地检索实践

更新:2026-10-13;状态:可用;验证范围:文档核对;操作需按所用环境验证。

用本地数据库理解查询、数据库构建与结果字段,并将正式结论建立在覆盖率和统计解释上。

适用范围与准备

先读BLAST原理与种子延伸和FASTA与FASTQ。这里介绍一般序列检索的软件接口;需要已安装 BLAST+、查询 FASTA 和来源明确的参考 FASTA,工具下载与完整手册在 NCBI BLAST Help。

本地接口示例

blastn -version
mkdir -p results blastdb
makeblastdb -in reference.fa -dbtype nucl -out blastdb/reference
blastn -query query.fa -db blastdb/reference \
  -outfmt '6 qseqid sseqid pident length qstart qend sstart send evalue bitscore qcovs' \
  -out results/hits.tsv

命令按官方接口核对,尚未在本机安装 BLAST+ 运行。reference.fa 和 query.fa 是需替换的输入,不是已经提供的数据。若查询和数据库是蛋白,应选择对应程序与 -dbtype prot;详见手册的程序选择表。

输出与检查点

格式 6 的该示例没有标题行,列顺序由 outfmt 定义。保存列名说明、程序版本、查询校验值和数据库构建日期。检查覆盖率、方向及多个局部命中;空输出也可能是阈值、程序或数据库不匹配,需要查看日志与输入。

最佳命中不是自动的功能注释或一对一同源证明。统计解释见比对相似性同源性与覆盖率,基因关系见基因标识符与同源基因。

深入入口:BLAST+ 应用手册、官方程序与文献索引。