<div dir="ltr"><div>I'm not sure what hit score is and whether it even applies to modern BLAST.<br></div><div><br></div><div>According to this
(<a href="http://etutorials.org/Misc/blast/Part+III+Practice/Chapter+7.+A+BLAST+Statistics+Tutorial/7.1+Basic+BLAST+Statistics/">http://etutorials.org/Misc/blast/Part+III+Practice/Chapter+7.+A+BLAST+Statistics+Tutorial/7.1+Basic+BLAST+Statistics/</a>),BLAST applies a lot of statistics to HSP scores/evalues, considering
database size, sequence composition and God knows what else. That book is
pretty old (2003), so maybe it's BLAST1?</div><div>Judging by the BLAST wikipedia page (<a href="https://en.wikipedia.org/wiki/BLAST#Algorithm">https://en.wikipedia.org/wiki/BLAST#Algorithm</a>), BLAST2 merges whatever it can into a single HSP and gives the score of that alignment as an HSP score.</div><div><br></div><div>And if HSPs cannot be merged, eg with repetitive sequences, the idea of a hit score doesn't make much sense anyway. I've just checked and with incompatible HSPs web-blast does not report any hit-wide statistics. The hits are just sorted by the best HSP, no matter whether the second one is about as good as the best or merely noticeable. You can see it with this sequence (<a href="https://gist.github.com/SynedraAcus/690870ff00bf00dd832a635fe0652f81">https://gist.github.com/SynedraAcus/690870ff00bf00dd832a635fe0652f81</a>) against protein nr. And it does allow pretty sizeable indels, on the order of hundreds of nucleotides at least.<br></div><div><br></div><div>So in most cases with BLAST2 you can just take the best HSP (most likely only one) and consider it done. If you expect alignments that can't be merged (large repeats/rearrangements or large indels/tracts of nonhomologous sequence), you probably need some specific statistics to cover that.<br></div><div><br></div><div>Also, out of pure curiosity: what do you use raw scores for? Some kind of distance calculations? I don't think I've ever seen them used.<br></div><div></div><div><br></div><br></div><div class="gmail_extra"><br><div class="gmail_quote">2018-01-11 20:15 GMT+08:00 Adam Sjøgren <span dir="ltr"><<a href="mailto:asjo@koldfront.dk" target="_blank">asjo@koldfront.dk</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"> Hi,<br>
<span class=""><br>
<br>
Peter writes:<br>
<br>
> Each of the HSPs will have its own raw score - it that not what you want?<br>
<br>
</span>To be honest: I wasn't sure; I am porting pieces of old code, trying to<br>
match it as much as possible, more or less with my head under my arm.<br>
<br>
I see now that we use hsp.score in other similar places, so I think I<br>
can use hsp.score and be done with it.<br>
<br>
(It is a different value, though, right?)<br>
<span class=""><br>
> Also you might find Bio.SearchIO more future proof over using Bio.Blast<br>
> (the later offers a more generic interface covering HMMER etc as well).<br>
<br>
</span>Thanks for the hint, I'll have to look into why my colleague chose<br>
Bio.Blast.<br>
<br>
<br>
Best regards,<br>
<br>
Adam<br>
<span class="HOEnZb"><font color="#888888"><br>
--<br>
"That's one of the remarkable things about life. It's Adam Sjøgren<br>
never so bad that it can't get worse." <a href="mailto:asjo@koldfront.dk">asjo@koldfront.dk</a><br>
</font></span><div class="HOEnZb"><div class="h5"><br>
______________________________<wbr>_________________<br>
Biopython mailing list - <a href="mailto:Biopython@mailman.open-bio.org">Biopython@mailman.open-bio.org</a><br>
<a href="http://mailman.open-bio.org/mailman/listinfo/biopython" rel="noreferrer" target="_blank">http://mailman.open-bio.org/<wbr>mailman/listinfo/biopython</a></div></div></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><font face="arial,helvetica,sans-serif">Alexey Morozov,<br>LIN SB RAS, bioinformatics group.<br>Irkutsk, Russia.<br></font></div>
</div>