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

clsdir=/data/mramseylab/classifiers

module load QIIME2/2019.7

qiime tools import \
  --input-path $clsdir/silva-taxonomy-AL/taxonomy-mod-AL.tsv \
  --type 'FeatureData[Taxonomy]' \
  --input-format TSVTaxonomyFormat \
  --output-path $clsdir/silva-taxonomy-AL/silva-mod-taxonomy-AL.qza
  
qiime metadata tabulate \
  --m-input-file $clsdir/silva-taxonomy-AL/silva-mod-taxonomy-AL.qza \
  --o-visualization $clsdir/silva-taxonomy-AL/silva-mod-taxonomy-AL.qzv