script for RNASeq analysis