We used the average strain expression signal detected by each probe set as the quantitative trait. QTL mapping was done using QTL Reaper.31, 32 A set of 3795 informative markers were used for mapping,28, 33 and up to a million permutations were performed to calculate the genome-wide empirical P-values and the associated likelihood ratio statistic score.