Fast and accurate DNA sequence alignment using Burrow-Wheeler Transform
Burrow-Wheeler Aligner for pairwise alignment of DNA
bwa$ bwa index reference.fa$ bwa mem reference.fa reads_R1.fastq reads_R2.fastq > alignment.sam$ bwa mem -t 8 reference.fa reads.fastq | samtools view -b -o output.bam