RNA-seq_code.sh