2AMF MiSeq Bioinformatics

salloc --account=microbiome -t 0-05:00

cd /gscratch/grandol1/2AMF/rawdata

unpigz --to-stdout /project/gtl/data/raw/2AMF/rawdata/2AMF_R1.fq.gz | split -l 1000000 -d --suffix-length=3 --additional-suffix=.fastq - 2AMF_R1_ ;
unpigz --to-stdout /project/gtl/data/raw/2AMF/rawdata/2AMF_R2.fq.gz | split -l 1000000 -d --suffix-length=3 --additional-suffix=.fastq - 2AMF_R2_

//project/gtl/data/raw/2AMF/rawdata/run_parse_count_onSplitInput.pl

above step running 12-30-22

cd /project/gtl/data/raw/2AMF/rawdata

./run_splitFastq_fwd.sh

./run_splitFastq_rev.sh

./run_aggregate.sh

cd /project/gtl/data/raw/2AMF/tfmergedreads

./run_slurm_mergereads.pl

cd /project/gtl/data/raw/2AMF/otu

./run_slurm_mkotu.pl

Â