I am trying to import Silva_132 DB but I cant find the right command for that. I found an answer that u gave before:
So, I downloaded the database and I am using this command:
qiime tools import --input-path /home/qiime2/Desktop/Datasets/Silva_132_release --output-path silva_132_db.qza --type FeatureData[Sequences]
But I am getting this error: Sequences is not a variant of FeatureData.field[‘type’]…
What am I doing wrong? I want to use this DB as BLAST+ taxonomy classifier.
After I will import the db ( according to your suggestion) where can I find a good example for a feature-classifier command? ( In qiime 1… you gave us few examples for each script, this was very useful for “slow” people like myself…)
Did you follow the link to another post at the bottom of the forum topic you listed? Reposting here:
That user shows how to import the necessary files and use them for the blast classifier, answering both of your questions. (the errors they report are related to the custom database they made; the commands they report there are technically correct).
My problem is that the database file that I downloaded composed of many files and not only the txt file as he mentioned.
I downloaded the file from Silva according to your link, unzip it and now I need to import it.
Shouldn’t I import all the files in it? according to Yer Lor question, he imports only one txt file.
Correct; collapse N times where N = the number of taxonomic levels you want to report.
You will need to export your feature table and the sequences and merge these outside of QIIME 2. It sounds like maybe you are planning to upload these to R for making plots? In which case phyloseq or whatever R package you are using should have methods for merging these data on the sample IDs.
The 7th level is species right? and as far as I understand one of the advantages of using DADA2 over qiime 1 procedure ( pick OTU, rep_seq , etc’ ) is that each Sequence Variants is now 100% “OTU” , so we don’t use anymore 97% clustering, right?
So in this case, can I get the sequences only for the 7th level phylogeny without export to R?
( Since I am going to analyze many samples from many runs, I am trying to build a quick and convenient pipeline so the less transformations/exports is better).
Thanks again for your answer, I am not taking it as granted at all !