@include "atac_rna/sc_atac_gex_counter_cs.mro"

call SC_ATAC_GEX_COUNTER_CS(
    sample_id                   = "pbmc_granulocyte_sorted_10k",
    sample_desc                 = "pbmc_granulocyte_sorted_10k",
    reference_path              = "/home/Genomes/cellranger/refdata-cellranger-arc-GRCh38-2020-A-2.0.0",
    sample_def                  = [
        {
            "fastq_id": null,
            "fastq_mode": "ILMN_BCL2FASTQ",
            "gem_group": null,
            "lanes": null,
            "library_type": "Gene Expression",
            "read_path": "/net/bmc-lab5/data/kellis/users/benjames/seq/pbmc_granulocyte_sorted_10k/gex",
            "sample_indices": ["any"],
            "sample_names": ["pbmc_granulocyte_sorted_10k"],
            "subsample_rate": null,
            "target_set": null,
            "target_set_name": null,
        },
        {
            "fastq_id": null,
            "fastq_mode": "ILMN_BCL2FASTQ",
            "gem_group": null,
            "lanes": null,
            "library_type": "Chromatin Accessibility",
            "read_path": "/net/bmc-lab5/data/kellis/users/benjames/seq/pbmc_granulocyte_sorted_10k/atac",
            "sample_indices": ["any"],
            "sample_names": ["pbmc_granulocyte_sorted_10k"],
            "subsample_rate": null,
            "target_set": null,
            "target_set_name": null,
        },
    ],
    force_cells                 = null,
    skip_compatibility_check    = false,
    rna_include_introns         = true,
    custom_peaks                = null,
    peak_qval                   = null,
    feature_linkage_max_dist_mb = null,
    k_means_max_clusters        = null,
    no_bam                      = false,
)
