Add selected results from doHMRF to the giotto object

addHMRF(
  gobject,
  spat_unit = NULL,
  feat_type = NULL,
  HMRFoutput,
  k = NULL,
  betas_to_add = NULL,
  hmrf_name = NULL
)

Arguments

gobject

giotto object

spat_unit

spatial unit

feat_type

feature type

HMRFoutput

HMRF output from doHMRF()

k

number of domains

betas_to_add

results from different betas that you want to add

hmrf_name

specify a custom name

Value

giotto object

Examples

g <- GiottoData::loadGiottoMini("visium")
#> 1. read Giotto object
#> 2. read Giotto feature information
#> 3. read Giotto spatial information
#> 3.1 read Giotto spatial shape information
#> 3.2 read Giotto spatial centroid information
#> 3.3 read Giotto spatial overlap information
#> 4. read Giotto image information
#> python already initialized in this session
#>  active environment : '/usr/bin/python3'
#>  python version : 3.10
#> checking default envname 'giotto_env'
#> a system default python environment was found
#> Using python path:
#>  "/usr/bin/python3"
spat_genes <- binSpect(g)
#> 
#> This is the single parameter version of binSpect
#> 
#> 1. matrix binarization complete
#> 
#> 2. spatial enrichment test completed
#> 
#> 3. (optional) average expression of high
#>  expressing cells calculated
#> 
#> 4. (optional) number of high expressing cells
#>  calculated

output_folder <- file.path(tempdir(), "HMRF")
if (!file.exists(output_folder)) dir.create(output_folder, recursive = TRUE)

out <- doHMRF(
    g,
    spatial_genes = spat_genes[seq_len(20)]$feats,
    expression_values = "scaled",
    spatial_network_name = "Delaunay_network",
    k = 6, betas = c(0, 10, 5),
    output_folder = output_folder
)
#> Error: package 'smfishHmrf' is not yet installed
#> 
#>  To install:
#> ## active python env: '/usr/bin/python3' 
#> ## python version: 3.10
#> ## restart session then use GiottoClass::set_giotto_python_path() if this is incorrect
#> reticulate::conda_install(envname = '/usr/bin/python3', packages = c('smfishHmrf'), pip = TRUE)

g <- addHMRF(
    gobject = g,
    HMRFoutput = out,
    k = 6,
    betas_to_add = 20,
    hmrf_name = "HMRF"
)
#> Error: object 'out' not found

spatPlot(
    gobject = g, cell_color = "HMRF_k6_b.20",
)
#> Error in (function (ggobject, instrs = NULL, sdimx = NULL, sdimy = NULL,     plot_method = c("ggplot", "scattermore"), cell_locations_metadata_selected,     cell_locations_metadata_other, cell_color = NULL, color_as_factor = TRUE,     cell_color_code = NULL, cell_color_gradient = NULL, gradient_midpoint = NULL,     gradient_style = "divergent", gradient_limits = NULL, select_cell_groups = NULL,     select_cells = NULL, point_size = 2, point_alpha = 1, point_border_col = "lightgrey",     point_border_stroke = 0.1, show_cluster_center = FALSE, show_center_label = TRUE,     center_point_size = 4, center_point_border_col = "black",     center_point_border_stroke = 0.1, label_size = 4, label_fontface = "bold",     show_other_cells = TRUE, other_cell_color = "lightgrey",     other_point_size = 1, show_legend = TRUE, ...) {    if (is.null(sdimx) || is.null(sdimy)) {        warning(wrap_txt("plot_method = ggplot,\n                        but spatial dimensions for sdimx and/or sdimy\n                        are not specified. It will default to the\n                        'sdimx' and 'sdimy' "))        sdimx <- "sdimx"        sdimy <- "sdimy"    }    plot_method <- match.arg(plot_method, c("ggplot", "scattermore"))    point_size <- as.numeric(point_size)    point_alpha <- as.numeric(point_alpha)    point_border_stroke <- as.numeric(point_border_stroke)    center_point_size <- as.numeric(center_point_size)    center_point_border_stroke <- as.numeric(center_point_border_stroke)    label_size <- as.numeric(label_size)    other_point_size <- as.numeric(other_point_size)    pl <- ggobject    if ((!is.null(select_cells) || !is.null(select_cell_groups)) &&         isTRUE(show_other_cells)) {        pl <- pl + ggplot2::geom_point(data = cell_locations_metadata_other,             mapping = aes_string(x = sdimx, sdimy), color = other_cell_color,             show.legend = FALSE, size = other_point_size, alpha = point_alpha)    }    if (is.null(cell_color)) {        cell_color <- "lightblue"        pl <- pl + giotto_point(plot_method = plot_method, data = cell_locations_metadata_selected,             mapping = aes_string(x = sdimx, y = sdimy), show.legend = show_legend,             shape = 21, fill = cell_color, size = point_size,             stroke = point_border_stroke, color = point_border_col,             alpha = point_alpha)    }    else if (length(cell_color) > 1L) {        if (is.numeric(cell_color) || is.factor(cell_color)) {            if (nrow(cell_locations_metadata_selected) != length(cell_color)) {                stop("\n vector needs to be the same lengths as\n                    number of cells \n")            }            cell_locations_metadata_selected[["temp_color"]] <- cell_color            pl <- pl + giotto_point(plot_method = plot_method,                 data = cell_locations_metadata_selected, mapping = aes_string2(x = sdimx,                   y = sdimy, fill = "temp_color"), show.legend = show_legend,                 shape = 21, size = point_size, color = point_border_col,                 stroke = point_border_stroke, alpha = point_alpha)        }        else if (is.character(cell_color)) {            if (!all(cell_color %in% grDevices::colors())) {                stop("cell_color is not numeric,\n                    a factor or vector of colors \n")            }            pl <- pl + giotto_point(plot_method = plot_method,                 data = cell_locations_metadata_selected, mapping = aes_string2(x = sdimx,                   y = sdimy), show.legend = show_legend, shape = 21,                 fill = cell_color, size = point_size, color = point_border_col,                 stroke = point_border_stroke, alpha = point_alpha)        }    }    else if (is.character(cell_color)) {        if (!cell_color %in% colnames(cell_locations_metadata_selected)) {            if (!cell_color %in% grDevices::colors()) {                stop(cell_color, " is not a color or a column name \n")            }            pl <- pl + giotto_point(plot_method = plot_method,                 data = cell_locations_metadata_selected, mapping = aes_string2(x = sdimx,                   y = sdimy), show.legend = show_legend, shape = 21,                 fill = cell_color, size = point_size, color = point_border_col,                 stroke = point_border_stroke, alpha = point_alpha)        }        else {            class_cell_color <- class(cell_locations_metadata_selected[[cell_color]])            if ((class_cell_color %in% c("integer", "numeric")) &&                 !isTRUE(color_as_factor)) {                if (!is.null(gradient_limits) && is.vector(gradient_limits) &&                   length(gradient_limits) == 2) {                  lower_lim <- gradient_limits[[1]]                  upper_lim <- gradient_limits[[2]]                  numeric_data <- cell_locations_metadata_selected[[cell_color]]                  limit_numeric_data <- ifelse(numeric_data >                     upper_lim, upper_lim, ifelse(numeric_data <                     lower_lim, lower_lim, numeric_data))                  cell_locations_metadata_selected[[cell_color]] <- limit_numeric_data                }                pl <- pl + giotto_point(plot_method = plot_method,                   data = cell_locations_metadata_selected, mapping = aes_string2(x = sdimx,                     y = sdimy, fill = cell_color), show.legend = show_legend,                   shape = 21, size = point_size, color = point_border_col,                   stroke = point_border_stroke, alpha = point_alpha)            }            else {                if (isTRUE(color_as_factor)) {                  factor_data <- factor(cell_locations_metadata_selected[[cell_color]])                  cell_locations_metadata_selected[[cell_color]] <- factor_data                }                if (isTRUE(show_cluster_center) || isTRUE(show_center_label)) {                  annotated_DT_centers <- cell_locations_metadata_selected[,                     .(center_1 = stats::median(get("sdimx")),                       center_2 = stats::median(get("sdimy"))),                     by = cell_color]                  factor_center_data <- factor(annotated_DT_centers[[cell_color]])                  annotated_DT_centers[[cell_color]] <- factor_center_data                }                pl <- pl + giotto_point(plot_method = plot_method,                   data = cell_locations_metadata_selected, mapping = aes_string2(x = sdimx,                     y = sdimy, fill = cell_color), show.legend = show_legend,                   shape = 21, size = point_size, color = point_border_col,                   stroke = point_border_stroke, alpha = point_alpha)                if (isTRUE(show_cluster_center) && (isTRUE(color_as_factor) ||                   class_cell_color %in% c("character", "factor"))) {                  pl <- pl + giotto_point(plot_method = plot_method,                     data = annotated_DT_centers, mapping = aes_string2(x = "center_1",                       y = "center_2", fill = cell_color), color = center_point_border_col,                     stroke = center_point_border_stroke, size = center_point_size,                     shape = 21, alpha = point_alpha)                }                if (isTRUE(show_center_label)) {                  pl <- pl + ggrepel::geom_text_repel(data = annotated_DT_centers,                     mapping = aes_string2(x = "center_1", y = "center_2",                       label = cell_color), size = label_size,                     fontface = label_fontface)                }            }            if (!is.null(cell_color_code)) {                pl <- pl + ggplot2::scale_fill_manual(values = cell_color_code)            }            else if (isTRUE(color_as_factor)) {                number_colors <- length(unique(factor_data))                cell_color_code <- set_default_color_discrete_cell(instrs = instrs)(n = number_colors)                names(cell_color_code) <- unique(factor_data)                pl <- pl + ggplot2::scale_fill_manual(values = cell_color_code)            }            else if (!isTRUE(color_as_factor)) {                if (is.null(gradient_midpoint)) {                  gradient_midpoint <- stats::median(cell_locations_metadata_selected[[cell_color]])                }                pl <- pl + set_default_color_continuous_cell(colors = cell_color_gradient,                   instrs = instrs, style = gradient_style, midpoint = gradient_midpoint,                   type = "fill", )            }        }    }    return(pl)})(ggobject = structure(list(data = structure(list(), class = "waiver"),     layers = list(), scales = <environment>, guides = <environment>,     mapping = structure(list(), names = character(0), class = "uneval"),     theme = structure(list(line = structure(list(colour = "black",         linewidth = 0.5, linetype = 1, lineend = "butt", arrow = FALSE,         inherit.blank = TRUE), class = c("element_line", "element"    )), rect = structure(list(fill = "white", colour = "black",         linewidth = 0.5, linetype = 1, inherit.blank = TRUE), class = c("element_rect",     "element")), text = structure(list(family = "", face = "plain",         colour = "black", size = 11, hjust = 0.5, vjust = 0.5,         angle = 0, lineheight = 0.9, margin = structure(c(0,         0, 0, 0), unit = 8L, class = c("margin", "simpleUnit",         "unit", "unit_v2")), debug = FALSE, inherit.blank = TRUE), class = c("element_text",     "element")), title = NULL, aspect.ratio = NULL, axis.title = NULL,         axis.title.x = structure(list(family = NULL, face = NULL,             colour = NULL, size = NULL, hjust = NULL, vjust = 1,             angle = NULL, lineheight = NULL, margin = structure(c(2.75,             0, 0, 0), unit = 8L, class = c("margin", "simpleUnit",             "unit", "unit_v2")), debug = NULL, inherit.blank = TRUE), class = c("element_text",         "element")), axis.title.x.top = structure(list(family = NULL,             face = NULL, colour = NULL, size = NULL, hjust = NULL,             vjust = 0, angle = NULL, lineheight = NULL, margin = structure(c(0,             0, 2.75, 0), unit = 8L, class = c("margin", "simpleUnit",             "unit", "unit_v2")), debug = NULL, inherit.blank = TRUE), class = c("element_text",         "element")), axis.title.x.bottom = NULL, axis.title.y = structure(list(            family = NULL, face = NULL, colour = NULL, size = NULL,             hjust = NULL, vjust = 1, angle = 90, lineheight = NULL,             margin = structure(c(0, 2.75, 0, 0), unit = 8L, class = c("margin",             "simpleUnit", "unit", "unit_v2")), debug = NULL,             inherit.blank = TRUE), class = c("element_text",         "element")), axis.title.y.left = NULL, axis.title.y.right = structure(list(            family = NULL, face = NULL, colour = NULL, size = NULL,             hjust = NULL, vjust = 1, angle = -90, lineheight = NULL,             margin = structure(c(0, 0, 0, 2.75), unit = 8L, class = c("margin",             "simpleUnit", "unit", "unit_v2")), debug = NULL,             inherit.blank = TRUE), class = c("element_text",         "element")), axis.text = structure(list(family = NULL,             face = NULL, colour = "grey30", size = structure(0.8, class = "rel"),             hjust = NULL, vjust = NULL, angle = NULL, lineheight = NULL,             margin = NULL, debug = NULL, inherit.blank = TRUE), class = c("element_text",         "element")), axis.text.x = structure(list(family = NULL,             face = NULL, colour = NULL, size = NULL, hjust = NULL,             vjust = 1, angle = NULL, lineheight = NULL, margin = structure(c(2.2,             0, 0, 0), unit = 8L, class = c("margin", "simpleUnit",             "unit", "unit_v2")), debug = NULL, inherit.blank = TRUE), class = c("element_text",         "element")), axis.text.x.top = structure(list(family = NULL,             face = NULL, colour = NULL, size = NULL, hjust = NULL,             vjust = 0, angle = NULL, lineheight = NULL, margin = structure(c(0,             0, 2.2, 0), unit = 8L, class = c("margin", "simpleUnit",             "unit", "unit_v2")), debug = NULL, inherit.blank = TRUE), class = c("element_text",         "element")), axis.text.x.bottom = NULL, axis.text.y = structure(list(            family = NULL, face = NULL, colour = NULL, size = NULL,             hjust = 1, vjust = NULL, angle = NULL, lineheight = NULL,             margin = structure(c(0, 2.2, 0, 0), unit = 8L, class = c("margin",             "simpleUnit", "unit", "unit_v2")), debug = NULL,             inherit.blank = TRUE), class = c("element_text",         "element")), axis.text.y.left = NULL, axis.text.y.right = structure(list(            family = NULL, face = NULL, colour = NULL, size = NULL,             hjust = 0, vjust = NULL, angle = NULL, lineheight = NULL,             margin = structure(c(0, 0, 0, 2.2), unit = 8L, class = c("margin",             "simpleUnit", "unit", "unit_v2")), debug = NULL,             inherit.blank = TRUE), class = c("element_text",         "element")), axis.text.theta = NULL, axis.text.r = structure(list(            family = NULL, face = NULL, colour = NULL, size = NULL,             hjust = 0.5, vjust = NULL, angle = NULL, lineheight = NULL,             margin = structure(c(0, 2.2, 0, 2.2), unit = 8L, class = c("margin",             "simpleUnit", "unit", "unit_v2")), debug = NULL,             inherit.blank = TRUE), class = c("element_text",         "element")), axis.ticks = structure(list(colour = "grey20",             linewidth = NULL, linetype = NULL, lineend = NULL,             arrow = FALSE, inherit.blank = TRUE), class = c("element_line",         "element")), axis.ticks.x = NULL, axis.ticks.x.top = NULL,         axis.ticks.x.bottom = NULL, axis.ticks.y = NULL, axis.ticks.y.left = NULL,         axis.ticks.y.right = NULL, axis.ticks.theta = NULL, axis.ticks.r = NULL,         axis.minor.ticks.x.top = NULL, axis.minor.ticks.x.bottom = NULL,         axis.minor.ticks.y.left = NULL, axis.minor.ticks.y.right = NULL,         axis.minor.ticks.theta = NULL, axis.minor.ticks.r = NULL,         axis.ticks.length = structure(2.75, unit = 8L, class = c("simpleUnit",         "unit", "unit_v2")), axis.ticks.length.x = NULL, axis.ticks.length.x.top = NULL,         axis.ticks.length.x.bottom = NULL, axis.ticks.length.y = NULL,         axis.ticks.length.y.left = NULL, axis.ticks.length.y.right = NULL,         axis.ticks.length.theta = NULL, axis.ticks.length.r = NULL,         axis.minor.ticks.length = structure(0.75, class = "rel"),         axis.minor.ticks.length.x = NULL, axis.minor.ticks.length.x.top = NULL,         axis.minor.ticks.length.x.bottom = NULL, axis.minor.ticks.length.y = NULL,         axis.minor.ticks.length.y.left = NULL, axis.minor.ticks.length.y.right = NULL,         axis.minor.ticks.length.theta = NULL, axis.minor.ticks.length.r = NULL,         axis.line = structure(list(), class = c("element_blank",         "element")), axis.line.x = NULL, axis.line.x.top = NULL,         axis.line.x.bottom = NULL, axis.line.y = NULL, axis.line.y.left = NULL,         axis.line.y.right = NULL, axis.line.theta = NULL, axis.line.r = NULL,         legend.background = structure(list(fill = NULL, colour = NA,             linewidth = NULL, linetype = NULL, inherit.blank = TRUE), class = c("element_rect",         "element")), legend.margin = structure(c(5.5, 5.5, 5.5,         5.5), unit = 8L, class = c("margin", "simpleUnit", "unit",         "unit_v2")), legend.spacing = structure(11, unit = 8L, class = c("simpleUnit",         "unit", "unit_v2")), legend.spacing.x = NULL, legend.spacing.y = NULL,         legend.key = NULL, legend.key.size = structure(1.2, unit = 3L, class = c("simpleUnit",         "unit", "unit_v2")), legend.key.height = NULL, legend.key.width = NULL,         legend.key.spacing = structure(5.5, unit = 8L, class = c("simpleUnit",         "unit", "unit_v2")), legend.key.spacing.x = NULL, legend.key.spacing.y = NULL,         legend.frame = NULL, legend.ticks = NULL, legend.ticks.length = structure(0.2, class = "rel"),         legend.axis.line = NULL, legend.text = structure(list(            family = NULL, face = NULL, colour = NULL, size = structure(0.8, class = "rel"),             hjust = NULL, vjust = NULL, angle = NULL, lineheight = NULL,             margin = NULL, debug = NULL, inherit.blank = TRUE), class = c("element_text",         "element")), legend.text.position = NULL, legend.title = structure(list(            family = NULL, face = NULL, colour = NULL, size = NULL,             hjust = 0, vjust = NULL, angle = NULL, lineheight = NULL,             margin = NULL, debug = NULL, inherit.blank = TRUE), class = c("element_text",         "element")), legend.title.position = NULL, legend.position = "right",         legend.position.inside = NULL, legend.direction = NULL,         legend.byrow = NULL, legend.justification = "center",         legend.justification.top = NULL, legend.justification.bottom = NULL,         legend.justification.left = NULL, legend.justification.right = NULL,         legend.justification.inside = NULL, legend.location = NULL,         legend.box = NULL, legend.box.just = NULL, legend.box.margin = structure(c(0,         0, 0, 0), unit = 1L, class = c("margin", "simpleUnit",         "unit", "unit_v2")), legend.box.background = structure(list(), class = c("element_blank",         "element")), legend.box.spacing = structure(11, unit = 8L, class = c("simpleUnit",         "unit", "unit_v2")), panel.background = structure(list(            fill = "white", colour = NA, linewidth = NULL, linetype = NULL,             inherit.blank = TRUE), class = c("element_rect",         "element")), panel.border = structure(list(fill = NA,             colour = "grey20", linewidth = NULL, linetype = NULL,             inherit.blank = TRUE), class = c("element_rect",         "element")), panel.spacing = structure(5.5, unit = 8L, class = c("simpleUnit",         "unit", "unit_v2")), panel.spacing.x = NULL, panel.spacing.y = NULL,         panel.grid = structure(list(colour = "grey92", linewidth = NULL,             linetype = NULL, lineend = NULL, arrow = FALSE, inherit.blank = TRUE), class = c("element_line",         "element")), panel.grid.major = NULL, panel.grid.minor = structure(list(            colour = NULL, linewidth = structure(0.5, class = "rel"),             linetype = NULL, lineend = NULL, arrow = FALSE, inherit.blank = TRUE), class = c("element_line",         "element")), panel.grid.major.x = NULL, panel.grid.major.y = NULL,         panel.grid.minor.x = NULL, panel.grid.minor.y = NULL,         panel.ontop = FALSE, plot.background = structure(list(            fill = NULL, colour = "white", linewidth = NULL,             linetype = NULL, inherit.blank = TRUE), class = c("element_rect",         "element")), plot.title = structure(list(family = NULL,             face = NULL, colour = NULL, size = structure(1.2, class = "rel"),             hjust = 0, vjust = 1, angle = NULL, lineheight = NULL,             margin = structure(c(0, 0, 5.5, 0), unit = 8L, class = c("margin",             "simpleUnit", "unit", "unit_v2")), debug = NULL,             inherit.blank = TRUE), class = c("element_text",         "element")), plot.title.position = "panel", plot.subtitle = structure(list(            family = NULL, face = NULL, colour = NULL, size = NULL,             hjust = 0, vjust = 1, angle = NULL, lineheight = NULL,             margin = structure(c(0, 0, 5.5, 0), unit = 8L, class = c("margin",             "simpleUnit", "unit", "unit_v2")), debug = NULL,             inherit.blank = TRUE), class = c("element_text",         "element")), plot.caption = structure(list(family = NULL,             face = NULL, colour = NULL, size = structure(0.8, class = "rel"),             hjust = 1, vjust = 1, angle = NULL, lineheight = NULL,             margin = structure(c(5.5, 0, 0, 0), unit = 8L, class = c("margin",             "simpleUnit", "unit", "unit_v2")), debug = NULL,             inherit.blank = TRUE), class = c("element_text",         "element")), plot.caption.position = "panel", plot.tag = structure(list(            family = NULL, face = NULL, colour = NULL, size = structure(1.2, class = "rel"),             hjust = 0.5, vjust = 0.5, angle = NULL, lineheight = NULL,             margin = NULL, debug = NULL, inherit.blank = TRUE), class = c("element_text",         "element")), plot.tag.position = "topleft", plot.tag.location = NULL,         plot.margin = structure(c(5.5, 5.5, 5.5, 5.5), unit = 8L, class = c("margin",         "simpleUnit", "unit", "unit_v2")), strip.background = structure(list(            fill = "grey85", colour = "grey20", linewidth = NULL,             linetype = NULL, inherit.blank = TRUE), class = c("element_rect",         "element")), strip.background.x = NULL, strip.background.y = NULL,         strip.clip = "inherit", strip.placement = "inside", strip.text = structure(list(            family = NULL, face = NULL, colour = "grey10", size = structure(0.8, class = "rel"),             hjust = NULL, vjust = NULL, angle = NULL, lineheight = NULL,             margin = structure(c(4.4, 4.4, 4.4, 4.4), unit = 8L, class = c("margin",             "simpleUnit", "unit", "unit_v2")), debug = NULL,             inherit.blank = TRUE), class = c("element_text",         "element")), strip.text.x = NULL, strip.text.x.bottom = NULL,         strip.text.x.top = NULL, strip.text.y = structure(list(            family = NULL, face = NULL, colour = NULL, size = NULL,             hjust = NULL, vjust = NULL, angle = -90, lineheight = NULL,             margin = NULL, debug = NULL, inherit.blank = TRUE), class = c("element_text",         "element")), strip.text.y.left = structure(list(family = NULL,             face = NULL, colour = NULL, size = NULL, hjust = NULL,             vjust = NULL, angle = 90, lineheight = NULL, margin = NULL,             debug = NULL, inherit.blank = TRUE), class = c("element_text",         "element")), strip.text.y.right = NULL, strip.switch.pad.grid = structure(2.75, unit = 8L, class = c("simpleUnit",         "unit", "unit_v2")), strip.switch.pad.wrap = structure(2.75, unit = 8L, class = c("simpleUnit",         "unit", "unit_v2"))), class = c("theme", "gg"), complete = TRUE, validate = TRUE),     coordinates = <environment>, facet = <environment>, plot_env = <environment>,     layout = <environment>, labels = structure(list(), names = character(0))), class = c("gg", "ggplot")), ext = new("SpatExtent", ptr = new("Rcpp_SpatExtent",     .xData = <environment>)), instrs = structure(list(python_path = "/usr/bin/python3",     show_plot = TRUE, return_plot = FALSE, save_plot = FALSE,     save_dir = NA, plot_format = "png", dpi = 300, units = "in",     height = 9, width = 9, is_docker = FALSE, active_spat_unit = "cell",     active_feat_type = "rna"), class = "giottoInstructions"),     sdimx = "sdimx", sdimy = "sdimy", plot_method = "ggplot",     cell_locations_metadata_selected = structure(list(cell_ID = c("AAAGGGATGTAGCAAG-1",     "AAATGGCATGTCTTGT-1", "AAATGGTCAATGTGCC-1", "AAATTAACGGGTAGCT-1",     "AACAACTGGTAGTTGC-1", "AACAGGAAATCGAATA-1", "AACAGGATGGGCCGCG-1",     "AACCATGGGATCGCTA-1", "AACCCAGAGACGGAGA-1", "AACCGAGCTTGGTCAT-1",     "AACCGTTGTGTTTGCT-1", "AACGATAGAAGGGCCG-1", "AACGATATGTCAACTG-1",     "AACGCGGTCTCCAGCC-1", "AACGTCAGACTAGTGG-1", "AACTCGATGGCGCAGT-1",     "AACTGATATTAGGCCT-1", "AACTGGGTCCCGACGT-1", "AACTTGCCCGTATGCA-1",     "AAGAGATGAATCGGTA-1", "AAGCTCGTGCCAAGTC-1", "AAGGAGCGGTTGGTGC-1",     "AAGGCGCGTAAAGCTT-1", "AAGGCTGTGCTCATCG-1", "AAGTAGAAGACCGGGT-1",     "AAGTAGTGACGCGAGG-1", "AATCCAAGGGCCTGAG-1", "AATGACTGTCAGCCGG-1",     "AATGCAACCGGGTACC-1", "AATGGTTCTCACAAGC-1", "AATTTGGTTCCAAAGA-1",     "ACAACAGCATGAGCTA-1", "ACAACGGTCCCTGCGA-1", "ACAAGGACAAGAGGTT-1",     "ACACAAAGACGGGTGG-1", "ACACACTTTCTACACG-1", "ACACCCAGCATGCAGC-1",     "ACACCCGAGAAATCCG-1", "ACACCTTAAGTAGGGC-1", "ACACGAGACTCCTTCT-1",     "ACACGGGAACTTAGGG-1", "ACATAATAAGGCGGTG-1", "ACCAAGAACGCGTGTC-1",     "ACCACAAGTTTCTATC-1", "ACCCTTCATCTGCGAA-1", "ACCGACTGAGTCCCAC-1",     "ACCTAATCGACTTCCT-1", "ACGATCATACATAGAG-1", "ACGCATTCGTGAGTAC-1",     "ACGCCAGATGATTTCT-1", "ACGCGGGCCAAGGACA-1", "ACGGAATTTAGCAAAT-1",     "ACGTAGATTGCTGATG-1", "ACGTGCGCCTCGTGCA-1", "ACGTTTAGTTGTGATC-1",     "ACTACCAGCTCTCTGG-1", "ACTATTCGTCCGTGGT-1", "ACTCAATAAAGGCACG-1",     "ACTGCCGTCGTAACTC-1", "ACTGTAGCACTTTGGA-1", "ACTGTCCAGGATTATA-1",     "ACTGTCTTCTTTAGAA-1", "ACTGTGCTAGTAGATC-1", "ACTTAGTACGACAAGA-1",     "ACTTCGCCATACGCAC-1", "ACTTGGGACCCGGTGG-1", "ACTTTCCTATAGCTTC-1",     "AGAACGTGGTACATTC-1", "AGACCATGGGATACAA-1", "AGACCCGCCCTCCTCG-1",     "AGACTAGCCTTCCAGA-1", "AGAGCGGGCTAATCAT-1", "AGAGCGTACAAGCTCG-1",     "AGATAACTTCAGGGCC-1", "AGATACGACTTCATAT-1", "AGATCTCAGGTGTGAT-1",     "AGATGACTCGCCCACG-1", "AGCACCAGTACTCACG-1", "AGCACTTAAGGACGCC-1",     "AGCATCGTCGATAATT-1", "AGCCCTAAGCGAAGTT-1", "AGCGACATCCCATTCA-1",     "AGCGACCAACGATATT-1", "AGCGGGAAGGGTCCAT-1", "AGCTATTTAATCCAAC-1",     "AGGACATCGCACGTCG-1", "AGGCCACCCGTTATGA-1", "AGGCCCATTGTACAGG-1",     "AGGCTATGGTTAGCTT-1", "AGGGCGTGATCGGCTA-1", "AGGGCTGCAGTTACAG-1",     "AGGGTTTAGTTCGGGA-1", "AGGTAACCTCCTATTC-1", "AGGTACGATATTGCCA-1",     "AGGTAGGTACAAAGCT-1", "AGGTTGAGGCACGCTT-1", "AGTAATGTCTTGCCGC-1",     "AGTAGGTAACATACAT-1", "AGTATAATACTAGGCA-1", "AGTCAGCCACCGCCTG-1",     "AGTCGGCCCAAACGAC-1", "AGTCTTCTCCTCAAAT-1", "AGTTCCTACAGAATTA-1",     "AGTTTCGCAGGTCGGA-1", "ATAAAGGCTCGGTCGT-1", "ATAACGGAGTCCAACG-1",     "ATAATAGTGTAGGGAC-1", "ATACCACGGGCAACTT-1", "ATACTAGCATGACCCT-1",     "ATAGACGAAGAGAAAG-1", "ATAGAGTACTGGGACA-1", "ATAGGGATATCCTTGA-1",     "ATAGGTTGGGCAGATG-1", "ATATCTTAGGGCCTTC-1", "ATATTCCCACAGGTCA-1",     "ATATTGCTGTCAAAGT-1", "ATATTTAACCCTCAAG-1", "ATCAATCTGGGCTGCA-1",     "ATCACTTCATCCTCGC-1", "ATCCAACGCAGTCATA-1", "ATCGACTCTTTCCGTT-1",     "ATCGCACGCCGGGAGA-1", "ATCTCCCTGCAATCTA-1", "ATCTGGTTAAGACTGT-1",     "ATCTGTAATTGTACCC-1", "ATCTTGACCTGCAACG-1", "ATGAAGCCAAGGAGCC-1",     "ATGACGCCGGCTCTAA-1", "ATGAGGAGTGTTAATC-1", "ATGATCGGGAATAGAC-1",     "ATGCACTACCGCATTG-1", "ATGCATGATCCAGGAT-1", "ATGCCGGTCTTGCATA-1",     "ATGCGACAGTCCCATT-1", "ATGCTCAGTGTTGCAT-1", "ATGGATTGACCAAACG-1",     "ATGTGCATCCGACGCA-1", "ATGTTACGAGCAATAC-1", "ATTAATACTACGCGGG-1",     "ATTAATTCGGTCACTC-1", "ATTACTAGCCTCTTGC-1", "ATTAGATTGATAGCGG-1",     "ATTATGCCATAGGGAG-1", "ATTCCTCCGCCAGTGC-1", "ATTCTCGTCTCTTTAG-1",     "ATTGAAGATCTTAGTG-1", "ATTGCTGCTCCTCCAT-1", "ATTGGGAATATCTTGG-1",     "ATTTAACTCGTATTAC-1", "ATTTACAGTTTACTGG-1", "ATTTGTTCCAGGGCTC-1",     "CAAACGGTCGCACTTT-1", "CAAACTATTGAGCTTC-1", "CAAATTGTCAGCAAGC-1",     "CAACGACCCGTTTACA-1", "CAACGGTTCTTGATAC-1", "CAACTCCTTGATCCCG-1",     "CAAGCAACGTCGGAGT-1", "CAAGCACCAAATGCCT-1", "CAATAAACCTTGGCCC-1",     "CAATTAAGGGTGATGA-1", "CACACAGGGATAGATT-1", "CACAGCACCCACGGCA-1",     "CACAGCTAGGGAGTGA-1", "CACCATCGGAGGAGAC-1", "CACCCACACGTCACCC-1",     "CACCGCCAGAAGGTTT-1", "CACCTAATCAGTTTAC-1", "CACCTCGATGGTGGAC-1",     "CACGAAAGTTAGTCCC-1", "CACGCACAGCGCAGCT-1", "CACGCAGCGAGGCTTT-1",     "CACTAAAGTTGCCTAT-1", "CACTCAGCTCTTGAGG-1", "CACTGTCCAAGTGAGA-1",     "CACTTAATCAGACGGA-1", "CACTTCGCCACAGGCT-1", "CAGAGCATGAGCTTGC-1",     "CAGATACTAACATAGT-1", "CAGCTCACTGAGACAT-1", "CAGCTCGTGCTTGTGT-1",     "CAGCTTAGTAGGTAGC-1", "CAGTGTCCGCAGAATG-1", "CAGTTCAAATTGACAC-1",     "CATACAAAGCCGAACC-1", "CATAGTAGCATAGTAG-1", "CATCATCTACCCGGAC-1",     "CATCGCCCGCGGCCAA-1", "CATGAACCGACATTTG-1", "CATGATGCACAATTCT-1",     "CATGCGTTAGACAGAA-1", "CATGCTGGCTCCAATT-1", "CATGGCAGGAAGATCG-1",     "CATGGTCTAGATACCG-1", "CATTACGTCGGCCCGT-1", "CATTATGCTTGTTGTG-1",     "CCAAACAGAACCCTCG-1", "CCAACGATGCACTGAT-1", "CCAAGACTTCTGCGAA-1",     "CCAAGCGTAACTCGTA-1", "CCAAGGAACAGAGAGG-1", "CCAATAGTGCCGTCGA-1",     "CCAATTGAATGTTAAT-1", "CCACGAGAAGAGAATC-1", "CCAGGGACGTGGCCTC-1",     "CCATAAACAACCCGAC-1", "CCATAGGTTGGCGTGG-1", "CCATCCATACCAAGTC-1",     "CCATGCCTGTTTAGTA-1", "CCATGCTCTGCAGGAA-1", "CCATTAGCGATAATCC-1",     "CCATTCCCTGCCCACA-1", "CCCGACCATAGTCCGC-1", "CCCGCCATGCTCCCGT-1",     "CCCGTCAGCGTCTGAC-1", "CCCTCATTCTGGAATT-1", "CCCTGAAATGAGTTGA-1",     "CCGAAAGTGGTGAGCA-1", "CCGCCGGTCAACACAC-1", "CCGGAATGGTTTCAGT-1",     "CCGGCGCATATTGGAT-1", "CCGGGACCCGCAGAGA-1", "CCGGGCGGTCTCGTCA-1",     "CCGTACCCAAGCGCCA-1", "CCGTGCCCATGACGGC-1", "CCTACGGCTCAGTCGA-1",     "CCTACTGCTTACACTT-1", "CCTAGTTAGTCGCATG-1", "CCTATGGGTTACCGTC-1",     "CCTCACCAATCTTGAC-1", "CCTCCGACAATTCAAG-1", "CCTCGAAGTGGACGGG-1",     "CCTCGCCAGCAAATTA-1", "CCTCTAATCTGCCAAG-1", "CCTCTCTCCCATCTAG-1",     "CCTCTGGCCTAGACGG-1", "CCTGAATATTTACATA-1", "CCTGCTATTTGAGAAG-1",     "CCTGGCTAGACCCGCC-1", "CCTGTCACCCGGGCTC-1", "CCTTCTTGATCCAGTG-1",     "CGAACGCCCAGTGCCG-1", "CGAAGTTGCTCTGTGT-1", "CGAGAGATGTGAACCT-1",     "CGAGCGTTGATCAGCC-1", "CGAGCTGGGCTTTAGG-1", "CGAGGCTAAATATGGC-1",     "CGCAAACACGAGTTAC-1", "CGCAATTACTTTCGGT-1", "CGCAATTAGGGTAATA-1",     "CGCAATTCTACAATAA-1", "CGCATTAGCTAATAGG-1", "CGCGCCCGACTTAATA-1",     "CGCGTTCATGAAATAC-1", "CGCTATACCGCCCACT-1", "CGCTATTCAATGTATG-1",     "CGCTCGACATAATGAT-1", "CGCTCTCCGTAGATTA-1", "CGCTGTGACGCCGCAC-1",     "CGCTTCCACTGAAATC-1", "CGCTTTCATACCGGTG-1", "CGGAGTTTGAGAGACA-1",     "CGGCACTCAAGAAAGT-1", "CGGCCACGCACAAAGT-1", "CGGGATCAATGTAAGA-1",     "CGGGCAGCTAAACCGC-1", "CGGTGTACTTGATCCC-1", "CGTACCTGATAGGCCT-1",     "CGTAGCGAATTGTCAG-1", "CGTATTAAGAGATCTA-1", "CGTCAGTGCGCACAAG-1",     "CGTCTGGAAGGGCCCG-1", "CGTGTCCCATTCGCGA-1", "CGTTGTTTCAATTCCC-1",     "CGTTTAAGCGGAGCAC-1", "CGTTTCGCTCATTACA-1", "CTAATTCGCACGCGCT-1",     "CTAATTTCAACAACAC-1", "CTACGCACGGAGTACC-1", "CTACTGCCACCTGACC-1",     "CTAGCATAGTATAATG-1", "CTAGTGAAGGACAGGA-1", "CTAGTTGGGCCCGGTA-1",     "CTATCGGGTCTCAACA-1", "CTATGTCACTAGCCCA-1", "CTATGTCTATTGAAAC-1",     "CTCAGGACTCACCTGT-1", "CTCATGGCTCACAATC-1", "CTCATTTGATGGGCGG-1",     "CTCCGGCCTAATATGC-1", "CTCCTAAGTTATGTCT-1", "CTCGAGGTCGAACAGT-1",     "CTCGCACCTATATAGT-1", "CTCGGTCCGTAGCCTG-1", "CTCTCACAATCGATGA-1",     "CTCTGCGAAGCAAGCA-1", "CTCTGGACGCCTGGTG-1", "CTCTGTTTGAGGATTC-1",     "CTGAAAGAGATCCGAC-1", "CTGCACCTGGAACCGC-1", "CTGGAAATGGATGCTT-1",     "CTGTATGGTGTAGAAA-1", "CTGTTCACTGCCTGTG-1", "CTGTTCATCTCACGGG-1",     "CTTAACTTCGAAGTAC-1", "CTTACACGGTATTCCA-1", "CTTACACTGGGAAATA-1",     "CTTATGTTGACTACCA-1", "CTTCATAGCTCAAGAA-1", "CTTCTTACGTCGTATA-1",     "CTTGTACTTGTTGACT-1", "CTTGTTTATGTAGCCA-1", "CTTTAACTTTCAAAGG-1",     "CTTTGGCTTTAGTAAA-1", "GAAACAGCCATGCAGT-1", "GAAACTCTAATGAAGG-1",     "GAAAGTGACTAACTGC-1", "GAAATCGCGCGCAACT-1", "GAAATGGCGGTGTTAG-1",     "GAAATTCACATCGCTG-1", "GAACACACATCAACCA-1", "GAACGACCGAATGATA-1",     "GAAGAACGGTGCAGGT-1", "GAAGCGTGAGGAATTT-1", "GAAGCTTGCTGACCGC-1",     "GAAGTGCTGGATAGCT-1", "GAATCGACATGGTCAC-1", "GAATCGCCGGACACGG-1",     "GAATGTTGGGTAATCT-1", "GAATTATAGTGAAAGG-1", "GAATTTCTCGCTGCAG-1",     "GACAACGCAGCTTACG-1", "GACACAAGGGAAGAAA-1", "GACCAGAGCCCTGTAG-1",     "GACGCCTGTTGCAGGG-1", "GACGTGTAGGGATTAT-1", "GACTAAGATCATGCAC-1",     "GACTAAGTAGGCTCAC-1", "GAGACTGATGGGTAGA-1", "GAGATCTGCTTGGCAT-1",     "GAGATCTGTCACTCCG-1", "GAGCATCATCCCTGGG-1", "GAGCGAGGGAGTACCG-1",     "GAGCTCTCGGACCTAA-1", "GAGGAATGGAGAGGTT-1", "GAGGCTATCAAAGTCG-1",     "GAGGGCATCGCGTATC-1", "GAGGTACATCCATCTT-1", "GAGTACGGGTATACAA-1",     "GATAAATCGGTGGATG-1", "GATATCTCATGCAATA-1", "GATCATTCCAAACATT-1",     "GATCCCTTTATACTGC-1", "GATCGACACTATCTGA-1", "GATCGGTGGCCATAAC-1",     "GATGCTACAAGCGCCT-1", "GATGGCGCACACATTA-1", "GATGTTTGTGCGAGAT-1",     "GATTCCCTTGTCGCAG-1", "GATTCCGCGTTTCCGT-1", "GCAAATATTACGCTTT-1",     "GCAAGTGCACAGAGAA-1", "GCACTAGTCGCGCTAT-1", "GCAGGACTATAGAATA-1",     "GCATCCCTAACTTTGA-1", "GCCCACCAAGGCTGTC-1", "GCCCAGTTGGTATGCC-1",     "GCCCGCGCGTTTGACA-1", "GCCGAAATTCCTACGT-1", "GCCTAGCGATCTGACC-1",     "GCCTATTCCGATATAG-1", "GCCTTCAGCCCTACCG-1", "GCGAAGCCATACCCGT-1",     "GCGAGTTCTGCAAAGA-1", "GCGCTAATTGAATAGA-1", "GCGGGAACCAGGCCCT-1",     "GCGGTCCCTAGACGCA-1", "GCGTCGAAATGTCGGT-1", "GCTAATACCGAATGCC-1",     "GCTAGTTTCATTGAGG-1", "GCTCTAAACCCTGACG-1", "GCTCTATGTTACGTGC-1",     "GCTCTCGGGTACCGAA-1", "GCTGCTCTCCGGACAC-1", "GCTGGCATATTCACCT-1",     "GCTTGATGATAATCAG-1", "GGAGAAGTCATTGGCA-1", "GGAGTGCCGCCCTGGA-1",     "GGATTCAGTACGGTGG-1", "GGCAATAGTCAATGAG-1", "GGCACTCCACTGGGCA-1",     "GGCATACAGGTAGCGG-1", "GGCGGAGTAATATTAG-1", "GGCGGTAGGATCATTG-1",     "GGCGTCCTATCCGCTG-1", "GGCTATTAAGTTGTAT-1", "GGCTCGTGCCACCAGC-1",     "GGCTCTGCTCCAACGC-1", "GGCTGGCTAGCTTAAA-1", "GGGACTGCATAGATAG-1",     "GGGAGAACTCACAGTA-1", "GGGAGTTAATGAGGCG-1", "GGGCGATATGTGTGAA-1",     "GGGCGGTCCTATTGTC-1", "GGGCGTACATTTATAT-1", "GGGCTGCCTAGGGCGA-1",     "GGGCTGGTTAGTCGCG-1", "GGGTATGTATGCACTT-1", "GGGTCACCGTGACGGT-1",     "GGGTGTTTCAGCTATG-1", "GGTAACTATGTATCTG-1", "GGTAGTGCTCGCACCA-1",     "GGTATTGCCGAGTTTA-1", "GGTCGGTAATTAGACA-1", "GGTGAGATGCAGATAA-1",     "GGTGCGGATAAGTGGC-1", "GGTTAGGCTTGGAGAA-1", "GGTTTACAATCTCAAT-1",     "GGTTTGAGTGCTGGAA-1", "GTAAGCGGGCAGTCAG-1", "GTAAGTAGGGTATACC-1",     "GTACGAGATTGCGACA-1", "GTACTAAGATTTGGAG-1", "GTACTGAGGTCGTAAC-1",     "GTACTGCATGAAGCGT-1", "GTACTTGGGCACTTCT-1", "GTAGACGTCGTTACAT-1",     "GTATAGGACTCAGTAG-1", "GTATCAGCTTGGGTTC-1", "GTCAAAGTTTACATAG-1",     "GTCACTCTCCAAATCT-1", "GTCAGAATAGTCTATG-1", "GTCATGCACCTCCGTT-1",     "GTCATGGACATGACTA-1", "GTCCCAACGTAAAGTA-1", "GTCCGGCTGAATTGCG-1",     "GTCGGATATCTCAGAC-1", "GTCGTCAATTATAAGG-1", "GTGAAACGTGCTCCAC-1",     "GTGCAACAAATGTGGC-1", "GTGCACGAAAGTGACT-1", "GTGCGGGTCTCCAAAT-1",     "GTGCTCAAGTACTGTC-1", "GTGGACGTGCTGAGAC-1", "GTGGGCTTAGACACAC-1",     "GTTCACAGGAGTCTAG-1", "GTTCCAGTCTGACCAT-1", "GTTCTTCCCTCGATGT-1",     "GTTGCACGGAGTTTCG-1", "GTTGTCGTGTTAGTTG-1", "GTTTCCTGGAGGGTGA-1",     "GTTTGACCAAATCCTA-1", "GTTTGGCCGCTCAGCG-1", "TAAATGAATCCGTTTC-1",     "TAACAAAGGGAGAAGC-1", "TAACGCTTTGAGAGCG-1", "TAACTATCGAAGGTCC-1",     "TAACTCCATGGAGGCT-1", "TAAGGAACTTGTGGGA-1", "TAAGGCAACATAAGAT-1",     "TAATACACAGTAGTAT-1", "TAATATTGAAATTCGC-1", "TACAAGTCTCGTGCAT-1",     "TACCAGCTAGGTTTAA-1", "TACCTCACGCTTGTAC-1", "TACGAACACGACTTCA-1",     "TACGAGAACTTCACGT-1", "TACGATGTTGATCATC-1", "TACGTGCACTATGCTG-1",     "TAGAGCTACGAAGAAC-1", "TAGAGTCTAAGCGAAC-1", "TAGCTAAGTCCGGGAG-1",     "TAGGGTGTTTCAAGAG-1", "TAGGTGAGCCCTACTC-1", "TAGGTTCGAGTTCGTC-1",     "TAGTGCCCTCCAGAGT-1", "TAGTTTATTCTTGCTT-1", "TATAAATCCACAAGCT-1",     "TATAAGTGAGGATAGC-1", "TATACACGCAAAGTAT-1", "TATCCATATCATGCGA-1",     "TATCCTGCATGGGAAT-1", "TATCGATGATTAAACG-1", "TATGCTCCCTACTTAC-1",     "TATGGATGTGCTACGC-1", "TATGGCCCGGCCTCGC-1", "TATGGGTACGTATCGT-1",     "TATTCAATTCTAATCC-1", "TATTCGTGCCAGAATA-1", "TATTTAGTCTAGATCG-1",     "TATTTGTTACCCTTTA-1", "TCAAACAACCGCGTCG-1", "TCAAACTTAGATTGTT-1",     "TCAACACATTGGGTAA-1", "TCAACCATGTTCGGGC-1", "TCACAGGTTATTGGGC-1",     "TCACGCATTGTAGATC-1", "TCACGGTCATCGCACA-1", "TCAGCAAATGCATCTC-1",     "TCAGGGTGTAACGTAA-1", "TCAGGTTCTTTGAGAA-1", "TCAGTACTGACCCGCG-1",     "TCATCACTCGAGCTCG-1", "TCATCCCAGAGGGTGG-1", "TCCACAATGGTTTACG-1",     "TCCAGGCGAGTACGGT-1", "TCCCAAAGACGAAGGA-1", "TCCCAGCTTTAGTCTG-1",     "TCCCAGGCTTAGCTAA-1", "TCCCGCGTACTCCTGG-1", "TCCCGGGTGTGCTGCT-1",     "TCCCGTCAGTCCCGCA-1", "TCCCTGGCTCGCTGGA-1", "TCCGATAATTGCCATA-1",     "TCCGATGACTGAGCTC-1", "TCCGATGGTGCGACAT-1", "TCCGATTACATTGCCG-1",     "TCCTCTACGAGATGGC-1", "TCCTTCAGTGGTCGAA-1", "TCCTTTCTTACGCTTA-1",     "TCGAAATTTAGGACCA-1", "TCGCATCCCTAAGTGT-1", "TCGCCGGTCGATCCGT-1",     "TCGCGTCCAGAAGGTC-1", "TCGCTCGATATATTCC-1", "TCGCTGGGCGGATTGT-1",     "TCGCTGTGCGTAAATC-1", "TCGGAATGCGCTCTGA-1", "TCGGAGTACATGAGTA-1",     "TCGTATAGTGCAATTA-1", "TCGTCCGCTGGCGTCT-1", "TCTACCGTCCACAAGC-1",     "TCTAGCATCTTCGATG-1", "TCTAGTTATCAGAAGA-1", "TCTATCATGCAGTTAC-1",     "TCTATCGGTCGCAACA-1", "TCTATTACGCTGGCGA-1", "TCTCTAATAGCTGGTA-1",     "TCTGAACCGGTCGGCT-1", "TCTTACAGAGGTACCG-1", "TCTTACGGCATCCGAC-1",     "TCTTCGATACCAATAA-1", "TCTTGATGCGTAGCGA-1", "TCTTGGTAACACCAAA-1",     "TGAAACTTATGCAAGC-1", "TGACCCACGTTAGACA-1", "TGAGAATGCTTTACCG-1",     "TGATCGGTTTGACCCT-1", "TGATCTCCGGCGCCAG-1", "TGATTCGTCTATCACT-1",     "TGATTTCCTCCTGACG-1", "TGCAAGAATGACGTAA-1", "TGCAGGATCGGCAAAG-1",     "TGCATGGATCGGATCT-1", "TGCCACCTGGCGAAAC-1", "TGCCTGATCAAACGAT-1",     "TGCGCAAAGCATTTGG-1", "TGCGCGATTAACGGAG-1", "TGCGGAGTAAAGGTGC-1",     "TGCGGCATAGTTCAAC-1", "TGCGGTGAAATTTCAT-1", "TGCTGGTTGGACAATT-1",     "TGGAAGACGAACACCA-1", "TGGAAGGATAAAGATG-1", "TGGAGTGATGCGATGA-1",     "TGGCCAAACTGAAGTA-1", "TGGCCGTATATTGACC-1", "TGGCTTTGGGTAGACA-1",     "TGGGAAATGCCTTTCC-1", "TGGGCACGTTCTATGG-1", "TGGGCCACAAGAGCGC-1",     "TGGTTCGTAGCAAAGG-1", "TGTAGTGATCTATAAT-1", "TGTCCCGACATAGCAC-1",     "TGTCCTAAGTCACCGC-1", "TGTCTACAGTTTCTGT-1", "TGTGACTACGCCAGTC-1",     "TGTGTCGCGAGTTGCA-1", "TGTTCTTCCATTGACT-1", "TGTTTCGGTACTTCTC-1",     "TTAACCAACCCTCCCT-1", "TTAAGCGCCTGACCCA-1", "TTAATCAGTACGTCAG-1",     "TTAATGTAGACCAGGT-1", "TTAATTTCAGACGCGG-1", "TTACAACTACGCATCC-1",     "TTACATCGTGGCCTGG-1", "TTACCATTGATTACCC-1", "TTAGCAACATGGATGT-1",     "TTATATTTGGCAATCC-1", "TTATCCAATCGAACTC-1", "TTATCTGTATCATAAC-1",     "TTATGACAAACTGGAT-1", "TTATTAGGGAAGCATC-1", "TTCAAAGTCTCTAGCC-1",     "TTCAACGACCCGACCG-1", "TTCAAGCCGAGCTGAG-1", "TTCATGGCGCAACAGG-1",     "TTCCTCGAGGGTGTCT-1", "TTCGACGGGAAGGGCG-1", "TTCGCACTCGCGTGCT-1",     "TTCGTACTCCAGAACG-1", "TTCTACTTGCGAGGGC-1", "TTCTAGGCCAATTGTG-1",     "TTCTTAGTGAACGGTG-1", "TTCTTAGTGGCTCAGA-1", "TTCTTGTAACCTAATG-1",     "TTGAAGAATTCCCAGG-1", "TTGAATATGGACTTTC-1", "TTGATCTAACTTTGTC-1",     "TTGATTATGCAGATGA-1", "TTGCACGGAGCAGCAC-1", "TTGCTCCCATACCGGA-1",     "TTGCTGAAGGAACCAC-1", "TTGCTGATCATGTTCG-1", "TTGGATTGGGTACCAC-1",     "TTGGCTCGCATGAGAC-1", "TTGGGACACTGCCCGC-1", "TTGGGCGGCGGTTGCC-1",     "TTGTAATCCGTACTCG-1", "TTGTATCACACAGAAT-1", "TTGTCGTTCAGTTACC-1",     "TTGTGGCCCTGACAGT-1", "TTGTTCAGTGTGCTAC-1"), in_tissue = c(1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), nr_feats = c(227L, 283L,     248L, 222L, 313L, 338L, 345L, 247L, 211L, 237L, 254L, 264L,     311L, 343L, 429L, 265L, 252L, 273L, 334L, 413L, 195L, 375L,     382L, 256L, 165L, 264L, 322L, 220L, 253L, 156L, 324L, 270L,     276L, 335L, 309L, 296L, 238L, 156L, 426L, 238L, 223L, 345L,     380L, 245L, 306L, 206L, 247L, 220L, 343L, 338L, 441L, 247L,     239L, 276L, 245L, 402L, 254L, 274L, 304L, 147L, 283L, 271L,     216L, 237L, 373L, 216L, 201L, 257L, 332L, 208L, 367L, 336L,     239L, 314L, 278L, 311L, 261L, 176L, 302L, 285L, 256L, 299L,     237L, 287L, 246L, 320L, 404L, 364L, 207L, 332L, 255L, 186L,     316L, 253L, 344L, 321L, 318L, 317L, 239L, 251L, 379L, 307L,     296L, 270L, 288L, 318L, 272L, 319L, 335L, 305L, 388L, 292L,     188L, 313L, 274L, 212L, 270L, 226L, 236L, 393L, 210L, 199L,     281L, 284L, 302L, 307L, 180L, 300L, 318L, 226L, 219L, 219L,     293L, 272L, 281L, 230L, 267L, 230L, 247L, 308L, 233L, 247L,     214L, 214L, 305L, 279L, 279L, 278L, 300L, 215L, 258L, 274L,     220L, 262L, 225L, 150L, 154L, 332L, 308L, 400L, 200L, 173L,     335L, 356L, 331L, 256L, 161L, 274L, 328L, 303L, 386L, 253L,     267L, 284L, 186L, 294L, 285L, 309L, 426L, 280L, 298L, 194L,     228L, 309L, 296L, 135L, 240L, 215L, 250L, 367L, 362L, 305L,     318L, 220L, 245L, 255L, 289L, 317L, 214L, 246L, 308L, 343L,     261L, 332L, 342L, 272L, 300L, 309L, 321L, 232L, 374L, 230L,     426L, 222L, 194L, 308L, 372L, 236L, 363L, 209L, 305L, 255L,     285L, 270L, 288L, 234L, 312L, 206L, 361L, 322L, 207L, 317L,     265L, 338L, 189L, 272L, 300L, 362L, 299L, 277L, 331L, 251L,     199L, 232L, 242L, 319L, 342L, 162L, 303L, 294L, 299L, 204L,     297L, 359L, 235L, 194L, 187L, 252L, 324L, 219L, 265L, 283L,     343L, 341L, 231L, 227L, 314L, 195L, 305L, 322L, 358L, 336L,     172L, 285L, 346L, 260L, 288L, 238L, 263L, 244L, 223L, 335L,     325L, 230L, 268L, 323L, 267L, 308L, 256L, 277L, 309L, 325L,     267L, 271L, 182L, 334L, 169L, 272L, 270L, 234L, 342L, 248L,     164L, 94L, 342L, 243L, 210L, 273L, 160L, 277L, 208L, 235L,     261L, 285L, 323L, 288L, 253L, 314L, 313L, 388L, 279L, 222L,     270L, 315L, 255L, 272L, 327L, 306L, 161L, 248L, 407L, 219L,     209L, 273L, 219L, 207L, 197L, 225L, 361L, 309L, 152L, 288L,     368L, 281L, 239L, 347L, 294L, 287L, 164L, 223L, 288L, 362L,     194L, 286L, 347L, 345L, 294L, 279L, 316L, 286L, 243L, 403L,     359L, 280L, 315L, 302L, 143L, 202L, 325L, 232L, 204L, 287L,     261L, 340L, 188L, 235L, 300L, 313L, 230L, 244L, 276L, 222L,     194L, 350L, 280L, 185L, 235L, 262L, 245L, 338L, 259L, 239L,     361L, 336L, 170L, 306L, 214L, 231L, 279L, 323L, 214L, 205L,     291L, 267L, 153L, 242L, 244L, 422L, 307L, 284L, 253L, 179L,     289L, 239L, 280L, 269L, 329L, 353L, 284L, 294L, 188L, 215L,     331L, 337L, 302L, 185L, 261L, 340L, 117L, 221L, 168L, 215L,     215L, 205L, 289L, 334L, 288L, 197L, 268L, 355L, 309L, 245L,     321L, 237L, 303L, 294L, 250L, 263L, 302L, 179L, 403L, 246L,     255L, 218L, 179L, 356L, 163L, 232L, 369L, 278L, 204L, 369L,     319L, 227L, 357L, 266L, 224L, 295L, 151L, 314L, 316L, 368L,     234L, 256L, 233L, 308L, 319L, 163L, 174L, 216L, 370L, 355L,     226L, 90L, 252L, 264L, 282L, 247L, 364L, 233L, 445L, 217L,     318L, 261L, 386L, 336L, 303L, 256L, 197L, 258L, 282L, 318L,     207L, 237L, 321L, 294L, 171L, 218L, 302L, 367L, 234L, 338L,     256L, 263L, 330L, 346L, 238L, 164L, 303L, 303L, 383L, 210L,     234L, 195L, 234L, 257L, 230L, 267L, 121L, 165L, 164L, 233L,     314L, 310L, 246L, 344L, 286L, 218L, 146L, 171L, 310L, 259L,     225L, 277L, 245L, 212L, 348L, 309L, 373L, 241L, 215L, 398L,     321L, 288L, 349L, 375L, 326L, 234L, 374L, 242L, 283L, 156L,     342L, 379L, 319L, 429L, 171L, 267L, 315L, 278L, 297L, 192L,     254L, 267L, 231L, 253L, 383L, 289L, 249L, 265L, 367L, 277L,     294L, 330L, 306L, 218L, 185L, 193L, 266L, 186L, 263L, 266L,     207L, 332L, 384L, 190L, 281L, 303L, 263L, 405L, 389L, 297L,     153L, 126L, 287L, 359L, 103L, 289L, 380L, 272L, 178L, 313L,     331L, 355L, 227L, 253L, 188L, 364L, 335L, 315L, 240L, 220L,     204L, 237L), perc_feats = c(35.8044164037855, 44.6372239747634,     39.1167192429022, 35.0157728706625, 49.3690851735016, 53.3123028391167,     54.416403785489, 38.9589905362776, 33.2807570977918, 37.3817034700315,     40.0630914826498, 41.6403785488959, 49.0536277602524, 54.1009463722397,     67.6656151419558, 41.7981072555205, 39.7476340694006, 43.0599369085174,     52.6813880126183, 65.1419558359621, 30.7570977917981, 59.1482649842271,     60.2523659305994, 40.3785488958991, 26.0252365930599, 41.6403785488959,     50.788643533123, 34.7003154574132, 39.9053627760252, 24.6056782334385,     51.1041009463722, 42.5867507886435, 43.5331230283912, 52.8391167192429,     48.7381703470032, 46.6876971608833, 37.5394321766561, 24.6056782334385,     67.192429022082, 37.5394321766561, 35.1735015772871, 54.416403785489,     59.9369085173502, 38.6435331230284, 48.2649842271293, 32.4921135646688,     38.9589905362776, 34.7003154574132, 54.1009463722397, 53.3123028391167,     69.5583596214511, 38.9589905362776, 37.6971608832808, 43.5331230283912,     38.6435331230284, 63.4069400630915, 40.0630914826498, 43.217665615142,     47.9495268138801, 23.186119873817, 44.6372239747634, 42.7444794952681,     34.0694006309148, 37.3817034700315, 58.8328075709779, 34.0694006309148,     31.7034700315457, 40.5362776025237, 52.3659305993691, 32.807570977918,     57.8864353312303, 52.9968454258675, 37.6971608832808, 49.5268138801262,     43.8485804416404, 49.0536277602524, 41.1671924290221, 27.7602523659306,     47.6340694006309, 44.9526813880126, 40.3785488958991, 47.1608832807571,     37.3817034700315, 45.2681388012618, 38.801261829653, 50.4731861198738,     63.7223974763407, 57.4132492113565, 32.6498422712934, 52.3659305993691,     40.2208201892744, 29.3375394321767, 49.8422712933754, 39.9053627760252,     54.2586750788643, 50.6309148264984, 50.1577287066246, 50,     37.6971608832808, 39.589905362776, 59.7791798107255, 48.4227129337539,     46.6876971608833, 42.5867507886435, 45.4258675078864, 50.1577287066246,     42.9022082018927, 50.3154574132492, 52.8391167192429, 48.1072555205047,     61.198738170347, 46.0567823343849, 29.6529968454259, 49.3690851735016,     43.217665615142, 33.4384858044164, 42.5867507886435, 35.6466876971609,     37.2239747634069, 61.98738170347, 33.1230283911672, 31.3880126182965,     44.3217665615142, 44.794952681388, 47.6340694006309, 48.4227129337539,     28.391167192429, 47.3186119873817, 50.1577287066246, 35.6466876971609,     34.5425867507886, 34.5425867507886, 46.2145110410095, 42.9022082018927,     44.3217665615142, 36.2776025236593, 42.1135646687697, 36.2776025236593,     38.9589905362776, 48.5804416403785, 36.7507886435331, 38.9589905362776,     33.7539432176656, 33.7539432176656, 48.1072555205047, 44.006309148265,     44.006309148265, 43.8485804416404, 47.3186119873817, 33.9116719242902,     40.6940063091483, 43.217665615142, 34.7003154574132, 41.3249211356467,     35.4889589905363, 23.6593059936909, 24.2902208201893, 52.3659305993691,     48.5804416403785, 63.0914826498423, 31.5457413249211, 27.2870662460568,     52.8391167192429, 56.1514195583596, 52.2082018927445, 40.3785488958991,     25.3943217665615, 43.217665615142, 51.7350157728707, 47.7917981072555,     60.8832807570978, 39.9053627760252, 42.1135646687697, 44.794952681388,     29.3375394321767, 46.3722397476341, 44.9526813880126, 48.7381703470032,     67.192429022082, 44.1640378548896, 47.0031545741325, 30.5993690851735,     35.9621451104101, 48.7381703470032, 46.6876971608833, 21.2933753943218,     37.8548895899054, 33.9116719242902, 39.4321766561514, 57.8864353312303,     57.0977917981073, 48.1072555205047, 50.1577287066246, 34.7003154574132,     38.6435331230284, 40.2208201892744, 45.583596214511, 50,     33.7539432176656, 38.801261829653, 48.5804416403785, 54.1009463722397,     41.1671924290221, 52.3659305993691, 53.9432176656151, 42.9022082018927,     47.3186119873817, 48.7381703470032, 50.6309148264984, 36.5930599369085,     58.9905362776025, 36.2776025236593, 67.192429022082, 35.0157728706625,     30.5993690851735, 48.5804416403785, 58.6750788643533, 37.2239747634069,     57.2555205047319, 32.9652996845426, 48.1072555205047, 40.2208201892744,     44.9526813880126, 42.5867507886435, 45.4258675078864, 36.9085173501577,     49.211356466877, 32.4921135646688, 56.9400630914826, 50.788643533123,     32.6498422712934, 50, 41.7981072555205, 53.3123028391167,     29.8107255520505, 42.9022082018927, 47.3186119873817, 57.0977917981073,     47.1608832807571, 43.6908517350158, 52.2082018927445, 39.589905362776,     31.3880126182965, 36.5930599369085, 38.1703470031546, 50.3154574132492,     53.9432176656151, 25.5520504731861, 47.7917981072555, 46.3722397476341,     47.1608832807571, 32.1766561514196, 46.8454258675079, 56.6246056782334,     37.0662460567823, 30.5993690851735, 29.4952681388013, 39.7476340694006,     51.1041009463722, 34.5425867507886, 41.7981072555205, 44.6372239747634,     54.1009463722397, 53.7854889589905, 36.4353312302839, 35.8044164037855,     49.5268138801262, 30.7570977917981, 48.1072555205047, 50.788643533123,     56.4668769716088, 52.9968454258675, 27.1293375394322, 44.9526813880126,     54.5741324921136, 41.0094637223975, 45.4258675078864, 37.5394321766561,     41.4826498422713, 38.4858044164038, 35.1735015772871, 52.8391167192429,     51.2618296529969, 36.2776025236593, 42.2712933753943, 50.9463722397476,     42.1135646687697, 48.5804416403785, 40.3785488958991, 43.6908517350158,     48.7381703470032, 51.2618296529969, 42.1135646687697, 42.7444794952681,     28.7066246056782, 52.6813880126183, 26.6561514195584, 42.9022082018927,     42.5867507886435, 36.9085173501577, 53.9432176656151, 39.1167192429022,     25.8675078864353, 14.8264984227129, 53.9432176656151, 38.3280757097792,     33.1230283911672, 43.0599369085174, 25.2365930599369, 43.6908517350158,     32.807570977918, 37.0662460567823, 41.1671924290221, 44.9526813880126,     50.9463722397476, 45.4258675078864, 39.9053627760252, 49.5268138801262,     49.3690851735016, 61.198738170347, 44.006309148265, 35.0157728706625,     42.5867507886435, 49.6845425867508, 40.2208201892744, 42.9022082018927,     51.5772870662461, 48.2649842271293, 25.3943217665615, 39.1167192429022,     64.1955835962145, 34.5425867507886, 32.9652996845426, 43.0599369085174,     34.5425867507886, 32.6498422712934, 31.0725552050473, 35.4889589905363,     56.9400630914826, 48.7381703470032, 23.9747634069401, 45.4258675078864,     58.0441640378549, 44.3217665615142, 37.6971608832808, 54.7318611987382,     46.3722397476341, 45.2681388012618, 25.8675078864353, 35.1735015772871,     45.4258675078864, 57.0977917981073, 30.5993690851735, 45.1104100946372,     54.7318611987382, 54.416403785489, 46.3722397476341, 44.006309148265,     49.8422712933754, 45.1104100946372, 38.3280757097792, 63.5646687697161,     56.6246056782334, 44.1640378548896, 49.6845425867508, 47.6340694006309,     22.5552050473186, 31.8611987381703, 51.2618296529969, 36.5930599369085,     32.1766561514196, 45.2681388012618, 41.1671924290221, 53.6277602523659,     29.6529968454259, 37.0662460567823, 47.3186119873817, 49.3690851735016,     36.2776025236593, 38.4858044164038, 43.5331230283912, 35.0157728706625,     30.5993690851735, 55.205047318612, 44.1640378548896, 29.1798107255521,     37.0662460567823, 41.3249211356467, 38.6435331230284, 53.3123028391167,     40.8517350157729, 37.6971608832808, 56.9400630914826, 52.9968454258675,     26.813880126183, 48.2649842271293, 33.7539432176656, 36.4353312302839,     44.006309148265, 50.9463722397476, 33.7539432176656, 32.3343848580442,     45.8990536277603, 42.1135646687697, 24.1324921135647, 38.1703470031546,     38.4858044164038, 66.5615141955836, 48.4227129337539, 44.794952681388,     39.9053627760252, 28.2334384858044, 45.583596214511, 37.6971608832808,     44.1640378548896, 42.4290220820189, 51.8927444794953, 55.6782334384858,     44.794952681388, 46.3722397476341, 29.6529968454259, 33.9116719242902,     52.2082018927445, 53.1545741324921, 47.6340694006309, 29.1798107255521,     41.1671924290221, 53.6277602523659, 18.4542586750789, 34.8580441640379,     26.4984227129338, 33.9116719242902, 33.9116719242902, 32.3343848580442,     45.583596214511, 52.6813880126183, 45.4258675078864, 31.0725552050473,     42.2712933753943, 55.993690851735, 48.7381703470032, 38.6435331230284,     50.6309148264984, 37.3817034700315, 47.7917981072555, 46.3722397476341,     39.4321766561514, 41.4826498422713, 47.6340694006309, 28.2334384858044,     63.5646687697161, 38.801261829653, 40.2208201892744, 34.384858044164,     28.2334384858044, 56.1514195583596, 25.7097791798107, 36.5930599369085,     58.2018927444795, 43.8485804416404, 32.1766561514196, 58.2018927444795,     50.3154574132492, 35.8044164037855, 56.3091482649842, 41.9558359621451,     35.3312302839117, 46.5299684542587, 23.8170347003155, 49.5268138801262,     49.8422712933754, 58.0441640378549, 36.9085173501577, 40.3785488958991,     36.7507886435331, 48.5804416403785, 50.3154574132492, 25.7097791798107,     27.4447949526814, 34.0694006309148, 58.3596214511041, 55.993690851735,     35.6466876971609, 14.1955835962145, 39.7476340694006, 41.6403785488959,     44.4794952681388, 38.9589905362776, 57.4132492113565, 36.7507886435331,     70.1892744479495, 34.2271293375394, 50.1577287066246, 41.1671924290221,     60.8832807570978, 52.9968454258675, 47.7917981072555, 40.3785488958991,     31.0725552050473, 40.6940063091483, 44.4794952681388, 50.1577287066246,     32.6498422712934, 37.3817034700315, 50.6309148264984, 46.3722397476341,     26.9716088328076, 34.384858044164, 47.6340694006309, 57.8864353312303,     36.9085173501577, 53.3123028391167, 40.3785488958991, 41.4826498422713,     52.0504731861199, 54.5741324921136, 37.5394321766561, 25.8675078864353,     47.7917981072555, 47.7917981072555, 60.410094637224, 33.1230283911672,     36.9085173501577, 30.7570977917981, 36.9085173501577, 40.5362776025237,     36.2776025236593, 42.1135646687697, 19.0851735015773, 26.0252365930599,     25.8675078864353, 36.7507886435331, 49.5268138801262, 48.8958990536278,     38.801261829653, 54.2586750788643, 45.1104100946372, 34.384858044164,     23.0283911671924, 26.9716088328076, 48.8958990536278, 40.8517350157729,     35.4889589905363, 43.6908517350158, 38.6435331230284, 33.4384858044164,     54.8895899053628, 48.7381703470032, 58.8328075709779, 38.01261829653,     33.9116719242902, 62.7760252365931, 50.6309148264984, 45.4258675078864,     55.0473186119874, 59.1482649842271, 51.4195583596215, 36.9085173501577,     58.9905362776025, 38.1703470031546, 44.6372239747634, 24.6056782334385,     53.9432176656151, 59.7791798107255, 50.3154574132492, 67.6656151419558,     26.9716088328076, 42.1135646687697, 49.6845425867508, 43.8485804416404,     46.8454258675079, 30.2839116719243, 40.0630914826498, 42.1135646687697,     36.4353312302839, 39.9053627760252, 60.410094637224, 45.583596214511,     39.2744479495268, 41.7981072555205, 57.8864353312303, 43.6908517350158,     46.3722397476341, 52.0504731861199, 48.2649842271293, 34.384858044164,     29.1798107255521, 30.4416403785489, 41.9558359621451, 29.3375394321767,     41.4826498422713, 41.9558359621451, 32.6498422712934, 52.3659305993691,     60.5678233438486, 29.9684542586751, 44.3217665615142, 47.7917981072555,     41.4826498422713, 63.8801261829653, 61.3564668769716, 46.8454258675079,     24.1324921135647, 19.8738170347003, 45.2681388012618, 56.6246056782334,     16.2460567823344, 45.583596214511, 59.9369085173502, 42.9022082018927,     28.0757097791798, 49.3690851735016, 52.2082018927445, 55.993690851735,     35.8044164037855, 39.9053627760252, 29.6529968454259, 57.4132492113565,     52.8391167192429, 49.6845425867508, 37.8548895899054, 34.7003154574132,     32.1766561514196, 37.3817034700315), total_expr = c(995.149096792437,     1123.87508033663, 1027.6321805555, 974.376567478206, 1174.48813214033,     1229.93527606038, 1270.28563239995, 911.358403783163, 951.278708856921,     1010.97394950848, 1031.55012182124, 976.976512221084, 1100.76640600405,     1088.97048616803, 1268.92451905763, 1057.93080366493, 1046.8524438131,     1068.14073423384, 1029.37419841557, 1105.3542575398, 756.067467900412,     1226.03463138342, 1155.08825338693, 1066.28854180742, 815.498835625545,     1056.53511523766, 1179.0790182727, 988.987323330922, 1089.74388658061,     793.354500897721, 1234.04065314962, 939.205312776025, 989.123814558527,     1256.25597644083, 1190.09953303933, 1122.35333613882, 1028.66590393627,     781.949571466476, 1339.25877446059, 967.449468638109, 953.496586802366,     1176.94349603553, 1131.60545555067, 1034.86098481125, 1177.43376394568,     930.407705288658, 932.59587539103, 990.257710633478, 1182.5976434172,     1171.79324237454, 1246.64446743499, 879.274403476489, 1009.20914139158,     1015.4549051031, 1029.86566230484, 1297.5859467352, 1053.69215832588,     1061.16875116526, 1140.99263931776, 762.425471437489, 1123.15058896042,     1073.44882400168, 910.352143707052, 964.932043522723, 1243.33446019895,     969.721356232917, 925.87200043943, 1060.2485945324, 1009.49692112379,     949.664355910048, 1250.53681267894, 1163.45688783831, 1038.46906418224,     1206.46750914705, 991.456279654484, 997.357131555138, 1091.84447670155,     842.928225081166, 1116.15772616475, 1040.45095086527, 1047.28964953029,     984.388443406971, 1037.06548787166, 1077.06196325346, 1037.68441755375,     1051.69796298035, 1327.30836251613, 1039.29757135186, 904.329956178986,     1091.65587964113, 981.038125067427, 824.276320730058, 1102.33125186862,     1028.35881178172, 1209.35223391061, 1148.32914591959, 1083.15493615302,     1146.2406548861, 1014.11466164354, 991.262995515024, 1220.06008923083,     1035.91661992888, 1118.10865459096, 1079.66071087084, 1102.18515789839,     1039.80144658185, 981.175500598067, 1166.45706090109, 1075.41582057356,     1167.15480280757, 1250.68584484001, 1059.94932393213, 874.975779798282,     1226.90950108715, 1030.72600386232, 915.995397257769, 1082.15889424651,     796.797152261997, 949.767098331134, 1251.26953154754, 946.327287119748,     902.47967297667, 1093.07423641624, 993.723176759765, 1086.96002286032,     1101.72369558932, 874.18669045772, 1058.74128516271, 1043.77200280931,     979.256753684475, 994.252510111631, 965.131609807606, 1051.11065199898,     911.68717695085, 1029.66398915018, 1012.21390996071, 1074.9051637742,     816.45032585686, 1041.5119999693, 1185.12339974744, 1027.0039502557,     1044.96962277558, 973.289959937676, 972.421777903221, 1104.53785971538,     1098.66323572614, 1128.04058035795, 1077.44118531241, 1112.47581432683,     967.213656150161, 1059.32497946659, 1107.22339286351, 978.298967013498,     1064.92393629462, 998.985204174203, 762.736582180235, 768.19107115664,     1249.53981791273, 1063.27908177887, 1202.63378098563, 899.94774392546,     831.877936592803, 1197.70658270915, 1094.41394652832, 1218.66288905589,     1026.48244857002, 802.000247573436, 1079.60961071875, 1167.82968317017,     1045.36469143145, 1256.67535272578, 1021.26155176075, 996.438613040115,     1109.41685925577, 881.136760565745, 1144.59915918244, 1074.26749751456,     1151.96347122388, 1263.57735936609, 1090.02426021202, 963.007769913904,     888.655889125914, 1009.54228175303, 1164.89581008832, 1018.559197982,     718.783693994834, 1010.07177391748, 725.149453415772, 1061.18518659809,     1249.60634396386, 1183.23244442763, 1150.64629611661, 1199.89998878901,     951.960642070437, 974.978142813579, 1072.11444016205, 1047.17760689003,     1183.63705593967, 961.81223692073, 827.478967927277, 982.566188568405,     1245.75136438842, 1092.90658864008, 1182.49033594431, 1092.54156165692,     1105.71239098483, 1108.22590525823, 1077.4877257274, 1146.78750735264,     1017.21585608919, 1281.53500730393, 1000.50322081535, 1140.00555363042,     969.383062179595, 897.610768041515, 1066.18757250121, 1331.85967041812,     920.746393913121, 1183.31269178798, 966.046657217733, 1114.29204658826,     790.034865817742, 974.644243700053, 1059.39994447989, 1116.37362691255,     864.811346199657, 1154.660976118, 941.054968199412, 1064.42130722015,     1047.71300246821, 943.227134512262, 1208.65301299366, 1082.26480812132,     1073.85033332098, 887.733713110524, 1078.25333746944, 1105.33441066995,     1183.05894170535, 1159.87658073433, 1048.87474981579, 1190.06507155507,     1056.14527958072, 910.686413887409, 976.646990022743, 986.520714576875,     1168.69541366701, 1169.35103006125, 808.938057330349, 783.495441198591,     1011.17894091539, 1157.09634472768, 914.755009467248, 1115.0769646888,     1087.0296586183, 932.924996683601, 912.065245677275, 765.244259771295,     1056.42425695882, 1242.45267172221, 983.093699078488, 1003.17064288329,     1133.9984313882, 1167.31197697986, 1253.85841977331, 1006.51378103346,     969.889826489821, 1068.52210868364, 897.465345402437, 1175.98623243131,     1199.69139626093, 1161.62239654015, 1206.59900637643, 824.10647622782,     1097.21257055963, 1193.49576339842, 1020.33407578144, 1016.20993616951,     1028.83243586842, 1093.70256385779, 1034.16320088611, 862.664184078829,     1225.41750603617, 1019.643317054, 887.738477818594, 1037.20609351039,     1239.8653801032, 994.248137991758, 1132.38729012393, 1033.18208806733,     1124.4934597053, 1118.78532067903, 1073.58636653789, 1040.50931670625,     1062.65436319572, 833.840834429282, 1145.35752804772, 839.880270676415,     1082.88127812057, 1085.7152585323, 986.199105223589, 1204.76267009876,     1010.97815876173, 813.837792380424, 537.345040594604, 1274.51043809521,     1005.43240774771, 951.2075962164, 1053.42500333923, 799.716564725356,     1053.69951440951, 960.009348355673, 862.504075516725, 1068.48522503977,     1128.21839539972, 1166.332620896, 1127.59983317868, 1064.99897931481,     1164.81656260949, 1156.92825690358, 1251.40886194961, 905.478299486095,     839.650460974747, 905.506954728022, 1204.18293937228, 1055.54324220363,     1076.89040667037, 1191.25260600872, 1140.55468230533, 792.313275182551,     1007.21466768834, 1079.391325878, 908.925155808981, 918.483300081163,     1102.54195190034, 964.929382967323, 843.56037114742, 915.929017571279,     1016.31598818642, 1099.09157056719, 1187.76678581518, 740.499824108306,     1110.31130985468, 1095.46166682945, 1022.16153973173, 922.663983867245,     1227.97832183403, 1142.76638220009, 1116.66917888116, 816.725892205788,     997.009072485337, 1142.55613683472, 1244.04227569505, 913.518215638606,     1112.68497726844, 1158.35178800196, 1087.24065765121, 1099.31107697547,     1094.66383786885, 993.334092333323, 1077.96681022923, 916.624103381498,     1239.55754411369, 1221.95260468809, 1114.7952012935, 1026.50029872949,     1061.07565480101, 737.847004777386, 934.572492553464, 1195.22327670301,     1013.86641796191, 950.65418604334, 1121.19086515392, 1049.28049129694,     1145.3084016355, 817.232105937776, 1011.62527071581, 1115.3919014715,     1165.4604071627, 993.017509138127, 949.439613887734, 1143.39816936969,     977.801139523569, 900.377271401379, 1048.72035277429, 1126.18725583115,     901.965833272067, 1004.18125408655, 1076.92277301814, 937.900704989224,     1102.5019886136, 1048.06892263228, 1012.92745524221, 1243.76663941583,     1109.04195902258, 829.806846955287, 1173.04947254143, 974.584517178106,     979.070156812593, 1064.7493483799, 1096.74586562509, 948.610609118514,     825.502438512828, 1131.43934399216, 995.525534917057, 776.515440845001,     978.531700696874, 1018.3612222727, 997.46421906072, 1148.02297529548,     1089.96394317106, 995.335936406541, 768.474865707, 1141.06750519548,     1005.50410423681, 1016.58176879078, 1070.41252733439, 1089.15667057228,     1143.09315966995, 1070.59509886759, 1157.16019334958, 876.905680269733,     970.85939387737, 1150.42038746921, 1181.19902565436, 1141.04191337342,     869.841203174466, 1062.09118733506, 1042.39255611508, 646.485932603773,     830.539929392575, 839.608056347894, 962.253703620869, 959.679777652939,     936.728390874946, 1042.30870964193, 1235.58491946341, 1135.16056375829,     925.636102229735, 1037.36917833249, 1190.05052710661, 1137.8258249324,     1028.27324898977, 1210.24141116166, 1003.70554343626, 1164.36194939349,     1142.27668553833, 1067.86906967352, 1039.76168942087, 1112.7628138853,     755.545931582921, 1269.75201822522, 1042.82660356759, 1030.11703395366,     977.385429900176, 869.276150198159, 1024.36221112439, 807.139680727052,     880.829823012884, 1070.35018777403, 1025.05807928026, 930.372904805134,     1065.02179067573, 1232.88828963314, 971.691543567229, 1220.36077601019,     1075.44208215308, 963.661613451544, 1170.97280330009, 773.416817667865,     1148.31656368132, 1172.16737015, 1088.63567790225, 791.360787896315,     1091.46995904337, 874.171641770601, 1136.750705269, 1176.78271533588,     759.175711977634, 833.015825765764, 974.451049760855, 1120.49398290782,     1230.11818704412, 1006.72358137007, 534.174213658673, 1069.5524682778,     1074.3703653083, 1025.84950615133, 921.826426331041, 1213.3739825443,     1007.63860950274, 1263.08154688765, 976.639081769557, 1003.64998144913,     1063.35173304575, 1195.25184747407, 1226.27460532745, 1150.2126685881,     1075.47157840457, 868.840285334027, 1060.18544407636, 1153.62603990586,     1153.76666784712, 856.927922199376, 997.528561347136, 1107.92272305668,     1137.26823747713, 844.941376479426, 956.236018798489, 1175.67768894008,     1095.70076163552, 1013.4656625271, 1194.38802909718, 1088.99443345778,     936.67022980206, 1176.93543473413, 1177.74723574121, 1021.64293042626,     802.713974099176, 1147.67871582472, 1131.01192302099, 1226.96256838525,     870.460785892088, 1022.22870709406, 896.795271403771, 996.2102724541,     1023.7651674657, 999.400651419049, 1077.23439007785, 609.822058312369,     799.892621773862, 803.519385118485, 992.90620935012, 1113.38678885843,     1148.7760839154, 1065.85921405778, 1084.6477018085, 979.82487122063,     992.516822490875, 758.82854410995, 837.209676696007, 1157.32310266157,     1050.26008916, 948.698388110398, 1101.76672867066, 1043.08842461029,     949.371092782475, 1072.44013370125, 1052.06463567415, 1116.35159775393,     1031.34610678917, 963.010326339346, 1072.16784833961, 990.116487658682,     1001.05172344948, 1234.29864311464, 1195.85293197176, 1024.11653330572,     1011.06024499406, 1221.37676532543, 1013.3624274076, 1123.15561440324,     778.323254312245, 1126.39926949168, 1272.46693250729, 1216.42004797225,     1361.8208304174, 841.747847838327, 1102.88307914987, 977.868009032042,     1122.16975079534, 1130.20864817178, 836.297449148413, 1050.90941494478,     1003.87747610452, 990.027544297921, 1049.97665609377, 1075.56023893662,     1075.64451501875, 1051.07682906703, 919.948881534829, 1183.25610320381,     1106.73404586847, 1120.56235656135, 1003.45069163949, 1097.55688659503,     822.603002712797, 879.414816012222, 884.498409861202, 1098.43925978458,     878.590133463639, 1027.29393858503, 1110.18874377088, 936.256932560715,     1190.07846348807, 916.592855409712, 894.260357637752, 1103.17140094924,     1101.23095357521, 1057.78357579323, 1226.20070515932, 1342.08106369467,     1149.18357250559, 789.595316615116, 525.519608557264, 942.686566618832,     1189.9151119911, 592.829149614415, 1119.45986306101, 912.305143466939,     940.494192466699, 805.428913159954, 1206.77802251838, 1083.72890979688,     1088.26723655449, 985.389052458918, 879.756775065605, 897.554063829254,     1115.16961458823, 1150.23011882017, 1184.01781830254, 999.641936650241,     994.664812118689, 930.548850034589, 1034.88927889182), leiden_clus = c(1,     2, 4, 4, 6, 2, 6, 5, 1, 5, 2, 5, 5, 3, 3, 2, 2, 5, 6, 7,     5, 4, 3, 4, 1, 4, 2, 1, 1, 1, 2, 1, 3, 6, 7, 2, 1, 1, 3,     1, 1, 3, 3, 4, 2, 1, 4, 1, 6, 2, 7, 5, 1, 5, 2, 3, 2, 3,     2, 1, 4, 1, 5, 5, 6, 4, 1, 5, 3, 1, 4, 2, 1, 4, 5, 5, 1,     1, 4, 1, 2, 5, 1, 2, 2, 3, 6, 5, 1, 3, 6, 1, 2, 2, 4, 2,     5, 2, 3, 5, 3, 3, 4, 2, 4, 5, 5, 4, 3, 4, 3, 5, 1, 6, 2,     1, 4, 1, 5, 3, 1, 1, 2, 5, 2, 3, 1, 5, 3, 1, 1, 3, 5, 5,     3, 1, 1, 5, 1, 2, 1, 1, 1, 3, 2, 4, 6, 2, 2, 4, 2, 4, 2,     2, 1, 1, 1, 4, 5, 6, 1, 1, 6, 3, 2, 3, 1, 2, 2, 6, 6, 4,     5, 2, 1, 4, 3, 3, 3, 4, 5, 6, 1, 4, 5, 1, 2, 5, 2, 6, 2,     2, 2, 1, 5, 3, 1, 4, 1, 1, 5, 4, 5, 2, 3, 2, 5, 5, 4, 1,     4, 4, 5, 1, 1, 3, 6, 5, 2, 3, 2, 1, 5, 2, 2, 5, 3, 2, 3,     3, 4, 2, 2, 5, 4, 2, 6, 3, 4, 3, 2, 2, 4, 1, 2, 2, 3, 1,     1, 6, 7, 1, 2, 3, 5, 1, 1, 4, 3, 1, 5, 6, 2, 4, 1, 4, 5,     1, 2, 4, 2, 2, 1, 4, 2, 4, 7, 1, 4, 1, 5, 6, 3, 5, 3, 2,     5, 2, 3, 1, 2, 3, 1, 3, 1, 2, 1, 2, 2, 5, 4, 3, 1, 1, 2,     4, 1, 2, 1, 2, 1, 5, 3, 4, 2, 2, 5, 2, 4, 2, 1, 1, 5, 6,     4, 2, 3, 2, 1, 5, 7, 1, 1, 4, 2, 7, 1, 1, 3, 6, 1, 2, 2,     5, 1, 6, 2, 2, 1, 1, 2, 4, 1, 3, 3, 3, 5, 3, 3, 2, 4, 3,     2, 4, 3, 2, 1, 1, 4, 1, 1, 1, 3, 5, 5, 1, 2, 2, 1, 3, 1,     3, 1, 3, 4, 1, 3, 1, 3, 3, 4, 4, 4, 3, 4, 4, 1, 1, 5, 6,     3, 5, 3, 5, 1, 1, 4, 7, 4, 2, 2, 5, 4, 2, 3, 3, 6, 6, 2,     3, 3, 1, 3, 2, 2, 1, 5, 3, 1, 5, 1, 4, 1, 1, 5, 3, 4, 2,     2, 2, 2, 1, 4, 3, 4, 4, 1, 5, 3, 1, 3, 3, 2, 3, 1, 7, 1,     5, 3, 1, 1, 3, 6, 2, 6, 4, 4, 4, 1, 2, 6, 3, 1, 6, 5, 4,     2, 1, 1, 1, 3, 2, 1, 1, 3, 6, 4, 5, 6, 1, 3, 5, 3, 2, 3,     4, 4, 3, 1, 1, 6, 3, 5, 2, 4, 4, 1, 4, 4, 3, 4, 2, 6, 5,     2, 3, 4, 1, 2, 3, 4, 5, 1, 6, 1, 5, 1, 6, 1, 2, 1, 1, 3,     4, 1, 3, 5, 1, 1, 1, 7, 3, 1, 1, 2, 4, 3, 3, 6, 4, 1, 3,     5, 5, 2, 5, 3, 1, 6, 2, 2, 1, 5, 6, 4, 6, 1, 1, 3, 2, 2,     1, 4, 2, 4, 2, 5, 3, 4, 1, 3, 4, 2, 3, 2, 1, 1, 4, 5, 1,     5, 1, 1, 2, 6, 1, 2, 2, 5, 3, 3, 6, 1, 1, 3, 5, 1, 4, 6,     4, 1, 3, 5, 3, 4, 1, 1, 3, 6, 4, 2, 1, 5, 1), custom_leiden = c(1,     3, 2, 2, 2, 3, 2, 4, 1, 1, 3, 4, 4, 5, 6, 4, 3, 1, 2, 1,     4, 2, 6, 2, 1, 2, 3, 1, 1, 1, 3, 1, 1, 2, 1, 3, 1, 1, 6,     1, 1, 5, 5, 2, 3, 1, 2, 1, 2, 3, 1, 4, 1, 4, 3, 6, 3, 1,     3, 1, 2, 1, 4, 4, 7, 2, 1, 1, 5, 1, 2, 3, 1, 2, 4, 4, 1,     1, 2, 1, 3, 1, 1, 3, 3, 5, 2, 4, 1, 5, 2, 1, 3, 3, 2, 3,     4, 3, 6, 4, 6, 5, 2, 3, 2, 7, 4, 2, 5, 2, 8, 4, 1, 2, 3,     1, 2, 1, 1, 6, 1, 1, 3, 4, 4, 6, 1, 4, 6, 1, 1, 1, 4, 4,     5, 1, 1, 4, 1, 3, 1, 1, 1, 6, 3, 2, 2, 3, 3, 2, 3, 2, 3,     3, 1, 1, 1, 2, 4, 7, 1, 1, 2, 5, 3, 5, 1, 3, 3, 2, 7, 2,     4, 3, 1, 2, 5, 1, 6, 2, 1, 2, 1, 2, 4, 1, 3, 4, 3, 7, 3,     3, 1, 1, 4, 8, 1, 2, 1, 1, 4, 2, 1, 3, 8, 3, 4, 4, 2, 1,     2, 2, 3, 1, 1, 5, 2, 4, 3, 8, 3, 1, 4, 3, 3, 4, 1, 3, 1,     5, 2, 3, 3, 4, 2, 3, 7, 8, 2, 1, 3, 3, 2, 1, 3, 3, 1, 1,     1, 2, 1, 1, 3, 8, 1, 1, 1, 2, 8, 1, 4, 2, 3, 2, 1, 2, 4,     1, 3, 2, 3, 3, 1, 2, 3, 2, 1, 1, 2, 1, 4, 2, 1, 4, 1, 3,     4, 3, 5, 1, 3, 8, 1, 5, 1, 3, 1, 3, 4, 1, 2, 5, 1, 1, 3,     2, 1, 3, 1, 3, 1, 4, 1, 2, 3, 3, 1, 3, 2, 3, 1, 1, 4, 2,     2, 3, 6, 3, 1, 4, 1, 1, 1, 2, 3, 1, 1, 1, 1, 2, 1, 3, 3,     4, 1, 2, 3, 3, 1, 1, 3, 2, 1, 6, 1, 5, 1, 5, 1, 3, 2, 6,     3, 2, 5, 3, 1, 1, 2, 1, 1, 1, 1, 4, 4, 1, 3, 3, 1, 5, 1,     6, 1, 1, 2, 1, 6, 1, 5, 1, 2, 2, 2, 1, 2, 2, 1, 1, 3, 2,     1, 4, 1, 4, 1, 1, 2, 1, 2, 3, 3, 4, 2, 3, 5, 5, 2, 7, 3,     1, 6, 1, 8, 3, 3, 1, 3, 5, 1, 4, 1, 2, 1, 1, 4, 1, 2, 3,     3, 3, 3, 1, 2, 5, 2, 2, 1, 4, 6, 1, 6, 1, 3, 1, 1, 1, 1,     4, 5, 3, 1, 5, 2, 3, 1, 2, 2, 2, 1, 3, 2, 1, 1, 2, 4, 2,     3, 1, 1, 1, 5, 3, 1, 1, 6, 2, 2, 4, 2, 1, 6, 1, 1, 3, 6,     2, 2, 1, 1, 1, 2, 6, 4, 3, 2, 2, 1, 2, 2, 1, 2, 4, 2, 4,     3, 6, 2, 1, 3, 1, 2, 4, 1, 2, 1, 4, 1, 2, 4, 3, 1, 1, 6,     2, 1, 5, 4, 1, 1, 1, 1, 1, 1, 1, 3, 2, 5, 8, 2, 2, 1, 1,     4, 4, 3, 4, 1, 1, 7, 3, 3, 1, 4, 4, 2, 7, 1, 1, 5, 3, 3,     1, 2, 3, 2, 3, 4, 6, 2, 1, 5, 2, 3, 1, 3, 1, 1, 2, 1, 1,     3, 1, 1, 3, 7, 1, 3, 3, 1, 6, 8, 2, 1, 1, 5, 4, 1, 2, 4,     2, 1, 1, 4, 5, 2, 3, 1, 1, 7, 2, 3, 1, 1, 1), sdimx = c(5477L,     5959L, 4720L, 5202L, 4101L, 5821L, 3344L, 3275L, 3895L, 6028L,     5615L, 4858L, 6440L, 5684L, 3964L, 3069L, 5752L, 5752L, 3069L,     3069L, 6234L, 4514L, 3826L, 4376L, 4445L, 3344L, 4445L, 4445L,     5271L, 3619L, 5821L, 6165L, 3275L, 3757L, 5064L, 5615L, 4652L,     3895L, 4170L, 6097L, 6096L, 5546L, 4583L, 4927L, 5684L, 5752L,     3276L, 5546L, 5408L, 6234L, 5477L, 3069L, 5821L, 5546L, 5890L,     4032L, 6303L, 5959L, 6028L, 3895L, 4996L, 5615L, 3344L, 3413L,     6372L, 4583L, 5477L, 6234L, 3138L, 5064L, 4376L, 5615L, 3138L,     4239L, 3619L, 6372L, 4996L, 5684L, 3620L, 4789L, 5684L, 5271L,     4996L, 3757L, 6303L, 3826L, 3826L, 4996L, 3276L, 4307L, 5684L,     5546L, 4032L, 6440L, 3344L, 6303L, 5340L, 4720L, 3482L, 3138L,     3895L, 5752L, 4652L, 5477L, 4652L, 6372L, 3275L, 3207L, 5408L,     4308L, 4514L, 6096L, 5615L, 5409L, 5271L, 3413L, 4996L, 5684L,     5821L, 3757L, 5408L, 5133L, 6165L, 6440L, 5133L, 4239L, 4376L,     5202L, 3620L, 4927L, 4308L, 5684L, 6372L, 6097L, 3207L, 5408L,     5133L, 3207L, 5752L, 5752L, 4239L, 4583L, 4858L, 3895L, 4720L,     4583L, 4170L, 4789L, 4789L, 3826L, 6096L, 3964L, 6372L, 5821L,     4652L, 3826L, 4376L, 3826L, 5752L, 6303L, 5684L, 3069L, 5271L,     5340L, 5408L, 4170L, 5271L, 6165L, 6028L, 3207L, 6165L, 4170L,     4720L, 6096L, 4170L, 4514L, 3344L, 6028L, 4376L, 5271L, 4101L,     3138L, 3757L, 4720L, 3619L, 4308L, 5959L, 4239L, 4170L, 5959L,     5546L, 5340L, 5133L, 5752L, 5684L, 4789L, 4514L, 3757L, 5477L,     6097L, 4514L, 4101L, 4376L, 6372L, 4858L, 6440L, 5477L, 3826L,     4858L, 4170L, 4101L, 3895L, 6303L, 4032L, 4720L, 4239L, 3895L,     3069L, 4995L, 3413L, 5064L, 4789L, 6441L, 3551L, 4514L, 3413L,     4583L, 6165L, 5959L, 3619L, 3482L, 6096L, 6303L, 6234L, 3688L,     6372L, 6096L, 4652L, 4170L, 5821L, 6234L, 6028L, 3620L, 6234L,     4376L, 5821L, 6096L, 3757L, 3482L, 4239L, 3138L, 5133L, 6372L,     4720L, 3688L, 3344L, 4308L, 3757L, 4445L, 4101L, 4651L, 3551L,     3138L, 5271L, 5271L, 5064L, 3344L, 5890L, 6096L, 4927L, 4583L,     5477L, 3964L, 3964L, 5271L, 4583L, 3551L, 3069L, 5133L, 5546L,     5753L, 3826L, 5890L, 5959L, 5752L, 4995L, 3275L, 5546L, 5477L,     5271L, 3275L, 5271L, 5890L, 4652L, 6028L, 4651L, 5477L, 6165L,     5959L, 4583L, 4858L, 3413L, 5064L, 3551L, 5546L, 4308L, 4514L,     3069L, 3757L, 5340L, 4720L, 3138L, 5821L, 3964L, 5408L, 6234L,     3344L, 5890L, 4789L, 4583L, 4652L, 5546L, 3963L, 5340L, 3620L,     5408L, 4101L, 4032L, 5959L, 5890L, 3276L, 4927L, 6165L, 4445L,     3344L, 3413L, 4858L, 4170L, 5064L, 4239L, 6303L, 3963L, 5202L,     3482L, 3344L, 3069L, 3482L, 5890L, 4170L, 4032L, 5202L, 5064L,     5340L, 4445L, 5890L, 4720L, 3207L, 3275L, 5821L, 6440L, 4032L,     4101L, 4170L, 3757L, 4101L, 3826L, 3826L, 3963L, 4308L, 4101L,     3275L, 4651L, 5408L, 5546L, 6028L, 4514L, 4307L, 5684L, 4445L,     3757L, 5202L, 3964L, 5202L, 5546L, 3344L, 3688L, 4032L, 4858L,     3482L, 6028L, 4376L, 3413L, 4376L, 5959L, 3482L, 4032L, 3826L,     6165L, 3207L, 5340L, 5615L, 6372L, 5202L, 3138L, 4927L, 5959L,     4720L, 3138L, 4445L, 5890L, 4239L, 6165L, 4858L, 3688L, 3688L,     4376L, 3688L, 6234L, 4101L, 5546L, 3688L, 5133L, 4996L, 4445L,     3895L, 3482L, 3895L, 5202L, 3551L, 3207L, 3344L, 3413L, 5408L,     5477L, 6234L, 4789L, 3688L, 4307L, 3757L, 5477L, 6096L, 5064L,     4789L, 3895L, 3895L, 3551L, 4858L, 4720L, 4514L, 4445L, 4583L,     5202L, 3413L, 3207L, 3619L, 3207L, 3551L, 3069L, 3963L, 6096L,     5615L, 4445L, 3069L, 5684L, 5064L, 4858L, 3895L, 3620L, 3482L,     4858L, 4996L, 5133L, 4927L, 3482L, 5821L, 4789L, 5133L, 4239L,     3895L, 3688L, 4789L, 4927L, 4583L, 5408L, 4927L, 3276L, 5202L,     6303L, 5615L, 5202L, 4376L, 5340L, 5684L, 3619L, 4720L, 4514L,     4583L, 4308L, 6234L, 6165L, 4101L, 3551L, 3207L, 6028L, 3688L,     5133L, 5340L, 5064L, 4514L, 5959L, 3551L, 6372L, 3069L, 5615L,     6440L, 4858L, 3413L, 3413L, 5959L, 5477L, 4308L, 3344L, 5408L,     3207L, 5546L, 6303L, 4032L, 3138L, 6303L, 5064L, 3964L, 5821L,     4170L, 4927L, 3963L, 5615L, 3551L, 4789L, 3757L, 4307L, 3207L,     5752L, 5890L, 5684L, 6234L, 4445L, 4927L, 4927L, 4652L, 4032L,     3551L, 5340L, 6303L, 5821L, 6440L, 6028L, 5821L, 4032L, 6028L,     4927L, 5752L, 4032L, 4376L, 5890L, 4032L, 6441L, 3826L, 4720L,     3069L, 6234L, 6165L, 4996L, 3482L, 5133L, 3276L, 4858L, 6234L,     3688L, 4170L, 5064L, 4996L, 4239L, 6440L, 5890L, 5959L, 4789L,     3688L, 3138L, 3138L, 5615L, 3482L, 5340L, 3620L, 5271L, 6372L,     4514L, 6372L, 6372L, 4583L, 4239L, 4308L, 3207L, 3413L, 4239L,     3551L, 5408L, 3620L, 4652L, 6441L, 4445L, 3482L, 5271L, 6165L,     5133L, 5133L, 4996L, 3757L, 6028L, 6441L, 4996L, 6303L, 5202L,     5340L, 5615L), sdimy = c(-4125L, -2808L, -5202L, -5322L,     -4604L, -3047L, -4724L, -2928L, -3047L, -4125L, -2688L, -3047L,     -4125L, -4005L, -4125L, -2569L, -3167L, -3646L, -4724L, -4484L,     -5442L, -4604L, -4125L, -4843L, -3526L, -5202L, -2808L, -4005L,     -4245L, -3526L, -2568L, -5082L, -3885L, -4723L, -4125L, -2928L,     -4125L, -3526L, -4245L, -4963L, -4484L, -3766L, -3287L, -4843L,     -2808L, -4604L, -5083L, -4244L, -5202L, -3287L, -5322L, -3047L,     -5202L, -3287L, -2928L, -4245L, -2928L, -4244L, -3167L, -4484L,     -5202L, -4843L, -2808L, -2928L, -4484L, -4723L, -4604L, -4484L,     -3886L, -4364L, -4604L, -3167L, -4125L, -5083L, -3047L, -4723L,     -3766L, -4484L, -4963L, -3885L, -3047L, -3287L, -4245L, -2808L,     -3167L, -3167L, -4604L, -3047L, -4364L, -3287L, -5442L, -5202L,     -2808L, -2688L, -5442L, -3885L, -3167L, -2808L, -4245L, -3167L,     -4245L, -4125L, -4843L, -2688L, -5322L, -4244L, -3407L, -5442L,     -3766L, -5202L, -4364L, -3766L, -4604L, -5442L, -2568L, -3407L,     -4963L, -5202L, -3766L, -4245L, -4005L, -4245L, -3167L, -4364L,     -3047L, -4364L, -3885L, -3167L, -4245L, -3167L, -4484L, -4244L,     -4963L, -5442L, -3766L, -4484L, -4005L, -3047L, -3885L, -2928L,     -3407L, -3526L, -3766L, -4005L, -2568L, -4963L, -5442L, -2688L,     -2928L, -4843L, -3526L, -4843L, -2568L, -2808L, -3646L, -3407L,     -3407L, -5083L, -3406L, -4364L, -4723L, -4245L, -4963L, -3646L,     -2568L, -3287L, -4484L, -3646L, -3406L, -5202L, -4125L, -4963L,     -5442L, -2568L, -3526L, -5322L, -3766L, -4364L, -4364L, -5202L,     -2928L, -5083L, -3047L, -4723L, -2808L, -3526L, -2568L, -2928L,     -2569L, -3766L, -3047L, -2688L, -3287L, -4364L, -3526L, -4364L,     -3885L, -4963L, -3646L, -5202L, -2928L, -5322L, -3646L, -3287L,     -4484L, -2928L, -3406L, -2928L, -5442L, -3766L, -5083L, -5202L,     -3646L, -3526L, -4005L, -3167L, -4963L, -3287L, -2808L, -4125L,     -2928L, -3167L, -5082L, -2688L, -2688L, -3167L, -4484L, -2688L,     -4723L, -3287L, -5202L, -3047L, -3406L, -4005L, -5322L, -3047L,     -4005L, -4364L, -4723L, -4723L, -3526L, -2928L, -5202L, -4723L,     -2688L, -3287L, -4723L, -4484L, -4005L, -4604L, -4364L, -4484L,     -3766L, -4484L, -2928L, -3526L, -3766L, -5202L, -4245L, -3407L,     -2928L, -4843L, -2688L, -5442L, -4005L, -5083L, -3287L, -3885L,     -2808L, -5083L, -2569L, -3167L, -4604L, -5322L, -2808L, -5442L,     -4604L, -3766L, -5442L, -4484L, -5322L, -5322L, -4364L, -5442L,     -4843L, -2568L, -3167L, -2808L, -3885L, -3526L, -2688L, -4723L,     -4125L, -3407L, -5082L, -2688L, -5083L, -2928L, -3526L, -3047L,     -4723L, -3646L, -3885L, -3885L, -2568L, -4963L, -3167L, -2808L,     -3526L, -2928L, -3766L, -3407L, -4005L, -5083L, -2808L, -2808L,     -4245L, -3406L, -4843L, -2808L, -3885L, -3526L, -2928L, -5322L,     -4724L, -3047L, -4364L, -2569L, -4005L, -5322L, -4604L, -3646L,     -4843L, -4963L, -2569L, -4604L, -4005L, -4484L, -4604L, -5322L,     -4843L, -2688L, -2928L, -3047L, -4484L, -5442L, -2569L, -2688L,     -3047L, -3047L, -2688L, -4843L, -4604L, -4484L, -4604L, -3287L,     -3526L, -3646L, -4244L, -3167L, -5442L, -4125L, -2808L, -5442L,     -3167L, -2688L, -3885L, -3646L, -5442L, -3646L, -4125L, -3167L,     -4963L, -5442L, -5322L, -3407L, -2808L, -3287L, -3766L, -3287L,     -4125L, -4364L, -4364L, -4963L, -4963L, -3885L, -4005L, -3526L,     -3526L, -4843L, -5322L, -5322L, -5083L, -4963L, -4963L, -5202L,     -3646L, -4364L, -2569L, -5083L, -5083L, -5202L, -4604L, -2928L,     -4125L, -5202L, -4963L, -4604L, -5442L, -3167L, -2688L, -5322L,     -5202L, -2688L, -3167L, -3167L, -4604L, -4244L, -2688L, -4005L,     -4364L, -3766L, -4484L, -2569L, -2569L, -3287L, -2808L, -3646L,     -3646L, -2808L, -4005L, -4843L, -4723L, -4364L, -5202L, -4604L,     -5083L, -2569L, -2569L, -2928L, -3287L, -3646L, -5322L, -3287L,     -4723L, -5083L, -3287L, -3047L, -4125L, -3047L, -4245L, -4843L,     -2688L, -4005L, -3766L, -4484L, -4125L, -3526L, -3167L, -4244L,     -3646L, -3287L, -4963L, -2568L, -3167L, -4963L, -5442L, -5442L,     -3766L, -2808L, -4723L, -4723L, -3885L, -4724L, -5442L, -5083L,     -2808L, -3646L, -3766L, -3407L, -3407L, -2928L, -4005L, -4244L,     -4364L, -4843L, -5083L, -5322L, -5322L, -3406L, -4125L, -4125L,     -4963L, -2569L, -4245L, -5083L, -5202L, -4005L, -4963L, -4604L,     -4843L, -4364L, -3287L, -2688L, -4843L, -4963L, -4364L, -5322L,     -4843L, -4484L, -5322L, -4005L, -5202L, -3406L, -3885L, -4245L,     -5083L, -4364L, -3047L, -4843L, -4723L, -3047L, -3526L, -4963L,     -4723L, -4604L, -3287L, -4843L, -5082L, -2688L, -3885L, -4484L,     -4005L, -5322L, -3407L, -3885L, -2928L, -3646L, -3766L, -3047L,     -4245L, -5082L, -5082L, -3766L, -2568L, -5202L, -3407L, -4604L,     -4604L, -4963L, -3167L, -4843L, -4125L, -3526L, -3646L, -3646L,     -4963L, -4484L, -3885L, -2688L, -2688L, -3766L, -2928L, -3646L,     -4723L, -4604L, -4364L, -3526L, -4005L, -3047L, -3406L, -4005L,     -5442L, -2568L, -5322L, -2568L, -3766L, -4125L, -5202L, -3406L,     -3526L, -4843L, -3406L, -4843L, -3287L, -4125L, -3646L, -5322L,     -3646L, -4364L, -2808L, -3885L, -4484L, -3047L, -5442L, -3646L,     -2808L, -3526L, -3766L, -4125L, -4245L, -4724L, -3885L, -3885L,     -3407L, -3287L, -4005L, -5083L, -5322L, -4723L, -4484L, -3766L,     -3885L, -3526L, -5202L, -3287L, -4005L, -4604L, -4843L, -5442L,     -2688L, -3885L, -3406L, -4125L)), class = c("data.table",     "data.frame"), row.names = c(NA, -624L), .internal.selfref = <pointer: 0x55dfcc1245c0>, sorted = "cell_ID"),     cell_locations_metadata_other = NULL, cell_color = "HMRF_k6_b.20",     color_as_factor = TRUE, cell_color_code = NULL, cell_color_gradient = NULL,     gradient_midpoint = NULL, gradient_style = c("divergent",     "sequential"), gradient_limits = NULL, select_cell_groups = NULL,     select_cells = NULL, point_size = 3, point_alpha = 1, show_cluster_center = FALSE,     show_center_label = FALSE, center_point_size = 4, label_size = 4,     label_fontface = "bold", show_other_cells = TRUE, other_cell_color = "lightgrey",     other_point_size = 1, show_legend = TRUE, point_border_stroke = 0.1,     point_border_col = "black", center_point_border_col = "black",     center_point_border_stroke = 0.1): HMRF_k6_b.20 is not a color or a column name