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

array=( A J An )

for i in "${array[@]}"
do
metadir=/data/mramseylab/metadata/plaque_meta_"$i".tsv
procdir=/data/mramseylab/raw_reads/trip"$i"/core-metric-results-"$i"

qiime diversity beta-correlation \
  --i-distance-matrix $procdir/bray_curtis_distance_matrix.qza \
  --m-metadata-file  $metadir \
  --m-metadata-column timepoint \
  --p-method pearson \
  --p-intersect-ids TRUE \ 
  --o-metadata-distance-matrix $procdir/time-matrix-"$i".qza \
  --o-mantel-scatter-visualization $procdir/mantel-scatter-"$i".qzv

  done