After mapping, the Python-based software HTseq-count was used to extract raw counts of annotated genes. Differentially expressed genes were identified using the R package edgeR, which uses counts per ...