#!/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"

rawdir=/data/mramseylab/16S/proc_reads/2022_AL_SUPP/2022/rep-seqs-single-2022.qza
clsdir=/data/mramseylab/16S/classifiers/homd-classifier
visdir=/data/mramseylab/16S/proc_reads/2022_AL_SUPP/2022
procdir=/data/mramseylab/16S/proc_reads/2022_AL_SUPP/2022sample-name/rarefied_table.qza
metadata=/data/mramseylab/16S/metadata/plaque_meta_2.tsv

#need to realign the representative sequences (yours) with classifer for taxonomic assignment

if
qiime feature-classifier classify-sklearn \
  --i-classifier $clsdir/HOMD-classifier.qza \
  --i-reads $rawdir \
  --o-classification $clsdir/subject-A-taxonomy-HOMD-2022.qza; then 
	echo "taxonomy aligned"
fi


if
qiime metadata tabulate \
  --m-input-file $clsdir/subject-A-taxonomy-HOMD-2022.qza \
  --o-visualization $clsdir/subject-A-taxonomy-HOMD-2022.qzv; then
	echo "metadata table made"
fi
#qzv can be used to see classification and sequence asignemnt for blast if wanted

if qiime taxa barplot \
  --i-table $procdir \
  --i-taxonomy $clsdir/subject-A-taxonomy-HOMD-2022.qza \
  --m-metadata-file $metadata \
  --o-visualization $visdir/taxabar-A-2022.qzv; then
    echo "barplot made!"
fi
echo " Script DONE!!!!"