md_analysis 0.1.0

molecular dynamics
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

use clap::{Args, Parser, Subcommand};

/// 分子动力学分析工具 — 读取 BAM 文件中的 dw/ar 字段并分 bin 计数
#[derive(Parser, Debug)]
#[command(version, about, long_about = None)]
pub struct CliArgs {
    #[command(subcommand)]
    pub commands: Commands,
    /// 输入的 BAM 文件路径
    #[arg(short = 'i')]
    pub bam_file: String,

    /// 输出目录（默认：当前目录）
    #[arg(default_value = ".")]
    pub output_dir: String,

    /// stats worker 数量（默认：可用核心数的一半）
    #[arg(long, short = 'w', default_value_t = 0)]
    pub workers: usize,
}

#[derive(Subcommand, Debug, Clone)]
pub enum Commands {
    RawBam(RawBamArgs),
    AlignedBam(AlignedBamArgs),
}

#[derive(Debug, Args, Clone)]
pub struct RawBamArgs {
    /// 仅处理物理长度 >= 此值的 reads（None = 不限制）
    #[arg(long = "length-min")]
    pub length_min: Option<usize>,

    /// 仅处理物理长度 <= 此值的 reads（None = 不限制）
    #[arg(long = "length-max")]
    pub length_max: Option<usize>,

    /// 仅收集前 n 个位置 + 后 n 个位置的 dw/ar 值（0 = 收集全部）
    #[arg(long = "slice-n", default_value_t = 0)]
    pub slice_n: usize,
}

#[derive(Debug, Args, Clone)]
pub struct AlignedBamArgs {
    #[arg(long = "ref-range")]
    pub ref_range: Option<String>,
}