metaseqR report

Summary

Analysis summary

Summary: The raw bam files, one for each RNA-Seq sample, were summarized to read counts table, using the Bioconductor package GenomicRanges. In the final read counts table, each row represented each column one RNA-Seq sample and each cell, the corresponding read counts associated with each row and column.The gene counts table was normalized for inherent systematic or experimental biases (e.g. sequencing depth, gene length, GC content bias etc.) using the Bioconductor package DESeq after removing genes that had zero counts over all the RNA-Seq samples (27776 genes). The output of the normalization algorithm was a table with normalized counts, which can be used for differential expression analysis with statistical algorithms developed specifically for count data. Prior to the statistical testing procedure, the gene read counts were filtered for possible artifacts that could affect the subsequent statistical testing procedures. Genes presenting any of the following were excluded from further analysis: i) genes with length less than 500 (3263 genes), ii) genes whose average reads per 100 bp was less than the 25^th quantile of the total normalized distribution of average reads per 100bp (0 genes with cutoff value 0 average reads per 100 bp), iii) genes with read counts below the median read counts of the total normalized count distribution (9068 genes with cutoff value 4 normalized read counts). The total number of genes excluded due to the application of gene filters was 3306. The total (unified) number of genes excluded due to the application of all filters was 38105. The resulting gene counts table was subjected to differential expression analysis for the contrasts Grp4aMinus versus Grp4aPlus, Grp4fMinus versus Grp4fPlus using the Bioconductor packages DESeq, edgeR, limma, NBPSeq. In order to combine the statistical significance from multiple algorithms and perform meta-analysis, the PANDORA weighted p-value across results method was applied. The final numbers of differentially expressed genes were (per contrast): for the contrast Grp4aMinus versus Grp4aPlus, 574 (16) statistically significant genes were found with a p-value (FDR or adjusted p-value) threshold of 0.05 and of these, 164 (10) were up-regulated, 211 (6) were down-regulated and 199 (0) were not differentially expressed according to an absolute fold change cutoff value of 1 in log₂ scale, for the contrast Grp4fMinus versus Grp4fPlus, 904 (16) statistically significant genes were found with a p-value (FDR or adjusted p-value) threshold of 0.05 and of these, 225 (4) were up-regulated, 337 (12) were down-regulated and 342 (0) were not differentially expressed according to an absolute fold change cutoff value of 1 in log₂ scale. Literature references for all the algorithms used can be found at the end of this report.

Input options

Read counts file: imported sam/bam/bed files
Conditions: Grp4aPlus, Grp4aMinus, Grp4fPlus, Grp4fMinus, Cplus1Plus
Samples included: YK3R1-5291-4aplus, YK3R2-5291-4aplus, YK3R3-5291-4aplus, YK3R4-5291-4aminus, YK3R5-5291-4aminus, YK3R6-5269-4fplus, YK3R7-5269-4fplus, YK3R8-5269-4fplus, YK3R9-5269-4fminus, YK3R10-5269-4fminus, YK3R11-Cplus-1plus
Samples excluded: none
Requested contrasts: Grp4aMinus_vs_Grp4aPlus, Grp4fMinus_vs_Grp4fPlus
Library sizes:

YK3R1-5291-4aplus: 2675935
YK3R2-5291-4aplus: 2224404
YK3R3-5291-4aplus: 1906090
YK3R4-5291-4aminus: 1858088
YK3R5-5291-4aminus: 3085115
YK3R6-5269-4fplus: 2871955
YK3R7-5269-4fplus: 2064189
YK3R8-5269-4fplus: 1637979
YK3R9-5269-4fminus: 3095887
YK3R10-5269-4fminus: 1966860
YK3R11-Cplus-1plus: 2468119

Annotation: download
Organism: mouse (Mus musculus), genome version alias mm10
Annotation source: Ensembl genomes
Count type: utr
Exon filters: none applied
Gene filters: length, avg.reads, expression, biotype

length
- length: 500

avg.reads
- average.per.bp: 100
- quantile: 0.25

expression
- median: TRUE
- mean: FALSE
- quantile: NA
- known: NA
- custom: NA

biotype
- pseudogene: FALSE
- snRNA: FALSE
- protein_coding: FALSE
- antisense: FALSE
- miRNA: FALSE
- snoRNA: FALSE
- lincRNA: FALSE
- processed_transcript: FALSE
- misc_RNA: FALSE
- rRNA: TRUE
- sense_intronic: FALSE
- sense_overlapping: FALSE
- polymorphic_pseudogene: FALSE
- IG_C_gene: FALSE
- IG_J_gene: FALSE
- IG_D_gene: FALSE
- IG_LV_gene: FALSE
- IG_V_gene: FALSE
- IG_V_pseudogene: TRUE
- TR_V_gene: FALSE
- TR_V_pseudogene: TRUE
- three_prime_overlapping_ncrna: FALSE

Filter application: after normalization
Normalization algorithm: DESeq
Normalization arguments: locfunc

[[list(locfunc = function (x, na.rm = FALSE) standardGeneric("median"))locfunc

Statistical algorithm(s): DESeq, edgeR, limma, NBPSeq
Statistical arguments for DESeq: method, sharingMode, fitType

method: blind
sharingMode: fit-only
fitType: local

Statistical arguments for edgeR: main.method, rowsum.filter, prior.df, trend, span, tag.method, grid.length, grid.range, offset, glm.method, subset, AveLogCPM, trend.method, dispersion, offset, weights, lib.size, prior.count, start, method, test, abundance.trend, robust, winsor.tail.p

main.method: classic
rowsum.filter: 5
prior.df: 10
trend: movingave
tag.method: grid
grid.length: 11
grid.range: -6, 6
glm.method: CoxReid
subset: 10000
trend.method: auto
prior.count: 0.125
method: auto
test: chisq
abundance.trend: TRUE
robust: FALSE
winsor.tail.p: 0.05, 0.1

Statistical arguments for limma: normalize.method

normalize.method: none

Statistical arguments for NBPSeq: main.method, model, tests, alternative

main.method: nbsmyth
model: log-linear-rel-mean, NBP
tests: HOA
alternative: two.sided

Meta-analysis method: PANDORA weighted p-value across results
Multiple testing correction: Benjamini-Hochberg FDR
p-value threshold: 0.05
Logarithmic tranformation offset: 1
Analysis preset: not available
Quality control plots: multidimensional scaling, biotype detection, biotype counts, sample and biotype saturation, filtered biotypes, correlation heatmap and correlogram, pairwise scatterplots between samples, boxplots, GC-content bias, transcript length bias, DEG biotype detection, volcano plot, DEG heatmap
Figure format: png, pdf
Output directory: /data/images/proton2/run342/www/metaseqr_pandora_run342a
Output data: Annotation, p-value, Adjusted p-value (FDR), Fold change, Statistics, Read counts
Output scale(s): log2 scale, Reads per Gene Model
Output values: Normalized values
Output statistics: Mean
Total run time: 16 minutes 24 seconds

Filtered genes

Number of filtered genes: 38105 which is the union of

Filtered because of zero reads: 27776
Filtered because of exon filters: 0
Filtered because of gene filters: 10329 which is the union of
- length: 3263 genes with filter cutoff value 500
- avg.reads: 0 genes with filter cutoff value 0
- expression: 9068 genes further decomposed to (filter name, filtered genes, filter cutoff):
- - median: 9068 genes with filter cutoff value 4
- biotype: 42 genes with filter cutoff value rRNA, IG_V_pseudogene, TR_V_pseudogene

Differentially expressed genes

Number of differentially expressed genes per contrast:

Grp4aMinus_vs_Grp4aPlus: 574 (16) statistically significant genes of which 164 (10) up regulated, 211 (6) down regulated and 199 (0) not differentially expressed according to a p-value (FDR or adjusted p-value) threshold of 0.05 and an absolute fold change cutoff value of 1 in log₂ scale. These numbers refer to the combined analysis performed by metaseqR. Per statistical algorithm, the differentially expressed genes are:
- DESeq: 493 (40) statistically significant genes of which 143 (18) up regulated, 163 (22) down regulated and 187 (0) not differentially expressed according to a p-value (FDR or adjusted p-value) threshold of 0.05 and an absolute fold change cutoff value of 1 in log₂ scale.
- edgeR: 566 (3) statistically significant genes of which 156 (2) up regulated, 194 (1) down regulated and 216 (0) not differentially expressed according to a p-value (FDR or adjusted p-value) threshold of 0.05 and an absolute fold change cutoff value of 1 in log₂ scale.
- limma: 587 (0) statistically significant genes of which 125 (0) up regulated, 188 (0) down regulated and 274 (0) not differentially expressed according to a p-value (FDR or adjusted p-value) threshold of 0.05 and an absolute fold change cutoff value of 1 in log₂ scale.
- NBPSeq: 800 (54) statistically significant genes of which 225 (26) up regulated, 383 (27) down regulated and 192 (1) not differentially expressed according to a p-value (FDR or adjusted p-value) threshold of 0.05 and an absolute fold change cutoff value of 1 in log₂ scale.
Grp4fMinus_vs_Grp4fPlus: 904 (16) statistically significant genes of which 225 (4) up regulated, 337 (12) down regulated and 342 (0) not differentially expressed according to a p-value (FDR or adjusted p-value) threshold of 0.05 and an absolute fold change cutoff value of 1 in log₂ scale. These numbers refer to the combined analysis performed by metaseqR. Per statistical algorithm, the differentially expressed genes are:
- DESeq: 567 (18) statistically significant genes of which 184 (2) up regulated, 257 (16) down regulated and 126 (0) not differentially expressed according to a p-value (FDR or adjusted p-value) threshold of 0.05 and an absolute fold change cutoff value of 1 in log₂ scale.
- edgeR: 1062 (16) statistically significant genes of which 253 (5) up regulated, 303 (11) down regulated and 506 (0) not differentially expressed according to a p-value (FDR or adjusted p-value) threshold of 0.05 and an absolute fold change cutoff value of 1 in log₂ scale.
- limma: 1031 (0) statistically significant genes of which 183 (0) up regulated, 309 (0) down regulated and 539 (0) not differentially expressed according to a p-value (FDR or adjusted p-value) threshold of 0.05 and an absolute fold change cutoff value of 1 in log₂ scale.
- NBPSeq: 1166 (70) statistically significant genes of which 305 (16) up regulated, 480 (51) down regulated and 381 (3) not differentially expressed according to a p-value (FDR or adjusted p-value) threshold of 0.05 and an absolute fold change cutoff value of 1 in log₂ scale.

Run log

[ 2017-06-22 17:24:25 ] INFO 2017-06-22 17:24:25: Data processing started...
[ 2017-06-22 17:24:25 ] INFO Read counts file: imported sam/bam/bed files
[ 2017-06-22 17:24:25 ] INFO Conditions: Grp4aPlus, Grp4aMinus, Grp4fPlus, Grp4fMinus, Cplus1Plus
[ 2017-06-22 17:24:25 ] INFO Samples to include: YK3R1-5291-4aplus, YK3R2-5291-4aplus, YK3R3-5291-4aplus, YK3R4-5291-4aminus, YK3R5-5291-4aminus, YK3R6-5269-4fplus, YK3R7-5269-4fplus, YK3R8-5269-4fplus, YK3R9-5269-4fminus, YK3R10-5269-4fminus, YK3R11-Cplus-1plus
[ 2017-06-22 17:24:25 ] INFO Samples to exclude: none
[ 2017-06-22 17:24:25 ] INFO Requested contrasts: Grp4aMinus_vs_Grp4aPlus, Grp4fMinus_vs_Grp4fPlus
[ 2017-06-22 17:24:25 ] INFO Annotation: download
[ 2017-06-22 17:24:25 ] INFO Organism: mm10
[ 2017-06-22 17:24:25 ] INFO Reference source: ensembl
[ 2017-06-22 17:24:25 ] INFO Count type: utr
[ 2017-06-22 17:24:25 ] INFO Exon filters: none applied
[ 2017-06-22 17:24:25 ] INFO Gene filters: length, avg.reads, expression, biotype
[ 2017-06-22 17:24:25 ] INFO length:
[ 2017-06-22 17:24:25 ] INFO length: 500
[ 2017-06-22 17:24:25 ] INFO avg.reads:
[ 2017-06-22 17:24:25 ] INFO average.per.bp: 100
[ 2017-06-22 17:24:25 ] INFO quantile: 0.25
[ 2017-06-22 17:24:25 ] INFO expression:
[ 2017-06-22 17:24:25 ] INFO median: TRUE
[ 2017-06-22 17:24:25 ] INFO mean: FALSE
[ 2017-06-22 17:24:25 ] INFO quantile: NA
[ 2017-06-22 17:24:25 ] INFO known: NA
[ 2017-06-22 17:24:25 ] INFO custom: NA
[ 2017-06-22 17:24:25 ] INFO biotype:
[ 2017-06-22 17:24:25 ] INFO pseudogene: FALSE
[ 2017-06-22 17:24:25 ] INFO snRNA: FALSE
[ 2017-06-22 17:24:25 ] INFO protein_coding: FALSE
[ 2017-06-22 17:24:25 ] INFO antisense: FALSE
[ 2017-06-22 17:24:25 ] INFO miRNA: FALSE
[ 2017-06-22 17:24:25 ] INFO snoRNA: FALSE
[ 2017-06-22 17:24:25 ] INFO lincRNA: FALSE
[ 2017-06-22 17:24:25 ] INFO processed_transcript: FALSE
[ 2017-06-22 17:24:25 ] INFO misc_RNA: FALSE
[ 2017-06-22 17:24:25 ] INFO rRNA: TRUE
[ 2017-06-22 17:24:25 ] INFO sense_intronic: FALSE
[ 2017-06-22 17:24:25 ] INFO sense_overlapping: FALSE
[ 2017-06-22 17:24:25 ] INFO polymorphic_pseudogene: FALSE
[ 2017-06-22 17:24:25 ] INFO IG_C_gene: FALSE
[ 2017-06-22 17:24:25 ] INFO IG_J_gene: FALSE
[ 2017-06-22 17:24:25 ] INFO IG_D_gene: FALSE
[ 2017-06-22 17:24:25 ] INFO IG_LV_gene: FALSE
[ 2017-06-22 17:24:25 ] INFO IG_V_gene: FALSE
[ 2017-06-22 17:24:25 ] INFO IG_V_pseudogene: TRUE
[ 2017-06-22 17:24:25 ] INFO TR_V_gene: FALSE
[ 2017-06-22 17:24:25 ] INFO TR_V_pseudogene: TRUE
[ 2017-06-22 17:24:25 ] INFO three_prime_overlapping_ncrna: FALSE
[ 2017-06-22 17:24:25 ] INFO Filter application: postnorm
[ 2017-06-22 17:24:25 ] INFO Normalization algorithm: deseq
[ 2017-06-22 17:24:25 ] INFO Normalization arguments:
[ 2017-06-22 17:24:25 ] INFO locfunc:
[ 2017-06-22 17:24:25 ] INFO [[
[ 2017-06-22 17:24:25 ] INFO list(locfunc = function (x, na.rm = FALSE) standardGeneric("median"))
[ 2017-06-22 17:24:25 ] INFO locfunc
[ 2017-06-22 17:24:25 ] INFO Statistical algorithm: deseq, edger, limma, nbpseq
[ 2017-06-22 17:24:25 ] INFO Statistical arguments:
[ 2017-06-22 17:24:25 ] INFO deseq: blind, fit-only, local
[ 2017-06-22 17:24:25 ] INFO edger: classic, 5, 10, movingave, NULL, grid, 11, c(-6, 6), NULL, CoxReid, 10000, NULL, auto, NULL, NULL, NULL, NULL, 0.125, NULL, auto, chisq, TRUE, FALSE, c(0.05, 0.1)
[ 2017-06-22 17:24:25 ] INFO limma: none
[ 2017-06-22 17:24:25 ] INFO nbpseq: nbsmyth, list(nbpseq = "log-linear-rel-mean", nbsmyth = "NBP"), HOA, two.sided
[ 2017-06-22 17:24:25 ] INFO Meta-analysis method: pandora
[ 2017-06-22 17:24:25 ] INFO Multiple testing correction: BH
[ 2017-06-22 17:24:25 ] INFO p-value threshold: 0.05
[ 2017-06-22 17:24:25 ] INFO Logarithmic transformation offset: 1
[ 2017-06-22 17:24:25 ] INFO Quality control plots: mds, biodetection, countsbio, saturation, readnoise, filtered, correl, pairwise, boxplot, gcbias, lengthbias, biodist, volcano, deheatmap
[ 2017-06-22 17:24:25 ] INFO Figure format: png, pdf
[ 2017-06-22 17:24:25 ] INFO Output directory: /data/images/proton2/run342/www/metaseqr_pandora_run342a
[ 2017-06-22 17:24:25 ] INFO Output data: annotation, p.value, adj.p.value, fold.change, meta.p.value, adj.meta.p.value, stats, counts, flags
[ 2017-06-22 17:24:25 ] INFO Output scale(s): log2, rpgm
[ 2017-06-22 17:24:25 ] INFO Output values: normalized
[ 2017-06-22 17:24:25 ] INFO Output statistics: mean
[ 2017-06-22 17:24:25 ] INFO Downloading gene annotation for mm10...
[ 2017-06-22 17:24:33 ] INFO Downloading transcript annotation for mm10...
[ 2017-06-22 17:27:05 ] INFO Converting annotation to GenomicRanges object...
[ 2017-06-22 17:27:05 ] INFO Merging transcript 3' UTRs to create unique gene/transcript models...
[ 2017-06-22 17:30:36 ] INFO Flanking merged transcript 3' UTRs per 500bp...
[ 2017-06-22 17:30:39 ] INFO Reading bam file YK3R1-5291-4aplus.bam for sample with name YK3R1-5291-4aplus. This might take some time...
[ 2017-06-22 17:30:39 ] INFO Reading bam file YK3R2-5291-4aplus.bam for sample with name YK3R2-5291-4aplus. This might take some time...
[ 2017-06-22 17:30:39 ] INFO Reading bam file YK3R3-5291-4aplus.bam for sample with name YK3R3-5291-4aplus. This might take some time...
[ 2017-06-22 17:30:39 ] INFO Reading bam file YK3R4-5291-4aminus.bam for sample with name YK3R4-5291-4aminus. This might take some time...
[ 2017-06-22 17:30:39 ] INFO Reading bam file YK3R5-5291-4aminus.bam for sample with name YK3R5-5291-4aminus. This might take some time...
[ 2017-06-22 17:30:39 ] INFO Reading bam file YK3R6-5269-4fplus.bam for sample with name YK3R6-5269-4fplus. This might take some time...
[ 2017-06-22 17:30:39 ] INFO Reading bam file YK3R7-5269-4fplus.bam for sample with name YK3R7-5269-4fplus. This might take some time...
[ 2017-06-22 17:30:39 ] INFO Reading bam file YK3R8-5269-4fplus.bam for sample with name YK3R8-5269-4fplus. This might take some time...
[ 2017-06-22 17:30:39 ] INFO Reading bam file YK3R9-5269-4fminus.bam for sample with name YK3R9-5269-4fminus. This might take some time...
[ 2017-06-22 17:30:40 ] INFO Reading bam file YK3R10-5269-4fminus.bam for sample with name YK3R10-5269-4fminus. This might take some time...
[ 2017-06-22 17:30:40 ] INFO Reading bam file YK3R11-Cplus-1plus.bam for sample with name YK3R11-Cplus-1plus. This might take some time...
[ 2017-06-22 17:30:42 ] INFO Counting reads overlapping with given annotation...
[ 2017-06-22 17:30:42 ] INFO ...for single-end reads...
[ 2017-06-22 17:30:42 ] INFO ...assuming forward sequenced reads...
[ 2017-06-22 17:30:43 ] INFO Counting reads overlapping with given annotation...
[ 2017-06-22 17:30:43 ] INFO ...for single-end reads...
[ 2017-06-22 17:30:43 ] INFO ...assuming forward sequenced reads...
[ 2017-06-22 17:30:43 ] INFO Counting reads overlapping with given annotation...
[ 2017-06-22 17:30:43 ] INFO ...for single-end reads...
[ 2017-06-22 17:30:43 ] INFO ...assuming forward sequenced reads...
[ 2017-06-22 17:30:44 ] INFO Counting reads overlapping with given annotation...
[ 2017-06-22 17:30:44 ] INFO ...for single-end reads...
[ 2017-06-22 17:30:44 ] INFO ...assuming forward sequenced reads...
[ 2017-06-22 17:30:44 ] INFO Counting reads overlapping with given annotation...
[ 2017-06-22 17:30:44 ] INFO ...for single-end reads...
[ 2017-06-22 17:30:44 ] INFO ...assuming forward sequenced reads...
[ 2017-06-22 17:30:44 ] INFO Counting reads overlapping with given annotation...
[ 2017-06-22 17:30:44 ] INFO ...for single-end reads...
[ 2017-06-22 17:30:44 ] INFO ...assuming forward sequenced reads...
[ 2017-06-22 17:30:45 ] INFO Counting reads overlapping with given annotation...
[ 2017-06-22 17:30:45 ] INFO ...for single-end reads...
[ 2017-06-22 17:30:45 ] INFO ...assuming forward sequenced reads...
[ 2017-06-22 17:30:45 ] INFO Counting reads overlapping with given annotation...
[ 2017-06-22 17:30:45 ] INFO ...for single-end reads...
[ 2017-06-22 17:30:45 ] INFO ...assuming forward sequenced reads...
[ 2017-06-22 17:30:45 ] INFO Counting reads overlapping with given annotation...
[ 2017-06-22 17:30:45 ] INFO ...for single-end reads...
[ 2017-06-22 17:30:45 ] INFO ...assuming forward sequenced reads...
[ 2017-06-22 17:30:46 ] INFO Counting reads overlapping with given annotation...
[ 2017-06-22 17:30:46 ] INFO ...for single-end reads...
[ 2017-06-22 17:30:46 ] INFO ...assuming forward sequenced reads...
[ 2017-06-22 17:30:46 ] INFO Counting reads overlapping with given annotation...
[ 2017-06-22 17:30:46 ] INFO ...for single-end reads...
[ 2017-06-22 17:30:46 ] INFO ...assuming forward sequenced reads...
[ 2017-06-22 17:31:05 ] INFO Checking chromosomes in transcript counts and gene annotation...
[ 2017-06-22 17:31:05 ] INFO Processing transcripts...
[ 2017-06-22 17:31:05 ] INFO Separating transcripts (UTR regions) per gene for YK3R1-5291-4aplus...
[ 2017-06-22 17:31:23 ] INFO Separating transcripts (UTR regions) per gene for YK3R2-5291-4aplus...
[ 2017-06-22 17:31:49 ] INFO Separating transcripts (UTR regions) per gene for YK3R3-5291-4aplus...
[ 2017-06-22 17:32:15 ] INFO Separating transcripts (UTR regions) per gene for YK3R4-5291-4aminus...
[ 2017-06-22 17:32:41 ] INFO Separating transcripts (UTR regions) per gene for YK3R5-5291-4aminus...
[ 2017-06-22 17:33:08 ] INFO Separating transcripts (UTR regions) per gene for YK3R6-5269-4fplus...
[ 2017-06-22 17:33:35 ] INFO Separating transcripts (UTR regions) per gene for YK3R7-5269-4fplus...
[ 2017-06-22 17:34:02 ] INFO Separating transcripts (UTR regions) per gene for YK3R8-5269-4fplus...
[ 2017-06-22 17:34:29 ] INFO Separating transcripts (UTR regions) per gene for YK3R9-5269-4fminus...
[ 2017-06-22 17:34:56 ] INFO Separating transcripts (UTR regions) per gene for YK3R10-5269-4fminus...
[ 2017-06-22 17:35:26 ] INFO Separating transcripts (UTR regions) per gene for YK3R11-Cplus-1plus...
[ 2017-06-22 17:35:55 ] INFO Saving gene model to /data/images/proton2/run342/www/metaseqr_pandora_run342a/data/gene_model.RData
[ 2017-06-22 17:35:58 ] INFO Summarizing count data...
[ 2017-06-22 17:36:22 ] INFO Removing genes with zero counts in all samples...
[ 2017-06-22 17:36:24 ] INFO Normalizing with: deseq
[ 2017-06-22 17:36:24 ] INFO Applying gene filter length...
[ 2017-06-22 17:36:24 ] INFO Threshold below which ignored: 500
[ 2017-06-22 17:36:24 ] INFO Applying gene filter avg.reads...
[ 2017-06-22 17:36:24 ] INFO Threshold below which ignored: 0
[ 2017-06-22 17:36:24 ] INFO Applying gene filter expression...
[ 2017-06-22 17:36:24 ] INFO Threshold below which ignored: 4
[ 2017-06-22 17:36:24 ] INFO Applying gene filter biotype...
[ 2017-06-22 17:36:24 ] INFO Biotypes ignored: rRNA, IG_V_pseudogene, TR_V_pseudogene
[ 2017-06-22 17:36:24 ] INFO 38105 genes filtered out
[ 2017-06-22 17:36:24 ] INFO 13684 genes remain after filtering
[ 2017-06-22 17:36:24 ] INFO Running statistical tests with: deseq
[ 2017-06-22 17:36:25 ] INFO Contrast: Grp4aMinus_vs_Grp4aPlus
[ 2017-06-22 17:36:31 ] INFO Contrast: Grp4fMinus_vs_Grp4fPlus
[ 2017-06-22 17:36:37 ] INFO Contrast Grp4aMinus_vs_Grp4aPlus: found 493 genes
[ 2017-06-22 17:36:37 ] INFO Contrast Grp4fMinus_vs_Grp4fPlus: found 567 genes
[ 2017-06-22 17:36:37 ] INFO Running statistical tests with: edger
[ 2017-06-22 17:36:39 ] INFO Contrast: Grp4aMinus_vs_Grp4aPlus
[ 2017-06-22 17:36:40 ] INFO Contrast: Grp4fMinus_vs_Grp4fPlus
[ 2017-06-22 17:36:41 ] INFO Contrast Grp4aMinus_vs_Grp4aPlus: found 566 genes
[ 2017-06-22 17:36:41 ] INFO Contrast Grp4fMinus_vs_Grp4fPlus: found 1062 genes
[ 2017-06-22 17:36:41 ] INFO Running statistical tests with: limma
[ 2017-06-22 17:36:41 ] INFO Contrast: Grp4aMinus_vs_Grp4aPlus
[ 2017-06-22 17:36:43 ] INFO Contrast: Grp4fMinus_vs_Grp4fPlus
[ 2017-06-22 17:36:45 ] INFO Contrast Grp4aMinus_vs_Grp4aPlus: found 587 genes
[ 2017-06-22 17:36:45 ] INFO Contrast Grp4fMinus_vs_Grp4fPlus: found 1031 genes
[ 2017-06-22 17:36:45 ] INFO Running statistical tests with: nbpseq
[ 2017-06-22 17:36:45 ] INFO Contrast: Grp4aMinus_vs_Grp4aPlus
[ 2017-06-22 17:36:52 ] INFO Contrast: Grp4fMinus_vs_Grp4fPlus
[ 2017-06-22 17:36:59 ] INFO Contrast Grp4aMinus_vs_Grp4aPlus: found 800 genes
[ 2017-06-22 17:36:59 ] INFO Contrast Grp4fMinus_vs_Grp4fPlus: found 1166 genes
[ 2017-06-22 17:36:59 ] INFO Performing meta-analysis with pandora
[ 2017-06-22 17:37:00 ] INFO Building output files...
[ 2017-06-22 17:37:03 ] INFO Contrast: Grp4aMinus_vs_Grp4aPlus
[ 2017-06-22 17:37:03 ] INFO Adding non-filtered data...
[ 2017-06-22 17:37:03 ] INFO binding annotation...
[ 2017-06-22 17:37:04 ] INFO binding p-values...
[ 2017-06-22 17:37:04 ] INFO binding FDRs...
[ 2017-06-22 17:37:04 ] INFO binding meta p-values...
[ 2017-06-22 17:37:05 ] INFO binding adjusted meta p-values...
[ 2017-06-22 17:37:05 ] INFO binding log2 normalized fold changes...
[ 2017-06-22 17:37:05 ] INFO binding normalized mean counts...
[ 2017-06-22 17:37:05 ] INFO binding normalized mean counts...
[ 2017-06-22 17:37:06 ] INFO binding all normalized counts for Grp4aMinus...
[ 2017-06-22 17:37:06 ] INFO binding all normalized counts for Grp4aPlus...
[ 2017-06-22 17:37:06 ] INFO binding filtering flags...
[ 2017-06-22 17:37:07 ] INFO Writing output...
[ 2017-06-22 17:37:07 ] INFO Adding filtered data...
[ 2017-06-22 17:37:07 ] INFO binding annotation...
[ 2017-06-22 17:37:07 ] INFO binding p-values...
[ 2017-06-22 17:37:07 ] INFO binding FDRs...
[ 2017-06-22 17:37:07 ] INFO binding meta p-values...
[ 2017-06-22 17:37:07 ] INFO binding adjusted meta p-values...
[ 2017-06-22 17:37:07 ] INFO binding log2 normalized fold changes...
[ 2017-06-22 17:37:07 ] INFO binding normalized mean counts...
[ 2017-06-22 17:37:08 ] INFO binding normalized mean counts...
[ 2017-06-22 17:37:09 ] INFO binding all normalized counts for Grp4aMinus...
[ 2017-06-22 17:37:09 ] INFO binding all normalized counts for Grp4aPlus...
[ 2017-06-22 17:37:09 ] INFO binding filtering flags...
[ 2017-06-22 17:37:10 ] INFO Writing output...
[ 2017-06-22 17:37:12 ] INFO Contrast: Grp4fMinus_vs_Grp4fPlus
[ 2017-06-22 17:37:12 ] INFO Adding non-filtered data...
[ 2017-06-22 17:37:12 ] INFO binding annotation...
[ 2017-06-22 17:37:13 ] INFO binding p-values...
[ 2017-06-22 17:37:13 ] INFO binding FDRs...
[ 2017-06-22 17:37:13 ] INFO binding meta p-values...
[ 2017-06-22 17:37:13 ] INFO binding adjusted meta p-values...
[ 2017-06-22 17:37:14 ] INFO binding log2 normalized fold changes...
[ 2017-06-22 17:37:14 ] INFO binding normalized mean counts...
[ 2017-06-22 17:37:14 ] INFO binding normalized mean counts...
[ 2017-06-22 17:37:15 ] INFO binding all normalized counts for Grp4fMinus...
[ 2017-06-22 17:37:15 ] INFO binding all normalized counts for Grp4fPlus...
[ 2017-06-22 17:37:16 ] INFO binding filtering flags...
[ 2017-06-22 17:37:16 ] INFO Writing output...
[ 2017-06-22 17:37:16 ] INFO Adding filtered data...
[ 2017-06-22 17:37:16 ] INFO binding annotation...
[ 2017-06-22 17:37:16 ] INFO binding p-values...
[ 2017-06-22 17:37:16 ] INFO binding FDRs...
[ 2017-06-22 17:37:16 ] INFO binding meta p-values...
[ 2017-06-22 17:37:16 ] INFO binding adjusted meta p-values...
[ 2017-06-22 17:37:17 ] INFO binding log2 normalized fold changes...
[ 2017-06-22 17:37:17 ] INFO binding normalized mean counts...
[ 2017-06-22 17:37:17 ] INFO binding normalized mean counts...
[ 2017-06-22 17:37:18 ] INFO binding all normalized counts for Grp4fMinus...
[ 2017-06-22 17:37:18 ] INFO binding all normalized counts for Grp4fPlus...
[ 2017-06-22 17:37:18 ] INFO binding filtering flags...
[ 2017-06-22 17:37:19 ] INFO Writing output...
[ 2017-06-22 17:37:21 ] INFO Creating quality control graphs...
[ 2017-06-22 17:37:21 ] INFO Plotting in png format...
[ 2017-06-22 17:37:21 ] INFO Plotting mds...
[ 2017-06-22 17:37:21 ] INFO Plotting biodetection...
[ 2017-06-22 17:37:23 ] INFO Plotting countsbio...
[ 2017-06-22 17:37:25 ] INFO Plotting saturation...
[ 2017-06-22 17:37:32 ] INFO Plotting readnoise...
[ 2017-06-22 17:37:34 ] INFO Plotting correl...
[ 2017-06-22 17:37:34 ] INFO Plotting pairwise...
[ 2017-06-22 17:38:10 ] INFO Plotting boxplot...
[ 2017-06-22 17:38:10 ] INFO Plotting gcbias...
[ 2017-06-22 17:38:11 ] INFO Plotting lengthbias...
[ 2017-06-22 17:38:11 ] INFO Plotting boxplot...
[ 2017-06-22 17:38:11 ] INFO Plotting gcbias...
[ 2017-06-22 17:38:12 ] INFO Plotting lengthbias...
[ 2017-06-22 17:38:13 ] INFO Plotting biodist...
[ 2017-06-22 17:38:13 ] INFO Contrast: Grp4aMinus_vs_Grp4aPlus
[ 2017-06-22 17:38:13 ] INFO Contrast: Grp4fMinus_vs_Grp4fPlus
[ 2017-06-22 17:38:14 ] INFO Plotting volcano...
[ 2017-06-22 17:38:14 ] INFO Contrast: Grp4aMinus_vs_Grp4aPlus
[ 2017-06-22 17:38:14 ] INFO Contrast: Grp4fMinus_vs_Grp4fPlus
[ 2017-06-22 17:38:14 ] INFO Plotting deheatmap...
[ 2017-06-22 17:38:14 ] INFO Contrast: Grp4aMinus_vs_Grp4aPlus
[ 2017-06-22 17:38:46 ] INFO Contrast: Grp4fMinus_vs_Grp4fPlus
[ 2017-06-22 17:39:15 ] INFO Plotting filtered...
[ 2017-06-22 17:39:16 ] INFO Plotting in pdf format...
[ 2017-06-22 17:39:16 ] INFO Plotting mds...
[ 2017-06-22 17:39:16 ] INFO Plotting biodetection...
[ 2017-06-22 17:39:16 ] INFO Plotting countsbio...
[ 2017-06-22 17:39:18 ] INFO Plotting saturation...
[ 2017-06-22 17:39:24 ] INFO Plotting readnoise...
[ 2017-06-22 17:39:26 ] INFO Plotting correl...
[ 2017-06-22 17:39:26 ] INFO Plotting pairwise...
[ 2017-06-22 17:39:45 ] INFO Plotting boxplot...
[ 2017-06-22 17:39:45 ] INFO Plotting gcbias...
[ 2017-06-22 17:39:46 ] INFO Plotting lengthbias...
[ 2017-06-22 17:39:47 ] INFO Plotting boxplot...
[ 2017-06-22 17:39:47 ] INFO Plotting gcbias...
[ 2017-06-22 17:39:48 ] INFO Plotting lengthbias...
[ 2017-06-22 17:39:49 ] INFO Plotting biodist...
[ 2017-06-22 17:39:49 ] INFO Contrast: Grp4aMinus_vs_Grp4aPlus
[ 2017-06-22 17:39:50 ] INFO Contrast: Grp4fMinus_vs_Grp4fPlus
[ 2017-06-22 17:39:50 ] INFO Plotting volcano...
[ 2017-06-22 17:39:50 ] INFO Contrast: Grp4aMinus_vs_Grp4aPlus
[ 2017-06-22 17:39:50 ] INFO Contrast: Grp4fMinus_vs_Grp4fPlus
[ 2017-06-22 17:39:50 ] INFO Plotting deheatmap...
[ 2017-06-22 17:39:50 ] INFO Contrast: Grp4aMinus_vs_Grp4aPlus
[ 2017-06-22 17:40:20 ] INFO Contrast: Grp4fMinus_vs_Grp4fPlus
[ 2017-06-22 17:40:48 ] INFO Plotting filtered...
[ 2017-06-22 17:40:48 ] INFO Creating HTML report...
[ 2017-06-22 17:40:48 ] INFO Compressing figures...

Figures

Multidimensional scaling

The Multi-Dimensional Scaling (MDS) plots comprise a means of visualizing the level of similarity of individual cases of a dataset. It is similar to Principal Component Analysis (PCA), but instead of using the covariance matrix to find similarities among cases, MDS uses absolute distance metrics such as the classical Euclidean distance. Because of the relative linear relations among sequencing samples, it provides a more realistic clustering among samples. MDS serves quality control and it can be interpreted as follows: when the distance among samples of the same biological condition in the MDS space is small, this is an indication of high correlation and reproducibility among them. When this distance is larger or heterogeneous (e.g. the 3rd sample of a triplicate set is further from the other 2), this constitutes an indication of low correlation and reproducibility among samples. It can help exclude poor samples from further analysis.

Biotype detection

The biotype detection bar diagrams are a set of quality control charts that show the percentage of each biotype in the genome (i.e. in the whole set of features provided, for example, protein coding genes, non coding RNAs or pseudogenes) in grey bars, which proportion has been detected in a sample before normalization and after a basic filtering by removing features with zero counts in red lined bars, and the percentage of each biotype within the sample in solid red bars. The difference between grey bars and solid red bars is that the grey bars show the percentage of a feature in the genome while the solid red bars show the percentage in the sample. Thus, the solid red bars may be sometimes higher than the grey bars because certain features (e.g. protein coding genes) may be detected within a sample with a higher proportion relatively to their presence in the genome, as compared with other features. For example, while the percentage of protein coding genes in the whole genome is already higher than other biotypes, this percentage is expected to be even higher in an RNA-Seq experiment where one expects protein-coding genes to exhibit greater abundance. The vertical green line separates the most abundant biotypes (on the left-hand side, corresponding to the left axis scale) from the rest (on the right-hand side, corresponding to the right axis scale). Otherwise, the lower abundance biotypes would be indistinguishable. Unexpected outcomes in this quality control chart (e.g. very low detection of protein coding genes) would signify possible low quality of a sample.

Biotype detection counts

The biotype detection counts boxplots are a set of quality control charts that depict both the biological classification for the detected features and the actual distribution of the read counts for each biological type. The boxplot comprises a means of summarizing the read counts distribution of a sample in the form of a bar with extending lines, as commonly used way of graphically presenting groups of numerical data. A boxplot also indicates which observations, if any, might be considered outliers and is able to visually show different types of populations, without making any assumptions of the underlying statistical distribution. The spacing between the different parts of the box help indicate variance, skewness and identify outliers. The thick bar inside the colored box is the median of the observations while the box extends over the Interquartile Range of the observations. The whiskers extend up (down) to +/-1.5xIQR. Unexpected outcomes (e.g. protein coding read count distribution similar to pseudogene read count distribution) indicates poor sample quality.

Read and biotype saturation

The read and biotype saturation plots are a set of quality control charts that depict the read count saturation levels at several sequencing depths. Thus, they comprise a means of assessing whether the sequencing depth of an RNA-Seq experiment is sufficient in order to detect the biological features under investigation. These quality control charts are separated in two subgroups: the first subgroup (read saturation per biotype for all samples) is a set of plots, one for each biological feature (e.g. protein coding, pseudogene, lincRNA, etc.), that depict the number of detected features in different sequencing depths and for all samples in the same plot. The second subgroup (read saturation per sample for all biotypes) is a set of plots similar to the above, but with, there is one pair of plots with two panels for each sample, presenting all biological features. The left panel depicts the saturation levels for the less abundatnt features, while the right panel, the saturation for the more abundant features, as placing them all together would make the less abundant features indistinguishable. All the saturation plots should be interpreted as follows: if the read counts for a biotype tend to be saturated, the respective curve should tend to reach a plateau at higher depths. Otherwise, more sequencing is needed for the specific biotype.

Read saturation per biotype for all samples

Read saturation per sample for all biotypes

RNA-Seq reads noise

The read noise plots depict the percentage of biological features detected when subsampling the total number of reads. Very steep curves in read noise plots indicate that although the sequencing depth reaches its maximum, a relatively small percentage of total features is detected, indicating that the level of background noise is relatively high. Less steep RNA composition curves, indicate less noise. When a sample's curve deviate from the rest, it could indicate lower or higher quality, depending on the curves of the rest of the samples.

Correlation plots

The sample correlation plots depict the accordance among the RNA-Seq samples, as this is manifested through the read counts table used with the metaseqr pipeline, with two representations that both use the correlation matrix (a matrix which depicts all the pairwise correlations between each pair of samples) of the read counts matrix. The first is a correlation clustered heatmap which depicts the correlations among samples as color-scaled image and the hierarchical clustering tree depicts the grouping of the samples according to their correlation. Samples from the same group that are not clustered together provides an indication that there might be a quality problem with the dataset. The second is a 'correlogram' plot, where again the samples are hierarchically clustered and grouped but this time correlations are presented as ellipses inside each cell. Each cell represents a pairwise comparison and each correlation coefficient is represented by an ellipse whose 'diameter', direction and color depict the accordance for that pair of samples. Highly correlated samples are depicted as ellipses with narrow diameter, while poorly correlated samples are depicted as ellipses with wide diameters. Also, highly correlated samples are depicted as ellipses with a left-to-right upwards direction while poorly correlated samples are depicted as ellipses with a right-to-left upwards direction.

Correlation heatmap

Data correlogram

Pairwise scatterplots

The pairwise comparison plots are split in three parts: the upper diagonal consists of simple scatterplots for all pairwise sample comparisons, together with their Pearson correlation coefficient. It is a simple measure of between sample correlation using all the available data points instead of only the correlation matrix. The lower diagonal consists of mean-difference plots for all pairwise sample comparisons. A mean-difference plot (or a Bland-Altman plots) is a method of data plotting used in analyzing the agreement between two different assays/variables. In this graphical method the differences (or alternatively the ratios) between the two variables are plotted against the averages of the two. Such a plot is useful, for example, for analyzing data with strong correlation between x and y axes, when the (x,y) dots on the plot are close to the diagonal x=y. In this case, the value of the transformed variable X is about the same as x and y and the variable Y shows the difference between x and y. In both represantations, irregular shapes of the red smoother lines are an indication of poor correlation between samples or of other systematic bias sources, which is usually corrected through data normalization.

Boxplots

The boxplot comprises a means of summarizing the read counts distribution of a sample in the form of a bar with extending lines, as a commonly used way of graphically presenting groups of numerical data. A boxplot also indicates which observations, if any, might be considered outliers and is able to visually show different types of populations, without making any assumptions about the underlying statistical distribution. The spacings between the different parts of the box help indicate variance, skewness and identify outliers. The thick bar inside the colored box is the median of the observations while the box extends over the Interquartile Range of the observations. The whiskers extend up (down) to +/-1.5xIQR. Boxplots at similar levels indicate good quality of the normalization. If boxplots remain at different levels after normalization, maybe another normalization algorithm may have to be examined. The un-normalized boxplots show the need for data normalization in order for the data from different samples to follow the same underlying distribution and statistical testing becoming possible.

Boxplot of un-normalized data

Boxplot of normalized data

GC content bias

The GC-content bias plot is a quality control chart that shows the possible dependence of the read counts (in log2 scale) under a gene to the GC content percentage of that gene. In order for the statistical tests to be able to detect statistical significance which occurs due to real biological effects and not by other systematic biases present in the data (e.g. a possible GC-content bias), the latter should be accounted for by the applied normalization algorithm. Although the tests are performed for each gene across biological conditions one could assume that the GC content does not represent a bias, as it is the same for the tested gene across samples and conditions. However, Risso et al. (2011) showed that the GC-content could have an impact in the statistical testing procedure. The GC-content bias plot depicts the dependence of the read counts to the GC content before and after normalization. The smoothing lines for each sample, should be as 'straight' as possible after normalization. In addition, if the smoothing lines differ significantly between biological conditions, this would constitute a possible quality warning.

GC content bias un-normalized

GC content bias normalized

Gene/transcript length bias

The gene/transcript length bias plot is a quality control chart that shows the possible dependence of the read counts (in log2 scale) under a gene to the length that gene (whole gene or sum of exons depending on the analysis). In order for the statistical tests to be able to detect statistical significance which occurs due to real biological effects and not by other systematic biases present in the data (e.g. a possible length bias), the latter should be accounted for by the applied normalization algorithm. Although the tests are performed for each gene across bioogical conditions, one could assume that the gene length does not represent a bias as it's the same for the tested gene across samples and conditions. However, it has been shown in several studies that the gene length could have an impact on the statistical testing procedure. The length bias plot depicts the dependence of the read counts to the gene/transcript length before and after normalization. The smoothing lines for each sample, should be as 'straight' as possible after normalization. In addition, if the smoothing lines differ significantly among biological conditions, this would constitute a possible quality warning.

Gene/transcript length bias un-normalized

Gene/transcript length bias normalized

Chromosome and biotype distribution of filtered genes

The chromosome and biotype distribution of filtered genes is a quality control chart with two rows and four panels: on the left panel of the first row, the bar chart depicts the numbers of filtered genes per chromosome (actual numbers shown above the bars). On the right panel of the first row, the bar chart depicts the numbers of filtered genes per biotype (actual numbers shown above the bars). On the left panel of the second row, the bar chart depicts the fraction of the filtered genes to the total genes per chromosome (actual percentages shown above the bars). On the right panel of the second row, the bar chart depicts the fraction of the filtered genes to the total genes per biotype (actual percentages shown above the bars). This plot should indicate possible quality problems when for example the filtered genes for a specific chromosome (or the fraction) is extremely higher than the rest. Generally, the fractions per chromosome should be uniform and the fractions per biotype should be proportional to the biotype fraction relative to the genome.

Volcano plots

A volcano plot is a scatterplot that is often used when analyzing high-throughput -omics data (e.g. microarray data, RNA-Seq data) to give an overview of interesting genes. The log2 fold change is plotted on the x-axis and the negative log10 p-value is plotted on the y-axis. A volcano plot combines the results of a statistical test (aka, p-values) with the magnitude of the change enabling quick visual identification of those genes that display large-magnitude changes that are also statistically significant. The horizontal dashed line sets the threshold for statistical significance, while the vertical dashed lines set the thresholds for biological significance. It should be noted that the volcano plots become harder to interpret when using more than one statistical algorithm and performing meta-analysis. This happens because the genes that have stronger evidence of being differentially expressed obtain lower p-values while the rest either remain at similar levels or obtain higher p-values. The result is a 'warped' volcano plot, with two main data clouds: one in the upper part of the plot, and one in the lower part of the plot. You can always zoom in when using interacting mode (the default).

Volcano plot for the contrast Grp4aMinus_vs_Grp4aPlus

Volcano plot for the contrast Grp4fMinus_vs_Grp4fPlus

DEG heatmaps

The Differentially Expressed Genes (DEGs) heatmaps depict how well samples from different conditions cluster together according to their expression values after normalization and statistical testing, for each requested statistical contrast. If samples from the same biological condition do not cluster together, this would constitute a warning sign regarding the quality of the samples. In addition, DEG heatmaps provide an initial view of possible clusters of co-expressed genes.

DEG heatmap for the contrast Grp4aMinus_vs_Grp4aPlus

DEG heatmap for the contrast Grp4fMinus_vs_Grp4fPlus

Chromosome and biotype distributions of DEGs

The chromosome and biotype distributions bar diagram for Differentially Expressed Genes (DEGs) is split in two panels: i) on the left panel DEGs are distributed per chromosome and the percentage of each chromosome in the genome is presented in grey bars, the percentage of DEGs in each chromosome is presented in red lined bars and the percentage of certain chromosomes in the distribution of DEGs is presented in solid red bars. ii) on the right panel, DEGs are distributed per biotype and the percentage of each biotype in the genome (i.e. in the whole set of features provided, for example, protein coding genes, non coding RNAs or pseudogenes) is presented in grey bars, the percentage of DEGs in each biotype is presented in blue lined bars and the percentage of each biotype in DEGs is presented in solid blue lines. The vertical green line separates the most abundant biotypes (on the left-hand side, corresponding to the left axis scale), from the rest (on the right-hand side, corresponding to the right axis scale). Otherwise, the lower abundance, biotypes would be indistinguishable.

Biotype distribution in DEG for the contrast Grp4aMinus_vs_Grp4aPlus

Biotype distribution in DEG for the contrast Grp4fMinus_vs_Grp4fPlus

Get all the figures in png, pdf format.

Results

DEG table for the contrast Grp4aMinus_vs_Grp4aPlus

The following table presents the top 10% statistically significant genes (use the download links below the table to retrieve the whole list) for the contrast Grp4aMinus_vs_Grp4aPlus. The fields of the table correspond to the requested features to be exported. The table can be searched using the search field on the top right.

chromosome	start	end	gene_id	gc_content	strand	gene_name	biotype	p-value_deseq	p-value_edger	p-value_limma	p-value_nbpseq	FDR_deseq	FDR_edger	FDR_limma	FDR_nbpseq	meta_p-value_Grp4aMinus_vs_Grp4aPlus	meta_FDR_Grp4aMinus_vs_Grp4aPlus	log2_normalized_fold_change_Grp4aMinus_vs_Grp4aPlus	log2_normalized_mean_counts_Grp4aMinus	rpgm_normalized_mean_counts_Grp4aMinus	log2_normalized_mean_counts_Grp4aPlus	rpgm_normalized_mean_counts_Grp4aPlus	log2_normalized_counts_YK3R4-5291-4aminus	log2_normalized_counts_YK3R5-5291-4aminus	rpgm_normalized_counts_YK3R4-5291-4aminus	rpgm_normalized_counts_YK3R5-5291-4aminus	log2_normalized_counts_YK3R1-5291-4aplus	log2_normalized_counts_YK3R2-5291-4aplus	log2_normalized_counts_YK3R3-5291-4aplus	rpgm_normalized_counts_YK3R1-5291-4aplus	rpgm_normalized_counts_YK3R2-5291-4aplus	rpgm_normalized_counts_YK3R3-5291-4aplus
chr14	68563388	68582086	ENSMUSG00000022057	0.3796	-	Adamdec1	protein_coding	3.14e-26	3.16e-05	2.63e-01	1.85e-28	4.30e-22	1.08e-01	9.46e-01	2.53e-24	2.64e-15	3.61e-11	-4.570868	6.14e+00	1.59e-01	3.83e+00	6.27e-03	2.81e+00	9.46e+00	2.69e-03	3.16e-01	4.46e+00	3.70e+00	3.32e+00	9.40e-03	5.37e-03	4.03e-03
chr15	65922443	65976804	ENSMUSG00000072511	0.4057	-	Hhla1	protein_coding	4.76e-12	2.88e-14	3.60e-04	4.29e-14	3.25e-08	3.94e-10	9.46e-01	1.96e-10	3.82e-11	2.61e-07	5.727920	0.00e+00	0.00e+00	5.69e+00	2.75e-02	0.00e+00	0.00e+00	0.00e+00	0.00e+00	5.95e+00	5.88e+00	5.25e+00	3.22e-02	3.06e-02	1.95e-02
chr7	46740497	46742980	ENSMUSG00000074115	0.5382	-	Saa1	protein_coding	4.23e-08	1.24e-04	3.65e-04	3.93e-14	7.23e-05	2.82e-01	9.46e-01	1.96e-10	9.31e-08	4.24e-04	2.578228	7.12e+00	1.14e-01	9.90e+00	6.83e-01	6.21e+00	8.02e+00	5.00e-02	1.78e-01	9.78e+00	9.44e+00	1.05e+01	6.02e-01	4.75e-01	9.71e-01
chr7	46751790	46754313	ENSMUSG00000057465	0.5452	+	Saa2	protein_coding	8.20e-08	3.04e-04	9.88e-05	1.17e-10	1.25e-04	4.18e-01	9.46e-01	2.28e-07	7.32e-07	2.14e-03	2.538612	5.78e+00	4.16e-02	8.46e+00	2.45e-01	5.09e+00	6.48e+00	2.27e-02	6.04e-02	8.15e+00	8.50e+00	8.73e+00	1.94e-01	2.49e-01	2.91e-01
chr7	30922372	30924681	ENSMUSG00000056978	0.5048	-	Hamp2	protein_coding	1.75e-10	7.82e-04	1.08e-02	2.51e-10	7.96e-07	5.35e-01	9.46e-01	4.29e-07	7.80e-07	2.14e-03	-2.195370	8.39e+00	3.08e-01	6.45e+00	6.68e-02	7.33e+00	9.45e+00	1.15e-01	5.01e-01	5.70e+00	6.66e+00	7.00e+00	3.67e-02	7.20e-02	9.15e-02
chr3	94693556	94704413	ENSMUSG00000068877	0.4995	+	Selenbp2	protein_coding	1.83e-06	2.28e-03	1.05e-03	5.11e-12	1.67e-03	6.07e-01	9.46e-01	1.40e-08	2.17e-06	4.95e-03	1.918881	9.76e+00	1.42e-01	1.17e+01	5.36e-01	1.02e+01	9.27e+00	1.88e-01	9.54e-02	1.19e+01	1.09e+01	1.21e+01	6.11e-01	2.99e-01	6.98e-01
chr2	35256380	35307892	ENSMUSG00000026879	0.5022	+	Gsn	protein_coding	5.60e-10	4.74e-04	7.23e-02	6.20e-09	1.92e-06	4.32e-01	9.46e-01	7.71e-06	3.30e-06	6.46e-03	-2.371789	6.59e+00	2.40e-02	4.72e+00	4.52e-03	5.04e+00	8.13e+00	4.93e-03	4.32e-02	5.04e+00	3.58e+00	5.52e+00	4.93e-03	1.70e-03	6.94e-03
chr9	123150946	123157432	ENSMUSG00000025784	0.5001	+	Clec3b	protein_coding	3.78e-09	1.56e-03	1.43e-01	4.75e-10	1.03e-05	5.75e-01	9.46e-01	7.22e-07	4.48e-06	7.66e-03	-3.954196	3.45e+00	1.95e-02	6.67e-01	6.51e-04	1.00e+00	5.91e+00	6.51e-04	3.84e-02	2.00e+00	0.00e+00	0.00e+00	1.95e-03	0.00e+00	0.00e+00
chr4	61435790	61439743	ENSMUSG00000096674	0.4087	-	Mup15	protein_coding	3.60e-05	2.64e-03	6.56e-03	2.52e-12	1.59e-02	6.29e-01	9.46e-01	8.62e-09	6.30e-06	9.45e-03	1.759479	1.10e+01	8.44e-01	1.26e+01	2.86e+00	1.15e+01	1.04e+01	1.15e+00	5.36e-01	1.26e+01	1.18e+01	1.35e+01	2.50e+00	1.40e+00	4.68e+00
chr4	60066469	60070475	ENSMUSG00000073842	0.4118	-	Mup7	protein_coding	1.00e-05	2.34e-03	7.08e-03	1.37e-11	6.52e-03	6.07e-01	9.46e-01	3.13e-08	6.90e-06	9.45e-03	1.826233	1.01e+01	4.26e-01	1.18e+01	1.51e+00	1.07e+01	9.50e+00	5.88e-01	2.63e-01	1.19e+01	1.09e+01	1.27e+01	1.42e+00	6.82e-01	2.43e+00
chr19	37697808	37701528	ENSMUSG00000024987	0.5278	+	Cyp26a1	protein_coding	9.21e-07	8.34e-06	4.42e-03	5.12e-07	9.70e-04	3.80e-02	9.46e-01	4.12e-04	1.15e-05	1.43e-02	-2.600559	5.37e+00	3.18e-02	2.72e+00	4.67e-03	4.64e+00	6.09e+00	1.68e-02	4.69e-02	1.58e+00	3.58e+00	3.00e+00	1.40e-03	7.70e-03	4.90e-03
chr14	55098578	55106593	ENSMUSG00000040701	0.5332	-	Ap1g2	protein_coding	9.20e-06	8.38e-07	1.64e-03	3.66e-06	6.30e-03	5.73e-03	9.46e-01	2.08e-03	1.47e-05	1.67e-02	3.906891	0.00e+00	0.00e+00	3.88e+00	7.96e-04	0.00e+00	0.00e+00	0.00e+00	0.00e+00	3.46e+00	4.00e+00	4.17e+00	5.68e-04	8.52e-04	9.66e-04
chr8	119992438	120052793	ENSMUSG00000031825	0.5098	+	Crispld2	protein_coding	2.43e-08	3.97e-03	6.62e-02	1.79e-08	4.75e-05	7.56e-01	9.46e-01	1.88e-05	1.84e-05	1.94e-02	-2.934112	5.02e+00	1.06e-02	2.77e+00	1.21e-03	3.46e+00	6.58e+00	2.02e-03	1.92e-02	2.81e+00	2.32e+00	3.17e+00	1.21e-03	8.06e-04	1.61e-03
chr4	61591929	61595871	ENSMUSG00000096688	0.4119	-	Mup17	protein_coding	1.59e-06	4.24e-03	4.17e-03	9.61e-09	1.56e-03	7.84e-01	9.46e-01	1.10e-05	2.28e-05	2.23e-02	1.835673	8.09e+00	1.24e-01	1.00e+01	4.43e-01	8.87e+00	7.31e+00	1.85e-01	6.26e-02	1.03e+01	9.17e+00	1.06e+01	5.03e-01	2.28e-01	5.98e-01
chr2	164777161	164779967	ENSMUSG00000017300	0.5426	-	Tnnc2	protein_coding	4.81e-09	5.63e-02	3.96e-01	3.98e-08	1.10e-05	1.00e+00	9.48e-01	3.89e-05	4.54e-05	4.15e-02	4.459432	0.00e+00	0.00e+00	2.00e+00	7.94e-03	0.00e+00	0.00e+00	0.00e+00	0.00e+00	6.00e+00	0.00e+00	0.00e+00	2.38e-02	0.00e+00	0.00e+00
chr4	60658466	60662411	ENSMUSG00000073834	0.4110	-	Mup11	protein_coding	1.39e-04	1.13e-02	1.03e-02	5.64e-10	4.75e-02	1.00e+00	9.46e-01	7.72e-07	5.49e-05	4.69e-02	1.538727	1.11e+01	5.86e-01	1.27e+01	1.70e+00	1.18e+01	1.03e+01	8.76e-01	2.96e-01	1.26e+01	1.20e+01	1.35e+01	1.43e+00	9.66e-01	2.71e+00
chr8	94172664	94173568	ENSMUSG00000031762	0.5569	+	Mt2	protein_coding	1.59e-07	6.71e-03	3.71e-02	2.49e-06	2.18e-04	8.90e-01	9.46e-01	1.48e-03	9.97e-05	7.94e-02	-1.682030	7.93e+00	2.09e-01	6.06e+00	6.46e-02	7.04e+00	8.82e+00	9.40e-02	3.24e-01	4.39e+00	6.41e+00	7.38e+00	1.44e-02	6.03e-02	1.19e-01
chr3	18071950	18243338	ENSMUSG00000039519	0.3883	-	Cyp7b1	protein_coding	2.51e-06	8.52e-03	7.62e-03	7.30e-07	2.14e-03	9.14e-01	9.46e-01	5.26e-04	1.04e-04	7.94e-02	1.640262	7.66e+00	1.44e-01	9.55e+00	4.50e-01	8.68e+00	6.64e+00	2.32e-01	5.63e-02	1.02e+01	9.29e+00	9.13e+00	6.77e-01	3.55e-01	3.19e-01
chr4	115299046	115332815	ENSMUSG00000066071	0.4031	+	Cyp4a12a	protein_coding	3.54e-04	4.57e-04	9.31e-04	2.08e-06	9.31e-02	4.32e-01	9.46e-01	1.30e-03	1.33e-04	9.23e-02	1.353231	9.67e+00	4.11e-01	1.11e+01	1.05e+00	1.01e+01	9.27e+00	5.24e-01	2.98e-01	1.13e+01	1.07e+01	1.12e+01	1.20e+00	8.32e-01	1.12e+00
chr7	123462291	123468004	ENSMUSG00000030762	0.5004	+	Aqp8	protein_coding	1.99e-07	2.61e-02	5.11e-02	1.25e-06	2.47e-04	1.00e+00	9.46e-01	8.54e-04	1.35e-04	9.23e-02	-1.653018	8.27e+00	2.23e-01	6.59e+00	7.04e-02	7.18e+00	9.36e+00	8.00e-02	3.65e-01	5.93e+00	5.78e+00	8.06e+00	3.33e-02	3.00e-02	1.48e-01
chr15	85205949	85286535	ENSMUSG00000006369	0.5036	+	Fbln1	protein_coding	3.83e-07	4.88e-02	5.20e-01	1.16e-07	4.37e-04	1.00e+00	9.55e-01	1.06e-04	1.83e-04	1.17e-01	-3.102811	2.98e+00	5.87e-03	1.46e+00	5.13e-04	0.00e+00	5.95e+00	0.00e+00	1.17e-02	2.81e+00	0.00e+00	1.58e+00	1.15e-03	0.00e+00	3.85e-04
chr4	58435255	58553898	ENSMUSG00000038668	0.4240	-	Lpar1	protein_coding	4.37e-06	1.87e-03	8.77e-02	1.90e-06	3.52e-03	6.05e-01	9.46e-01	1.24e-03	1.92e-04	1.17e-01	-2.731804	4.20e+00	6.42e-03	1.64e+00	7.85e-04	2.58e+00	5.81e+00	1.07e-03	1.18e-02	1.58e+00	3.32e+00	0.00e+00	4.28e-04	1.93e-03	0.00e+00
chr10	99263231	99267489	ENSMUSG00000019960	0.4750	+	Dusp6	protein_coding	2.88e-04	1.07e-04	6.91e-04	6.95e-05	8.57e-02	2.82e-01	9.46e-01	2.38e-02	1.96e-04	1.17e-01	1.416168	7.31e+00	4.29e-02	8.71e+00	1.15e-01	7.48e+00	7.13e+00	4.82e-02	3.76e-02	8.86e+00	8.90e+00	8.38e+00	1.26e-01	1.29e-01	8.96e-02
chr15	6445327	6498751	ENSMUSG00000022149	0.4003	+	C9	protein_coding	9.39e-04	9.53e-04	2.12e-02	2.60e-07	2.07e-01	5.67e-01	9.46e-01	2.22e-04	2.65e-04	1.51e-01	1.441347	9.91e+00	1.00e-01	1.12e+01	2.72e-01	9.75e+00	1.01e+01	8.95e-02	1.11e-01	1.03e+01	1.21e+01	1.11e+01	1.31e-01	4.52e-01	2.34e-01
chr17	75005568	75392512	ENSMUSG00000001870	0.4266	+	Ltbp1	protein_coding	2.94e-03	1.72e-04	3.19e-03	2.28e-05	4.73e-01	3.36e-01	9.46e-01	1.11e-02	4.38e-04	2.35e-01	-2.857981	3.84e+00	1.81e-03	8.62e-01	1.34e-04	3.58e+00	4.09e+00	1.47e-03	2.14e-03	1.58e+00	0.00e+00	1.00e+00	2.68e-04	0.00e+00	1.34e-04
chr5	137061504	137072268	ENSMUSG00000037411	0.5016	-	Serpine1	protein_coding	4.82e-05	2.79e-03	1.41e-02	2.09e-05	2.00e-02	6.29e-01	9.46e-01	1.06e-02	4.46e-04	2.35e-01	3.222392	1.00e+00	2.40e-04	4.05e+00	4.23e-03	1.00e+00	1.00e+00	2.40e-04	2.40e-04	4.09e+00	4.91e+00	3.17e+00	3.83e-03	6.95e-03	1.92e-03
chr4	115486200	115496142	ENSMUSG00000028715	0.4146	-	Cyp4a14	protein_coding	9.72e-05	4.21e-02	1.94e-02	7.21e-07	3.84e-02	1.00e+00	9.46e-01	5.26e-04	4.89e-04	2.46e-01	-1.568915	9.33e+00	1.94e-01	7.01e+00	6.53e-02	9.74e+00	8.92e+00	2.48e-01	1.41e-01	6.60e+00	9.08e+00	5.36e+00	2.79e-02	1.56e-01	1.16e-02
chr12	110690605	110702728	ENSMUSG00000021270	0.4790	-	Hsp90aa1	protein_coding	3.33e-05	3.28e-03	3.63e-02	1.61e-05	1.52e-02	6.71e-01	9.46e-01	8.67e-03	5.03e-04	2.46e-01	1.567989	7.05e+00	1.56e-02	8.40e+00	4.65e-02	7.18e+00	6.92e+00	1.70e-02	1.42e-02	9.36e+00	8.51e+00	7.31e+00	7.78e-02	4.30e-02	1.87e-02
chr6	55170626	55175043	ENSMUSG00000003477	0.5124	-	Inmt	protein_coding	2.48e-03	1.60e-03	4.09e-04	9.52e-05	4.19e-01	5.75e-01	9.46e-01	3.18e-02	6.26e-04	2.96e-01	-1.114222	1.03e+01	4.68e-01	9.20e+00	2.16e-01	1.03e+01	1.04e+01	4.48e-01	4.88e-01	9.11e+00	9.26e+00	9.24e+00	2.03e-01	2.24e-01	2.21e-01
chr8	94179082	94180327	ENSMUSG00000031765	0.5425	+	Mt1	protein_coding	1.53e-05	7.87e-03	4.82e-02	3.90e-05	8.98e-03	9.04e-01	9.46e-01	1.64e-02	6.90e-04	3.15e-01	-1.336142	8.78e+00	2.36e-01	7.40e+00	9.33e-02	7.97e+00	9.59e+00	1.16e-01	3.57e-01	6.41e+00	7.29e+00	8.51e+00	3.90e-02	7.24e-02	1.68e-01
chr16	4941415	4964237	ENSMUSG00000022515	0.4782	-	Anks3	protein_coding	5.34e-04	3.66e-04	1.88e-03	1.01e-03	1.33e-01	4.18e-01	9.46e-01	1.65e-01	7.80e-04	3.44e-01	1.747613	3.73e+00	4.72e-03	5.49e+00	1.68e-02	4.00e+00	3.46e+00	5.67e-03	3.78e-03	5.70e+00	5.49e+00	5.29e+00	1.93e-02	1.66e-02	1.44e-02
chr15	74994877	74998031	ENSMUSG00000075602	0.4887	-	Ly6a	protein_coding	8.38e-04	4.42e-04	1.13e-02	1.11e-04	1.94e-01	4.32e-01	9.46e-01	3.45e-02	8.27e-04	3.53e-01	1.545185	6.10e+00	2.44e-02	7.59e+00	7.19e-02	5.67e+00	6.52e+00	1.73e-02	3.15e-02	7.98e+00	6.71e+00	8.07e+00	8.69e-02	3.60e-02	9.28e-02
chr5	35652023	35679782	ENSMUSG00000029096	0.5294	-	Htra3	protein_coding	1.69e-05	1.47e-02	2.44e-01	1.65e-05	9.27e-03	1.00e+00	9.46e-01	8.67e-03	1.00e-03	4.15e-01	-2.689299	3.18e+00	1.85e-03	1.53e+00	2.10e-04	1.00e+00	5.36e+00	9.02e-05	3.61e-03	1.00e+00	1.00e+00	2.58e+00	9.02e-05	9.02e-05	4.51e-04
chr15	54571366	54602846	ENSMUSG00000024479	0.4201	+	Mal2	protein_coding	3.33e-04	2.20e-03	1.09e-02	2.73e-04	9.10e-02	6.07e-01	9.46e-01	6.43e-02	1.22e-03	4.89e-01	1.359573	6.85e+00	3.43e-02	8.12e+00	8.84e-02	7.01e+00	6.69e+00	3.81e-02	3.04e-02	8.78e+00	8.10e+00	7.48e+00	1.31e-01	8.13e-02	5.30e-02
chr1	90765923	90843971	ENSMUSG00000048126	0.4734	-	Col6a3	protein_coding	2.60e-05	8.90e-03	1.17e-01	1.20e-04	1.23e-02	9.26e-01	9.46e-01	3.57e-02	1.34e-03	5.19e-01	-1.839776	5.16e+00	6.25e-03	3.84e+00	1.66e-03	3.81e+00	6.51e+00	1.58e-03	1.09e-02	3.46e+00	4.25e+00	3.81e+00	1.21e-03	2.18e-03	1.58e-03
chr19	11536801	11558467	ENSMUSG00000024678	0.3944	+	Ms4a4d	protein_coding	1.61e-03	3.65e-04	5.03e-03	1.18e-03	3.08e-01	4.18e-01	9.46e-01	1.70e-01	1.37e-03	5.19e-01	-1.736966	4.85e+00	1.70e-02	3.16e+00	4.70e-03	4.46e+00	5.25e+00	1.23e-02	2.17e-02	3.17e+00	3.32e+00	3.00e+00	4.70e-03	5.29e-03	4.11e-03
chr5	107908053	107987085	ENSMUSG00000029270	0.4204	-	Fam69a	protein_coding	8.59e-04	9.04e-04	5.03e-03	1.10e-03	1.96e-01	5.62e-01	9.46e-01	1.70e-01	1.44e-03	5.22e-01	1.748021	3.53e+00	2.35e-03	5.42e+00	8.35e-03	4.25e+00	2.81e+00	3.52e-03	1.17e-03	5.58e+00	5.04e+00	5.64e+00	9.20e-03	6.27e-03	9.59e-03
chr5	72647795	72671078	ENSMUSG00000067219	0.4310	+	Nipal1	protein_coding	2.47e-05	1.21e-02	2.74e-01	6.14e-05	1.23e-02	1.00e+00	9.46e-01	2.15e-02	1.50e-03	5.22e-01	-2.192645	4.10e+00	4.68e-03	2.64e+00	9.07e-04	2.32e+00	5.88e+00	6.04e-04	8.76e-03	1.58e+00	3.17e+00	3.17e+00	3.02e-04	1.21e-03	1.21e-03
chr2	84936579	84958509	ENSMUSG00000027074	0.4679	+	Slc43a3	protein_coding	1.23e-03	1.05e-02	7.10e-04	6.35e-04	2.59e-01	9.73e-01	9.46e-01	1.22e-01	1.55e-03	5.22e-01	-1.260180	7.65e+00	2.82e-02	6.38e+00	1.17e-02	7.77e+00	7.54e+00	3.04e-02	2.59e-02	6.09e+00	6.57e+00	6.49e+00	9.39e-03	1.32e-02	1.25e-02
chr15	91145884	91191807	ENSMUSG00000055782	0.3861	-	Abcd2	protein_coding	1.64e-04	1.48e-02	1.59e-02	2.14e-04	5.49e-02	1.00e+00	9.46e-01	5.32e-02	1.70e-03	5.22e-01	-2.070389	4.60e+00	6.20e-03	2.71e+00	1.30e-03	3.81e+00	5.39e+00	2.99e-03	9.42e-03	2.81e+00	2.32e+00	3.00e+00	1.38e-03	9.19e-04	1.61e-03
chr8	123891769	123894751	ENSMUSG00000031972	0.5675	-	Acta1	protein_coding	6.39e-06	8.00e-02	3.95e-01	4.54e-05	4.86e-03	1.00e+00	9.48e-01	1.83e-02	1.74e-03	5.22e-01	3.502500	0.00e+00	0.00e+00	1.67e+00	3.77e-03	0.00e+00	0.00e+00	0.00e+00	0.00e+00	5.00e+00	0.00e+00	0.00e+00	1.13e-02	0.00e+00	0.00e+00
chr17	25718926	25727419	ENSMUSG00000019214	0.5615	-	Chtf18	protein_coding	1.67e-02	1.28e-03	3.51e-03	1.28e-04	1.00e+00	5.75e-01	9.46e-01	3.64e-02	1.76e-03	5.22e-01	-2.906891	2.82e+00	3.60e-04	0.00e+00	0.00e+00	3.32e+00	2.32e+00	4.99e-04	2.22e-04	0.00e+00	0.00e+00	0.00e+00	0.00e+00	0.00e+00	0.00e+00
chr1	164796644	164824266	ENSMUSG00000026574	0.4370	+	Dpt	protein_coding	1.50e-05	2.77e-02	6.94e-01	3.53e-05	8.98e-03	1.00e+00	9.71e-01	1.61e-02	1.78e-03	5.22e-01	-2.189825	3.57e+00	1.56e-02	2.99e+00	3.07e-03	1.00e+00	6.15e+00	4.38e-04	3.07e-02	3.17e+00	3.00e+00	2.81e+00	3.51e-03	3.07e-03	2.63e-03
chr8	70072924	70080066	ENSMUSG00000036151	0.5149	+	Tm6sf2	protein_coding	2.12e-03	3.25e-04	4.86e-03	3.25e-03	3.72e-01	4.18e-01	9.46e-01	3.18e-01	1.82e-03	5.22e-01	-1.325875	5.89e+00	1.14e-02	4.60e+00	4.45e-03	5.52e+00	6.25e+00	8.58e-03	1.43e-02	4.46e+00	4.58e+00	4.75e+00	4.00e-03	4.38e-03	4.96e-03
chr7	128523616	128546981	ENSMUSG00000030847	0.4872	+	Bag3	protein_coding	1.15e-04	6.46e-03	1.40e-01	1.09e-04	4.18e-02	8.87e-01	9.46e-01	3.45e-02	1.83e-03	5.22e-01	1.395342	7.18e+00	7.97e-02	8.26e+00	2.11e-01	7.43e+00	6.94e+00	9.31e-02	6.64e-02	9.40e+00	8.55e+00	6.85e+00	3.67e-01	2.03e-01	6.21e-02
chr7	110768206	110812405	ENSMUSG00000005686	0.5018	+	Ampd3	protein_coding	2.27e-02	1.87e-03	2.11e-03	1.28e-04	1.00e+00	6.05e-01	9.46e-01	3.64e-02	1.84e-03	5.22e-01	-2.906891	2.90e+00	9.81e-04	0.00e+00	0.00e+00	2.81e+00	3.00e+00	9.06e-04	1.06e-03	0.00e+00	0.00e+00	0.00e+00	0.00e+00	0.00e+00	0.00e+00
chr6	11900292	11907497	ENSMUSG00000029632	0.3973	-	Ndufa4	protein_coding	1.68e-03	2.35e-03	1.96e-03	1.60e-03	3.08e-01	6.07e-01	9.46e-01	2.03e-01	1.88e-03	5.22e-01	1.189767	6.81e+00	7.63e-02	8.03e+00	1.75e-01	7.14e+00	6.48e+00	9.36e-02	5.89e-02	8.23e+00	7.95e+00	7.91e+00	2.00e-01	1.65e-01	1.60e-01
chr3	90603771	90606045	ENSMUSG00000001020	0.4853	+	S100a4	protein_coding	2.57e-05	7.79e-02	1.05e-01	6.13e-05	1.23e-02	1.00e+00	9.46e-01	2.15e-02	1.89e-03	5.22e-01	2.309855	3.00e+00	4.78e-03	4.81e+00	2.64e-02	3.00e+00	3.00e+00	4.78e-03	4.78e-03	3.46e+00	6.39e+00	4.58e+00	6.83e-03	5.67e-02	1.57e-02
chr3	90601904	90605885	ENSMUSG00000105518	0.5080	+	Gm42674	protein_coding	2.57e-05	7.85e-02	1.05e-01	6.13e-05	1.23e-02	1.00e+00	9.46e-01	2.15e-02	1.90e-03	5.22e-01	2.309855	3.00e+00	2.80e-03	4.81e+00	1.54e-02	3.00e+00	3.00e+00	2.80e-03	2.80e-03	3.46e+00	6.39e+00	4.58e+00	3.99e-03	3.31e-02	9.19e-03
chr15	10177238	10349180	ENSMUSG00000005268	0.4018	+	Prlr	protein_coding	5.75e-04	8.48e-04	1.89e-02	1.44e-03	1.40e-01	5.52e-01	9.46e-01	1.86e-01	1.91e-03	5.22e-01	-1.183712	7.41e+00	9.34e-03	6.30e+00	4.08e-03	6.85e+00	7.98e+00	5.82e-03	1.29e-02	5.81e+00	6.44e+00	6.64e+00	2.81e-03	4.39e-03	5.05e-03
chr15	66929321	66969640	ENSMUSG00000005125	0.4895	-	Ndrg1	protein_coding	3.95e-03	2.97e-04	3.36e-03	4.61e-03	5.31e-01	4.18e-01	9.46e-01	3.77e-01	2.06e-03	5.33e-01	-1.249778	6.02e+00	5.11e-03	4.72e+00	2.10e-03	5.95e+00	6.09e+00	4.87e-03	5.35e-03	4.39e+00	4.52e+00	5.25e+00	1.60e-03	1.76e-03	2.96e-03
chr4	62147932	62150863	ENSMUSG00000066153	0.3997	-	Mup21	protein_coding	4.50e-04	2.41e-02	9.02e-03	1.95e-04	1.16e-01	1.00e+00	9.46e-01	4.94e-02	2.09e-03	5.33e-01	1.213008	8.19e+00	2.50e-01	9.49e+00	5.81e-01	8.81e+00	7.57e+00	3.52e-01	1.49e-01	9.98e+00	9.28e+00	9.20e+00	7.93e-01	4.88e-01	4.63e-01
chr11	54902453	54910377	ENSMUSG00000018339	0.5380	+	Gpx3	protein_coding	4.05e-05	4.15e-02	2.13e-01	5.43e-05	1.73e-02	1.00e+00	9.46e-01	2.06e-02	2.10e-03	5.33e-01	-2.169925	4.25e+00	1.53e-02	2.87e+00	3.02e-03	2.58e+00	5.91e+00	2.39e-03	2.82e-02	3.00e+00	2.81e+00	2.81e+00	3.34e-03	2.86e-03	2.86e-03
chr6	136872435	136875823	ENSMUSG00000030218	0.4402	-	Mgp	protein_coding	7.58e-06	6.82e-02	2.61e-01	1.45e-04	5.46e-03	1.00e+00	9.46e-01	3.89e-02	2.10e-03	5.33e-01	-1.782902	5.03e+00	2.12e-02	3.71e+00	5.90e-03	3.32e+00	6.74e+00	3.32e-03	3.91e-02	3.00e+00	5.13e+00	3.00e+00	2.58e-03	1.25e-02	2.58e-03
chr15	34306680	34436242	ENSMUSG00000022324	0.4365	+	Matn2	protein_coding	2.31e-04	7.74e-03	2.33e-01	5.04e-05	7.01e-02	9.04e-01	9.46e-01	1.97e-02	2.14e-03	5.33e-01	-2.629357	2.98e+00	9.08e-03	1.33e+00	9.76e-04	1.00e+00	4.95e+00	5.86e-04	1.76e-02	2.00e+00	1.00e+00	1.00e+00	1.76e-03	5.86e-04	5.86e-04
chr3	66219910	66225805	ENSMUSG00000027832	0.4383	+	Ptx3	protein_coding	1.94e-04	3.09e-02	1.30e-01	3.39e-05	6.02e-02	1.00e+00	9.46e-01	1.60e-02	2.27e-03	5.54e-01	3.544321	0.00e+00	0.00e+00	2.50e+00	3.44e-03	0.00e+00	0.00e+00	0.00e+00	0.00e+00	1.00e+00	4.91e+00	1.58e+00	3.22e-04	9.34e-03	6.44e-04
chr2	119972699	120027408	ENSMUSG00000033902	0.4729	+	Mapkbp1	protein_coding	3.10e-03	1.40e-03	6.60e-03	1.18e-03	4.79e-01	5.75e-01	9.46e-01	1.70e-01	2.41e-03	5.79e-01	-1.817623	4.54e+00	3.18e-03	2.43e+00	8.00e-04	4.75e+00	4.32e+00	3.67e-03	2.68e-03	2.00e+00	1.58e+00	3.70e+00	4.23e-04	2.82e-04	1.69e-03
chr8	85083270	85098282	ENSMUSG00000005142	0.5026	+	Man2b1	protein_coding	1.28e-03	1.55e-03	8.39e-03	3.02e-03	2.65e-01	5.75e-01	9.46e-01	3.06e-01	2.66e-03	6.28e-01	-1.314605	5.97e+00	5.75e-03	4.68e+00	2.26e-03	5.55e+00	6.39e+00	4.10e-03	7.40e-03	4.25e+00	4.70e+00	5.09e+00	1.61e-03	2.23e-03	2.94e-03

Download the DEG result list for Grp4aMinus_vs_Grp4aPlus.
Download the whole result list for Grp4aMinus_vs_Grp4aPlus.

DEG table for the contrast Grp4fMinus_vs_Grp4fPlus

The following table presents the top 10% statistically significant genes (use the download links below the table to retrieve the whole list) for the contrast Grp4fMinus_vs_Grp4fPlus. The fields of the table correspond to the requested features to be exported. The table can be searched using the search field on the top right.

chromosome	start	end	gene_id	gc_content	strand	gene_name	biotype	p-value_deseq	p-value_edger	p-value_limma	p-value_nbpseq	FDR_deseq	FDR_edger	FDR_limma	FDR_nbpseq	meta_p-value_Grp4fMinus_vs_Grp4fPlus	meta_FDR_Grp4fMinus_vs_Grp4fPlus	log2_normalized_fold_change_Grp4fMinus_vs_Grp4fPlus	log2_normalized_mean_counts_Grp4fMinus	rpgm_normalized_mean_counts_Grp4fMinus	log2_normalized_mean_counts_Grp4fPlus	rpgm_normalized_mean_counts_Grp4fPlus	log2_normalized_counts_YK3R9-5269-4fminus	log2_normalized_counts_YK3R10-5269-4fminus	rpgm_normalized_counts_YK3R9-5269-4fminus	rpgm_normalized_counts_YK3R10-5269-4fminus	log2_normalized_counts_YK3R6-5269-4fplus	log2_normalized_counts_YK3R7-5269-4fplus	log2_normalized_counts_YK3R8-5269-4fplus	rpgm_normalized_counts_YK3R6-5269-4fplus	rpgm_normalized_counts_YK3R7-5269-4fplus	rpgm_normalized_counts_YK3R8-5269-4fplus
chr7	131032053	131121630	ENSMUSG00000047517	0.4550	+	Dmbt1	protein_coding	9.02e-19	9.64e-02	5.47e-01	1.93e-17	1.23e-14	7.38e-01	8.73e-01	1.73e-13	9.79e-10	1.34e-05	-3.636625	4.32e+00	4.83e-02	2.16e+00	3.66e-03	8.63e+00	0.00e+00	9.65e-02	0.00e+00	1.00e+00	0.00e+00	5.49e+00	2.44e-04	0.00e+00	1.07e-02
chr3	90612882	90624181	ENSMUSG00000001025	0.4742	+	S100a6	protein_coding	9.28e-10	2.39e-05	3.39e-02	2.53e-17	2.54e-06	3.37e-02	6.54e-01	1.73e-13	1.18e-08	8.03e-05	-2.596740	9.66e+00	2.86e-01	7.52e+00	4.70e-02	8.24e+00	1.11e+01	7.02e-02	5.01e-01	7.51e+00	8.34e+00	6.70e+00	4.21e-02	7.51e-02	2.39e-02
chr3	90603771	90606045	ENSMUSG00000001020	0.4853	+	S100a4	protein_coding	1.11e-12	2.97e-03	2.82e-01	1.20e-16	6.73e-09	2.75e-01	7.62e-01	5.46e-13	1.82e-08	8.03e-05	-3.346423	6.05e+00	1.68e-01	3.69e+00	1.59e-02	3.17e+00	8.92e+00	5.46e-03	3.31e-01	2.81e+00	5.93e+00	2.32e+00	4.10e-03	4.10e-02	2.73e-03
chr3	90601904	90605885	ENSMUSG00000105518	0.5080	+	Gm42674	protein_coding	1.47e-12	3.18e-03	2.83e-01	2.28e-16	6.73e-09	2.78e-01	7.62e-01	7.82e-13	2.35e-08	8.03e-05	-3.307429	6.05e+00	9.84e-02	3.70e+00	9.58e-03	3.17e+00	8.92e+00	3.19e-03	1.94e-01	2.81e+00	5.98e+00	2.32e+00	2.40e-03	2.48e-02	1.60e-03
chr18	52516069	52529867	ENSMUSG00000024529	0.4156	-	Lox	protein_coding	4.40e-10	7.03e-04	4.50e-01	3.56e-14	1.51e-06	1.51e-01	8.34e-01	9.74e-11	2.65e-07	7.26e-04	-4.405141	3.40e+00	1.26e-02	1.19e+00	3.78e-04	0.00e+00	6.81e+00	0.00e+00	2.52e-02	1.58e+00	2.00e+00	0.00e+00	4.54e-04	6.80e-04	0.00e+00
chr16	20668313	20692884	ENSMUSG00000045983	0.4855	+	Eif4g1	protein_coding	5.98e-06	3.98e-06	9.28e-04	9.68e-12	6.82e-03	1.21e-02	6.54e-01	2.21e-08	6.80e-07	1.55e-03	-1.886460	1.10e+01	1.12e-01	9.10e+00	3.03e-02	1.07e+01	1.14e+01	8.56e-02	1.39e-01	9.68e+00	8.89e+00	8.73e+00	4.35e-02	2.50e-02	2.24e-02
chr4	61857321	61857876	ENSMUSG00000081053	0.4442	-	Mup-ps17	unprocessed_pseudogene	2.46e-07	2.27e-07	1.96e-03	1.97e-07	4.81e-04	1.56e-03	6.54e-01	1.68e-04	2.16e-06	4.21e-03	-2.362570	6.30e+00	6.63e-02	3.90e+00	1.22e-02	6.82e+00	5.78e+00	8.95e-02	4.31e-02	4.25e+00	4.46e+00	3.00e+00	1.44e-02	1.68e-02	5.59e-03
chr10	126893113	126893908	ENSMUSG00000112423	0.5201	+	AC122018.2	processed_pseudogene	6.69e-07	5.27e-07	2.89e-03	9.43e-08	1.02e-03	2.40e-03	6.54e-01	9.93e-05	3.13e-06	5.19e-03	-3.273018	4.80e+00	2.24e-02	9.36e-01	1.60e-03	5.21e+00	4.39e+00	2.88e-02	1.60e-02	2.81e+00	0.00e+00	0.00e+00	4.79e-03	0.00e+00	0.00e+00
chr10	59949669	59951834	ENSMUSG00000020108	0.5803	-	Ddit4	protein_coding	3.63e-09	5.21e-04	2.80e-02	2.56e-09	8.29e-06	1.46e-01	6.54e-01	3.50e-06	3.41e-06	5.19e-03	-2.758992	5.96e+00	4.11e-02	3.64e+00	5.66e-03	7.23e+00	4.70e+00	7.03e-02	1.18e-02	4.00e+00	3.00e+00	3.91e+00	7.08e-03	3.30e-03	6.61e-03
chr7	46751790	46754313	ENSMUSG00000057465	0.5452	+	Saa2	protein_coding	4.83e-07	2.46e-04	6.87e-02	1.02e-09	8.27e-04	1.08e-01	6.54e-01	1.75e-06	9.56e-06	1.30e-02	2.554102	5.34e+00	2.82e-02	7.39e+00	1.69e-01	5.73e+00	4.95e+00	3.57e-02	2.06e-02	5.64e+00	7.55e+00	8.97e+00	3.37e-02	1.28e-01	3.45e-01
chr5	139150223	139186510	ENSMUSG00000025857	0.5066	+	Dnaaf5	protein_coding	3.36e-05	1.17e-07	4.94e-04	6.10e-06	3.06e-02	1.56e-03	6.54e-01	3.34e-03	1.04e-05	1.30e-02	-2.030762	6.30e+00	2.32e-02	4.21e+00	5.44e-03	6.21e+00	6.39e+00	2.17e-02	2.46e-02	4.52e+00	4.52e+00	3.58e+00	6.53e-03	6.53e-03	3.27e-03
chr19	36111965	36119844	ENSMUSG00000024803	0.4397	-	Ankrd1	protein_coding	1.16e-06	7.28e-04	2.81e-01	2.09e-09	1.59e-03	1.51e-01	7.62e-01	3.18e-06	2.65e-05	3.03e-02	-4.108524	2.75e+00	1.10e-02	3.33e-01	1.67e-04	0.00e+00	5.49e+00	0.00e+00	2.21e-02	1.00e+00	0.00e+00	0.00e+00	5.02e-04	0.00e+00	0.00e+00
chr19	41482645	41562246	ENSMUSG00000025019	0.3869	+	Lcor	protein_coding	2.63e-05	2.22e-05	4.38e-03	4.24e-07	2.57e-02	3.37e-02	6.54e-01	3.05e-04	3.23e-05	3.26e-02	-1.833210	8.12e+00	6.65e-02	6.10e+00	1.85e-02	7.91e+00	8.34e+00	5.67e-02	7.63e-02	6.97e+00	6.30e+00	5.04e+00	2.94e-02	1.85e-02	7.59e-03
chr17	33996017	34000333	ENSMUSG00000061232	0.5383	-	H2-K1	protein_coding	8.62e-04	2.78e-05	6.11e-03	8.42e-09	2.68e-01	3.37e-02	6.54e-01	1.05e-05	3.33e-05	3.26e-02	1.539726	1.04e+01	4.20e-01	1.18e+01	1.22e+00	1.06e+01	1.02e+01	4.74e-01	3.67e-01	1.11e+01	1.20e+01	1.24e+01	6.73e-01	1.27e+00	1.72e+00
chr7	30942368	30944032	ENSMUSG00000050440	0.4913	-	Hamp	protein_coding	7.96e-05	1.82e-04	1.24e-03	3.48e-07	5.73e-02	9.23e-02	6.54e-01	2.65e-04	5.00e-05	4.56e-02	1.723977	7.73e+00	1.55e-01	9.39e+00	5.14e-01	7.78e+00	7.67e+00	1.61e-01	1.49e-01	8.75e+00	9.65e+00	9.76e+00	3.17e-01	5.90e-01	6.35e-01
chr4	61670177	61674136	ENSMUSG00000078674	0.4104	-	Mup18	protein_coding	1.90e-04	2.89e-03	3.98e-02	4.98e-10	1.04e-01	2.73e-01	6.54e-01	9.74e-07	5.75e-05	4.91e-02	1.821052	9.08e+00	2.46e-01	1.08e+01	8.70e-01	9.85e+00	8.32e+00	3.66e-01	1.26e-01	9.43e+00	1.15e+01	1.15e+01	2.74e-01	1.19e+00	1.15e+00
chr6	84571414	84593908	ENSMUSG00000063415	0.5464	-	Cyp26b1	protein_coding	5.84e-05	5.57e-04	1.07e-04	6.43e-06	4.44e-02	1.46e-01	6.54e-01	3.39e-03	6.87e-05	5.53e-02	-1.537417	8.60e+00	9.43e-02	7.06e+00	3.23e-02	8.65e+00	8.55e+00	9.77e-02	9.09e-02	7.07e+00	6.95e+00	7.16e+00	3.24e-02	3.00e-02	3.46e-02
chr17	34112023	34122634	ENSMUSG00000024335	0.5520	-	Brd2	protein_coding	6.25e-04	4.44e-06	1.32e-03	3.20e-05	2.25e-01	1.21e-02	6.54e-01	1.18e-02	1.04e-04	7.66e-02	-1.383257	8.77e+00	2.94e-02	7.34e+00	1.12e-02	8.74e+00	8.80e+00	2.88e-02	2.99e-02	7.83e+00	7.06e+00	7.14e+00	1.53e-02	8.90e-03	9.44e-03
chr3	152153447	152154116	ENSMUSG00000084343	0.5194	+	Gm16198	processed_pseudogene	4.29e-05	1.79e-05	8.50e-03	2.12e-05	3.67e-02	3.37e-02	6.54e-01	9.05e-03	1.08e-04	7.66e-02	-2.963474	3.53e+00	9.58e-03	5.28e-01	5.32e-04	4.25e+00	2.81e+00	1.44e-02	4.79e-03	1.58e+00	0.00e+00	0.00e+00	1.60e-03	0.00e+00	0.00e+00
chr10	91082940	91102607	ENSMUSG00000019975	0.4197	+	Ikbip	protein_coding	2.86e-06	1.89e-03	2.44e-01	1.19e-07	3.56e-03	2.31e-01	7.40e-01	1.09e-04	1.12e-04	7.66e-02	3.017074	2.75e+00	7.14e-04	4.82e+00	6.62e-03	3.17e+00	2.32e+00	9.52e-04	4.76e-04	2.81e+00	7.13e+00	4.52e+00	7.14e-04	1.65e-02	2.62e-03
chr9	74848463	74891476	ENSMUSG00000093662	0.4696	+	Gm20649	lincRNA	1.85e-05	5.97e-04	2.04e-02	2.27e-06	1.94e-02	1.46e-01	6.54e-01	1.55e-03	1.50e-04	9.80e-02	-2.423682	5.40e+00	3.63e-02	2.36e+00	6.12e-03	4.75e+00	6.04e+00	2.08e-02	5.19e-02	3.00e+00	4.09e+00	0.00e+00	5.59e-03	1.28e-02	0.00e+00
chr1	153740349	153745468	ENSMUSG00000026475	0.5045	+	Rgs16	protein_coding	1.82e-04	1.47e-03	3.50e-03	4.01e-06	1.03e-01	2.19e-01	6.54e-01	2.29e-03	2.47e-04	1.54e-01	-1.653801	8.08e+00	9.70e-02	6.37e+00	3.06e-02	7.72e+00	8.44e+00	7.33e-02	1.21e-01	7.11e+00	6.04e+00	5.95e+00	4.78e-02	2.27e-02	2.13e-02
chr14	78492246	78536860	ENSMUSG00000022016	0.3949	-	Akap11	protein_coding	4.61e-04	1.25e-04	3.43e-03	2.36e-05	1.97e-01	9.22e-02	6.54e-01	9.49e-03	2.61e-04	1.55e-01	-1.510962	8.10e+00	5.74e-02	6.48e+00	2.00e-02	8.01e+00	8.20e+00	5.35e-02	6.13e-02	7.24e+00	6.32e+00	5.88e+00	3.14e-02	1.65e-02	1.21e-02
chr12	109541001	109571726	ENSMUSG00000021268	0.5022	+	Meg3	lincRNA	5.38e-05	4.57e-04	2.61e-03	1.03e-04	4.33e-02	1.32e-01	6.54e-01	2.92e-02	2.85e-04	1.57e-01	-1.988318	5.35e+00	3.99e-03	3.14e+00	9.31e-04	5.25e+00	5.46e+00	3.69e-03	4.29e-03	2.58e+00	2.58e+00	4.25e+00	4.99e-04	4.99e-04	1.80e-03
chr12	49776115	49776775	ENSMUSG00000112605	0.5023	-	CR974589.1	processed_pseudogene	2.47e-04	1.96e-04	4.87e-02	2.85e-06	1.25e-01	9.59e-02	6.54e-01	1.77e-03	2.86e-04	1.57e-01	-2.754888	4.33e+00	2.08e-02	1.88e+00	2.40e-03	3.17e+00	5.49e+00	6.39e-03	3.51e-02	2.32e+00	1.00e+00	2.32e+00	3.19e-03	7.99e-04	3.19e-03
chr15	81400138	81456807	ENSMUSG00000022401	0.4163	+	Xpnpep3	protein_coding	1.70e-03	9.63e-06	6.80e-04	9.02e-04	3.65e-01	2.20e-02	6.54e-01	1.10e-01	3.17e-04	1.66e-01	-1.431157	6.32e+00	4.58e-03	4.89e+00	1.66e-03	6.29e+00	6.36e+00	4.46e-03	4.69e-03	4.86e+00	5.04e+00	4.75e+00	1.62e-03	1.85e-03	1.51e-03
chr11	70651850	70654644	ENSMUSG00000018293	0.5571	-	Pfn1	protein_coding	1.57e-03	3.63e-05	1.30e-03	1.57e-04	3.60e-01	3.82e-02	6.54e-01	3.76e-02	3.28e-04	1.66e-01	1.388229	6.99e+00	3.07e-02	8.38e+00	8.08e-02	7.25e+00	6.73e+00	3.62e-02	2.52e-02	8.03e+00	8.65e+00	8.44e+00	6.26e-02	9.64e-02	8.33e-02
chr3	121787500	121788572	ENSMUSG00000105597	0.3523	-	4633401B06Rik	TEC	2.13e-03	3.58e-04	2.50e-03	1.51e-05	4.34e-01	1.22e-01	6.54e-01	7.12e-03	4.11e-04	1.96e-01	-3.247928	3.14e+00	6.79e-03	0.00e+00	0.00e+00	2.58e+00	3.70e+00	3.99e-03	9.58e-03	0.00e+00	0.00e+00	0.00e+00	0.00e+00	0.00e+00	0.00e+00
chr7	126363824	126369570	ENSMUSG00000030742	0.5069	-	Lat	protein_coding	1.35e-04	5.58e-03	2.03e-01	2.11e-07	8.38e-02	3.53e-01	7.17e-01	1.70e-04	4.24e-04	1.96e-01	-3.807355	2.38e+00	2.09e-03	0.00e+00	0.00e+00	0.00e+00	4.75e+00	0.00e+00	4.17e-03	0.00e+00	0.00e+00	0.00e+00	0.00e+00	0.00e+00	0.00e+00
chr4	61591929	61595871	ENSMUSG00000096688	0.4119	-	Mup17	protein_coding	1.17e-03	1.20e-02	8.28e-02	3.22e-08	3.16e-01	4.33e-01	6.55e-01	3.67e-05	4.40e-04	1.96e-01	1.576975	9.51e+00	3.19e-01	1.09e+01	9.51e-01	1.02e+01	8.85e+00	4.54e-01	1.83e-01	9.48e+00	1.19e+01	1.14e+01	2.82e-01	1.52e+00	1.05e+00
chr9	35204206	35211055	ENSMUSG00000039048	0.4778	-	Foxred1	protein_coding	1.47e-03	2.66e-05	1.32e-03	7.55e-04	3.53e-01	3.37e-02	6.54e-01	9.94e-02	4.45e-04	1.96e-01	-1.334984	7.01e+00	1.18e-02	5.66e+00	4.63e-03	6.94e+00	7.09e+00	1.12e-02	1.24e-02	5.81e+00	5.29e+00	5.88e+00	5.06e-03	3.50e-03	5.34e-03
chr9	78230656	78242684	ENSMUSG00000074183	0.4591	+	Gsta1	protein_coding	1.20e-03	2.95e-05	4.21e-03	4.15e-04	3.16e-01	3.37e-02	6.54e-01	6.76e-02	4.98e-04	2.08e-01	2.647698	1.29e+00	1.10e-03	3.96e+00	1.08e-02	1.58e+00	1.00e+00	1.47e-03	7.36e-04	4.17e+00	3.70e+00	4.00e+00	1.25e-02	8.84e-03	1.10e-02
chr6	11900292	11907497	ENSMUSG00000029632	0.3973	-	Ndufa4	protein_coding	6.75e-04	1.76e-04	9.42e-03	5.68e-05	2.37e-01	9.23e-02	6.54e-01	1.81e-02	5.02e-04	2.08e-01	1.456638	7.14e+00	9.36e-02	8.49e+00	2.58e-01	7.16e+00	7.12e+00	9.50e-02	9.23e-02	7.88e+00	9.22e+00	8.37e+00	1.57e-01	3.97e-01	2.21e-01
chr12	40201567	40248504	ENSMUSG00000001627	0.4154	-	Ifrd1	protein_coding	3.38e-03	1.34e-04	6.34e-04	2.68e-04	5.38e-01	9.22e-02	6.54e-01	5.16e-02	5.27e-04	2.12e-01	-1.114194	9.55e+00	4.96e-02	8.43e+00	2.29e-02	9.56e+00	9.54e+00	4.99e-02	4.94e-02	8.48e+00	8.48e+00	8.34e+00	2.37e-02	2.36e-02	2.14e-02
chr7	101263034	101302251	ENSMUSG00000047767	0.4544	-	Atg16l2	protein_coding	1.01e-04	8.83e-03	4.51e-03	2.56e-05	6.89e-02	4.03e-01	6.54e-01	9.74e-03	5.66e-04	2.17e-01	-1.584963	7.59e+00	7.96e-03	5.93e+00	2.63e-03	7.88e+00	7.30e+00	9.55e-03	6.38e-03	6.39e+00	6.30e+00	5.09e+00	3.37e-03	3.17e-03	1.34e-03
chr13	119690462	119708086	ENSMUSG00000093930	0.4159	+	Hmgcs1	protein_coding	2.80e-03	5.51e-03	2.72e-03	2.75e-06	5.07e-01	3.51e-01	6.54e-01	1.77e-03	5.83e-04	2.17e-01	-1.251504	1.12e+01	8.06e-01	9.88e+00	3.38e-01	1.10e+01	1.13e+01	7.17e-01	8.94e-01	1.04e+01	9.66e+00	9.61e+00	4.59e-01	2.82e-01	2.74e-01
chr3	144570304	144597680	ENSMUSG00000037072	0.4220	+	Selenof	protein_coding	2.29e-03	1.14e-03	9.88e-04	4.57e-05	4.54e-01	1.88e-01	6.54e-01	1.60e-02	5.86e-04	2.17e-01	1.295598	8.54e+00	5.09e-02	9.83e+00	1.25e-01	8.72e+00	8.37e+00	5.72e-02	4.47e-02	9.55e+00	1.01e+01	9.81e+00	1.02e-01	1.52e-01	1.22e-01
chr11	119077573	119086221	ENSMUSG00000039989	0.5687	-	Cbx4	protein_coding	1.53e-03	5.00e-05	4.67e-03	4.94e-04	3.60e-01	4.54e-02	6.54e-01	7.43e-02	6.48e-04	2.27e-01	1.649503	4.62e+00	4.21e-03	6.26e+00	1.36e-02	5.00e+00	4.25e+00	5.32e-03	3.09e-03	5.91e+00	6.07e+00	6.82e+00	1.01e-02	1.13e-02	1.92e-02
chr1	57388237	57407101	ENSMUSG00000048495	0.3937	-	Tyw5	protein_coding	1.70e-03	4.66e-05	9.31e-03	2.53e-04	3.65e-01	4.54e-02	6.54e-01	4.95e-02	6.57e-04	2.27e-01	2.514573	1.29e+00	3.18e-04	4.31e+00	2.42e-03	2.58e+00	0.00e+00	6.36e-04	0.00e+00	4.17e+00	4.58e+00	4.17e+00	2.16e-03	2.93e-03	2.16e-03
chr6	122577792	122602444	ENSMUSG00000040613	0.4301	-	Apobec1	protein_coding	1.92e-03	1.02e-04	1.66e-03	6.36e-04	3.99e-01	8.22e-02	6.54e-01	8.70e-02	6.75e-04	2.27e-01	-1.645083	5.60e+00	6.41e-03	3.93e+00	1.96e-03	5.43e+00	5.78e+00	5.61e-03	7.22e-03	3.81e+00	4.39e+00	3.58e+00	1.74e-03	2.67e-03	1.47e-03
chr8	85492576	85527560	ENSMUSG00000031700	0.4852	+	Gpt2	protein_coding	5.76e-03	3.66e-04	2.25e-03	4.71e-05	6.45e-01	1.22e-01	6.54e-01	1.61e-02	6.88e-04	2.27e-01	-1.112773	1.08e+01	1.16e-01	9.68e+00	5.35e-02	1.07e+01	1.09e+01	1.10e-01	1.22e-01	9.98e+00	9.71e+00	9.36e+00	6.48e-02	5.38e-02	4.21e-02
chr17	29090979	29100722	ENSMUSG00000023067	0.5389	+	Cdkn1a	protein_coding	7.88e-04	4.54e-04	8.99e-03	7.79e-05	2.58e-01	1.32e-01	6.54e-01	2.32e-02	7.08e-04	2.27e-01	-2.125531	4.87e+00	1.20e-02	2.77e+00	2.45e-03	4.25e+00	5.49e+00	6.97e-03	1.70e-02	2.00e+00	3.32e+00	3.00e+00	1.16e-03	3.49e-03	2.71e-03
chr19	5795690	5802672	ENSMUSG00000092341	0.4249	-	Malat1	lincRNA	3.70e-04	1.54e-03	2.44e-02	1.85e-05	1.67e-01	2.19e-01	6.54e-01	8.17e-03	7.12e-04	2.27e-01	-1.439652	8.59e+00	1.62e-01	7.02e+00	5.93e-02	9.10e+00	8.09e+00	2.16e-01	1.07e-01	8.03e+00	7.00e+00	6.02e+00	1.03e-01	5.01e-02	2.52e-02
chr1	71557150	71579631	ENSMUSG00000026192	0.4420	+	Atic	protein_coding	3.28e-03	5.31e-05	2.30e-03	8.39e-04	5.29e-01	4.54e-02	6.54e-01	1.06e-01	7.61e-04	2.37e-01	-1.279670	7.28e+00	2.16e-02	5.99e+00	8.82e-03	7.10e+00	7.47e+00	1.88e-02	2.44e-02	6.27e+00	6.07e+00	5.64e+00	1.05e-02	9.14e-03	6.79e-03
chr5	21543559	21575904	ENSMUSG00000039883	0.4349	+	Lrrc17	protein_coding	1.24e-04	5.25e-02	5.70e-01	1.11e-07	8.09e-02	6.51e-01	8.82e-01	1.08e-04	8.01e-04	2.44e-01	-1.666967	8.09e+00	2.43e-01	7.24e+00	7.64e-02	5.70e+00	1.05e+01	1.67e-02	4.70e-01	6.00e+00	9.07e+00	6.66e+00	2.07e-02	1.76e-01	3.28e-02
chr10	128624534	128626747	ENSMUSG00000025362	0.4783	-	Rps26	protein_coding	2.82e-03	1.72e-04	1.26e-02	8.66e-05	5.07e-01	9.23e-02	6.54e-01	2.52e-02	8.52e-04	2.51e-01	1.383235	7.41e+00	6.65e-02	8.71e+00	1.74e-01	7.64e+00	7.17e+00	7.74e-02	5.56e-02	7.98e+00	8.80e+00	9.34e+00	9.76e-02	1.73e-01	2.52e-01
chr16	84954440	85173707	ENSMUSG00000022892	0.4279	-	App	protein_coding	3.19e-03	2.14e-04	3.31e-03	2.47e-04	5.29e-01	9.79e-02	6.54e-01	4.95e-02	8.64e-04	2.51e-01	-1.212840	8.73e+00	1.99e-01	7.48e+00	8.54e-02	8.59e+00	8.87e+00	1.79e-01	2.18e-01	7.94e+00	7.10e+00	7.41e+00	1.14e-01	6.35e-02	7.89e-02
chr15	85734983	85802819	ENSMUSG00000022383	0.4601	+	Ppara	protein_coding	9.42e-04	7.98e-04	7.46e-03	1.10e-04	2.78e-01	1.52e-01	6.54e-01	3.04e-02	8.86e-04	2.53e-01	-1.638261	6.59e+00	2.27e-02	4.69e+00	7.14e-03	6.55e+00	6.63e+00	2.21e-02	2.33e-02	5.78e+00	4.70e+00	3.58e+00	1.29e-02	5.95e-03	2.62e-03
chr5	115241412	115272898	ENSMUSG00000041740	0.4778	-	Rnf10	protein_coding	2.49e-03	2.15e-04	5.08e-03	2.85e-04	4.80e-01	9.79e-02	6.54e-01	5.35e-02	9.38e-04	2.62e-01	-1.316710	7.84e+00	2.45e-02	6.45e+00	9.78e-03	7.74e+00	7.94e+00	2.28e-02	2.62e-02	7.03e+00	6.48e+00	5.83e+00	1.39e-02	9.42e-03	6.00e-03
chr5	72809968	72812074	ENSMUSG00000106739	0.4338	-	Gm43499	TEC	7.88e-04	5.47e-03	1.67e-02	1.41e-05	2.58e-01	3.51e-01	6.54e-01	6.89e-03	1.00e-03	2.75e-01	-3.017922	3.49e+00	9.98e-03	6.67e-01	5.32e-04	4.39e+00	2.58e+00	1.60e-02	3.99e-03	1.00e+00	1.00e+00	0.00e+00	7.99e-04	7.99e-04	0.00e+00
chr12	51133551	51134702	ENSMUSG00000112926	0.5208	-	CT571252.1	processed_pseudogene	5.70e-04	3.26e-04	2.02e-02	3.20e-04	2.20e-01	1.22e-01	6.54e-01	5.91e-02	1.05e-03	2.81e-01	-1.974005	4.61e+00	2.12e-02	2.57e+00	4.79e-03	5.32e+00	3.91e+00	3.12e-02	1.12e-02	3.70e+00	2.00e+00	2.00e+00	9.58e-03	2.40e-03	2.40e-03
chr1	120067377	120074074	ENSMUSG00000050777	0.5361	-	Tmem37	protein_coding	1.20e-03	8.04e-04	7.17e-03	2.13e-04	3.16e-01	1.52e-01	6.54e-01	4.59e-02	1.10e-03	2.84e-01	1.218907	8.21e+00	1.86e-01	9.37e+00	4.34e-01	8.26e+00	8.16e+00	1.92e-01	1.80e-01	9.94e+00	9.03e+00	9.13e+00	6.20e-01	3.29e-01	3.53e-01
chrX	155323918	155340754	ENSMUSG00000025289	0.3851	-	Prdx4	protein_coding	2.12e-04	7.13e-04	2.25e-02	4.58e-04	1.12e-01	1.51e-01	6.54e-01	7.06e-02	1.12e-03	2.84e-01	1.385431	6.30e+00	2.07e-02	7.58e+00	5.45e-02	6.04e+00	6.55e+00	1.70e-02	2.44e-02	8.41e+00	7.16e+00	7.18e+00	8.86e-02	3.72e-02	3.77e-02
chr7	35119293	35121928	ENSMUSG00000034957	0.6259	+	Cebpa	protein_coding	8.22e-03	2.53e-04	2.10e-03	3.90e-04	7.65e-01	1.08e-01	6.54e-01	6.59e-02	1.14e-03	2.84e-01	-1.042479	9.94e+00	3.52e-01	8.90e+00	1.71e-01	9.77e+00	1.01e+01	3.10e-01	3.93e-01	9.03e+00	8.74e+00	8.94e+00	1.85e-01	1.52e-01	1.75e-01
chr9	74861921	74891781	ENSMUSG00000043013	0.4587	+	Onecut1	protein_coding	3.78e-04	1.85e-03	1.94e-02	1.34e-04	1.67e-01	2.29e-01	6.54e-01	3.40e-02	1.16e-03	2.84e-01	-1.714246	5.98e+00	2.01e-02	4.22e+00	5.92e-03	5.32e+00	6.64e+00	1.14e-02	2.88e-02	3.70e+00	5.25e+00	3.70e+00	3.49e-03	1.08e-02	3.49e-03
chr4	116551633	116555936	ENSMUSG00000055900	0.4064	-	Tmem69	protein_coding	4.30e-03	1.75e-04	2.12e-03	1.26e-03	6.05e-01	9.23e-02	6.54e-01	1.30e-01	1.19e-03	2.84e-01	-1.241889	7.21e+00	3.43e-02	5.96e+00	1.44e-02	7.06e+00	7.37e+00	3.06e-02	3.80e-02	6.29e+00	5.73e+00	5.86e+00	1.78e-02	1.20e-02	1.32e-02
chr8	94179082	94180327	ENSMUSG00000031765	0.5425	+	Mt1	protein_coding	1.58e-03	1.37e-02	6.02e-03	1.66e-05	3.60e-01	4.51e-01	6.54e-01	7.55e-03	1.21e-03	2.84e-01	1.368664	8.49e+00	1.70e-01	9.82e+00	4.40e-01	8.77e+00	8.21e+00	2.03e-01	1.37e-01	9.35e+00	1.04e+01	9.68e+00	3.02e-01	6.38e-01	3.81e-01
chr15	81318921	81360765	ENSMUSG00000022404	0.4325	-	Slc25a17	protein_coding	1.70e-03	1.95e-03	1.28e-02	5.32e-05	3.65e-01	2.33e-01	6.54e-01	1.75e-02	1.23e-03	2.84e-01	-1.376148	8.51e+00	1.86e-01	6.96e+00	7.14e-02	8.40e+00	8.63e+00	1.71e-01	2.01e-01	7.79e+00	7.11e+00	5.98e+00	1.13e-01	6.99e-02	3.16e-02
chr4	155601416	155623340	ENSMUSG00000042202	0.4552	+	Slc35e2	protein_coding	1.42e-03	7.20e-04	5.11e-02	4.41e-05	3.51e-01	1.51e-01	6.54e-01	1.59e-02	1.23e-03	2.84e-01	-1.465883	7.77e+00	4.46e-02	6.52e+00	1.60e-02	6.88e+00	8.66e+00	2.00e-02	6.91e-02	6.77e+00	6.78e+00	6.02e+00	1.85e-02	1.87e-02	1.10e-02
chr11	75461715	75470899	ENSMUSG00000038217	0.5005	+	Tlcd2	protein_coding	7.28e-03	1.35e-04	1.47e-03	1.67e-03	7.22e-01	9.22e-02	6.54e-01	1.45e-01	1.25e-03	2.84e-01	1.130086	7.33e+00	3.87e-02	8.47e+00	8.50e-02	7.52e+00	7.14e+00	4.38e-02	3.35e-02	8.33e+00	8.67e+00	8.40e+00	7.69e-02	9.75e-02	8.05e-02
chr8	25601601	25719667	ENSMUSG00000054823	0.4305	+	Nsd3	protein_coding	4.55e-03	1.82e-04	5.07e-03	6.35e-04	6.05e-01	9.23e-02	6.54e-01	8.70e-02	1.28e-03	2.86e-01	-1.183057	8.22e+00	1.92e-02	7.04e+00	8.43e-03	7.93e+00	8.50e+00	1.54e-02	2.30e-02	7.14e+00	7.31e+00	6.66e+00	8.89e-03	1.00e-02	6.35e-03
chr6	56777524	56797813	ENSMUSG00000059486	0.4121	-	Kbtbd2	protein_coding	2.77e-03	8.13e-04	7.32e-03	2.45e-04	5.07e-01	1.52e-01	6.54e-01	4.95e-02	1.42e-03	3.12e-01	-1.305503	8.03e+00	6.52e-02	6.62e+00	2.62e-02	7.95e+00	8.11e+00	6.15e-02	6.89e-02	7.31e+00	6.61e+00	5.93e+00	3.95e-02	2.42e-02	1.50e-02
chr8	11654924	11701063	ENSMUSG00000031509	0.4681	-	1700016D06Rik	protein_coding	9.81e-04	7.51e-04	1.18e-02	5.06e-04	2.80e-01	1.52e-01	6.54e-01	7.45e-02	1.45e-03	3.12e-01	-1.459432	6.48e+00	4.79e-02	4.95e+00	1.71e-02	6.93e+00	6.02e+00	6.26e-02	3.31e-02	5.78e+00	4.32e+00	4.75e+00	2.79e-02	9.83e-03	1.35e-02
chr9	31386192	31421333	ENSMUSG00000042185	0.4474	+	Nfrkb	protein_coding	2.90e-03	1.44e-03	3.86e-03	2.82e-04	5.16e-01	2.19e-01	6.54e-01	5.35e-02	1.46e-03	3.12e-01	-2.672425	3.08e+00	2.54e-04	3.33e-01	1.13e-05	3.00e+00	3.17e+00	2.37e-04	2.71e-04	0.00e+00	1.00e+00	0.00e+00	0.00e+00	3.38e-05	0.00e+00
chr5	137528127	137533510	ENSMUSG00000029713	0.6009	-	Gnb2	protein_coding	3.01e-03	1.46e-03	7.24e-03	2.04e-04	5.22e-01	2.19e-01	6.54e-01	4.50e-02	1.60e-03	3.35e-01	-1.245774	8.53e+00	1.72e-01	7.36e+00	7.23e-02	8.06e+00	9.01e+00	1.18e-01	2.27e-01	7.51e+00	7.30e+00	7.27e+00	8.00e-02	6.94e-02	6.76e-02
chr2	153616933	153632781	ENSMUSG00000056941	0.4698	-	Commd7	protein_coding	5.46e-03	1.49e-04	5.40e-03	1.55e-03	6.38e-01	9.23e-02	6.54e-01	1.43e-01	1.62e-03	3.35e-01	-1.240177	7.00e+00	3.48e-02	5.78e+00	1.46e-02	6.64e+00	7.35e+00	2.64e-02	4.31e-02	6.07e+00	5.58e+00	5.70e+00	1.76e-02	1.25e-02	1.36e-02
chr5	33897916	33936413	ENSMUSG00000029111	0.4385	-	Nelfa	protein_coding	3.59e-03	3.42e-04	3.38e-03	1.75e-03	5.56e-01	1.22e-01	6.54e-01	1.45e-01	1.64e-03	3.35e-01	-1.237864	6.94e+00	3.74e-02	5.71e+00	1.57e-02	6.67e+00	7.21e+00	3.05e-02	4.44e-02	5.67e+00	5.49e+00	5.98e+00	1.51e-02	1.33e-02	1.87e-02
chr6	71388634	71440637	ENSMUSG00000002222	0.4099	-	Rmnd5a	protein_coding	1.21e-02	5.87e-04	7.21e-03	1.51e-04	8.54e-01	1.46e-01	6.54e-01	3.68e-02	1.67e-03	3.35e-01	-1.059921	1.05e+01	2.12e-01	9.43e+00	1.02e-01	1.04e+01	1.07e+01	1.87e-01	2.37e-01	9.88e+00	9.39e+00	9.02e+00	1.34e-01	9.61e-02	7.42e-02
chr7	19344778	19356524	ENSMUSG00000003549	0.5267	+	Ercc1	protein_coding	4.85e-04	9.52e-04	5.59e-02	3.26e-04	2.01e-01	1.67e-01	6.54e-01	5.95e-02	1.70e-03	3.38e-01	3.115477	0.00e+00	0.00e+00	2.71e+00	6.61e-04	0.00e+00	0.00e+00	0.00e+00	0.00e+00	3.32e+00	3.81e+00	1.00e+00	7.76e-04	1.12e-03	8.62e-05
chr10	81384428	81388352	ENSMUSG00000078440	0.5804	+	Dohh	protein_coding	6.47e-03	3.11e-04	4.37e-03	1.17e-03	6.81e-01	1.22e-01	6.54e-01	1.27e-01	1.79e-03	3.44e-01	-1.191784	7.63e+00	4.49e-02	6.41e+00	1.95e-02	7.42e+00	7.83e+00	3.84e-02	5.13e-02	6.86e+00	6.19e+00	6.19e+00	2.60e-02	1.63e-02	1.63e-02
chr7	104409025	104465193	ENSMUSG00000030921	0.4104	-	Trim30a	protein_coding	5.96e-04	1.03e-03	1.47e-02	1.18e-03	2.20e-01	1.75e-01	6.54e-01	1.27e-01	1.80e-03	3.44e-01	-1.817623	4.48e+00	3.87e-03	2.43e+00	9.76e-04	4.95e+00	4.00e+00	5.17e-03	2.58e-03	1.58e+00	3.70e+00	2.00e+00	3.44e-04	2.07e-03	5.17e-04
chr14	34310727	34345144	ENSMUSG00000021794	0.4228	+	Glud1	protein_coding	2.96e-02	2.61e-03	1.14e-02	1.21e-05	1.00e+00	2.61e-01	6.54e-01	6.14e-03	1.81e-03	3.44e-01	-0.970598	1.32e+01	1.46e+00	1.22e+01	7.43e-01	1.31e+01	1.33e+01	1.34e+00	1.57e+00	1.26e+01	1.24e+01	1.17e+01	9.42e-01	7.89e-01	4.98e-01
chr10	95480805	95501940	ENSMUSG00000062981	0.4277	-	Mrpl42	protein_coding	6.25e-03	1.63e-04	3.44e-03	3.33e-03	6.64e-01	9.23e-02	6.54e-01	2.03e-01	1.85e-03	3.46e-01	1.309544	5.27e+00	9.89e-03	6.57e+00	2.49e-02	5.49e+00	5.04e+00	1.14e-02	8.33e-03	6.23e+00	6.87e+00	6.61e+00	1.93e-02	3.02e-02	2.52e-02
chr18	20944625	20983848	ENSMUSG00000033107	0.4397	+	Rnf125	protein_coding	3.05e-03	4.29e-04	9.85e-03	9.87e-04	5.22e-01	1.32e-01	6.54e-01	1.17e-01	1.89e-03	3.49e-01	-1.268778	7.24e+00	7.83e-02	5.86e+00	3.22e-02	7.23e+00	7.25e+00	7.78e-02	7.89e-02	6.29e+00	6.29e+00	5.00e+00	4.02e-02	4.02e-02	1.62e-02
chr1	118458652	118459267	ENSMUSG00000107722	0.4237	-	2900060B14Rik	antisense	9.13e-04	6.78e-03	5.42e-03	4.25e-04	2.78e-01	3.79e-01	6.54e-01	6.79e-02	1.94e-03	3.52e-01	-1.546488	6.20e+00	5.75e-02	4.53e+00	1.92e-02	6.32e+00	6.09e+00	6.22e-02	5.28e-02	5.00e+00	5.00e+00	3.58e+00	2.44e-02	2.44e-02	8.66e-03
chrX	96096042	96168552	ENSMUSG00000031207	0.4125	+	Msn	protein_coding	5.72e-03	4.26e-04	2.43e-03	2.79e-03	6.45e-01	1.32e-01	6.54e-01	1.88e-01	2.02e-03	3.52e-01	-1.223167	6.67e+00	1.81e-02	5.47e+00	7.64e-03	6.43e+00	6.92e+00	1.50e-02	2.11e-02	5.55e+00	5.39e+00	5.46e+00	8.11e-03	7.23e-03	7.58e-03
chr19	5490455	5495201	ENSMUSG00000056201	0.5264	+	Cfl1	protein_coding	3.93e-03	6.13e-04	4.30e-03	1.61e-03	5.78e-01	1.46e-01	6.54e-01	1.44e-01	2.02e-03	3.52e-01	1.305413	5.94e+00	2.04e-02	7.21e+00	5.10e-02	6.07e+00	5.81e+00	2.23e-02	1.86e-02	6.79e+00	7.65e+00	7.17e+00	3.72e-02	6.76e-02	4.83e-02
chr13	101680563	101768217	ENSMUSG00000041417	0.4414	-	Pik3r1	protein_coding	8.96e-03	2.69e-04	2.49e-03	2.83e-03	7.90e-01	1.11e-01	6.54e-01	1.88e-01	2.03e-03	3.52e-01	-1.071681	7.79e+00	1.88e-02	6.73e+00	8.90e-03	7.60e+00	7.98e+00	1.63e-02	2.13e-02	6.81e+00	6.77e+00	6.61e+00	9.38e-03	9.13e-03	8.20e-03
chr13	12395027	12440289	ENSMUSG00000050244	0.4288	+	Heatr1	protein_coding	5.39e-03	5.95e-04	3.13e-03	1.70e-03	6.38e-01	1.46e-01	6.54e-01	1.45e-01	2.03e-03	3.52e-01	-1.277873	6.73e+00	1.01e-02	5.45e+00	4.10e-03	6.49e+00	6.97e+00	8.43e-03	1.17e-02	5.78e+00	5.36e+00	5.21e+00	5.11e-03	3.79e-03	3.41e-03
chr12	83950608	83952951	ENSMUSG00000046791	0.6139	+	Riox1	protein_coding	5.73e-03	3.76e-04	6.67e-03	1.37e-03	6.45e-01	1.22e-01	6.54e-01	1.32e-01	2.11e-03	3.60e-01	-1.416208	6.04e+00	4.66e-02	4.66e+00	1.70e-02	5.64e+00	6.44e+00	3.39e-02	5.94e-02	4.95e+00	4.70e+00	4.32e+00	2.07e-02	1.73e-02	1.31e-02
chr3	125865271	125938619	ENSMUSG00000032854	0.3789	-	Ugt8a	protein_coding	7.87e-03	2.02e-03	3.37e-03	4.00e-04	7.43e-01	2.33e-01	6.54e-01	6.59e-02	2.15e-03	3.62e-01	-2.700440	2.70e+00	1.24e-03	0.00e+00	0.00e+00	2.58e+00	2.81e+00	1.12e-03	1.35e-03	0.00e+00	0.00e+00	0.00e+00	0.00e+00	0.00e+00	0.00e+00
chr18	64340364	64398488	ENSMUSG00000045991	0.4448	+	Onecut2	protein_coding	6.19e-03	3.09e-03	1.84e-03	6.32e-04	6.64e-01	2.75e-01	6.54e-01	8.70e-02	2.17e-03	3.62e-01	-0.981708	1.02e+01	8.86e-02	9.20e+00	4.48e-02	1.03e+01	1.00e+01	9.68e-02	8.04e-02	9.03e+00	9.29e+00	9.27e+00	3.97e-02	4.78e-02	4.71e-02
chr1	9547948	9580673	ENSMUSG00000025911	0.4316	+	Adhfe1	protein_coding	4.54e-03	8.02e-04	4.34e-03	1.62e-03	6.05e-01	1.52e-01	6.54e-01	1.44e-01	2.25e-03	3.63e-01	1.340307	5.70e+00	3.15e-03	7.00e+00	8.06e-03	5.78e+00	5.61e+00	3.33e-03	2.96e-03	6.55e+00	6.99e+00	7.44e+00	5.74e-03	7.77e-03	1.07e-02
chr14	52006077	52020733	ENSMUSG00000049295	0.5403	-	Zfp219	protein_coding	8.03e-03	3.17e-04	3.47e-03	2.92e-03	7.53e-01	1.22e-01	6.54e-01	1.92e-01	2.25e-03	3.63e-01	-1.403516	5.60e+00	2.77e-02	4.14e+00	1.01e-02	5.58e+00	5.61e+00	2.74e-02	2.80e-02	4.70e+00	3.81e+00	3.91e+00	1.46e-02	7.58e-03	8.16e-03
chr17	83350931	83480361	ENSMUSG00000032624	0.4079	+	Eml4	protein_coding	1.44e-03	8.99e-04	8.60e-03	2.41e-03	3.51e-01	1.62e-01	6.54e-01	1.77e-01	2.27e-03	3.63e-01	-1.448901	5.45e+00	1.29e-02	4.01e+00	4.55e-03	5.86e+00	5.04e+00	1.65e-02	9.29e-03	3.58e+00	4.46e+00	4.00e+00	3.19e-03	6.09e-03	4.35e-03
chr10	128900989	128912822	ENSMUSG00000025351	0.4956	+	Cd63	protein_coding	2.22e-03	4.12e-03	2.63e-02	1.17e-04	4.46e-01	3.18e-01	6.54e-01	3.10e-02	2.30e-03	3.63e-01	-1.314423	8.46e+00	1.27e-01	6.89e+00	5.10e-02	8.50e+00	8.41e+00	1.32e-01	1.23e-01	7.77e+00	7.30e+00	5.58e+00	7.89e-02	5.71e-02	1.71e-02
chr7	100996568	101012866	ENSMUSG00000032860	0.5171	-	P2ry2	protein_coding	5.30e-03	6.59e-04	6.44e-03	1.26e-03	6.38e-01	1.49e-01	6.54e-01	1.30e-01	2.31e-03	3.63e-01	-1.430792	6.05e+00	1.27e-02	4.49e+00	4.58e-03	5.98e+00	6.13e+00	1.20e-02	1.34e-02	5.25e+00	4.52e+00	3.70e+00	7.17e-03	4.26e-03	2.32e-03
chr5	31202664	31220687	ENSMUSG00000029147	0.4314	-	Ppm1g	protein_coding	4.87e-03	2.25e-03	1.42e-02	2.40e-04	6.22e-01	2.50e-01	6.54e-01	4.95e-02	2.47e-03	3.85e-01	-1.266646	8.30e+00	7.55e-02	6.91e+00	3.12e-02	8.08e+00	8.53e+00	6.39e-02	8.71e-02	7.71e+00	6.89e+00	6.13e+00	4.95e-02	2.79e-02	1.63e-02
chr6	125460266	125494791	ENSMUSG00000030342	0.5016	-	Cd9	protein_coding	5.03e-04	5.20e-02	7.33e-02	2.19e-05	2.02e-01	6.51e-01	6.54e-01	9.10e-03	2.55e-03	3.89e-01	-1.501131	7.89e+00	1.01e-01	5.64e+00	3.53e-02	8.85e+00	6.93e+00	1.59e-01	4.17e-02	8.02e+00	5.21e+00	3.70e+00	8.93e-02	1.24e-02	4.14e-03
chr5	136136146	136141615	ENSMUSG00000039754	0.5333	+	Alkbh4	protein_coding	1.06e-02	4.13e-03	2.60e-03	4.00e-04	8.30e-01	3.18e-01	6.54e-01	6.59e-02	2.60e-03	3.89e-01	-2.700440	2.70e+00	1.42e-03	0.00e+00	0.00e+00	2.58e+00	2.81e+00	1.29e-03	1.54e-03	0.00e+00	0.00e+00	0.00e+00	0.00e+00	0.00e+00	0.00e+00
chrX	73086293	73097095	ENSMUSG00000057836	0.4542	-	Xlr3a	protein_coding	1.79e-03	2.53e-03	1.08e-02	1.06e-03	3.77e-01	2.57e-01	6.54e-01	1.21e-01	2.68e-03	3.89e-01	1.559084	4.64e+00	2.13e-03	6.12e+00	6.46e-03	4.52e+00	4.75e+00	1.95e-03	2.31e-03	5.43e+00	6.67e+00	6.25e+00	3.73e-03	8.97e-03	6.66e-03

Download the DEG result list for Grp4fMinus_vs_Grp4fPlus.
Download the whole result list for Grp4fMinus_vs_Grp4fPlus.

References

Statham, A.L., Strbenac, D., Coolen, M.W., Stirzaker, C., Clark, S.J., Robinson, M.D. (2010). Repitools: an R package for the analysis of enrichment-based epigenomic data. Bioinformatics 26(13), 1662-1663.
Anders, S., and Huber, W. (2010). Differential expression analysis for sequence count data. Genome Biol 11, R106.
Robinson, M.D., McCarthy, D.J., and Smyth, G.K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139-140.
Smyth, G. (2005). Limma: linear models for microarray data. In Bioinformatics and Computational Biology Solutions using R and Bioconductor, G. R., C. V., D. S., I. R., and H. W., eds. (New York, Springer), pp. 397-420.
Di, Y, Schafer, D., Cumbie, J.S., and Chang, J.H. (2011). The NBP Negative Binomial Model for Assessing Differential Gene Expression from RNA-Seq. Statistical Applications in Genetics and Molecular Biology 10(1), 1-28.
Planet, E., Attolini, C.S., Reina, O., Flores, O., and Rossell, D. (2012). htSeqTools: high-throughput sequencing quality control, processing and visualization in R. Bioinformatics 28, 589-590.
Tarazona, S., Garcia-Alcalde, F., Dopazo, J., Ferrer, A., and Conesa, A. (2011). Differential expression in RNA-seq: a matter of depth. Genome Res 21, 2213-2223.
Risso, D., Schwartz, K., Sherlock, G., and Dudoit, S. (2011). GC-content normalization for RNA-Seq data. BMC Bioinformatics 12, 480.
Benjamini, Y., and Hochberg, Y. (1995). Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B (Methodological) 57, 289-300.