Rapidly evolving taxa are excellent models for understanding the mechanisms that give rise to biodiversity. However, developing an accurate historical framework for comparative analysis of such lineages remains a challenge due to ubiquitous incomplete lineage sorting and introgression. Here, we use a whole-genome alignment, multiple locus-sampling strategies, and locus-based and SNP-based species-tree methods to infer a species tree for eastern North American Neodiprion species, a clade of pine-feeding sawflies (Order: Hymenopteran; Family: Diprionidae). We recovered a well-supported species tree that—except for three uncertain relationships—is robust to different strategies for analyzing whole-genome data. Despite this consistency, underlying gene-tree discordance is high. To understand this discordance, we use multiple regression to model topological discordance as a function of several genomic features. We find that gene-tree discordance tends to be higher in regions of the genome that may be more prone to gene-tree estimation error, as indicated by a lower density of parsimony-informative sites, a higher density of genes, a higher average pairwise genetic distance, and gene trees with lower average bootstrap support. Also, contrary to the expectation that discordance via incomplete lineage sorting is reduced in low-recombination regions of the genome, we find a negative correlation between recombination rate and topological discordance. We offer potential explanations for this pattern and hypothesize that it may be unique to lineages that have diverged with gene flow. Our analysis also reveals an unexpected discordance hotspot on Chromosome 1, which contains several genes potentially involved in mitochondrial-nuclear interactions and produces a gene-tree that resembles a highly discordant mitochondrial tree. Based on these observations, we hypothesize that our genome-wide scan for topological discordance has identified a nuclear locus involved in a mito-nuclear incompatibility. Together, these results demonstrate how phylogenomic analysis coupled with high-quality, annotated genomes can generate novel hypotheses about the mechanisms that drive divergence and produce variable genealogical histories across genomes.