[MRG] Adding cANI Rust implementation to PR #2943

[mr-eyes]

Rust implementation for ANI estimation

Usage example:

use crate::c_ani::containment_to_distance;

fn main() {
    let containment = 0.1;
    let ksize = 31;
    let scaled = 1000; // u64
    let n_unique_kmers = 1000;
    let sequence_len_bp = n_unique_kmers * scaled;
    let confidence = Some(0.95);
    let estimate_ci = Some(true);
    let prob_threshold = Some(1e-3);

    let result = containment_to_distance(
        containment,
        ksize,
        scaled,
        Some(n_unique_kmers),
        Some(sequence_len_bp),
        confidence,
        estimate_ci,
        prob_threshold,
    );


    match result {
        Ok(ci_result) => {
            let ani_result = 1.0 - ci_result.point_estimate;
            println!("ANI: {}", ani_result);
        }

        Err(e) => println!("Error occurred: {:?}", e),
    }
}

[codecov]

Codecov Report

Attention: 56 lines in your changes are missing coverage. Please review.

Comparison is base (01adde3) 92.85% compared to head (fb03dc6) 58.29%.

❗ Current head fb03dc6 differs from pull request most recent head eb8731a. Consider uploading reports for the commit eb8731a to get more accurate results

Files	Patch %	Lines
src/core/src/c_ani.rs	0.00%	56 Missing ⚠️

Additional details and impacted files

@@                  Coverage Diff                  @@
##           upd-gatherresult    #3005       +/-   ##
=====================================================
- Coverage             92.85%   58.29%   -34.57%     
=====================================================
  Files                   104       32       -72     
  Lines                 12527     2918     -9609     
  Branches               2624        0     -2624     
=====================================================
- Hits                  11632     1701     -9931     
- Misses                  594     1217      +623     
+ Partials                301        0      -301     

Flag	Coverage Δ
hypothesis-py	?
python	?
rust	58.29% <0.00%> (?)

Flags with carried forward coverage won't be shown. Click here to find out more.

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

[bluegenes]

hi @mr-eyes if you change the branch to 'use-upd-gatheresult' i can go ahead and merge into my PR.

cheers for setting up - this ensures you get credit for your contributions :)

[mr-eyes]

hi @mr-eyes if you change the branch to 'use-upd-gatheresult' i can go ahead and merge into my PR.

cheers for setting up - this ensures you get credit for your contributions :)

upd-gatherresult is now the base. Thank you!

[bluegenes]

thx!

[mr-eyes]

thx!

You are most welcome!