pub fn load_feature_bc_matrix(
    outs: impl AsRef<Path>,
    features: &mut Vec<String>,
    barcodes: &mut Vec<String>,
    gex_sparse_matrix: &mut Vec<Vec<(i32, i32)>>
)