Map pair-end reads (sequences) to the indexed genome using 32 [t]hreads and compress the result to save space
bwa mem -t 32 path/to/reference.fa path/to/read_pair_end_1.fq.gz path/to/read_pair_end_2.fq.gz | gzip > path/to/alignment_pair_end.sam.gz