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

echo "emails sent"

module load QIIME2/2024.2

echo "module loaded"

procdir=/data/mramseylab/16S/proc_reads/2024_AL_SUPP_6000/filteredsamples_2024-2/
clsdir=/data/mramseylab/16S/classifiers/homd-classifier/exported-taxa/
metadata=/data/mramseylab/16S/metadata/meta_BSFS.tsv
visdir=/data/mramseylab/16S/visualizations/AL_2024/

if qiime taxa barplot \
  --i-table $procdir/rarefied_table.qza \
  --i-taxonomy $clsdir/modifiedHOMDforbars.qza \
  --m-metadata-file $metadata \
  --o-visualization $procdir/HOMD-filter-bar.qzv; then
    echo "barplot made!"
fi
echo " Script DONE!!!!"