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

// Basic AGC Archive Extraction Example
//
// This example demonstrates:
// - Opening an AGC archive
// - Listing all samples
// - Extracting a single sample
// - Accessing contig names and sequences

use anyhow::Result;
use ragc_core::{Decompressor, DecompressorConfig};

fn main() -> Result<()> {
    // Open the AGC archive
    let archive_path = std::env::args()
        .nth(1)
        .unwrap_or_else(|| "data.agc".to_string());

    println!("Opening archive: {archive_path}");

    let mut decompressor = Decompressor::open(&archive_path, DecompressorConfig::default())?;

    // List all samples in the archive
    let samples = decompressor.list_samples();
    println!("\nFound {} samples:", samples.len());
    for (i, sample) in samples.iter().enumerate() {
        println!("  {}: {}", i + 1, sample);
        if i >= 9 {
            println!("  ... ({} more)", samples.len() - 10);
            break;
        }
    }

    // Extract the first sample (if available)
    if let Some(first_sample) = samples.first() {
        println!("\nExtracting sample: {first_sample}");

        let contigs = decompressor.get_sample(first_sample)?;

        println!("  Found {} contigs:", contigs.len());
        for (contig_name, sequence) in &contigs {
            println!("    {}: {} bp", contig_name, sequence.len());
        }

        // Display first 100bp of first contig
        if let Some((contig_name, sequence)) = contigs.first() {
            let preview = String::from_utf8_lossy(&sequence[..sequence.len().min(100)]);
            println!("\nFirst 100bp of {contig_name}:");
            println!("{preview}");
        }
    }

    // Close the decompressor
    decompressor.close()?;

    Ok(())
}