Known Ensembl gene levels were quantified by HTSeq version 0.6.0 in intersection-strict mode (the BAM file was streamed to HTSeq through novosort version 1.0.1, as HTSeq accepts read-name-sorted alignments). This provides an integral count of reads for each gene in each sample to be used in downstream analyses (a sample-by-gene “read count matrix”).