/project/microbiome/data_queue/seq/LowReadII/rawdata
salloc --account=microbiome -t 0-06:00
mkdir -p /gscratch/grandol1/LowReadII/rawdata
cd /gscratch/grandol1/LowReadII/rawdata
unpigz --to-stdout /project/microbiome/data_queue/seq/LowReadII/rawdata/Low-Read-II_S1_L001_R1_001.fastq | split -l 1000000 -d --suffix-length=3 --additional-suffix=.fastq - LowReadII_R1_ ;
unpigz --to-stdout /project/microbiome/data_queue/seq/LowReadII/rawdata/Low-Read-II_S1_L001_R2_001.fastq | split -l 1000000 -d --suffix-length=3 --additional-suffix=.fastq - LowReadII_R2_
//project/microbiome/data_queue/seq/LowReadII/rawdata/run_parse_count_onSplitInput.pl
cd /project/microbiome/data_queue/seq/LowReadII/rawdata
./run_splitFastq_fwd.sh
./run_splitFastq_rev.sh
./run_aggregate.sh
cd /project/microbiome/data_queue/seq/LowReadII/tfmergedreads
./run_slurm_mergereads.pl
cd /project/microbiome/data_queue/seq/LowReadII/otu
./run_slurm_mkotu.pl
vsearch -sintax OTUFILE -db REFERENCEDATABASE -tabbedout OUTPUT -sintax_cutoff 0.8 -strand both -threads 32