ParBLiSS / FastANI

Fast Whole-Genome Similarity (ANI) Estimation
Apache License 2.0
374 stars 67 forks source link

Matrix output does not equal default output #60

Closed mdiorio371 closed 4 years ago

mdiorio371 commented 4 years ago

The values presented in the matrix format are not the same as the default output.

Using the command: fastANI --rl test_refs.txt --ql test_refs.txt -o out_test --matrix -t 4

I get out_test: CP003047.1.txt CP003047.1.txt 100 1543 1545 CP003047.1.txt CP006876.1.txt 98.4822 1455 1545 CP003047.1.txt CP004027.1.txt 98.2008 1410 1545 CP004027.1.txt CP004027.1.txt 100 1541 1544 CP004027.1.txt CP006876.1.txt 98.33 1417 1544 CP004027.1.txt CP003047.1.txt 98.1959 1414 1544 CP006876.1.txt CP006876.1.txt 100 1587 1587 CP006876.1.txt CP003047.1.txt 98.4951 1464 1587 CP006876.1.txt CP004027.1.txt 98.2722 1435 1587

and out_test.matrix: CP003047.1.txt CP004027.1.txt 98.198334 CP006876.1.txt 98.488670 98.301117

The values are close but disagree in the 10^-3 range, and the difference does not seem to come up from rounding. Should one value be considered more reliable?

cjain7 commented 4 years ago

Please see #36

mdiorio371 commented 4 years ago

Oh I see, it's the mean value. Thank you!