Read mapping and variant calling