SNPsplit workflow
The general workflow for working with SNPsplit is outlined below.
Note
Steps 1-3 are only required as a one-off
-
Obtain the SNP annotation file from the Mouse Genomes Project. The latest file is v8 (12/2021).
-
Download the Mouse reference genome from Ensembl (GRCm39)
-
Create an N-masked reference genome using the VCF file from 1.) (either single strain or dual hybrid). This step creates the SNP file required for the
SNPsplit
step (6.))
Note
Steps 4-6 need to be carried out for each sample
-
Use aligner of choice to index the new N-masked reference genome
-
Run alignments against the N-masked reference using your favourite aligner
-
Run
SNPsplit
on the resulting BAM file, using the SNP file generated in 3.)
Once this workflow has completed, proceed to downstream analysis of your choice (outside the scope of SNPsplit).
If you are still using the older mouse genome build GRCm38 you may find some information here.