The sequencing data were first assessed for quality control using FastQC (v.0.11.5, Babraham Bioinformatics, Cambridge, UK). Then all sequenced libraries were mapped to the human genome (GRCh38/hg38) using the STAR RNA-seq aligner (v.2.5) (Dobin et al., 2013) with the following parameter: “--outSAMmapqUnique 60”. The reads distribution across the genome was assessed using bamutils (from ngsutils v.0.5.9) (Breese & Liu, 2013). Uniquely mapped sequencing reads were assigned to GRCh38/hg38 refGene genes using featureCounts (subread v.1.5.1) (Liao, Smyth, Shi, 2014) with the following parameters: “−s 2 –p –Q 10”.
