#!/bin/bash
#SBATCH -t 48:00:00
#SBATCH --nodes=1 --ntasks-per-node=14
#SBATCH --mem=120Gb
#SBATCH --export=NONE
#SBATCH --mail-user=thpalma@uri.edu
#SBATCH --mail-type=END,FAIL

# Purge all previously loaded modules
module Purge

# Load SAMtools modules
module load SAMtools/1.9-foss-2018b

# Array of file paths
files=(TP4_trimmed/TP4_trimmed.sam TP5_trimmed/TP5_trimmed.sam TP6_trimmed/TP6_trimmed.sam)

for file in "${files[@]}"; do
  # Extract file name and folder name
  filename="$(basename "$file")"
  folder="$(dirname "$file")"

  # Convert .sam file to .bam file
  bam_file="$folder/${filename%.*}.bam"
  samtools view -bS "$file" > "$bam_file"

  # Sort .bam file
  sorted_bam_file="$folder/${filename%.*}.sorted.bam"
  samtools sort "$bam_file" -o "$sorted_bam_file"

  # Create bam index file
  samtools index "$sorted_bam_file"

  # Create IGV folder and move generated files inside
  mkdir "$folder/IGV"
  mv "$bam_file" "$folder/IGV"
  mv "$sorted_bam_file" "$folder/IGV"
  mv "$sorted_bam_file".bai "$folder/IGV"
done