Length of RNAseq reads
Hi,
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?
Thanks
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
TPM: counts for transcripts non included in ENCODE?
written by R •
-
Best aligner to identify single point mutation
written by silas008
-
STAR error: EXITING because of FATAL ERROR in reads input: short read sequence line: 0
written by silas008
-
What type of GTF file should I use for miRNAseq analysis?
written by silas008
-
Bowtie2 mapped just 50% of the reads
written by silas008
-
Extracting reads belonging to different RNA families from Rfam
written by Tm
-
Mapping miRNA with Bowtie & Bowtie2
written by gkuffel22 •
-
De Novo assembly/mapping
written by Juan Ledesma •
-
Problem with bad quality bases in RNAseq
written by silas008
-
Assessing The Quality Of Mirna Rnaseq Data
written by Ashutosh Pandey
I suggest you to first perform adapter trimming as they are the most likely source of the additional length. As for read length filtering, a simple perl script or awk script can easily achieve what you would like to do.
Actually, I would be more flexible about read length. Previous posts (How can I accurately identify miRNA sequences from small RNA seq results) suggest to keep reads between 18 and 30.