Biostar Beta. Not for public use.
Exonerate, How To Search A Whole Genome For A Gene, Including Mismatches?
Entering edit mode
9.3 years ago
John • 750

I used to be able to do this:

exonerate whole _genome.fa gene_ to_search.fa

And it would return exact matches or close matches and show alignments. I tried this again recently, and it only returns exact matches. Did the syntax change?

Thanks, John

dna alignment • 1.9k views
Entering edit mode
9.3 years ago

Exonerate still is pretty much the same. But, you're using your whole genome as a query. This is useful only when you have a LOT of thing to search and wants to reduce memory IO. Exonerate always return the best match unless you explicitly specify a "suboptimal threshold" (e. g. using --percent).

Anyway, there's a possibility that your case to be an output artifact related to the geneseed threshold, fixed in the 2.0.0 -> 2.1.0 transition. Just check the changelog at Exonerate's source.

Entering edit mode
9.3 years ago
ALchEmiXt ♦ 1.9k

Google returned me the man page here. You can/should specify -q and -t. But why use this tool? When fastas grow big I am almost sure blast+ to tabular output or limiting its output to the number of hits you expect will do a faster multi-thraded job directly on fastas as well.

My 2ct


Login before adding your answer.

Similar Posts
Loading Similar Posts
Help About
Access RSS

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

Powered by the version 2.3.3