# Index reference genome
srun -c 1 -p medium --mem 10G -J wgsIndexGenome -o log_%j 01_scripts/06_gatk_dictionnary_reference.sh

# Find regions to re-align
srun -c 1 -p medium --mem 30G -J wgsFindIndels -o log_%j 01_scripts/07_gatk_realign_targets.sh

# Re-align around indels