Disentangling the effects of demography and selection has remained a focal point of population genetic analysis. Knowledge about mutation and recombination is essential in this endeavour; however, despite clear evidence that both mutation and recombination rates vary across genomes, it is common practice to model both rates as fixed. In this study, we quantify how this unaccounted for rate heterogeneity may impact inference using common approaches for inferring selection (DFE-alpha, Grapes, and polyDFE) and/or demography (fastsimcoal2 andδaδi). We demonstrate that, if not properly modelled, this heterogeneity can increase uncertainty in the estimation of demographic and selective parameters and in some scenarios may result in mis-leading inference. These results highlight the importance of quantifying the fundamental evolutionary parameters of mutation and recombination prior to utilizing population genomic data to quantify the effects of genetic drift (i.e., as modulated by demographic history) and selection; or, at the least, that the effects of uncertainty in these parameters can and should be directly modelled in downstream inference.Significance StatementDespite evidence from numerous species that mutation and recombination rates vary along the genome, both rates tend to be modelled as fixed when performing population genetic inference. The impact of failing to account for this rate heterogeneity on the estimation of demographic and selective parameters has yet to be well quantified; thus, we here study this effect by comparing inference under both fixed and variable rate scenarios. Our results demonstrate that unaccounted for rate heterogeneity can increase uncertainty and lead to mis-inference in certain scenarios. We highlight the importance of utilizing mutation and recombination rate maps where possible, and of modelling the uncertainty underlying these estimates.