Entering edit mode
6.0 years ago
glady
▴
320
Hello all, I had some miRNA samples, I performed the mapping of the miRNA reads to miRBase mature file by BWA. And the raw counts were generated by the following command:
samtools view -F 260 sample_1.sam | \
cut -f 3 | sort | uniq -c | \
awk '{printf("%s\t%s\n", $2, $1)}' > sample_1-counts.txt
How can I obtain the RPKM values from these raw counts? I mean, how can I normalize these raw counts? Or is there any other way like some tools which can be used for performing quantification on the mirna.sam file?
Thank you.
I think you are mistakenly mentioned as mRNA samples. I hope that is miRNA, am i corrrect? If i am, here is the code for RPKM.
Yes, you are correct. It is miRNA. Thank you for the reply.
I'm getting an error in the above code:
This normalisation which i have mentioned is fully intrasample based. So, you should calculate the tota readcountl of perticular miRNA from across all the samples. Ex: If you have 5 samples, you should calculate the total readcount of perticular miRNA from all 5 samples. Do like this for all miRNAs and add in the separate column titled as Total in your sample_1-counts file.
Is it correct that you have mRNA samples but quantify miRNAs? Does the library prep or rather RNA extraction properly capture the miRNAs? For RPKM, please check this video. Tools that can do that can be found on google.
Its miRNA and not MRNA, my typing mistake. Sorry, I have edited the post. Thank you