构建qiime2 blast genbank数据库

anneng

https://bmcgenomdata.biomedcentral.com/articles/10.1186/s12863-022-01067-5

anneng

https://forum.qiime2.org/t/processing-filtering-and-evaluating-the-silva-database-and-other-reference-sequence-data-with-rescript/15494

qiime2有个工具可以直接做
https://github.com/bokulich-lab/RESCRIPt#using-rescript

anneng

rescript 验证
conda create -y -n rescript

conda activate rescript

conda install -c conda-forge -c bioconda -c qiime2 -c https://packages.qiime2.org/qiime2/2023.5/tested/ -c defaults qiime2 q2cli q2templates q2-types q2-longitudinal q2-feature-classifier 'q2-types-genomics>2023.2' "pandas>=0.25.3" xmltodict ncbi-datasets-pylib

pip install git+https://ghproxy.com//https://github.com/bokulich-lab/RESCRIPt.git

创建COI库
qiime rescript get-ncbi-data --p-query '(cytochrome c oxidase subunit 1[Title] OR cytochrome c oxidase subunit I[Title] OR cytochrome oxidase subunit 1[Title] OR cytochrome oxidase subunit I[Title] OR COX1[Title] OR CO1[Title] OR COI[Title]) AND BARCODE[KYWD])' --verbose --output-dir NCBIdata_BOLDonly2 --o-sequences coi_genbank.qza --o-taxonomy coi_tax.qza

创建16s库