@include "atac/_sc_atac_analyzer_stages.mro"

call PERFORM_DIFFERENTIAL_ANALYSIS as COMPUTE_DIFFERENTIAL_ATAC(
    peaks                   = "/gpfs/data/sequence/results/aifantislab/2021-runs/2021-07-27/cellranger/count-NYU070/SC_ATAC_GEX_COUNTER_CS/SC_ATAC_GEX_COUNTER/_ATAC_MATRIX_COMPUTER/_PEAK_CALLER/DETECT_PEAKS/fork0/join-u45960ddab2/files/peaks.bed",
    reference_path          = "/gpfs/data/sequence/cellranger-refdata/refdata-cellranger-arc-GRCh38-2020-A",
    filtered_peak_bc_matrix = "/gpfs/data/sequence/results/aifantislab/2021-runs/2021-07-27/cellranger/count-NYU070/SC_ATAC_GEX_COUNTER_CS/SC_ATAC_GEX_COUNTER/_JOINT_CELL_DETECTOR/MERGE_ATAC_RNA_MATRICES/fork0/join-u45960deb66/files/atac.h5",
    filtered_tf_bc_matrix   = "/gpfs/data/sequence/results/aifantislab/2021-runs/2021-07-27/cellranger/count-NYU070/SC_ATAC_GEX_COUNTER_CS/SC_ATAC_GEX_COUNTER/_SC_ATAC_GEX_ANALYZER/_PEAK_ANNOTATOR/GENERATE_TF_MATRIX/fork0/join-u45960df02d/files/filtered_tf_bc_matrix.h5",
    factorization           = ["lsa"],
    clustering              = "/gpfs/data/sequence/results/aifantislab/2021-runs/2021-07-27/cellranger/count-NYU070/SC_ATAC_GEX_COUNTER_CS/SC_ATAC_GEX_COUNTER/_SC_ATAC_GEX_ANALYZER/_ATAC_CLUSTERING_COMPUTER/COMBINE_CLUSTERING/fork0/chnk0-u45960df3fa/files/clustering",
)
