介绍一种从基因组中提取片段的方法,很简单!
流程
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
| # 安装blast wget https://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/LATEST/ncbi-blast-2.11.0+-x64-linux.tar.gz # 解压 tar -zxvf ncbi-blast-2.11.0+-x64-linux.tar.gz # 解压文件 echo 'PATH=$PATH:~/ncbi-blast-2.11.0+/' >> ~/.bashrc # 将路径加入到环境变量 source ~/.bashrc # 重新编译路径 # 还要安装一个ncbi-blast sudo apt install ncbi-blast+-legacy
# 建库 formatdb -i scaffolds.fasta -p F -o F # 比对 blastall -i query.fasta -d scaffolds.fasta -o G19.blast.out -p blastn # -i:比对的序列文件 # -d:参考基因组库 # -o:输出文件 # -p:比对方式
|
参考资料:
Blast - 简书 (jianshu.com)
[Linux 安装Blast本地版 - 知乎 (zhihu.com)](