Hello friends,

Based on this post Question: Differential gene expression based on read counts using DESeq package.

Using htseq_count, I generated the raw read count for each gene for my cancer and normal sample. Then I used DESeq and DESeq2 package to identify differential gene expression.

1) First 4 rows are from DESeq output

2) Next 4 rows are from DESeq 2 output

3) Then followed by HTseq_count output

Just to understand, I manually computed the log2(fold_change) for both DESeq and DESeq2 R package output.

• For DESeq output, I first divided (value2/value1) and then took log(value2/value1,2). I received the same value as computed by DESeq.

• But for DESeq2 output, I divided (val2/val1) and then took log(val2/val1,2). The DESeq2 computed log2 fold change is different from mine. Am I missing something?