Herpes simplex virus type 1 causes mucocutaneous lesions, and is the leading cause of infectious blindness in the United States. Animal studies have shown that the severity of HSV-1 ocular disease is influenced by three main factors; innate immunity, host immune response and viral strain. We previously showed that mixed infection with two avirulent HSV-1 strains (OD4 and CJ994) resulted in recombinants that exhibit a range of disease phenotypes from severe to avirulent, suggesting epistatic interactions were involved. The goal of this study was to develop a quantitative trait locus (QTL) analysis of HSV-1 ocular virulence determinants and to identify virulence associated SNPs. Blepharitis and stromal keratitis quantitative scores were characterized for 40 OD4:CJ994 recombinants. Viral titers in the eye were also measured. Virulence quantitative trait locus mapping (vQTLmap) was performed using the Lasso, Random Forest, and Ridge regression methods to identify significant phenotypically meaningful regions for each ocular disease parameter. The most predictive Ridge regression model identified several phenotypically meaningful SNPs for blepharitis and stromal keratitis. Notably, phenotypically meaningful nonsynonymous variations were detected in the UL24, UL29 (ICP8), UL41 (VHS), UL53 (gK), UL54 (ICP27), UL56, ICP4, US1 (ICP22), US3 and gG genes. Network analysis revealed that many of these variations were in HSV-1 regulatory networks and viral genes that affect innate immunity. Several genes previously implicated in virulence were identified, validating this approach, while other genes were novel. Several novel polymorphisms were also identified in these genes. This approach provides a framework that will be useful for identifying virulence genes in other pathogenic viruses, as well as epistatic effects that affect HSV-1 ocular virulence.