A computational procedure to identify and annotate somatic mutations in next-generation sequencing data