# Do you use SNAP for short-read mapping?

I am calling SNPs from WGS samples produced at my lab. I am currently using bwa-mem for mapping Illumina reads as it is recommended by GATK best practice. However, bwa is a bit slow. I heard from my colleague that SNAP is much faster than bwa. I tried it on a small set of reads and it is indeed faster. However, I am not sure how it works with downstream SNP callers, so here are my questions: have you used SNAP for short-read mapping? What is your experience? Does SNAP work well with SNP callers like GATK and freebayes? Thanks!

• If you have the option of using a paid solution, I would strongly urge you to try Sentieon's implementation of the bwa algorithm. I've been using it for a few months and it is several times faster than vanilla bwa. Jun 3 '17 at 13:52

As such you'd be better off using BWA MEM because I assume accurately called variation is always better than fast and incorrectly called variation. The main issue you'll have is ensuring shorter split hits and mapping quality are reported in the same way as bwa MEM -M which GATK/Picard is expecting. Ultimately however you'd be better off posting this question on the GATK forum.