Exonerate raw score to percent identity
1
0
Entering edit mode
2.5 years ago

Does anyone know what exactly does line “raw score” in exonerate output mean?

For example I attach here a screenshot of a full hit, where the length of a query is 74bp, yet the raw score is 370, which is much greater than 74. The original command was:

exonerate --verbose 0 -E --model affine:bestfit --percent 20 --singlepass 1 --showtargetgff 1 query.fasta target.fasta > exonerate_out.out

Can you explain what exactly raw score means, and how do I transform it to percent identity?

screenshot of an output

exonerate SmithWaterman alignment • 917 views
ADD COMMENT
2
Entering edit mode
2.5 years ago
Carambakaracho ★ 2.9k

From exonerate's user manual:

Alignment Scores

Currently, only the raw alignment scores are displayed. This score just is the sum of transistion scores used in the dynamic programming. For example, in the case of a Smith-Waterman alignment, this will be the sum of the substitution matrix scores and the gap penalties.

Conversion of the raw score is unnecessarily complicated, but with --ryo roll your own option you can somehow use

--ryo
%p[is]
Percent {id,similarity} over the equivalenced portions of the alignment. (ie. %pi == 100*(%ei / %et))

It's been ages since I used exonerate, last memory is me having troubles compiling the 2.4.0 binary (which might be completely wrong)

ADD COMMENT
0
Entering edit mode

Thank you, that helps - i forgot about ryo option because i was using biopython's wrapper to parse the output. Since we are here - do you know if i am able to access ryo format options from biopython?

ADD REPLY

Login before adding your answer.

Traffic: 1192 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