ragc-core 0.1.1

Core compression and decompression algorithms for the AGC genome compression format
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94

#![allow(clippy::all)]
// Integration test for splitter identification
// Should produce output identical to C++ test_splitters

use ahash::AHashSet;
use ragc_core::determine_splitters;

fn print_splitters(splitters: &AHashSet<u64>, max_count: usize) {
    println!("Count: {}", splitters.len());
    print!("Splitters: ");

    // Convert to sorted vector for consistent output
    let mut sorted: Vec<u64> = splitters.iter().copied().collect();
    sorted.sort_unstable();

    for (i, &s) in sorted.iter().enumerate().take(max_count) {
        if i > 0 {
            print!(" ");
        }
        print!("{s:x}");
    }
    println!();
}

fn main() {
    // Test 1: Two contigs with truly different k-mers
    println!("# Test 1: Simple case (AAAC, GGGT)");
    {
        let contigs = vec![
            vec![0, 0, 0, 1], // AAAC
            vec![2, 2, 2, 3], // GGGT
        ];

        let (splitters, _singletons, _duplicates) = determine_splitters(&contigs, 3, 100);
        print_splitters(&splitters, 50);
    }

    // Test 2: Identical contigs (no singletons)
    println!("\n# Test 2: No singletons (identical contigs)");
    {
        let contigs = vec![vec![0, 1, 2, 3], vec![0, 1, 2, 3]];

        let (splitters, _singletons, _duplicates) = determine_splitters(&contigs, 3, 100);
        print_splitters(&splitters, 50);
    }

    // Test 3: All unique k-mers
    println!("\n# Test 3: All unique (AAAA, CCCC)");
    {
        let contigs = vec![
            vec![0, 0, 0, 0], // AAAA
            vec![1, 1, 1, 1], // CCCC
        ];

        let (splitters, _singletons, _duplicates) = determine_splitters(&contigs, 3, 100);
        print_splitters(&splitters, 50);
    }

    // Test 4: Mixed scenario
    println!("\n# Test 4: Mixed scenario");
    {
        let contigs = vec![
            vec![0, 1, 2, 3, 0, 1, 2, 3, 0, 1, 2, 3], // ACGTACGTACGT (repeating)
            vec![0, 0, 0, 0, 1, 1, 1, 1, 2, 2, 2, 2], // AAAACCCCGGGG (unique patterns)
        ];

        let (splitters, _singletons, _duplicates) = determine_splitters(&contigs, 3, 100);
        print_splitters(&splitters, 50);
    }

    // Test 5: Longer realistic example
    println!("\n# Test 5: Longer sequences");
    {
        let contigs = vec![
            vec![0, 1, 2, 3, 0, 1, 2, 3, 0, 1, 2, 3, 0, 1, 2, 3, 0, 1, 2, 3], // ACGTACGTACGTACGTACGT
            vec![3, 2, 1, 0, 3, 2, 1, 0, 3, 2, 1, 0, 3, 2, 1, 0, 3, 2, 1, 0], // TGCATGCATGCATGCATGCA
            vec![0, 0, 1, 1, 2, 2, 3, 3, 0, 0, 1, 1, 2, 2, 3, 3, 0, 0, 1, 1], // AACCGGTTAACCGGTTAACC
        ];

        let (splitters, _singletons, _duplicates) = determine_splitters(&contigs, 5, 100);
        print_splitters(&splitters, 100);
    }

    // Test 6: Single contig (all singletons)
    println!("\n# Test 6: Single contig");
    {
        let contigs = vec![
            vec![0, 1, 2, 3, 1, 0, 3, 2, 0, 2, 1, 3], // ACGTCATGACGT (mostly unique)
        ];

        let (splitters, _singletons, _duplicates) = determine_splitters(&contigs, 3, 100);
        print_splitters(&splitters, 100);
    }
}