metaseqR report

Summary

Analysis summary

Summary: The raw bam files, one for each RNA-Seq sample, were summarized to a gene read counts table, using the Bioconductor package GenomicRanges. In the final read counts table, each row represented one gene, 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 edgeR after removing genes that had zero counts over all the RNA-Seq samples (1201 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 (232 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.13038 average reads per 100 bp), iii) genes with read counts below the median read counts of the total normalized count distribution (9161 genes with cutoff value 549 normalized read counts). The total number of genes excluded due to the application of gene filters was 4039. The total (unified) number of genes excluded due to the application of all filters was 10393. The resulting gene counts table was subjected to differential expression analysis for the contrasts WT versus HET using the Bioconductor package edgeR. The final numbers of differentially expressed genes were (per contrast): for the contrast WT versus HET, 414 (91) statistically significant genes were found with a p-value (FDR or adjusted p-value) threshold of 0.05 and of these, 53 (9) were up-regulated, 200 (82) were down-regulated and 161 (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: WT, HET
Samples: SCR20, SCR21, SCR22, SCR23
Requested contrasts: WT_vs_HET
Library sizes:

SCR20: 44715186
SCR21: 36665351
SCR22: 34275109
SCR23: 45423883

Annotation: download
Organism: mouse (Mus musculus), genome version alias mm9
Annotation source: RefSeq database
Count type: gene
Exon filters: min.active.exons

min.active.exons
- exons.per.gene: 5
- min.exons: 2
- frac: 0.2

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
- lincRNA: FALSE
- snoRNA: FALSE
- processed_transcript: FALSE
- misc_RNA: FALSE
- rRNA: TRUE
- sense_overlapping: FALSE
- sense_intronic: FALSE
- polymorphic_pseudogene: FALSE
- non_coding: FALSE
- three_prime_overlapping_ncrna: FALSE
- IG_C_gene: FALSE
- IG_J_gene: FALSE
- IG_D_gene: FALSE
- IG_V_gene: FALSE
- ncrna_host: FALSE

Filter application: after normalization
Normalization algorithm: edgeR
Normalization arguments: method, refColumn, logratioTrim, sumTrim, doWeighting, Acutoff, p

method: TMM
logratioTrim: 0.3
sumTrim: 0.05
doWeighting: TRUE
Acutoff: -1e+10
p: 0.75

Statistical algorithm(s): edgeR
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

Meta-analysis method: no meta-analysis, reported p-values from the first supplied statistical algorithm
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, correlation heatmap and correlogram, pairwise scatterplots between samples, boxplots, volcano plot, DEG biotype detection, DEG heatmap
Figure format: png, pdf
Output directory: /data/images/proton/Cereghini/run142/metaseqr_wt_vs_het_mm9_test2
Output data: Annotation, p-value, Adjusted p-value (FDR), Fold change, Statistics, Read counts
Output scale(s): log2 scale
Output values: Normalized values
Output statistics: Mean, Median, Standard deviation, Median Absolute Deviation (MAD), Coefficient of Variation, Robust Coefficient of Variation
Total run time: 18 minutes 40 seconds

Filtered genes

Number of filtered genes: 10393 which is the union of

Filtered because of zero reads: 1201
Filtered because of exon filters: 0
Filtered because of gene filters: 9192 which is the union of
- length: 232 genes with filter cutoff value 500
- avg.reads: 3804 genes with filter cutoff value 0.1303781
- expression: 9161 genes further decomposed to (filter name, filtered genes, filter cutoff):
- - median: 9161 genes with filter cutoff value 549
- biotype: 2 genes with filter cutoff value rRNA

Differentially expressed genes

Number of differentially expressed genes per contrast:

WT_vs_HET: 414 (91) statistically significant genes of which 53 (9) up regulated, 200 (82) down regulated and 161 (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.

Run log

[ 2015-02-24 16:12:43 ] INFO 2015-02-24 16:12:43: Data processing started...
[ 2015-02-24 16:12:43 ] INFO Read counts file: imported sam/bam/bed files
[ 2015-02-24 16:12:43 ] INFO Conditions: WT, HET
[ 2015-02-24 16:12:43 ] INFO Samples: SCR20, SCR21, SCR22, SCR23
[ 2015-02-24 16:12:43 ] INFO Requested contrasts: WT_vs_HET
[ 2015-02-24 16:12:43 ] INFO Annotation: download
[ 2015-02-24 16:12:43 ] INFO Organism: mm9
[ 2015-02-24 16:12:43 ] INFO Reference source: refseq
[ 2015-02-24 16:12:43 ] INFO Count type: gene
[ 2015-02-24 16:12:43 ] INFO Exon filters: min.active.exons
[ 2015-02-24 16:12:43 ] INFO min.active.exons:
[ 2015-02-24 16:12:43 ] INFO exons.per.gene: 5
[ 2015-02-24 16:12:43 ] INFO min.exons: 2
[ 2015-02-24 16:12:43 ] INFO frac: 0.2
[ 2015-02-24 16:12:43 ] INFO Gene filters: length, avg.reads, expression, biotype
[ 2015-02-24 16:12:43 ] INFO length:
[ 2015-02-24 16:12:43 ] INFO length: 500
[ 2015-02-24 16:12:43 ] INFO avg.reads:
[ 2015-02-24 16:12:43 ] INFO average.per.bp: 100
[ 2015-02-24 16:12:43 ] INFO quantile: 0.25
[ 2015-02-24 16:12:43 ] INFO expression:
[ 2015-02-24 16:12:43 ] INFO median: TRUE
[ 2015-02-24 16:12:43 ] INFO mean: FALSE
[ 2015-02-24 16:12:43 ] INFO quantile: NA
[ 2015-02-24 16:12:43 ] INFO known: NA
[ 2015-02-24 16:12:43 ] INFO custom: NA
[ 2015-02-24 16:12:43 ] INFO biotype:
[ 2015-02-24 16:12:43 ] INFO pseudogene: FALSE
[ 2015-02-24 16:12:43 ] INFO snRNA: FALSE
[ 2015-02-24 16:12:43 ] INFO protein_coding: FALSE
[ 2015-02-24 16:12:43 ] INFO antisense: FALSE
[ 2015-02-24 16:12:43 ] INFO miRNA: FALSE
[ 2015-02-24 16:12:43 ] INFO lincRNA: FALSE
[ 2015-02-24 16:12:43 ] INFO snoRNA: FALSE
[ 2015-02-24 16:12:43 ] INFO processed_transcript: FALSE
[ 2015-02-24 16:12:43 ] INFO misc_RNA: FALSE
[ 2015-02-24 16:12:43 ] INFO rRNA: TRUE
[ 2015-02-24 16:12:43 ] INFO sense_overlapping: FALSE
[ 2015-02-24 16:12:43 ] INFO sense_intronic: FALSE
[ 2015-02-24 16:12:43 ] INFO polymorphic_pseudogene: FALSE
[ 2015-02-24 16:12:43 ] INFO non_coding: FALSE
[ 2015-02-24 16:12:43 ] INFO three_prime_overlapping_ncrna: FALSE
[ 2015-02-24 16:12:43 ] INFO IG_C_gene: FALSE
[ 2015-02-24 16:12:43 ] INFO IG_J_gene: FALSE
[ 2015-02-24 16:12:43 ] INFO IG_D_gene: FALSE
[ 2015-02-24 16:12:43 ] INFO IG_V_gene: FALSE
[ 2015-02-24 16:12:43 ] INFO ncrna_host: FALSE
[ 2015-02-24 16:12:43 ] INFO Filter application: postnorm
[ 2015-02-24 16:12:43 ] INFO Normalization algorithm: edger
[ 2015-02-24 16:12:43 ] INFO Normalization arguments:
[ 2015-02-24 16:12:43 ] INFO method: TMM
[ 2015-02-24 16:12:43 ] INFO logratioTrim: 0.3
[ 2015-02-24 16:12:43 ] INFO sumTrim: 0.05
[ 2015-02-24 16:12:43 ] INFO doWeighting: TRUE
[ 2015-02-24 16:12:43 ] INFO Acutoff: -1e+10
[ 2015-02-24 16:12:43 ] INFO p: 0.75
[ 2015-02-24 16:12:43 ] INFO Statistical algorithm: edger
[ 2015-02-24 16:12:43 ] INFO Statistical arguments:
[ 2015-02-24 16:12:43 ] 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)
[ 2015-02-24 16:12:43 ] INFO Meta-analysis method: none
[ 2015-02-24 16:12:43 ] INFO Multiple testing correction: BH
[ 2015-02-24 16:12:43 ] INFO p-value threshold: 0.05
[ 2015-02-24 16:12:43 ] INFO Logarithmic transformation offset: 1
[ 2015-02-24 16:12:43 ] INFO Quality control plots: mds, biodetection, countsbio, saturation, readnoise, correl, pairwise, boxplot, volcano, biodist, deheatmap
[ 2015-02-24 16:12:43 ] INFO Figure format: png, pdf
[ 2015-02-24 16:12:43 ] INFO Output directory: /data/images/proton/Cereghini/run142/metaseqr_wt_vs_het_mm9_test2
[ 2015-02-24 16:12:43 ] INFO Output data: annotation, p.value, adj.p.value, fold.change, stats, counts
[ 2015-02-24 16:12:43 ] INFO Output scale(s): log2
[ 2015-02-24 16:12:43 ] INFO Output values: normalized
[ 2015-02-24 16:12:43 ] INFO Output statistics: mean, median, sd, mad, cv, rcv
[ 2015-02-24 16:12:43 ] INFO Downloading gene annotation for mm9...
[ 2015-02-24 16:12:57 ] INFO Converting annotation to GenomicRanges object...
[ 2015-02-24 16:12:57 ] INFO Getting DNA sequences...
[ 2015-02-24 16:13:04 ] INFO Getting GC content...
[ 2015-02-24 16:13:05 ] INFO Converting annotation to GenomicRanges object...
[ 2015-02-24 16:13:05 ] INFO Reading bam file sort_uniq.bam for sample with name SCR20. This might take some time...
[ 2015-02-24 16:13:05 ] INFO Reading bam file sort_uniq.bam for sample with name SCR21. This might take some time...
[ 2015-02-24 16:13:05 ] INFO Reading bam file sort_uniq.bam for sample with name SCR22. This might take some time...
[ 2015-02-24 16:13:05 ] INFO Reading bam file sort_uniq.bam for sample with name SCR23. This might take some time...
[ 2015-02-24 16:13:51 ] INFO Counting reads overlapping with given annotation...
[ 2015-02-24 16:13:51 ] INFO ...for single-end reads...
[ 2015-02-24 16:13:51 ] INFO ...using strandedness...
[ 2015-02-24 16:14:16 ] INFO Counting reads overlapping with given annotation...
[ 2015-02-24 16:14:16 ] INFO ...for single-end reads...
[ 2015-02-24 16:14:16 ] INFO ...using strandedness...
[ 2015-02-24 16:38:55 ] INFO 2015-02-24 16:38:55: Data processing started...
[ 2015-02-24 16:38:55 ] INFO Read counts file: imported sam/bam/bed files
[ 2015-02-24 16:38:55 ] INFO Conditions: WT, HET
[ 2015-02-24 16:38:55 ] INFO Samples: SCR20, SCR21, SCR22, SCR23
[ 2015-02-24 16:38:55 ] INFO Requested contrasts: WT_vs_HET
[ 2015-02-24 16:38:55 ] INFO Annotation: download
[ 2015-02-24 16:38:55 ] INFO Organism: mm9
[ 2015-02-24 16:38:55 ] INFO Reference source: refseq
[ 2015-02-24 16:38:55 ] INFO Count type: gene
[ 2015-02-24 16:38:55 ] INFO Exon filters: min.active.exons
[ 2015-02-24 16:38:55 ] INFO min.active.exons:
[ 2015-02-24 16:38:55 ] INFO exons.per.gene: 5
[ 2015-02-24 16:38:55 ] INFO min.exons: 2
[ 2015-02-24 16:38:55 ] INFO frac: 0.2
[ 2015-02-24 16:38:55 ] INFO Gene filters: length, avg.reads, expression, biotype
[ 2015-02-24 16:38:55 ] INFO length:
[ 2015-02-24 16:38:55 ] INFO length: 500
[ 2015-02-24 16:38:55 ] INFO avg.reads:
[ 2015-02-24 16:38:55 ] INFO average.per.bp: 100
[ 2015-02-24 16:38:55 ] INFO quantile: 0.25
[ 2015-02-24 16:38:55 ] INFO expression:
[ 2015-02-24 16:38:55 ] INFO median: TRUE
[ 2015-02-24 16:38:55 ] INFO mean: FALSE
[ 2015-02-24 16:38:55 ] INFO quantile: NA
[ 2015-02-24 16:38:55 ] INFO known: NA
[ 2015-02-24 16:38:55 ] INFO custom: NA
[ 2015-02-24 16:38:55 ] INFO biotype:
[ 2015-02-24 16:38:55 ] INFO pseudogene: FALSE
[ 2015-02-24 16:38:55 ] INFO snRNA: FALSE
[ 2015-02-24 16:38:55 ] INFO protein_coding: FALSE
[ 2015-02-24 16:38:55 ] INFO antisense: FALSE
[ 2015-02-24 16:38:55 ] INFO miRNA: FALSE
[ 2015-02-24 16:38:55 ] INFO lincRNA: FALSE
[ 2015-02-24 16:38:55 ] INFO snoRNA: FALSE
[ 2015-02-24 16:38:55 ] INFO processed_transcript: FALSE
[ 2015-02-24 16:38:55 ] INFO misc_RNA: FALSE
[ 2015-02-24 16:38:55 ] INFO rRNA: TRUE
[ 2015-02-24 16:38:55 ] INFO sense_overlapping: FALSE
[ 2015-02-24 16:38:55 ] INFO sense_intronic: FALSE
[ 2015-02-24 16:38:55 ] INFO polymorphic_pseudogene: FALSE
[ 2015-02-24 16:38:55 ] INFO non_coding: FALSE
[ 2015-02-24 16:38:55 ] INFO three_prime_overlapping_ncrna: FALSE
[ 2015-02-24 16:38:55 ] INFO IG_C_gene: FALSE
[ 2015-02-24 16:38:55 ] INFO IG_J_gene: FALSE
[ 2015-02-24 16:38:55 ] INFO IG_D_gene: FALSE
[ 2015-02-24 16:38:55 ] INFO IG_V_gene: FALSE
[ 2015-02-24 16:38:55 ] INFO ncrna_host: FALSE
[ 2015-02-24 16:38:55 ] INFO Filter application: postnorm
[ 2015-02-24 16:38:55 ] INFO Normalization algorithm: edger
[ 2015-02-24 16:38:55 ] INFO Normalization arguments:
[ 2015-02-24 16:38:55 ] INFO method: TMM
[ 2015-02-24 16:38:55 ] INFO logratioTrim: 0.3
[ 2015-02-24 16:38:55 ] INFO sumTrim: 0.05
[ 2015-02-24 16:38:55 ] INFO doWeighting: TRUE
[ 2015-02-24 16:38:55 ] INFO Acutoff: -1e+10
[ 2015-02-24 16:38:55 ] INFO p: 0.75
[ 2015-02-24 16:38:55 ] INFO Statistical algorithm: edger
[ 2015-02-24 16:38:55 ] INFO Statistical arguments:
[ 2015-02-24 16:38:55 ] 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)
[ 2015-02-24 16:38:55 ] INFO Meta-analysis method: none
[ 2015-02-24 16:38:55 ] INFO Multiple testing correction: BH
[ 2015-02-24 16:38:55 ] INFO p-value threshold: 0.05
[ 2015-02-24 16:38:55 ] INFO Logarithmic transformation offset: 1
[ 2015-02-24 16:38:55 ] INFO Quality control plots: mds, biodetection, countsbio, saturation, readnoise, correl, pairwise, boxplot, volcano, biodist, deheatmap
[ 2015-02-24 16:38:55 ] INFO Figure format: png, pdf
[ 2015-02-24 16:38:55 ] INFO Output directory: /data/images/proton/Cereghini/run142/metaseqr_wt_vs_het_mm9_test2
[ 2015-02-24 16:38:55 ] INFO Output data: annotation, p.value, adj.p.value, fold.change, stats, counts
[ 2015-02-24 16:38:55 ] INFO Output scale(s): log2
[ 2015-02-24 16:38:55 ] INFO Output values: normalized
[ 2015-02-24 16:38:55 ] INFO Output statistics: mean, median, sd, mad, cv, rcv
[ 2015-02-24 16:38:55 ] INFO Downloading gene annotation for mm9...
[ 2015-02-24 16:39:10 ] INFO Converting annotation to GenomicRanges object...
[ 2015-02-24 16:39:10 ] INFO Getting DNA sequences...
[ 2015-02-24 16:39:16 ] INFO Getting GC content...
[ 2015-02-24 16:39:18 ] INFO Converting annotation to GenomicRanges object...
[ 2015-02-24 16:39:18 ] INFO Reading bam file sort_uniq.bam for sample with name SCR20. This might take some time...
[ 2015-02-24 16:39:18 ] INFO Reading bam file sort_uniq.bam for sample with name SCR21. This might take some time...
[ 2015-02-24 16:39:18 ] INFO Reading bam file sort_uniq.bam for sample with name SCR22. This might take some time...
[ 2015-02-24 16:39:18 ] INFO Reading bam file sort_uniq.bam for sample with name SCR23. This might take some time...
[ 2015-02-24 16:40:09 ] INFO Counting reads overlapping with given annotation...
[ 2015-02-24 16:40:09 ] INFO ...for single-end reads...
[ 2015-02-24 16:40:09 ] INFO ...using strandedness...
[ 2015-02-24 16:40:36 ] INFO Counting reads overlapping with given annotation...
[ 2015-02-24 16:40:36 ] INFO ...for single-end reads...
[ 2015-02-24 16:40:36 ] INFO ...using strandedness...
[ 2015-02-24 16:41:50 ] INFO Counting reads overlapping with given annotation...
[ 2015-02-24 16:41:50 ] INFO ...for single-end reads...
[ 2015-02-24 16:41:50 ] INFO ...using strandedness...
[ 2015-02-24 16:42:25 ] INFO Counting reads overlapping with given annotation...
[ 2015-02-24 16:42:25 ] INFO ...for single-end reads...
[ 2015-02-24 16:42:25 ] INFO ...using strandedness...
[ 2015-02-24 16:55:37 ] INFO Exporting raw read counts table to /data/images/proton/Cereghini/run142/metaseqr_wt_vs_het_mm9_test2/lists/raw_counts_table.txt.gz
[ 2015-02-24 16:55:37 ] INFO Saving gene model to /data/images/proton/Cereghini/run142/metaseqr_wt_vs_het_mm9_test2/data/gene_model.RData
[ 2015-02-24 16:55:37 ] INFO Removing genes with zero counts in all samples...
[ 2015-02-24 16:55:38 ] INFO Normalizing with: edger
[ 2015-02-24 16:55:38 ] INFO Applying gene filter length...
[ 2015-02-24 16:55:38 ] INFO Threshold below which ignored: 500
[ 2015-02-24 16:55:38 ] INFO Applying gene filter avg.reads...
[ 2015-02-24 16:55:38 ] INFO Threshold below which ignored: 0.130378096479791
[ 2015-02-24 16:55:38 ] INFO Applying gene filter expression...
[ 2015-02-24 16:55:38 ] INFO Threshold below which ignored: 549
[ 2015-02-24 16:55:38 ] INFO Applying gene filter biotype...
[ 2015-02-24 16:55:38 ] INFO Biotypes ignored: rRNA
[ 2015-02-24 16:55:38 ] INFO 10393 genes filtered out
[ 2015-02-24 16:55:38 ] INFO 10701 genes remain after filtering
[ 2015-02-24 16:55:38 ] INFO Running statistical tests with: edger
[ 2015-02-24 16:55:39 ] INFO Contrast: WT_vs_HET
[ 2015-02-24 16:55:40 ] INFO Contrast WT_vs_HET: found 414 genes
[ 2015-02-24 16:55:40 ] INFO Exporting and compressing normalized read counts table to /data/images/proton/Cereghini/run142/metaseqr_wt_vs_het_mm9_test2/lists/normalized_counts_table.txt
[ 2015-02-24 16:55:40 ] INFO Building output files...
[ 2015-02-24 16:55:40 ] INFO Contrast: WT_vs_HET
[ 2015-02-24 16:55:40 ] INFO Adding non-filtered data...
[ 2015-02-24 16:55:40 ] INFO binding annotation...
[ 2015-02-24 16:55:41 ] INFO binding p-values...
[ 2015-02-24 16:55:41 ] INFO binding FDRs...
[ 2015-02-24 16:55:41 ] INFO binding log2 normalized fold changes...
[ 2015-02-24 16:55:41 ] INFO binding normalized mean counts...
[ 2015-02-24 16:55:42 ] INFO binding normalized median counts...
[ 2015-02-24 16:55:43 ] INFO binding normalized count sds...
[ 2015-02-24 16:55:43 ] INFO binding normalized count MADs...
[ 2015-02-24 16:55:45 ] INFO binding normalized count CVs...
[ 2015-02-24 16:55:46 ] INFO binding normalized counts RCVs...
[ 2015-02-24 16:55:48 ] INFO binding normalized mean counts...
[ 2015-02-24 16:55:49 ] INFO binding normalized median counts...
[ 2015-02-24 16:55:50 ] INFO binding normalized count sds...
[ 2015-02-24 16:55:50 ] INFO binding normalized count MADs...
[ 2015-02-24 16:55:52 ] INFO binding normalized count CVs...
[ 2015-02-24 16:55:53 ] INFO binding normalized counts RCVs...
[ 2015-02-24 16:55:56 ] INFO binding all normalized counts for WT...
[ 2015-02-24 16:55:56 ] INFO binding all normalized counts for HET...
[ 2015-02-24 16:55:56 ] INFO Writing output...
[ 2015-02-24 16:55:56 ] INFO Adding filtered data...
[ 2015-02-24 16:55:56 ] INFO binding annotation...
[ 2015-02-24 16:55:56 ] INFO binding p-values...
[ 2015-02-24 16:55:56 ] INFO binding FDRs...
[ 2015-02-24 16:55:56 ] INFO binding log2 normalized fold changes...
[ 2015-02-24 16:55:56 ] INFO binding normalized mean counts...
[ 2015-02-24 16:55:56 ] INFO binding normalized median counts...
[ 2015-02-24 16:55:57 ] INFO binding normalized count sds...
[ 2015-02-24 16:55:58 ] INFO binding normalized count MADs...
[ 2015-02-24 16:55:59 ] INFO binding normalized count CVs...
[ 2015-02-24 16:56:00 ] INFO binding normalized counts RCVs...
[ 2015-02-24 16:56:02 ] INFO binding normalized mean counts...
[ 2015-02-24 16:56:03 ] INFO binding normalized median counts...
[ 2015-02-24 16:56:04 ] INFO binding normalized count sds...
[ 2015-02-24 16:56:04 ] INFO binding normalized count MADs...
[ 2015-02-24 16:56:06 ] INFO binding normalized count CVs...
[ 2015-02-24 16:56:06 ] INFO binding normalized counts RCVs...
[ 2015-02-24 16:56:09 ] INFO binding all normalized counts for WT...
[ 2015-02-24 16:56:09 ] INFO binding all normalized counts for HET...
[ 2015-02-24 16:56:09 ] INFO Writing output...
[ 2015-02-24 16:56:10 ] INFO Creating quality control graphs...
[ 2015-02-24 16:56:10 ] INFO Plotting in png format...
[ 2015-02-24 16:56:10 ] INFO Plotting mds...
[ 2015-02-24 16:56:11 ] INFO Plotting biodetection...
[ 2015-02-24 16:56:12 ] INFO Plotting countsbio...
[ 2015-02-24 16:56:13 ] INFO Plotting saturation...
[ 2015-02-24 16:56:17 ] INFO Plotting readnoise...
[ 2015-02-24 16:56:17 ] INFO Plotting correl...
[ 2015-02-24 16:56:17 ] INFO Plotting pairwise...
[ 2015-02-24 16:56:32 ] INFO Plotting boxplot...
[ 2015-02-24 16:56:33 ] INFO Plotting boxplot...
[ 2015-02-24 16:56:33 ] INFO Plotting volcano...
[ 2015-02-24 16:56:33 ] INFO Contrast: WT_vs_HET
[ 2015-02-24 16:56:34 ] INFO Plotting biodist...
[ 2015-02-24 16:56:34 ] INFO Contrast: WT_vs_HET
[ 2015-02-24 16:56:34 ] INFO Plotting deheatmap...
[ 2015-02-24 16:56:34 ] INFO Contrast: WT_vs_HET
[ 2015-02-24 16:57:01 ] INFO Plotting in pdf format...
[ 2015-02-24 16:57:01 ] INFO Plotting mds...
[ 2015-02-24 16:57:01 ] INFO Plotting biodetection...
[ 2015-02-24 16:57:02 ] INFO Plotting countsbio...
[ 2015-02-24 16:57:03 ] INFO Plotting saturation...
[ 2015-02-24 16:57:06 ] INFO Plotting readnoise...
[ 2015-02-24 16:57:06 ] INFO Plotting correl...
[ 2015-02-24 16:57:06 ] INFO Plotting pairwise...
[ 2015-02-24 16:57:09 ] INFO Plotting boxplot...
[ 2015-02-24 16:57:09 ] INFO Plotting boxplot...
[ 2015-02-24 16:57:09 ] INFO Plotting volcano...
[ 2015-02-24 16:57:09 ] INFO Contrast: WT_vs_HET
[ 2015-02-24 16:57:09 ] INFO Plotting biodist...
[ 2015-02-24 16:57:09 ] INFO Contrast: WT_vs_HET
[ 2015-02-24 16:57:09 ] INFO Plotting deheatmap...
[ 2015-02-24 16:57:09 ] INFO Contrast: WT_vs_HET
[ 2015-02-24 16:57:35 ] INFO Creating HTML report...
[ 2015-02-24 16:57:35 ] 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

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 WT_vs_HET

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 WT_vs_HET

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 WT_vs_HET

Get all the figures in png, pdf format.

Results

DEG table for the contrast WT_vs_HET

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 WT_vs_HET. 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_edger	FDR_edger	log2_normalized_fold_change_WT_vs_HET	log2_normalized_mean_counts_WT	log2_normalized_median_counts_WT	log2_normalized_sd_counts_WT	log2_normalized_mad_counts_WT	log2_normalized_cv_counts_WT	log2_normalized_rcv_counts_WT	log2_normalized_mean_counts_HET	log2_normalized_median_counts_HET	log2_normalized_sd_counts_HET	log2_normalized_mad_counts_HET	log2_normalized_cv_counts_HET	log2_normalized_rcv_counts_HET	log2_normalized_counts_SCR20	log2_normalized_counts_SCR21	log2_normalized_counts_SCR22	log2_normalized_counts_SCR23
chr4	115158804	115168746	NM_007822	0.4146	-	Cyp4a14	protein_coding	1.49e-11	1.02e-07	-2.575494	10.77	10.77	0.158520	1.66e-01	1.47e-02	1.54e-02	8.19	8.19	0.229980	0.241101	2.81e-02	2.94e-02	10.65	10.88	8.03	8.35
chr10	86984539	87046882	NM_008777	0.4290	+	Pah	protein_coding	1.91e-11	1.02e-07	-2.406964	11.46	11.46	0.191030	2.00e-01	1.67e-02	1.75e-02	9.05	9.05	0.157255	0.164860	1.74e-02	1.82e-02	11.59	11.32	8.94	9.16
chr1	164805177	164828843	NM_018881	0.4155	-	Fmo2	protein_coding	1.74e-10	6.22e-07	-2.362501	11.55	11.55	0.012948	1.36e-02	1.12e-03	1.18e-03	9.16	9.16	0.407802	0.427521	4.45e-02	4.67e-02	11.56	11.54	8.87	9.45
chr16	21761357	21787907	NM_023737	0.4381	-	Ehhadh	protein_coding	2.79e-10	7.46e-07	-2.335077	12.46	12.46	0.287033	3.01e-01	2.30e-02	2.42e-02	10.13	10.13	0.183502	0.192376	1.81e-02	1.90e-02	12.66	12.25	10.00	10.26
chr13	23931302	23952583	NM_134069	0.4132	+	Slc17a3	protein_coding	4.66e-10	9.97e-07	-2.569366	10.69	10.69	0.062817	6.59e-02	5.87e-03	6.16e-03	8.07	8.07	0.583790	0.612020	7.23e-02	7.58e-02	10.74	10.65	7.66	8.48
chr15	89164903	89167438	NM_019977	0.5599	+	Miox	protein_coding	4.98e-09	8.88e-06	-2.450686	12.90	12.90	0.391692	4.11e-01	3.04e-02	3.18e-02	10.43	10.43	0.502605	0.526909	4.82e-02	5.05e-02	13.17	12.62	10.79	10.08
chr18	21230844	21258700	NM_008586	0.4271	+	Mep1b	protein_coding	6.40e-09	9.79e-06	-2.518780	11.02	11.02	0.652251	6.84e-01	5.92e-02	6.20e-02	8.53	8.53	0.501221	0.525459	5.87e-02	6.16e-02	11.48	10.56	8.18	8.89
chr6	24038282	24118092	NM_019481	0.3808	-	Slc13a1	protein_coding	7.65e-09	1.02e-05	-2.893350	10.18	10.18	0.276929	2.90e-01	2.72e-02	2.85e-02	7.22	7.22	0.697437	0.731163	9.66e-02	1.01e-01	9.98	10.37	6.73	7.71
chr6	85780380	85781884	NM_023455	0.5070	-	Nat8	protein_coding	9.69e-09	1.15e-05	-2.988087	9.81	9.81	0.225103	2.36e-01	2.29e-02	2.40e-02	6.76	6.76	0.700461	0.734333	1.04e-01	1.09e-01	9.65	9.97	7.26	6.27
chr4	129886352	129899539	NM_001034872	0.4952	-	Gm853	protein_coding	1.57e-08	1.68e-05	-2.661441	10.10	10.10	0.188871	1.98e-01	1.87e-02	1.96e-02	7.42	7.42	0.467207	0.489799	6.30e-02	6.60e-02	10.24	9.97	7.75	7.09
chr7	51480806	51484987	NM_010639	0.5182	+	Klk1	protein_coding	2.41e-08	2.35e-05	-2.171737	10.45	10.45	0.034361	3.60e-02	3.29e-03	3.45e-03	8.26	8.26	0.355088	0.372259	4.30e-02	4.51e-02	10.42	10.47	8.01	8.51
chr9	123502937	123541683	NM_011731	0.4345	-	Slc6a20b	protein_coding	4.08e-08	3.44e-05	-2.156374	10.43	10.43	0.104436	1.09e-01	1.00e-02	1.05e-02	8.25	8.25	0.379464	0.397813	4.60e-02	4.82e-02	10.50	10.35	7.98	8.52
chr4	19745873	19769360	NM_028746	0.3922	+	Slc7a13	protein_coding	4.18e-08	3.44e-05	-2.106708	10.39	10.39	0.127254	1.33e-01	1.23e-02	1.28e-02	8.28	8.28	0.229845	0.240959	2.78e-02	2.91e-02	10.30	10.48	8.44	8.11
chr1	129625989	129664141	NM_001033041	0.4242	+	Acmsd	protein_coding	5.40e-08	4.12e-05	-2.700908	9.59	9.59	0.027770	2.91e-02	2.90e-03	3.04e-03	6.87	6.87	0.442168	0.463550	6.44e-02	6.75e-02	9.61	9.57	7.18	6.55
chr7	126705221	126744208	NM_146197	0.4246	+	Acsm2	protein_coding	5.77e-08	4.12e-05	-2.123691	13.69	13.69	0.286018	3.00e-01	2.09e-02	2.19e-02	11.58	11.58	0.190874	0.200104	1.65e-02	1.73e-02	13.49	13.89	11.44	11.71
chr19	56362427	56394579	NM_146101	0.4806	+	Habp2	protein_coding	1.06e-07	7.07e-05	-1.867103	11.08	11.08	0.235213	2.47e-01	2.12e-02	2.23e-02	9.23	9.23	0.057920	0.060721	6.28e-03	6.58e-03	11.25	10.92	9.27	9.18
chr8	95725116	95753518	NM_133660	0.4213	-	Ces1e	protein_coding	1.32e-07	8.34e-05	-2.824108	9.14	9.14	0.043417	4.55e-02	4.75e-03	4.98e-03	6.24	6.24	0.726175	0.761291	1.16e-01	1.22e-01	9.11	9.17	5.73	6.75
chr10	127286983	127298232	NM_147222	0.4302	+	Rdh19	protein_coding	1.64e-07	9.74e-05	-2.081582	10.91	10.91	0.097381	1.02e-01	8.92e-03	9.36e-03	8.77	8.77	0.629297	0.659727	7.18e-02	7.53e-02	10.98	10.84	8.32	9.21
chr3	129539656	129578246	NM_007686	0.4327	+	Cfi	protein_coding	3.27e-07	1.84e-04	-1.992356	10.25	10.25	0.154419	1.62e-01	1.51e-02	1.58e-02	8.27	8.27	0.042994	0.045073	5.20e-03	5.45e-03	10.36	10.14	8.30	8.24
chr19	8410112	8479595	NM_177002	0.3896	-	Slc22a30	protein_coding	4.43e-07	2.37e-04	-2.668947	11.08	11.08	0.231647	2.43e-01	2.09e-02	2.19e-02	8.14	8.14	1.316639	1.380307	1.62e-01	1.70e-01	10.91	11.24	7.21	9.07
chr9	119057160	119066211	NM_146230	0.5579	-	Acaa1b	protein_coding	6.56e-07	3.20e-04	-1.791768	11.61	11.61	0.002607	2.73e-03	2.24e-04	2.35e-04	9.80	9.80	0.357565	0.374856	3.65e-02	3.83e-02	11.61	11.61	10.05	9.55
chr7	31278676	31298421	NM_019546	0.4983	+	Prodh2	protein_coding	6.57e-07	3.20e-04	-1.777856	12.39	12.39	0.039587	4.15e-02	3.19e-03	3.35e-03	10.61	10.61	0.248548	0.260567	2.34e-02	2.46e-02	12.37	12.42	10.78	10.43
chr6	133799872	133803685	NM_010594	0.4156	-	Kap	protein_coding	7.21e-07	3.35e-04	-2.396087	14.73	14.73	0.115785	1.21e-01	7.86e-03	8.24e-03	12.18	12.18	0.969823	1.016720	7.96e-02	8.35e-02	14.65	14.81	11.50	12.87
chr13	20885981	21116123	NM_012054	0.4086	+	Aoah	protein_coding	7.98e-07	3.56e-04	-1.948578	10.30	10.30	0.351026	3.68e-01	3.41e-02	3.57e-02	8.37	8.37	0.215239	0.225647	2.57e-02	2.70e-02	10.06	10.55	8.22	8.52
chr11	61486283	61534301	NM_001033227	0.5335	-	Slc5a10	protein_coding	1.04e-06	4.47e-04	-1.744532	10.93	10.93	0.024659	2.59e-02	2.26e-03	2.37e-03	9.17	9.17	0.315018	0.330251	3.44e-02	3.60e-02	10.91	10.94	9.39	8.94
chr17	34993318	34999459	NM_008198	0.5199	-	Cfb	protein_coding	1.23e-06	5.08e-04	-1.677997	10.91	10.91	0.076978	8.07e-02	7.06e-03	7.40e-03	9.23	9.23	0.042406	0.044457	4.59e-03	4.82e-03	10.96	10.85	9.26	9.20
chr2	165230794	165298697	NM_054055	0.5036	-	Slc13a3	protein_coding	1.67e-06	6.63e-04	-1.755167	11.33	11.33	0.094173	9.87e-02	8.31e-03	8.71e-03	9.55	9.55	0.438697	0.459911	4.60e-02	4.82e-02	11.27	11.40	9.86	9.24
chr19	7747565	7785800	NM_144785	0.3942	-	Slc22a19	protein_coding	2.14e-06	8.19e-04	-1.893913	11.65	11.65	0.211206	2.21e-01	1.81e-02	1.90e-02	9.70	9.70	0.645419	0.676630	6.66e-02	6.98e-02	11.80	11.50	9.24	10.15
chr10	61533769	61535319	NR_002687	0.5455	+	Gm5424	retrotransposed	2.73e-06	1.01e-03	-1.704801	11.74	11.74	0.197912	2.07e-01	1.69e-02	1.77e-02	10.03	10.03	0.296859	0.311214	2.96e-02	3.10e-02	11.60	11.88	10.24	9.82
chr7	135409210	135415947	NM_133254	0.5347	+	Slc5a2	protein_coding	3.85e-06	1.37e-03	-1.826139	10.88	10.88	0.569040	5.97e-01	5.23e-02	5.49e-02	9.09	9.09	0.315090	0.330326	3.47e-02	3.63e-02	10.47	11.28	8.87	9.31
chr17	45767152	45796258	NM_001013767	0.4818	-	Capn11	protein_coding	4.10e-06	1.41e-03	4.428360	4.05	4.05	0.489231	5.13e-01	1.21e-01	1.27e-01	7.47	7.47	2.624582	2.751498	3.51e-01	3.68e-01	3.70	4.39	9.33	5.61
chr19	42800347	42827265	NM_029011	0.5115	-	Pyroxd2	protein_coding	4.31e-06	1.44e-03	-1.561990	11.44	11.44	0.043041	4.51e-02	3.76e-03	3.95e-03	9.88	9.88	0.013023	0.013653	1.32e-03	1.38e-03	11.47	11.41	9.89	9.87
chr15	9265352	9300625	NM_144845	0.3925	+	Ugt3a2	protein_coding	5.21e-06	1.69e-03	-1.833567	10.38	10.38	0.097960	1.03e-01	9.44e-03	9.89e-03	8.50	8.50	0.528374	0.553925	6.21e-02	6.51e-02	10.45	10.31	8.13	8.88
chr6	128612402	128627392	NM_053165	0.3971	+	Clec2h	protein_coding	5.40e-06	1.70e-03	-2.274033	9.21	9.21	0.077385	8.11e-02	8.40e-03	8.81e-03	6.91	6.91	0.462053	0.484397	6.69e-02	7.01e-02	9.16	9.27	6.58	7.24
chr13	23962141	23987599	NM_001170638	0.4083	+	Slc17a1	protein_coding	6.16e-06	1.88e-03	-1.773370	12.49	12.49	0.231490	2.43e-01	1.85e-02	1.94e-02	10.67	10.67	0.559412	0.586464	5.24e-02	5.50e-02	12.32	12.65	10.27	11.07
chr9	44607482	44611034	NM_138951	0.5283	-	Ttc36	protein_coding	6.35e-06	1.89e-03	-1.665374	12.04	12.04	0.013853	1.45e-02	1.15e-03	1.21e-03	10.34	10.34	0.428420	0.449137	4.14e-02	4.34e-02	12.05	12.03	10.64	10.04
chr7	126761341	126806029	NM_054094	0.4148	+	Acsm1	protein_coding	7.06e-06	2.04e-03	-1.531271	11.43	11.43	0.092568	9.70e-02	8.10e-03	8.49e-03	9.90	9.90	0.078058	0.081833	7.89e-03	8.27e-03	11.49	11.36	9.95	9.84
chrY	234230	280254	NM_011419	0.3756	+	Kdm5d	protein_coding	7.94e-06	2.24e-03	6.081262	3.68	3.68	0.316401	3.32e-01	8.59e-02	9.01e-02	7.59	7.59	4.330249	4.539644	5.71e-01	5.98e-01	3.91	3.46	10.65	4.52
chr9	119102095	119118223	NM_133980	0.4909	-	Slc22a13	protein_coding	8.32e-06	2.24e-03	-1.620393	10.79	10.79	0.141152	1.48e-01	1.31e-02	1.37e-02	9.15	9.15	0.347793	0.364611	3.80e-02	3.98e-02	10.69	10.89	9.40	8.91
chr17	34999546	35019045	NM_013484	0.4733	-	C2	protein_coding	8.38e-06	2.24e-03	-2.117734	9.55	9.55	0.064164	6.73e-02	6.72e-03	7.05e-03	7.42	7.42	0.296916	0.311274	4.00e-02	4.20e-02	9.50	9.59	7.21	7.63
chr13	93521540	93564773	NM_011582	0.4601	-	Thbs4	protein_coding	9.04e-06	2.36e-03	4.004106	8.23	8.23	0.697528	7.31e-01	8.47e-02	8.88e-02	11.26	11.26	2.776207	2.910455	2.47e-01	2.58e-01	7.74	8.73	13.23	9.30

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

Download the raw read counts table for the experiment.
Download the normalized read counts table for the experiment.

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.
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.
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.
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.