I need to do the mapping of RNAseq reads on mouse genome reference to find the miRNAs count. My data contains reads length 20-50nt, but miRNAs has about 22nt, so I wanna map just reads about this length.
What is the best:
- To filter (extract of the data) just reads 22nt
- To cut reads bigger than 22nt to obtain just reads with this length?