How to run quantifier from miRDeep2 (miRNA seqs from MirGeneDB)
0
0
Entering edit mode
16 months ago
Apex92 ▴ 320

Dear all,

I want to perform miRNA quantification using quantifier.pl from miRDeep2. As default miRDeep2 is compatible with miRBase but I have to perform the quantification with mature and precursor sequences (fasta format) of miRNAs from MirGeneDB. I ran the tool as below:

quantifier.pl -p pre_miRNAs.fa -m mature_miRNAs.fa -r reads.fa

I do not get the expected output. I get these messages:

0 mature mappings to precursors

and

Warning: 0 mature sequences mapped to any of your given precursor sequences

The output file miRNAs_expressed_all_samples_1691413314.csv is empty considering that I have the mapping of my reads to provided miRNA sequences. I also do find mature sequences in precursor sequences when I search them using the grep function in bash.

Could someone please help me to solve the problem? Thank you.

quantifier alignment mapping mirdeep2 • 629 views
ADD COMMENT
0
Entering edit mode

I solved this issue. The reason was that the precursor IDs did not match the mature IDs so I removed the "_pre" and the script worked as expected.

ADD REPLY

Login before adding your answer.

Traffic: 1870 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6