#!/bin/bash
#SBATCH -t 4:00:00
#SBATCH --nodes=1 --ntasks-per-node=1
#SBATCH --mem=24g
#SBATCH --export=NONE
#SBATCH --mail-user=alex_labossiere@uri.edu
#SBATCH --mail-type=BEGIN,END,FAIL


module load QIIME2/2019.7

rawdir=/data/mramseylab/proc_reads/fil-nocontrolv1v2/filtered_table_NoControls.qza
clsdir=/data/mramseylab/classifiers/silva-taxonomy-AL/silva-taxonomy-AL-v1v2.qza
metadata=/data/mramseylab/metadata/plaque_meta_3.5.tsv
visdir=/data/mramseylab/visualizations/AL_2022

qiime taxa barplot \
  --i-table $rawdir \
  --i-taxonomy $clsdir \
  --m-metadata-file $metadata \
  --o-visualization $visdir/taxa-bar-plots-v1v2.qzv

# takes the output file from classifier and uses it as imput, uses metadata from user to place taxa within samples, visualization shows % abundance of organisms
# within given samples. Can use dropdown from view.qiime2.org to seperate samples based on subject and days, you can now see how the difference in organisms %
# abundaces changes over time of the expirment
# recall this is what we need but there is bias within our samples, we need to put this data through a series of filters in order to take out
# any remainder "noise" or bias within samples (this can be bias from gDNA found in water, saliva, mucin, etc.)