deepbiop-bam 0.1.11

Deep Learning Processing Library for Bam 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

use std::path::PathBuf;

use ahash::HashMap;
use anyhow::Result;
use noodles::sam::record::Cigar;
use pyo3::prelude::*;

use crate::chimeric;
use crate::cigar::calc_softclips;

use pyo3_stub_gen::derive::*;

/// Calculate the number of chimeric reads in a BAM file.
#[gen_stub_pyfunction(module = "deepbiop.bam")]
#[pyfunction]
fn count_chimeric_reads_for_path(bam: PathBuf, threads: Option<usize>) -> Result<usize> {
    chimeric::count_chimeric_reads_for_path(bam, threads)
}

/// Calculate the number of chimeric reads in multiple BAM files.
#[gen_stub_pyfunction(module = "deepbiop.bam")]
#[pyfunction]
fn count_chimeric_reads_for_paths(
    bams: Vec<PathBuf>,
    threads: Option<usize>,
) -> Result<HashMap<PathBuf, usize>> {
    Ok(chimeric::count_chimeric_reads_for_paths(&bams, threads))
}

/// Calculate left and right soft clips from a cigar string.
#[gen_stub_pyfunction(module = "deepbiop.bam")]
#[pyfunction]
fn left_right_soft_clip(cigar_string: &str) -> Result<(usize, usize)> {
    let cigar = Cigar::new(cigar_string.as_bytes());
    calc_softclips(&cigar)
}

// register bam sub module
pub fn register_bam_module(parent_module: &Bound<'_, PyModule>) -> PyResult<()> {
    let sub_module_name = "bam";
    let child_module = PyModule::new_bound(parent_module.py(), sub_module_name)?;

    child_module.add_function(wrap_pyfunction!(left_right_soft_clip, &child_module)?)?;
    child_module.add_function(wrap_pyfunction!(
        count_chimeric_reads_for_path,
        &child_module
    )?)?;
    child_module.add_function(wrap_pyfunction!(
        count_chimeric_reads_for_paths,
        &child_module
    )?)?;

    parent_module.add_submodule(&child_module)?;
    Ok(())
}