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

WORK_DIR=/mnt/scratch/ws/psbelokopytova/202208081122polya_data/4c_plasmid/2022-06-20/
THR=5
source activate nn_Polina_tf2
cd $WORK_DIR;
[ -d stat ] || mkdir hic_stat;
declare -a files=("C2" "D1203_E1")

for file in "${files[@]}"
do
python /mnt/scratch/ws/psbelokopytova/202208081122polya_data/4c_plasmid/hic_statistics2.py ${file} $WORK_DIR/align/${file}/${file}_pr_plasmid1_cutTaiI_supplement.sorted.sam $WORK_DIR
done

#declare -a files=("K21" "K22")

#for file in "${files[@]}"
#do
#python /mnt/scratch/ws/psbelokopytova/202208081122polya_data/4c_plasmid/hic_statistics2.py ${file} $WORK_DIR/align/${file}/${file}_pr_plasmid1_cmv_to_cmv_cutTaiI_supplement.sorted.sam
#python /mnt/scratch/ws/psbelokopytova/202208081122polya_data/4c_plasmid/hic_statistics2.py ${file} $WORK_DIR/align/${file}/${file}_pr_plasmid1_empty_to_empty_cutTaiI_supplement.sorted.sam
#done

