RNA-sequencing analysis excels in its ability to infer transcriptomic differences between discrete genetic states. This advantage is often leveraged to explore the consequences of single mutations causing disease, where prior knowledge of expected outcomes may not be established. Successful interpretation of such RNA-seq studies reveals the direct impacts of the mutation, as well as the homeostatic responses of the biological system. Recent studies have highlighted that, when homozygous mutations are studied in non-isogenic backgrounds, genes from the same chromosome as the mutation tend to be over-represented among differentially expressed (DE) genes. One hypothesis suggests that DE genes chromosomally linked to the mutation may not be true biological responses to the disruption of the mutation but, instead, result from differences in the representation of expression quantitative trait loci (eQTLs) that differ between the sample groups being compared. This can be problematic since including spurious DE genes in a functional enrichment study may result in incorrect inferences of mutation effect. Here we show thatchromosomally co-located differentially expressed genes(CC-DEGs) can also be observed in analyses of dominant mutations in heterozygotes. We define a method and a metric to quantify, in RNA-seq data, localised differential allelic representation (DAR) between groups of samples subject to DE analysis. We show how the DAR metric can predict regions prone to eQTL-driven differential expression, and how it can improve functional enrichment analyses through gene exclusion or weighting of gene-level rankings. Advantageously, this improved ability to identify probable eQTLs also reveals examples of CC-DEGs thatarelikely to be functionally related to a mutant phenotype. This situation was predicted by R.A. Fisher in 1930 as due to selection for advantageous linkage disequilibrium after chromosomal rearrangements. By comparing the genomes of zebrafish (Danio rerio) and medaka (Oryzias latipes), a teleost with a conserved ancestral karyotype, we find possible examples of chromosomal aggregation of CC-DEGs during evolution of the zebrafish lineage. The ability to identify and exclude eQTL artefacts from true transcriptomic responses to mutation provides exciting opportunities for developing new approaches for analysing RNA-seq data. Our DAR metric provides a solid foundation for addressing the eQTL issue in new and existing datasets because it relies solely on RNA-seq data, which will thus improve our understanding of the mechanisms of disease-causing mutations.