Batch verify hairpin structure
I downloaded hairpin and mature sequence from miRBase and predicted the possible locations of some miRNAs in my genome. Now I want to verify the hairpin structure to filter these locations. RNAfold can give the results, but it can only verify them one by one. I want to ask if there is any software that can batch filter these locations? Thank you! enter image description here
• 1,629 views
•
link
written
by
153348734 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Filter out miRNA from ncRNA dataset
written by binaryCode •
-
Differential expression pattern of known miRNAs
written by Rey •
-
Several questions about miRDeep2
written by 2822462298 •
-
smallRNA sequencing workflow and Database usage
written by Konstantinos Yeles
-
Retrieve genomic coordinates (locus start, locus end) from blastn hits?
written by sombrajo92 •
-
mirDeep2 using bowtie vs. bwa - why do more aligned reads yield less miRNA
written by Michael Dondrup
What do you mean by "verify" and by that the
RNAfoldcan only verify then "one by one". TheRNAfoldwill happily predict thousands of structures for thousands of your input sequences at once.What is acceptable for you as "hairpins" is entirely different problem though. I'm not aware of any software that would classify structures as "hairpin" or "other", although, there is set of tools called
RNAshapesand one of the functions (I thinkRapidShapes) is to extract abstract shapes (for example((()))would be[]if I remember correctly.Depending on what is acceptable for you as hairpin, you might be better of using custom script, that e.g. checks whether there is any
(character occurring after the first)character in the dot bracket notation string (or vice versa). Beware, that this may hold true for mostly unpaired structures (even with only one or two basepairs), which is presumably not what you want, so setting some minimum number of base-pairs may be in order.Hope this helps
Thank you! I'm a novice and don't know much about miRNA.