在群体遗传学或基因家族分析中,我们常常需要从多个基因组中提取某个基因及其上下游调控区的同源序列。然而,直接根据参考基因组坐标截取往往不可靠(存在结构变异),而 BLAST 流程又容易在集群共享存储上遇到各种“Bus error”。本文记录了一套稳健的批量提取流程,最终成功从 20 个玉米自交系中提取了目标区域,序列长度一致。一、软件安装(无 root 权限,使用 Conda)1. 安装 Minicondacd/your/software/path wget https://repo.anaconda.com/miniconda/Miniconda3-latest-Linux-x86_64.sh bash Miniconda3-latest-Linux-x86_64.sh -b -p ./miniconda3 eval"$(./miniconda3/bin/conda shell.bash hook)" conda init source~/.bashrc2. 创建环境并安装BLAST、bedtools、samtoolscondacreate -n blast_env -c bioconda blast bedtools samtools conda activate blast_env验证安装是否成功:blastn