Differentially gene expression analysis with RNA-seq data is quite common nowadays, and there are pretty good Bioconductor packages for that: limma::voom, DESeq2 …
The code for that part is quite simple, being super quick to get a list of de-regulated genes.