kraken建库错误

kraken 出现错误

运行命令kraken-build --standard --threads 24 --db MY_KRAKEN_DATABASE，返回错误：

Found jellyfish v1.1.12
--2019-04-26 16:43:23--  ftp://ftp.ncbi.nlm.nih.gov/pub/taxonomy/accession2taxid/nucl_est.accession2taxid.gz
           => “nucl_est.accession2taxid.gz”
正在解析主机 ftp.ncbi.nlm.nih.gov (ftp.ncbi.nlm.nih.gov)... 130.14.250.13, 2607:f220:41e:250::11
正在连接 ftp.ncbi.nlm.nih.gov (ftp.ncbi.nlm.nih.gov)|130.14.250.13|:21... 已连接。
正在以 anonymous 登录 ... 登录成功！
==> SYST ... 完成。   ==> PWD ... 完成。
==> TYPE I ... 完成。 ==> CWD (1) /pub/taxonomy/accession2taxid ... 完成。
==> SIZE nucl_est.accession2taxid.gz ... 完成。
==> PASV ... 完成。   ==> RETR nucl_est.accession2taxid.gz ...
文件 “nucl_est.accession2taxid.gz” 不存在。

查看NCBI ftp数据库后发现，est 和gss文件ncbi不再提供，因此需要修改kraken相关配置文件：

搜索并修改download_taxonomy.sh文件，注释掉第22、24行

Ref：https://github.com/DerrickWood/kraken/commit/69a682775754676a44139a1826a63f0dd9674f4c