#PBS -q bl2x220g7q
#PBS -l select=1:ncpus=5:mem=10g
#PBS -l walltime=01:00:00
#PBS -l cput=01:00:00

source ~/programms/ame/bin/activate
WORK_DIR=/mnt/storage/home/psbelokopytova/scratch/202112281307data_Polya/4c_plasmid/2021-11-22/
cd $WORK_DIR;
[ -d stat ] || mkdir cutadapt;
# Declare a string array with type
declare -a files=("P4C5" "P4C25")

for file in "${files[@]}"
do
cutadapt -a AGATCGGAAGAGCACACGTCTGAACTCCAGTCAC -A AGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGTAGATCTCGGTGGTCGCCGTATCATT -o $WORK_DIR/cutadapt/cutIL_${file}_1.fq.gz -p $WORK_DIR/cutadapt/cutIL_${file}_2.fq.gz $WORK_DIR/raw/211102_M067_V350034057_L04_${file}_1.fq.gz $WORK_DIR/raw/211102_M067_V350034057_L04_${file}_2.fq.gz -m 8;
wait;
cutadapt -g pr_plasmid1=GAATCAGGGGATAACGCAGGAAAGAACATG -y '{name}' -g pr_plasmid2=AACAGCTATGACCATTACGCCAAGCTTGGCCTGCAGGTCGACTCTAGAGGATCCCCGGGTACCGAGCTCGAATTCACTGGCCGTCGTTTTACAACGT -y '{name}' -G pr_plasmid1=GAATCAGGGGATAACGCAGGAAAGAACATG -y '{name}' -G pr_plasmid2=AACAGCTATGACCATTACGCCAAGCTTGGCCTGCAGGTCGACTCTAGAGGATCCCCGGGTACCGAGCTCGAATTCACTGGCCGTCGTTTTACAACGT -y '{name}' -o $WORK_DIR/cutadapt/cut_barc_praimers_cutIL_${file}_1.fq.gz -p $WORK_DIR/cutadapt/cut_barc_praimers_cutIL_${file}_2.fq.gz $WORK_DIR/cutadapt/cutIL_${file}_1.fq.gz $WORK_DIR/cutadapt/cutIL_${file}_2.fq.gz --overlap 10;
done
