align-hg38-single-end.sh

#!/bin/bash
#SBATCH -c 12                               # Request one core
#SBATCH -N 1                               # Request one node (if you request more than one core with -c, also using
                                           # -N 1 means all cores will be on the same node)
#SBATCH -t 1-00:00                         # Runtime in D-HH:MM format
#SBATCH -p medium                           # Partition to run in
#SBATCH --mem=24G                          # Memory total in MB (for all cores)
#SBATCH -o hostname_%j.out                 # File to which STDOUT will be written, including job ID
#SBATCH -e hostname_%j.err                 # File to which STDERR will be written, including job ID
#SBATCH --mail-type=ALL                    # Type of email notification- BEGIN,END,FAIL,ALL

module load gcc/6.2.0
module load bowtie2/2.3.4.3

for library in SRR*
do
	echo $library

	if test -f "${library}/${library}.bam"; then
		echo "$library already processed."
		continue
	fi

	cd $library

	# alignment against hg38
	bowtie2 -x /n/data2/dfci/medonc/cwu/livius/GRCh38/GCA_000001405.15_GRCh38_no_alt_analysis_set.fna.bowtie_index --very-sensitive -p 12 -U ${library}_1.fastq.gz -S tmp.sam

	# sort aligned reads
	samtools sort tmp.sam -@ 12 > tmp.sorted.sam

	# remove duplicates
	gatk MarkDuplicates -I tmp.sorted.sam -O ${library}.bam --REMOVE_DUPLICATES true --REMOVE_SEQUENCING_DUPLICATES true -M MarkDuplicates.metrics.txt

	# index file 
	samtools index ${library}.bam

	# clean up
	rm tmp*

	cd ..

done