rule align_reads:
	input:
		genome = multiext("01_raw_data/GENOME", ".amb", ".ann", ".bwt", ".pac", ".sa"),
		r1 = "01_raw_data/{sample}_1.fastq.gz",
		r2 = "01_raw_data/{sample}_2.fastq.gz"
	output: "03_sam_files/{sample}.sam"
	log: "00_logs/{sample}_align_reads_err.log"
	shell:
		"bwa mem -t THREADS '01_raw_data/GENOME' '{input.r1}' '{input.r2}' > '{output}' "
		"2> '{log}'"
