bam_tide 1.2.0

A fast and memory-efficient BAM processing toolkit for coverage calculation and quantification, designed as a scalable alternative to deeptools bamCoverage for large sequencing datasets. And additional BAM tools.
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

// src/bin/bam_coverage.rs
use bam_tide::bed_data::BedData;
use bam_tide::cli::CoverageCli;
use clap::{Parser};

fn main() {
    let opts = CoverageCli::parse();

    let bed = match BedData::from_bam_with_policy(&opts) {
        Ok(b) => b,
        Err(e) => {
            eprintln!("Error: {e}");
            std::process::exit(1);
        }
    };

    if let Err(e) = write_output(&bed, &opts.outfile) {
        eprintln!("Error: {e}");
        std::process::exit(1);
    }
}

pub fn write_output(bed: &BedData, outfile: &str) -> Result<(), String> {
    if outfile.ends_with(".bw") || outfile.ends_with(".bigwig") {
        bed.write_bigwig(outfile)
    } else {
        bed.write_bedgraph(outfile).map_err(|e| e.to_string())
    }
}