RNA-seq data processing and quality check
Raw data (raw reads) of fastq format were firstly processed through in-house perl scripts. At the same time, Q20, Q30, and GC content the clean data were calculated. All the downstream analyses were based on the clean data with high quality. Reference genome and gene model annotation files were downloaded from genome website. Hisat2 v2.0.5 was used to build the index of the reference genome and to align paired-end clean reads to the reference genome.