There are a few ways to do this. If you will eventually look at more SNP sites I would recommend you do a full variant calling pipeline. For example: 1. align the fastq raw reads to reference genome (you can use BWA-mem) 2. Remove duplications and sort using samtools and/or picard. 3. Call variants (e.g. using GATK) 4. Annotation and see if the expected SNP(s) is presented.
However, are you only interested in one single SNP? If so, you can just do step one - alignment, and download a genome visualization software such as IGV, and just zoom in to your SNP of interest and take a quick look.