script for working out directionality in relation to RNAseq reads