HI again Reza, I have a few suggestions for enhancements of the output of miRador:
The file 'precursors.fa' matches the files 'preAnnotatedCandidates.csv' and 'preAnnotatedCandidates.fa'. However, in the README, it states that the files 'finalAnnotatedCandidates.csv' and 'finalAnnotatedCandidates.fa' are the definitive output (after annotation). The names of the entries can change if one compares the 'finalAnnotatedCandidates' relative to the 'preAnnotatedCandidates'. Also, the number of entries diverges from the number found in the 'preAnnotatedCandidates'. It would be very useful if a "final" file of precursor sequences were output. At present, there is no file of precursor FASTA sequences that matches the 'final' output.
The file 'precursors.fa' always shows the top (+) genomic strand, even when the MIRNA locus is on the (-) genomic strand. Because MIRNA precursors are single-stranded, I suggest doing the reverse-comps for loci that are on the (-) genomic strand.
Related to the above, it would be helpful in many circumstances if miRador would list the genomic coordinates of the precursors. At present, only the genomic coordinates of the mature miRNA and the miRNA* are shown (in the 'finalAnnotatedCandidates.csv' file).
HI again Reza, I have a few suggestions for enhancements of the output of miRador:
I hope these suggestions are helpful!
Mike