Anyone know what I did wrong with my archae analysis, using the SILVA database for archaea?
Because when validating the sequences obtained in qiime these turn out to be bacteria and not archaea. The validation was carried out in NCBI BLAST
qiime tools import --type ‘SampleData[SequencesWithQuality]’ --input-path secuencias --input-format CasavaOneEightSingleLanePerSampleDirFmt --output-path demux-single-end-todas16s.qza
qiime tools validate demux-single-end-todas16s.qza
qiime demux summarize --i-data demux-single-end-todas16s.qza --o-visualization demux-single-end-todas16s.qzv
qiime dada2 denoise-single –i- demultiplexed-seqs demux-single-end-todas16s.qza –p-trunc-len 440 –p-n-threads 0 –o-table table-dada2.qza –o-representative-sequences rep-seq-dada2-todas16s.qza –o-denoising-stats stats-dada2.qza
(qiime2-2019.4) [email protected]:~/16s$ qiime tools import --type ‘FeatureData[Sequence]’ --input-path arqueas.fna --output-path reference_arqueas.qza
(qiime2-2019.4) [email protected]:~/16s$ qiime feature-classifier extract-reads --i-sequences reference_arqueas.qza --p-f-primer CCTACGGGNGGCWGCAG --p-r-primer GACTACHVGGGTATCTAATCC --p-trunc-len 440 --o-reads ref-seqs-arqueas.qza
(qiime2-2019.4) [email protected]:~/16s$ qiime feature-classifier fit-classifier-naive-bayes --i-reference-reads ref-seqs-arqueas.qza --i-reference-taxonomy archaea_ncbi_rrna_tax.qza --o-classifier classifier_arqueas.qza
(qiime2-2019.4) [email protected]:~/16s$ qiime feature-classifier classify-sklearn --i-classifier classifier_arqueas.qza --i-reads rep-seq-dada2-todas16s.qza --o-classification taxonomy_arqueas.qza
(qiime2-2019.4) [email protected]:~/16s$ qiime taxa barplot --i-table table-dada2-todas16s.qza --i-taxonomy taxonomy_arqueas.qza --m-metadata-file sample-metadata.tsv --o-visualization taxa-bar-plots_arqueas.qzv