Building full-length rRNA database using RESCRIPt

Hi @SoilRotifer ,

Thanks for your further explanation. I managed to download all bacterial genomes from NCBI RefSeq during the weekend, and have extracted more than 20000 full length rRNA sequences from those genomes. I agree with you that the actual hard work is to assign taxonomy to those sequences. Taxonkit seems to be really handy (thank you so much for recommending it!). My current plan is to create a separate tsv file that contains the taxonomy info of all those sequences, and then import them to qiime2, as suggested in this post.

Creating a custom reference database - User Support - QIIME 2 Forum

Thanks again for your help.

Tongda

1 Like