I am looking to identify a stem-loop DNA secondary structure in genomic DNA (not RNA) sequences.
This structure is usually ~ 10nt from the 3' terminus, and it is usually ~15-20nt long.
There seem to be > a dozen options to predict secondary structure in DNA (like at this page: https://omictools.com/dna-structures-category)
So I'm seeking advice forum members' advice for what would satisfy the following criteria:
command line use,
fast & high throughput to process entire genomes,
prefer multi-threaded option to speed up runs,
will return palindrome and sequence of stem-loop,
prefer it will allow palindrome length as pre-filter for the run, (though post-processing of results is OK, but not preferred),
prefer it will allow palindrome position with respect to DNA sequence 3' end, (though post-processing of results is OK, but not preferred), and
returns score of predicted secondary structure, in case I need / want to use that info later.
And anything else I might missed out, which you think would be important for my goals.
I've never predicted DNA secondary structure, much less used software like UNAfold or Mfold, or EMBOSS's 'einverted' or 'palindrome'. So, I look forward to your responses. Thanks!