added read RNAseq function