Bug: Bio::SearchDist (or the underlying
Bio::Ext::Align) reports incorrect E-values.
BioPerl version: 1.4
Description of the problem:
Bio::SearchDist does not provide meaningful E-values when trying to fit an extreme value distribution (EVD) to a set of integer scores.
Test script and files:
The test input file below contains a sorted set of 20724 integer scores. The scores are distributed around a median of zero with a range between -204 and +1073.
'test_evd.pl' is a simple test script to fit EVD to these scores and calculate E-value for each score. The fit_evd() method seems to complete successfully (it returns 1). However, in the output, the first 29 "E-values" are 1075530112, the 30th is 1, and all others are zero.
Links to script and input/output files: