for i in {58..64}; do bwa mem GAPDH_reference/GAPDH.fasta 2021-11-23-data/${i}_R1.fq.gz 2021-11-23-data/${i}_R2.fq.gz | samtools view -bS | samtools sort > 2021-11-23-alignments/${i}.sorted.bam; done
