I have a project using CO1 BOLD database to detect invertbrates.
I used DADA2 to create table, rep-seq and stats.
I used feature-classifier classify-consensus-blast to create a taxa table
I use barplot to create taxa table
My question is: How can I get the taxa ID and rep seq for each taxon in each sample?
I did use this to create a feature table which contains ID, sequence, taxon, confidence. However, I also need to add the sequence ID i.e. 1174152 in the following example. How can I extract the sequence ID in the BOLD database along with the sequences.
hmm... this is a feature ID in a separate file (the reference database), so there is not really a good way to link this to the feature table and taxonomy that you currently have (as the feature IDs are distinct). It sounds like you want to display the ID of the closest hit in the reference for each query. For this, you could (outside of QIIME 2): run blast, extract the query ID and target ID columns, then merge this with the other files using qiime metadata tabulate.
if confidence information is present in the taxonomy file it will be included in the merged visualization output by qiime metadata tabulate
I created the taxa table using barplot pluggin. I found that each taxon may associated with multiple sequences with different confidence. Is there anyway that I can associate the taxon per sample to the appropriate sequence?
I created the taxa table using barplot pluggin. I found that each taxon may associated with multiple sequences with different confidence. Is there anyway that I can associate the taxon per sample to the appropriate sequence?
If the barplot can't associate the taxon with the rep-sequence, can you recommend a taxonomy plugin can do this?
My goal is to associate the rep-sequence with taxon per sample.
Yes see my advice above — metadata tabulate is what you want. The tutorial at the link above describes exactly how to do this. If you are expecting something different, I suggest writing out a small example here to describe the format that you are looking for.
Sorry to bother you again. Here are the script that I used to create barplot and taxa associate with the rep seq. However, I need per sample taxa table created in the barplot and add the rep-seq to the taxa table.
I am new to the Qiime. I would appreciate you more help. Thanks, Jin