sample_sheet <- glue::glue("sample_sheets/tmrc2_samples_20220106.xlsx")

1 Introduction

This is mostly just a run of this worksheet to reacquaint myself with it.

This document is intended to provide a general overview of the TMRC2 samples which have thus far been sequenced. In some cases, this includes only those samples starting in 2019; in other instances I am including our previous (2015-2016) samples.

In all cases the processing performed was:

  1. Default trimming was performed.
  2. Hisat2 was used to map the remaining reads against the Leishmania panamensis genome revision 36.
  3. The alignments from hisat2 were used to count reads/gene against the revision 36 annotations with htseq.
  4. These alignments were also passed to the pileup functionality of samtools and the vcf/bcf utilities in order to make a matrix of all observed differences between each sample with respect to the reference.

The analyses in this document use the matrices of counts/gene from #3 and variants/position from #4 in order to provide some images and metrics describing the samples we have sequenced so far.

2 Annotations

Everything which follows depends on the Existing TriTrypDB annotations revision 46, circa 2019. The following block loads a database of these annotations and turns it into a matrix where the rows are genes and columns are all the annotation types provided by TriTrypDB.

The same database was used to create a matrix of orthologous genes between L.panamensis and all of the other species in the TriTrypDB.

tt <- sm(library(EuPathDB))
tt <- sm(library(org.Lpanamensis.MHOMCOL81L13.v46.eg.db))
pan_db <- org.Lpanamensis.MHOMCOL81L13.v46.eg.db
all_fields <- columns(pan_db)

all_lp_annot <- sm(load_orgdb_annotations(
    pan_db,
    keytype = "gid",
    fields = c("annot_gene_entrez_id", "annot_gene_name",
               "annot_strand", "annot_chromosome", "annot_cds_length",
               "annot_gene_product")))$genes

lp_go <- sm(load_orgdb_go(pan_db))
lp_lengths <- all_lp_annot[, c("gid", "annot_cds_length")]
colnames(lp_lengths)  <- c("ID", "length")
all_lp_annot[["annot_gene_product"]] <- tolower(all_lp_annot[["annot_gene_product"]])
orthos <- sm(EuPathDB::extract_eupath_orthologs(db = pan_db))

hisat_annot <- all_lp_annot
## rownames(hisat_annot) <- paste0("exon_", rownames(hisat_annot), ".E1")

3 Load a genome

meta <- EuPathDB::download_eupath_metadata(webservice="tritrypdb")
## Unable to find species names for 2 species.
## Leishmania sp. Ghana 2012 LV757 MHOMGH2012GH5, Leishmania sp. Namibia MPRO/NA/1975/252/LV425
## Appending to an existing file: EuPathDB/metadata/biocv3.14_tritrypdbv55_metadata.csv
## Appending to an existing file: EuPathDB/metadata/GRanges_biocv3.14_tritrypdbv55_metadata.csv
## Appending to an existing file: EuPathDB/metadata/OrgDb_biocv3.14_tritrypdbv55_metadata.csv
## Appending to an existing file: EuPathDB/metadata/TxDb_biocv3.14_tritrypdbv55_metadata.csv
## Appending to an existing file: EuPathDB/metadata/OrganismDbi_biocv3.14_tritrypdbv55_metadata.csv
## Appending to an existing file: EuPathDB/metadata/BSgenome_biocv3.14_tritrypdbv55_metadata.csv
## Appending to an existing file: EuPathDB/metadata/biocv3.14_tritrypdbv55_invalid_metadata.csv
## Appending to an existing file: EuPathDB/metadata/GRanges_biocv3.14_tritrypdbv55_invalid_metadata.csv
## Appending to an existing file: EuPathDB/metadata/OrgDb_biocv3.14_tritrypdbv55_invalid_metadata.csv
## Appending to an existing file: EuPathDB/metadata/TxDb_biocv3.14_tritrypdbv55_invalid_metadata.csv
## Appending to an existing file: EuPathDB/metadata/OrganismDbi_biocv3.14_tritrypdbv55_invalid_metadata.csv
## Appending to an existing file: EuPathDB/metadata/BSgenome_biocv3.14_tritrypdbv55_invalid_metadata.csv
lp_entry <- EuPathDB::get_eupath_entry(species="Leishmania panamensis", metadata=meta)
## Found the following hits: Leishmania panamensis MHOM/COL/81/L13, Leishmania panamensis strain MHOM/PA/94/PSC-1, choosing the first.
## Using: Leishmania panamensis MHOM/COL/81/L13.
colnames(lp_entry)
##  [1] "AnnotationVersion"  "AnnotationSource"   "BiocVersion"       
##  [4] "DataProvider"       "Genome"             "GenomeSource"      
##  [7] "GenomeVersion"      "NumArrayGene"       "NumChipChipGene"   
## [10] "NumChromosome"      "NumCodingGene"      "NumCommunity"      
## [13] "NumContig"          "NumEC"              "NumEST"            
## [16] "NumGene"            "NumGO"              "NumOrtholog"       
## [19] "NumOtherGene"       "NumPopSet"          "NumProteomics"     
## [22] "NumPseudogene"      "NumRNASeq"          "NumRTPCR"          
## [25] "NumSNP"             "NumTFBS"            "Organellar"        
## [28] "ReferenceStrain"    "MegaBP"             "PrimaryKey"        
## [31] "ProjectID"          "RecordClassName"    "SourceID"          
## [34] "SourceVersion"      "TaxonomyID"         "TaxonomyName"      
## [37] "URLGenome"          "URLGFF"             "URLProtein"        
## [40] "Coordinate_1_based" "Maintainer"         "SourceUrl"         
## [43] "Tags"               "BsgenomePkg"        "GrangesPkg"        
## [46] "OrganismdbiPkg"     "OrgdbPkg"           "TxdbPkg"           
## [49] "Taxon"              "Genus"              "Species"           
## [52] "Strain"             "BsgenomeFile"       "GrangesFile"       
## [55] "OrganismdbiFile"    "OrgdbFile"          "TxdbFile"          
## [58] "GenusSpecies"       "TaxonUnmodified"    "TaxonCanonical"    
## [61] "TaxonXref"
testing_panamensis <- "BSGenome.Leishmania.panamensis.MHOMCOL81L13.v53"
## testing_panamensis <- EuPathDB::make_eupath_bsgenome(entry=lp_entry, eu_version="v46")
library(as.character(testing_panamensis), character.only=TRUE)
## Loading required package: BSgenome
## Loading required package: Biostrings
## Loading required package: XVector
## 
## Attaching package: 'Biostrings'
## The following object is masked from 'package:base':
## 
##     strsplit
## Loading required package: rtracklayer
genome <- get0(as.character(testing_panamensis))

4 TODO:

Resequence samples: TMRC20002, TMRC20006, TMRC20004 (maybe TMRC20008 and TMRC20029)

5 Generate Expressionsets and Sample Estimation

The process of sample estimation takes two primary inputs:

  1. The sample sheet, which contains all the metadata we currently have on hand, including filenames for the outputs of #3 and #4 above.
  2. The gene annotations.

An expressionset is a data structure used in R to examine RNASeq data. It is comprised of annotations, metadata, and expression data. In the case of our processing pipeline, the location of the expression data is provided by the filenames in the metadata.

The first lines of the following block create the Expressionset. All of the following lines perform various normalizations and generate plots from it.

5.1 Notes

The following samples are much lower coverage:

  • TMRC20002
  • TMRC20006
  • TMRC20007
  • TMRC20008

20210610: I made some manual changes to the sample sheet which I downloaded, filling in some zymodeme with ‘unknown’

5.2 TODO:

  1. Do the multi-gene family removal right here instead of way down at the bottom
  2. Add zymodeme snps to the annotation later.
  3. Start phylogenetic analysis of variant table.
sanitize_columns <- c("passagenumber", "clinicalresponse", "clinicalcategorical",
                      "zymodemecategorical", "zymodemecategorical")
lp_expt <- sm(create_expt(sample_sheet,
                          gene_info = hisat_annot,
                          id_column = "hpglidentifier",
                          file_column = "lpanamensisv36hisatfile")) %>%
  set_expt_conditions(fact = "zymodemecategorical") %>%
  subset_expt(nonzero = 8550) %>%
  subset_expt(coverage = 5000000) %>%
  semantic_expt_filter(semantic = c("amastin", "gp63", "leishmanolysin"),
                       semantic_column = "annot_gene_product") %>%
  sanitize_expt_metadata(columns = sanitize_columns) %>%
  set_expt_factors(columns = sanitize_columns, class = "factor")
## The samples (and read coverage) removed when filtering 8550 non-zero genes are:
## TMRC20002 TMRC20006 
##  11681227   6670348
## subset_expt(): There were 75, now there are 73 samples.
## The samples removed (and read coverage) when filtering samples with less than 5e+06 reads are:
## TMRC20004 TMRC20029 
##    564812   1658096
## subset_expt(): There were 73, now there are 71 samples.
## semantic_expt_filter(): Removed 68 genes.
libsizes <- plot_libsize(lp_expt)
pp(file = "images/lp_expt_libsizes.png", image = libsizes$plot, width = 14, height = 9)

## I think samples 7,10 should be removed at minimum, probably also 9,11
nonzero <- plot_nonzero(lp_expt)
pp(file = "images/lp_nonzero.png", image = nonzero$plot, width = 9, height = 9)
## Warning: ggrepel: 50 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
## Warning: ggrepel: 53 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps

lp_box <- plot_boxplot(lp_expt)
## 5364 entries are 0.  We are on a log scale, adding 1 to the data.
pp(file = "images/lp_expt_boxplot.png", image = lp_box, width = 12, height = 9)

filter_plot <- plot_libsize_prepost(lp_expt)
filter_plot$lowgene_plot
## Warning: Using alpha for a discrete variable is not advised.

filter_plot$count_plot

5.3 Distribution Visualization

Najib’s favorite plots are of course the PCA/TNSE. These are nice to look at in order to get a sense of the relationships between samples. They also provide a good opportunity to see what happens when one applies different normalizations, surrogate analyses, filters, etc. In addition, one may set different experimental factors as the primary ‘condition’ (usually the color of plots) and surrogate ‘batches’.

5.4 By Susceptilibity

Column ‘Q’ in the sample sheet, make a categorical version of it with these parameters:

  • 0 <= x <= 35 is resistant
  • 36 <= x <= 48 is ambiguous
  • 49 <= x is sensitive
starting <- as.numeric(pData(lp_expt)[["susceptibilityinfectionreduction32ugmlsbvhistoricaldata"]])
## Warning: NAs introduced by coercion
sus_categorical <- starting
na_idx <- is.na(starting)
sus_categorical[na_idx] <- "unknown"

resist_idx <- starting <= 0.35
sus_categorical[resist_idx] <- "resistant"
indeterminant_idx <- starting >= 0.36 & starting <= 0.48
sus_categorical[indeterminant_idx] <- "ambiguous"
susceptible_idx <- starting >= 0.49
sus_categorical[susceptible_idx] <- "sensitive"

pData(lp_expt$expressionset)[["sus_category"]] <- sus_categorical
clinical_colors <- list(
    ##    "z2.1" = "#0000cc",
    ##    "z2.3" = "#874400",
    ##    "z2.2" = "#df7000",
    ##    "z2.4" = "#cc0000",
    "z2.1" = "#874400",
    "z2.2" = "#0000cc",
    "z2.3" = "#cc0000",
    "z2.4" = "#df7000",
    "unknown" = "#cbcbcb",
    "null" = "#000000")
clinical_samples <- lp_expt %>%
  set_expt_batches(fact = sus_categorical) %>%
  set_expt_colors(clinical_colors)

clinical_norm <- sm(normalize_expt(clinical_samples, norm = "quant", transform = "log2",
                                   convert = "cpm", batch = FALSE, filter = TRUE))
zymo_pca <- plot_pca(clinical_norm, plot_title = "PCA of parasite expression values",
                     plot_labels = FALSE)
pp(file = "images/zymo_pca_sus_shape.png", image = zymo_pca$plot)

only_two_types <- subset_expt(clinical_samples, subset = "condition=='z2.3'|condition=='z2.2'")
## subset_expt(): There were 71, now there are 59 samples.
only_two_norm <- sm(normalize_expt(only_two_types, norm = "quant", transform = "log2",
                                   convert = "cpm", batch = FALSE, filter = TRUE))
onlytwo_pca <- plot_pca(only_two_norm, plot_title = "PCA of z2.2 and z2.3 parasite expression values",
                     plot_labels = FALSE)
pp(file = "images/zymo_z2.2_z2.3_pca_sus_shape.png", image = onlytwo_pca$plot)

zymo_3dpca <- plot_3d_pca(zymo_pca)
zymo_3dpca$plot
clinical_n <- sm(normalize_expt(clinical_samples, transform = "log2",
                                convert = "cpm", batch = FALSE, filter = TRUE))
zymo_tsne <- plot_tsne(clinical_n, plot_title = "TSNE of parasite expression values")
zymo_tsne$plot

clinical_nb <- normalize_expt(clinical_samples, convert = "cpm", transform = "log2",
                         filter = TRUE, batch = "svaseq")
## Removing 142 low-count genes (8568 remaining).
## batch_counts: Before batch/surrogate estimation, 1008 entries are x==0: 0%.
## batch_counts: Before batch/surrogate estimation, 3380 entries are 0<x<1: 1%.
## Setting 370 low elements to zero.
## transform_counts: Found 370 values equal to 0, adding 1 to the matrix.
clinical_nb_pca <- plot_pca(clinical_nb, plot_title = "PCA of parasite expression values",
                            plot_labels = FALSE)
pp(file = "images/clinical_nb_pca_sus_shape.png", image = clinical_nb_pca$plot)

clinical_nb_tsne <- plot_tsne(clinical_nb, plot_title = "TSNE of parasite expression values")
clinical_nb_tsne$plot
## Warning: ggrepel: 68 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps

corheat <- plot_corheat(clinical_norm, plot_title = "Correlation heatmap of parasite
                 expression values
")
corheat$plot

plot_sm(clinical_norm)$plot
## Performing correlation.

5.5 By Cure/Fail status

cf_colors <- list(
    "cure" = "#006f00",
    "fail" = "#9dffa0",
    "unknown" = "#cbcbcb",
    "notapplicable" = "#000000")
cf_expt <- set_expt_conditions(lp_expt, fact = "clinicalcategorical") %>%
  set_expt_batches(fact = sus_categorical) %>%
  set_expt_colors(cf_colors)
## Warning in set_expt_colors(., cf_colors): Colors for the following categories
## are not being used: notapplicable.
cf_norm <- normalize_expt(cf_expt, convert = "cpm", transform = "log2",
                          norm = "quant", filter = TRUE)
## Removing 142 low-count genes (8568 remaining).
## transform_counts: Found 2 values equal to 0, adding 1 to the matrix.
start_cf <- plot_pca(cf_norm, plot_title = "PCA of parasite expression values",
                     plot_labels = FALSE)
pp(file = "images/cf_sus_shape.png", image = start_cf$plot)

cf_nb <- normalize_expt(cf_expt, convert = "cpm", transform = "log2",
                        norm = "quant", filter = TRUE, batch = "svaseq")
## Warning in normalize_expt(cf_expt, convert = "cpm", transform = "log2", :
## Quantile normalization and sva do not always play well together.
## Removing 142 low-count genes (8568 remaining).
## batch_counts: Before batch/surrogate estimation, 2 entries are x==0: 0%.
## batch_counts: Before batch/surrogate estimation, 4130 entries are 0<x<1: 1%.
## Setting 154 low elements to zero.
## transform_counts: Found 154 values equal to 0, adding 1 to the matrix.
cf_nb_pca <- plot_pca(cf_nb, plot_title = "PCA of parasite expression values",
                      plot_labels = FALSE)
pp(file = "images/cf_sus_share_nb.png", image = cf_nb_pca$plot)

cf_norm <- normalize_expt(cf_expt, transform = "log2", convert = "cpm",
                          filter = TRUE, norm = "quant")
## Removing 142 low-count genes (8568 remaining).
## transform_counts: Found 2 values equal to 0, adding 1 to the matrix.
test <- pca_information(cf_norm,
                        expt_factors = c("clinicalcategorical", "zymodemecategorical",
                                         "pathogenstrain", "passagenumber"),
                        num_components = 6, plot_pcas = TRUE)
test$anova_p
##                           PC1       PC2     PC3       PC4       PC5    PC6
## clinicalcategorical 2.850e-01 0.3311343 0.47800 1.728e-03 0.6934731 0.4209
## zymodemecategorical 4.608e-08 0.0009469 0.33206 1.487e-02 0.0009543 0.2642
## pathogenstrain      7.092e-01 0.7763033 0.84356 4.512e-06 0.0181051 0.5619
## passagenumber       8.896e-01 0.2377294 0.04096 2.795e-02 0.2229372 0.4130
test$cor_heatmap

sus_colors <- list(
    "resistant" = "#8563a7",
    "sensitive" = "#8d0000",
    "ambiguous" = "#cbcbcb",
    "unknown" = "#000000")
sus_expt <- set_expt_conditions(lp_expt, fact = "sus_category") %>%
  set_expt_batches(fact = "zymodemecategorical") %>%
  set_expt_colors(colors = sus_colors) %>%
  subset_expt(subset = "batch!='z24'") %>%
  subset_expt(subset = "batch!='z21'")
## subset_expt(): There were 71, now there are 70 samples.
## subset_expt(): There were 70, now there are 67 samples.
sus_norm <- normalize_expt(sus_expt, transform = "log2", convert = "cpm",
                           norm = "quant", filter = TRUE)
## Removing 143 low-count genes (8567 remaining).
## transform_counts: Found 2 values equal to 0, adding 1 to the matrix.
sus_pca <- plot_pca(sus_norm, plot_title = "PCA of parasite expression values",
                    plot_labels = FALSE)
pp(file = "images/sus_norm_pca.png", image = sus_pca[["plot"]])

sus_nb <- normalize_expt(sus_expt, transform = "log2", convert = "cpm",
                         batch = "svaseq", filter = TRUE)
## Removing 143 low-count genes (8567 remaining).
## batch_counts: Before batch/surrogate estimation, 939 entries are x==0: 0%.
## batch_counts: Before batch/surrogate estimation, 3146 entries are 0<x<1: 1%.
## Setting 208 low elements to zero.
## transform_counts: Found 208 values equal to 0, adding 1 to the matrix.
sus_nb_pca <- plot_pca(sus_nb, plot_title = "PCA of parasite expression values",
                       plot_labels = FALSE)
pp(file = "images/sus_nb_pca.png", image = sus_nb_pca[["plot"]])

6 Zymodeme analyses

The following sections perform a series of analyses which seek to elucidate differences between the zymodemes 2.2 and 2.3 either through differential expression or variant profiles.

6.1 Differential expression

6.1.1 With respect to zymodeme attribution

TODO: Do this with and without sva and compare the results.

zy_expt <- subset_expt(lp_expt, subset = "condition=='z2.2'|condition=='z2.3'")
## subset_expt(): There were 71, now there are 59 samples.
zy_norm <- normalize_expt(zy_expt, filter = TRUE, convert = "cpm", norm = "quant")
## Removing 159 low-count genes (8551 remaining).
zy_de_nobatch <- sm(all_pairwise(zy_expt, filter = TRUE, model_batch = "svaseq"))
zy_de <- sm(all_pairwise(zy_expt, filter = TRUE, model_batch = "svaseq"))
zy_table <- sm(combine_de_tables(zy_de, excel = glue::glue("excel/zy_tables-v{ver}.xlsx")))
zy_sig <- sm(extract_significant_genes(zy_table, excel = glue::glue("excel/zy_sig-v{ver}.xlsx")))

6.1.2 Images of zymodeme DE

pp(file = "images/zymo_ma.png", image = zy_table[["plots"]][["z23_vs_z22"]][["deseq_ma_plots"]][["plot"]])

6.2 With respect to cure/failure

In contrast, we can search for genes which are differentially expressed with respect to cure/failure status.

cf_de <- sm(all_pairwise(cf_expt, filter = TRUE, model_batch = "svaseq"))
cf_table <- sm(combine_de_tables(cf_de, excel = glue::glue("excel/cf_tables-v{ver}.xlsx")))
cf_sig <- sm(extract_significant_genes(cf_table, excel = glue::glue("excel/cf_sig-v{ver}.xlsx")))

6.3 With respect to susceptibility

Finally, we can use our category of susceptibility and look for genes which change from sensitive to resistant. Keep in mind, though, that for the moment we have a lot of ambiguous and unknown strains.

sus_de <- sm(all_pairwise(sus_expt, filter = TRUE, model_batch = "svaseq"))
sus_table <- sm(combine_de_tables(sus_de, excel = glue::glue("excel/sus_tables-v{ver}.xlsx")))
sus_sig <- sm(extract_significant_genes(sus_table, excel = glue::glue("excel/sus_sig-v{ver}.xlsx")))
knitr::kable(head(sus_sig$deseq$ups$sensitive_vs_resistant, n = 20))
gid annotgeneproduct annotgenetype chromosome start end strand annotgeneentrezid annotgenename annotstrand annotchromosome annotcdslength length deseq_logfc deseq_adjp edger_logfc edger_adjp limma_logfc limma_adjp basic_nummed basic_denmed basic_numvar basic_denvar basic_logfc basic_t basic_p basic_adjp deseq_basemean deseq_lfcse deseq_stat deseq_p ebseq_fc ebseq_logfc ebseq_c1mean ebseq_c2mean ebseq_mean ebseq_var ebseq_postfc ebseq_ppee ebseq_ppde ebseq_adjp edger_logcpm edger_lr edger_p limma_ave limma_t limma_b limma_p limma_adjp_ihw deseq_adjp_ihw edger_adjp_ihw ebseq_adjp_ihw basic_adjp_ihw lfc_meta lfc_var lfc_varbymed p_meta p_var
LPAL13_000035800 LPAL13_000035800 hypothetical protein protein coding LPAL13_SCAF000500 737 1006 - reverse Not Assigned 270.0 269 15.350 0 14.480 0e+00 10.660 0.2306 4.9910 -3.9520 17.904 0.3803 8.943 10.920 0e+00 0.0000 2911.00 1.1280 13.600 0 30205.04 14.883 0.1400 4530.73 3171.56 1.650e+07 1223.564 0.00 0.00 0.00 6.7310 83.11 0 2.1490 1.635 -4.2650 0.1070 2.807e-01 2.877e-38 6.208e-16 0.000e+00 6.641e-09 13.500 0.000e+00 0.000e+00 3.567e-02 3.816e-03
LPAL13_320026300 LPAL13_320026300 hypothetical protein, conserved protein coding LpaL13_32 754268 755485 - reverse 32 1218.0 1217 12.980 0 13.470 0e+00 9.953 0.2542 4.6270 -4.0140 18.533 0.6270 8.642 10.230 0e+00 0.0000 1527.00 1.0900 11.910 0 17204.67 14.070 0.1308 2423.22 1696.30 2.403e+06 629.468 0.00 0.00 0.00 5.8010 73.69 0 1.9250 1.562 -4.3400 0.1233 3.157e-01 6.162e-29 3.920e-14 0.000e+00 9.705e-09 11.980 1.188e-01 9.915e-03 4.110e-02 5.068e-03
LPAL13_000044900 LPAL13_000044900 actin-related protein 2, putative protein coding LPAL13_SCAF000645 507 1685 - reverse Not Assigned 1179.0 1178 12.190 0 13.610 0e+00 9.092 0.1953 3.7750 -4.1700 17.373 0.1233 7.944 10.000 0e+00 0.0000 849.30 1.1340 10.750 0 128040.28 16.966 0.0000 1280.39 896.27 6.350e+05 343.849 1.00 0.00 1.00 4.9560 69.01 0 1.2680 1.757 -4.1320 0.0837 2.367e-01 2.157e-23 3.415e-13 0.000e+00 3.138e-08 11.780 1.188e-01 1.008e-02 2.790e-02 2.335e-03
LPAL13_000053200 LPAL13_000053200 hypothetical protein protein coding LPAL13_SCAF000804 5037 5249 - reverse Not Assigned 213.0 212 8.925 0 10.330 0e+00 6.031 0.0185 0.8570 -4.1700 9.600 0.1233 5.027 8.459 0e+00 0.0000 77.23 1.0440 8.550 0 13478.05 13.718 0.0000 134.77 94.34 9.270e+03 36.806 1.00 0.00 1.00 1.5470 58.05 0 -0.8949 3.084 -2.2650 0.0030 2.308e-02 1.372e-14 1.926e-11 0.000e+00 3.411e-07 8.406 3.528e-02 4.197e-03 1.003e-03 3.018e-06
LPAL13_000051300 LPAL13_000051300 hypothetical protein, conserved protein coding LPAL13_SCAF000772 11 2344 + forward Not Assigned 2334.0 2333 8.786 0 9.749 0e+00 4.677 0.1557 0.4463 -3.9870 11.397 0.5398 4.433 6.594 0e+00 0.0000 146.30 1.2060 7.288 0 2017.81 10.979 0.0962 214.30 150.04 8.669e+04 62.049 0.00 0.00 0.00 2.4650 36.51 0 -1.0200 1.909 -3.9620 0.0608 1.903e-01 1.304e-10 1.564e-07 0.000e+00 7.692e-06 7.477 1.518e+00 2.031e-01 2.025e-02 1.230e-03
LPAL13_000040700 LPAL13_000040700 hypothetical protein, conserved protein coding LPAL13_SCAF000598 54 1067 + forward Not Assigned 1014.0 1013 6.648 0 7.096 0e+00 3.375 0.0531 -1.0270 -3.9700 6.795 0.5482 2.943 5.481 0e+00 0.0001 21.87 0.9994 6.652 0 227.13 7.827 0.1428 34.70 24.33 9.986e+02 9.762 0.00 0.00 0.00 -0.0905 32.62 0 -2.2010 2.554 -3.1400 0.0130 6.566e-02 5.267e-09 7.944e-07 0.000e+00 8.347e-05 5.792 7.091e-01 1.224e-01 4.337e-03 5.642e-05
LPAL13_000017600 LPAL13_000017600 hypothetical protein, conserved protein coding LPAL13_SCAF000146 359 586 + forward Not Assigned 228.0 227 6.564 0 6.548 0e+00 6.287 0.0429 4.1180 -1.1470 5.268 2.5720 5.266 8.300 0e+00 0.0000 613.10 0.6613 9.927 0 72.70 6.184 13.2097 961.12 676.75 4.516e+05 58.224 0.00 1.00 0.00 4.4880 59.91 0 2.2600 2.665 -2.6780 0.0097 5.346e-02 6.951e-20 9.352e-12 8.934e-01 3.762e-07 6.754 2.432e+00 3.600e-01 3.238e-03 3.146e-05
LPAL13_300029400 LPAL13_300029400 hypothetical protein, conserved protein coding LpaL13_30 853953 854150 - reverse 30 198.0 197 6.390 0 6.323 0e+00 5.128 0.0011 1.5080 -2.5360 1.854 1.9106 4.045 8.519 0e+00 0.0000 84.63 0.7095 9.007 0 59.75 5.901 2.1189 127.20 89.67 1.056e+04 22.395 0.00 0.00 0.00 1.6260 58.02 0 -0.0869 4.290 0.4427 0.0001 1.350e-03 2.648e-16 1.926e-11 0.000e+00 2.140e-06 5.953 3.312e-01 5.564e-02 2.036e-05 1.244e-09
LPAL13_000011700 LPAL13_000011700 hypothetical protein protein coding LPAL13_SCAF000076 101 364 - reverse Not Assigned 264.0 263 6.200 0 7.639 0e+00 3.206 0.0185 -1.4720 -4.1700 6.800 0.1233 2.698 5.362 0e+00 0.0002 14.87 1.1620 5.337 0 2690.50 11.394 0.0000 26.89 18.83 5.831e+02 7.965 1.00 0.00 1.00 -0.5642 26.82 0 -2.8600 3.083 -2.3420 0.0030 2.262e-02 5.009e-06 9.045e-06 0.000e+00 1.999e-04 5.654 9.394e-01 1.661e-01 1.003e-03 3.018e-06
LPAL13_080010600 LPAL13_080010600 hypothetical protein, conserved protein coding LpaL13_08 195555 195749 - reverse 8 195.0 194 6.002 0 7.324 0e+00 2.837 0.0176 -1.9210 -4.1700 4.935 0.1233 2.249 5.207 0e+00 0.0003 11.70 1.0980 5.466 0 2013.65 10.976 0.0000 20.13 14.09 5.608e+02 6.299 1.00 0.00 1.00 -0.9047 27.91 0 -3.0720 3.108 -2.3430 0.0028 2.206e-02 2.256e-06 5.812e-06 0.000e+00 2.196e-04 5.091 1.949e+00 3.828e-01 9.341e-04 2.617e-06
LPAL13_040019400 LPAL13_040019400 hypothetical protein protein coding LpaL13_04 440768 441127 - reverse 4 360.0 359 5.679 0 5.727 0e+00 3.639 0.0059 -0.4572 -3.3380 1.891 1.1727 2.881 7.087 0e+00 0.0000 28.83 0.8258 6.877 0 43.77 5.452 0.8450 37.41 26.44 2.057e+03 8.868 0.00 0.00 0.00 0.1553 36.77 0 -1.8060 3.603 -1.3140 0.0006 7.378e-03 1.373e-09 1.607e-07 0.000e+00 6.646e-06 5.066 6.063e-02 1.197e-02 2.041e-04 1.250e-07
LPAL13_200050100 LPAL13_200050100 hypothetical protein protein coding LpaL13_20.1 1627529 1627717 + forward 20.1 189.0 188 5.525 0 5.473 0e+00 4.659 0.0009 2.3090 -1.8430 1.742 2.3634 4.152 8.155 0e+00 0.0000 121.50 0.5792 9.539 0 24.99 4.643 8.3178 208.07 148.14 2.549e+04 17.788 0.00 1.00 0.00 2.1720 62.54 0 0.7885 4.367 0.9945 0.0000 1.122e-03 2.478e-18 4.481e-12 9.791e-01 8.497e-06 5.208 3.918e-01 7.523e-02 1.551e-05 7.217e-10
LPAL13_080010800 LPAL13_080010800 hypothetical protein protein coding LpaL13_08 199409 199792 - reverse 8 384.0 383 5.448 0 6.865 0e+00 2.069 0.1208 -2.1170 -4.1700 4.788 0.1233 2.053 4.821 0e+00 0.0006 11.10 1.0070 5.407 0 1600.96 10.645 0.0000 16.00 11.20 4.009e+02 5.222 1.00 0.00 1.00 -0.8462 30.73 0 -3.0170 2.073 -3.7940 0.0422 1.489e-01 3.581e-06 2.225e-06 0.000e+00 4.671e-04 4.368 3.348e+00 7.664e-01 1.405e-02 5.922e-04
LPAL13_350011800 LPAL13_350011800 hypothetical protein, conserved protein coding LpaL13_35 171009 171242 + forward 35 234.0 233 5.289 0 5.278 0e+00 4.720 0.0033 2.6120 -1.1520 2.957 1.1319 3.764 8.418 0e+00 0.0000 178.20 0.5669 9.329 0 31.86 4.994 9.5389 304.18 215.79 6.840e+04 23.769 0.00 1.00 0.00 2.6990 59.94 0 1.0920 3.836 -0.0887 0.0003 4.173e-03 1.649e-17 9.352e-12 9.791e-01 1.578e-07 5.163 8.144e-01 1.577e-01 9.547e-05 2.734e-08
LPAL13_170014500 LPAL13_170014500 hypothetical protein, conserved protein coding LpaL13_17 361708 362040 + forward 17 333.0 332 5.267 0 5.133 1e-04 3.102 0.0171 -0.8103 -3.2750 6.452 1.5010 2.465 4.134 2e-04 0.0020 22.51 1.0160 5.184 0 42.99 5.426 1.0338 44.87 31.72 1.873e+03 10.109 0.00 0.00 0.00 -0.2222 20.80 0 -2.3740 3.121 -2.1820 0.0027 2.110e-02 8.269e-06 1.362e-04 0.000e+00 1.681e-03 4.355 1.129e-01 2.591e-02 9.014e-04 2.423e-06
LPAL13_000011800 LPAL13_000011800 hypothetical protein, conserved protein coding LPAL13_SCAF000076 446 640 - reverse Not Assigned 195.0 194 5.215 0 5.876 0e+00 1.508 0.2266 -2.2090 -3.9770 4.642 0.5524 1.768 3.842 5e-04 0.0039 11.97 0.9815 5.313 0 100.45 6.650 0.1540 16.46 11.57 6.253e+02 5.470 0.42 0.58 0.42 -0.8171 26.37 0 -2.9410 1.650 -4.2470 0.1039 2.807e-01 4.541e-06 1.267e-05 6.120e-01 3.178e-03 3.750 3.144e+00 8.385e-01 3.463e-02 3.598e-03
LPAL13_000035500 LPAL13_000035500 hypothetical protein, conserved protein coding LPAL13_SCAF000492 7045 7410 + forward Not Assigned 366.0 365 4.506 0 4.504 0e+00 4.085 0.0261 4.2390 0.8008 4.573 0.5655 3.439 7.496 0e+00 0.0000 529.10 0.5721 7.875 0 18.87 4.238 51.0721 963.80 689.98 4.309e+05 18.173 0.00 1.00 0.00 4.2760 43.65 0 2.7960 2.914 -2.1340 0.0049 3.282e-02 1.711e-12 7.000e-09 8.934e-01 5.314e-07 4.559 7.298e-01 1.601e-01 1.633e-03 7.997e-06
LPAL13_000026500 LPAL13_000026500 hypothetical protein protein coding LPAL13_SCAF000301 144 494 - reverse Not Assigned 351.0 350 4.424 0 4.383 0e+00 3.024 0.0953 0.2999 -2.3730 5.646 2.0132 2.673 4.398 1e-04 0.0012 46.88 0.7218 6.129 0 25.13 4.651 3.0463 76.80 54.67 5.564e+03 12.336 0.00 1.00 0.00 0.9356 28.87 0 -0.7437 2.214 -3.5240 0.0304 1.179e-01 1.053e-07 4.659e-06 9.682e-01 1.058e-03 3.944 1.908e-03 4.838e-04 1.013e-02 3.081e-04
LPAL13_000014000 LPAL13_000014000 hypothetical protein protein coding LPAL13_SCAF000119 655 942 + forward Not Assigned 288.0 287 4.336 0 4.322 0e+00 4.000 0.0071 2.3130 -1.0970 1.833 1.9220 3.410 7.178 0e+00 0.0000 129.60 0.5067 8.557 0 16.74 4.065 12.1331 203.23 145.90 1.498e+04 13.093 0.00 1.00 0.00 2.2670 54.00 0 1.0460 3.522 -0.8407 0.0008 8.896e-03 1.192e-14 1.052e-10 9.791e-01 1.809e-05 4.374 6.912e-01 1.580e-01 2.644e-04 2.098e-07
LPAL13_220019500 LPAL13_220019500 hypothetical protein protein coding LpaL13_22 578260 578538 + forward 22 279.0 278 3.869 0 3.866 0e+00 3.186 0.0322 3.3350 0.3472 2.695 0.7936 2.988 7.415 0e+00 0.0000 287.50 0.4908 7.882 0 13.63 3.768 32.7940 446.95 322.70 8.664e+04 12.589 0.00 1.00 0.00 3.4060 45.98 0 2.3050 2.810 -2.3530 0.0065 3.956e-02 1.711e-12 2.772e-09 9.791e-01 9.724e-07 3.704 9.768e-02 2.637e-02 2.181e-03 1.427e-05 
knitr::kable(head(sus_sig$deseq$downs$sensitive_vs_resistant, n = 20))
gid annotgeneproduct annotgenetype chromosome start end strand annotgeneentrezid annotgenename annotstrand annotchromosome annotcdslength length deseq_logfc deseq_adjp edger_logfc edger_adjp limma_logfc limma_adjp basic_nummed basic_denmed basic_numvar basic_denvar basic_logfc basic_t basic_p basic_adjp deseq_basemean deseq_lfcse deseq_stat deseq_p ebseq_fc ebseq_logfc ebseq_c1mean ebseq_c2mean ebseq_mean ebseq_var ebseq_postfc ebseq_ppee ebseq_ppde ebseq_adjp edger_logcpm edger_lr edger_p limma_ave limma_t limma_b limma_p limma_adjp_ihw deseq_adjp_ihw edger_adjp_ihw ebseq_adjp_ihw basic_adjp_ihw lfc_meta lfc_var lfc_varbymed p_meta p_var
LPAL13_000033300 LPAL13_000033300 hypothetical protein, conserved protein coding LPAL13_SCAF000463 551 811 + forward Not Assigned 261.0 260 -4.464 0.0026 -4.408 0.0030 -7.262 0e+00 -3.3430 3.4550 12.7819 0.0608 -6.797 -10.000 0 0e+00 134.700 1.2270 -3.637 0.0003 0.1605 -2.6393 339.45 54.476 139.969 2.623e+04 0.1688 0.0000 0.0000 0.0000 2.2950 12.97 0.0003 -0.9429 -8.222 13.700 0 4.930e-08 3.156e-03 3.526e-03 0.000e+00 3.138e-08 -5.378 0.000e+00 0.000e+00 1.971e-04 2.955e-08
LPAL13_000038400 LPAL13_000038400 expression-site associated gene (esag3), putative protein coding LPAL13_SCAF000573 101 1360 + forward Not Assigned 1260.0 1259 -2.815 0.0000 -2.813 0.0000 -3.784 0e+00 4.6300 8.2100 3.5430 0.0361 -3.581 -9.948 0 0e+00 3720.000 0.5328 -5.283 0.0000 0.1947 -2.3605 9563.26 1862.157 4172.487 1.879e+07 0.1996 0.0000 0.0000 0.0000 7.0800 31.12 0.0000 5.7700 -6.294 8.671 0 4.873e-06 6.859e-06 1.468e-06 0.000e+00 4.396e-08 -3.241 3.350e-02 -1.034e-02 6.075e-08 3.322e-15
LPAL13_350063000 LPAL13_350063000 hypothetical protein protein coding LpaL13_35 1964328 1964543 - reverse 35 216.0 215 -2.757 0.0000 -2.737 0.0000 -3.704 0e+00 -2.1430 1.1750 2.1229 0.2166 -3.318 -10.830 0 0e+00 21.340 0.4698 -5.868 0.0000 0.1753 -2.5125 58.91 10.316 24.895 7.197e+02 0.1934 0.0000 1.0000 0.0000 -0.3625 34.37 0.0000 -1.4530 -8.082 12.190 0 4.930e-08 3.292e-07 3.841e-07 9.791e-01 1.500e-09 -3.037 1.963e-02 -6.464e-03 2.994e-09 6.631e-18
LPAL13_140019300 LPAL13_140019300 bt1 family, putative protein coding LpaL13_14 530784 531350 + forward 14 567.0 566 -2.653 0.0000 -2.651 0.0000 -2.671 0e+00 4.7250 7.1000 0.6657 1.0271 -2.375 -7.180 0 0e+00 1940.000 0.3605 -7.359 0.0000 0.1910 -2.3882 5126.37 979.260 2223.392 6.182e+06 0.1961 0.0000 1.0000 0.0000 6.1420 64.15 0.0000 5.4100 -7.759 14.430 0 8.114e-08 4.900e-11 2.142e-12 9.560e-01 4.518e-05 -2.777 1.910e-01 -6.879e-02 2.735e-11 2.228e-21
LPAL13_310035500 LPAL13_310035500 hypothetical protein protein coding LpaL13_31 1198439 1198957 - reverse 31 519.0 518 -2.525 0.0033 -2.426 0.0013 -3.477 0e+00 -4.1400 -0.4281 4.1920 0.4519 -3.712 -8.576 0 0e+00 7.035 0.7117 -3.548 0.0004 0.2972 -1.7503 19.35 5.744 9.826 3.556e+02 0.3460 0.0000 0.0000 0.0000 -1.9660 14.90 0.0001 -3.1990 -7.723 8.857 0 1.220e-07 3.383e-03 1.342e-03 0.000e+00 5.704e-08 -2.836 1.131e-02 -3.987e-03 1.671e-04 3.986e-08
LPAL13_000012000 LPAL13_000012000 hypothetical protein protein coding LPAL13_SCAF000080 710 1159 - reverse Not Assigned 450.0 449 -2.495 0.0010 -2.491 0.0003 -4.411 0e+00 0.1042 3.9370 8.7982 0.1681 -3.832 -6.689 0 0e+00 212.100 0.6323 -3.946 0.0001 0.2657 -1.9123 490.41 130.281 238.320 6.077e+04 0.2767 0.0243 0.9757 0.0243 2.9530 18.17 0.0000 1.3280 -5.852 6.927 0 2.056e-05 1.155e-03 4.073e-04 9.680e-01 8.497e-06 -3.046 6.800e-02 -2.232e-02 3.334e-05 1.707e-09
LPAL13_310039200 LPAL13_310039200 hypothetical protein protein coding LpaL13_31 1301745 1301972 - reverse 31 228.0 227 -2.348 0.0000 -2.345 0.0000 -2.711 0e+00 1.4340 3.7450 1.6362 0.2106 -2.311 -8.385 0 0e+00 198.300 0.3819 -6.150 0.0000 0.3053 -1.7118 426.88 130.311 219.282 4.004e+04 0.3156 0.1979 0.8021 0.1979 2.8580 43.82 0.0000 2.0540 -6.751 10.360 0 1.849e-06 8.256e-08 7.702e-09 7.442e-01 8.324e-08 -2.621 8.058e-02 -3.074e-02 1.905e-09 6.879e-18
LPAL13_340039600 LPAL13_340039600 hypothetical protein protein coding LpaL13_34 1247554 1247757 - reverse 34 204.0 203 -2.336 0.0001 -2.334 0.0000 -3.414 0e+00 1.1870 4.2690 4.1252 0.1033 -3.081 -7.803 0 0e+00 232.400 0.4936 -4.733 0.0000 0.2185 -2.1941 605.18 132.244 274.125 6.529e+04 0.2234 0.0000 1.0000 0.0000 3.0740 26.01 0.0000 1.9780 -5.646 6.302 0 3.888e-05 5.505e-05 1.474e-05 8.934e-01 7.620e-07 -2.690 1.959e-03 -7.283e-04 9.841e-07 1.141e-12
LPAL13_310031300 LPAL13_310031300 hypothetical protein, conserved protein coding LpaL13_31 1084772 1085059 - reverse 31 288.0 287 -2.255 0.0014 -2.248 0.0008 -3.859 0e+00 -0.9701 2.1590 4.1051 0.7645 -3.129 -6.824 0 0e+00 65.480 0.5902 -3.821 0.0001 0.2715 -1.8808 150.60 40.886 73.800 7.795e+03 0.2884 0.1534 0.8466 0.1534 1.2550 16.24 0.0001 -0.2054 -6.412 8.198 0 4.634e-06 1.439e-03 9.162e-04 8.673e-01 2.258e-06 -2.799 7.637e-02 -2.729e-02 6.295e-05 4.451e-09
LPAL13_310031000 LPAL13_310031000 hypothetical protein, conserved protein coding LpaL13_31 1075172 1075459 - reverse 31 288.0 287 -2.217 0.0000 -2.201 0.0000 -3.375 0e+00 -1.7270 0.9896 3.2962 0.5154 -2.717 -6.778 0 0e+00 27.090 0.4444 -4.989 0.0000 0.3096 -1.6916 59.78 18.499 30.883 1.138e+03 0.3299 0.3062 0.6938 0.3062 0.0538 26.16 0.0000 -1.1310 -7.852 11.790 0 8.114e-08 2.381e-05 1.185e-05 7.316e-01 2.436e-06 -2.598 0.000e+00 0.000e+00 3.075e-07 9.255e-14
LPAL13_000038500 LPAL13_000038500 hypothetical protein protein coding LPAL13_SCAF000575 39 251 + forward Not Assigned 213.0 212 -2.142 0.0030 -2.130 0.0032 -3.662 0e+00 -1.9850 1.3650 4.9749 0.6608 -3.350 -6.945 0 0e+00 32.660 0.5980 -3.582 0.0003 0.2950 -1.7611 81.03 23.898 41.037 2.717e+03 0.3194 0.0219 0.9781 0.0219 0.2414 12.84 0.0003 -1.2750 -6.274 7.073 0 6.329e-06 3.021e-03 3.685e-03 9.680e-01 2.258e-06 -2.540 4.506e-02 -1.774e-02 2.268e-04 3.858e-08
LPAL13_140019100 LPAL13_140019100 bt1 family, putative protein coding LpaL13_14 525164 525514 + forward 14 351.0 350 -2.085 0.0000 -2.084 0.0000 -2.270 0e+00 3.9610 6.0230 0.5051 0.5074 -2.062 -8.394 0 0e+00 900.900 0.2929 -7.119 0.0000 0.2498 -2.0009 2157.85 539.110 1024.734 8.522e+05 0.2536 0.0000 1.0000 0.0000 5.0350 60.97 0.0000 4.5830 -8.750 18.360 0 1.675e-08 2.969e-10 8.398e-12 8.934e-01 2.258e-06 -2.201 2.558e-02 -1.162e-02 8.466e-13 5.626e-25
LPAL13_000012100 LPAL13_000012100 hypothetical protein protein coding LPAL13_SCAF000080 1637 1894 - reverse Not Assigned 258.0 257 -2.077 0.0107 -2.069 0.0076 -3.736 0e+00 -2.2070 1.1620 6.5920 0.6332 -3.370 -6.277 0 0e+00 31.370 0.6642 -3.128 0.0018 0.3004 -1.7351 69.77 20.952 35.598 2.066e+03 0.3283 0.0073 0.9927 0.0073 0.2324 10.72 0.0011 -1.3790 -5.597 4.952 0 4.490e-05 1.086e-02 9.088e-03 9.780e-01 1.036e-05 -2.487 1.387e-01 -5.575e-02 9.408e-04 7.864e-07
LPAL13_050005000 LPAL13_050005000 hypothetical protein protein coding LpaL13_05 3394 3612 - reverse 5 219.0 218 -2.022 0.0038 -2.019 0.0018 -2.980 0e+00 -0.0244 2.6770 3.6342 0.1568 -2.701 -7.147 0 0e+00 92.660 0.5776 -3.500 0.0005 0.2979 -1.7473 195.27 58.154 99.289 7.853e+03 0.3064 0.0000 1.0000 0.0000 1.7490 14.27 0.0002 0.4891 -5.487 5.509 0 6.053e-05 3.865e-03 2.037e-03 8.934e-01 2.298e-06 -2.259 1.872e-02 -8.287e-03 2.080e-04 5.564e-08
LPAL13_340039700 LPAL13_340039700 snare domain containing protein, putative protein coding LpaL13_34 1248192 1248947 - reverse 34 756.0 755 -1.890 0.0000 -1.889 0.0000 -2.149 0e+00 4.6280 6.7120 0.7919 0.0938 -2.084 -10.970 0 0e+00 1426.000 0.3010 -6.278 0.0000 0.2803 -1.8347 3280.43 919.659 1627.891 1.661e+06 0.2841 0.0000 1.0000 0.0000 5.6970 47.03 0.0000 5.2740 -7.285 12.540 0 2.533e-07 5.661e-08 2.434e-09 9.698e-01 1.500e-09 -1.969 6.690e-03 -3.397e-03 3.045e-10 7.844e-20
LPAL13_170015400 LPAL13_170015400 hypothetical protein, conserved protein coding LpaL13_17 395975 396307 + forward 17 333.0 332 -1.739 0.0000 -1.736 0.0000 -2.077 0e+00 1.3560 3.2620 1.0759 0.1243 -1.906 -8.631 0 0e+00 150.600 0.2986 -5.825 0.0000 0.3481 -1.5226 296.39 103.157 161.127 1.504e+04 0.3526 0.0000 1.0000 0.0000 2.4540 35.86 0.0000 1.9690 -6.586 9.689 0 2.759e-06 4.103e-07 2.448e-07 9.791e-01 4.953e-08 -1.803 1.522e-02 -8.438e-03 5.796e-09 1.373e-17
LPAL13_330021800 LPAL13_330021800 lipase precursor-like protein protein coding LpaL13_33 571895 572905 + forward 33 1011.0 1010 -1.576 0.0000 -1.576 0.0000 -1.595 0e+00 6.4530 8.1050 0.6133 0.5172 -1.652 -6.480 0 0e+00 4358.000 0.2767 -5.695 0.0000 0.3083 -1.6977 9733.17 3000.563 5020.346 1.650e+07 0.3126 0.0000 1.0000 0.0000 7.3110 39.18 0.0000 6.9130 -5.664 6.219 0 3.616e-05 6.451e-07 4.515e-08 1.000e+00 4.736e-05 -1.620 6.154e-02 -3.799e-02 1.274e-07 4.399e-14
LPAL13_140019200 LPAL13_140019200 inositol-3-phosphate synthase protein coding LpaL13_14 527711 529291 + INO1 forward 14 1581.0 1580 -1.465 0.0000 -1.463 0.0000 -1.540 0e+00 8.8290 10.4100 0.2203 0.4390 -1.579 -7.488 0 0e+00 19840.000 0.2244 -6.528 0.0000 0.3520 -1.5065 40764.45 14347.428 22272.534 2.188e+08 0.3546 0.0000 1.0000 0.0000 9.4970 50.86 0.0000 9.2710 -6.797 10.650 0 1.585e-06 9.542e-09 3.150e-10 1.000e+00 4.504e-05 -1.544 5.157e-02 -3.340e-02 1.385e-09 5.474e-18
LPAL13_330021900 LPAL13_330021900 hypothetical protein, conserved protein coding LpaL13_33 574061 574276 + forward 33 216.0 215 -1.454 0.0000 -1.452 0.0000 -1.764 0e+00 2.2980 4.0580 0.6084 0.8373 -1.759 -5.816 0 3e-04 260.300 0.2674 -5.439 0.0000 0.2851 -1.8106 620.52 176.883 309.975 1.118e+05 0.2952 0.0007 0.9993 0.0007 3.2420 30.71 0.0000 2.8150 -6.296 8.738 0 5.888e-06 3.031e-06 1.890e-06 9.791e-01 2.490e-04 -1.614 6.533e-03 -4.047e-03 3.803e-08 1.813e-16
LPAL13_210015500 LPAL13_210015500 core histone h2a/h2b/h3/h4, putative protein coding LpaL13_21 326108 326506 + forward 21 399.0 398 -1.410 0.0012 -1.410 0.0003 -1.831 3e-04 4.9850 6.6460 2.6364 0.2834 -1.661 -4.840 0 4e-04 2218.000 0.3640 -3.875 0.0001 0.5020 -0.9943 3301.46 1657.261 2150.521 3.170e+06 0.5213 0.9841 0.0159 0.9841 6.3340 18.83 0.0000 5.6570 -4.790 3.040 0 2.888e-04 1.217e-03 2.228e-04 4.526e-02 4.717e-04 -1.600 2.080e-03 -1.300e-03 4.373e-05 2.979e-09 
sus_ma <- sus_table[["plots"]][["sensitive_vs_resistant"]][["deseq_ma_plots"]][["plot"]]
pp(file = "images/sus_ma.png", image = sus_ma)

## test <- ggplt(sus_ma)

6.4 Ontology searches

Now let us look for ontology categories which are increased in the 2.3 samples followed by the 2.2 samples.

## Gene categories more represented in the 2.3 group.
zy_go_up <- sm(simple_goseq(sig_genes = zy_sig[["deseq"]][["ups"]][[1]],
                            go_db = lp_go, length_db = lp_lengths))

## Gene categories more represented in the 2.2 group.
zy_go_down <- sm(simple_goseq(sig_genes = zy_sig[["deseq"]][["downs"]][[1]],
                              go_db = lp_go, length_db = lp_lengths))

6.4.1 A couple plots from the differential expression

6.4.1.1 Number of genes in agreement among DE methods, 2.3 more than 2.2

In the function ‘combined_de_tables()’ above, one of the tasks performed is to look at the agreement among DESeq2, limma, and edgeR. The following show a couple of these for the set of genes observed with a fold-change >= |2| and adjusted p-value <= 0.05.

zy_table[["venns"]][[1]][["p_lfc1"]][["up_noweight"]]

6.4.1.2 Number of genes in agreement among DE methods, 2.2 more than 2.3

zy_table[["venns"]][[1]][["p_lfc1"]][["down_noweight"]]

6.4.1.3 goseq ontology plots of groups of genes, 2.3 more than 2.2

zy_go_up$pvalue_plots$bpp_plot_over

6.4.1.4 goseq ontology plots of groups of genes, 2.2 more than 2.3

zy_go_down$pvalue_plots$bpp_plot_over

6.5 Look for agreement between sensitivity and zymodemes

Remind myself, the data structures are (zy|sus)_(de|table|sig).

zy_df <- zy_table[["data"]][["z23_vs_z22"]]
sus_df <- sus_table[["data"]][["sensitive_vs_resistant"]]

both_df <- merge(zy_df, sus_df, by = "row.names")
plot_df <- both_df[, c("deseq_logfc.x", "deseq_logfc.y")]
rownames(plot_df) <- both_df[["Row.names"]]
colnames(plot_df) <- c("z23_vs_z22", "sensitive_vs_resistant")

compare <- plot_linear_scatter(plot_df)
## Warning in plot_multihistogram(df): NAs introduced by coercion
pp(file = "images/compare_sus_zy.png", image = compare$scatter)

compare$cor
## 
##  Pearson's product-moment correlation
## 
## data:  df[, 1] and df[, 2]
## t = -244, df = 8549, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.9378 -0.9325
## sample estimates:
##     cor 
## -0.9352

6.6 Zymodeme enzyme gene IDs

Najib read me an email listing off the gene names associated with the zymodeme classification. I took those names and cross referenced them against the Leishmania panamensis gene annotations and found the following:

They are:

  1. ALAT: LPAL13_120010900 – alanine aminotransferase
  2. ASAT: LPAL13_340013000 – aspartate aminotransferase
  3. G6PD: LPAL13_000054100 – glucase-6-phosphate 1-dehydrogenase
  4. NH: LPAL13_14006100, LPAL13_180018500 – inosine-guanine nucleoside hydrolase
  5. MPI: LPAL13_320022300 (maybe) – mannose phosphate isomerase (I chose phosphomannose isomerase)

Given these 6 gene IDs (NH has two gene IDs associated with it), I can do some looking for specific differences among the various samples.

6.6.1 Expression levels of zymodeme genes

The following creates a colorspace (red to green) heatmap showing the observed expression of these genes in every sample.

my_genes <- c("LPAL13_120010900", "LPAL13_340013000", "LPAL13_000054100",
              "LPAL13_140006100", "LPAL13_180018500", "LPAL13_320022300",
              "other")
my_names <- c("ALAT", "ASAT", "G6PD", "NHv1", "NHv2", "MPI", "other")

zymo_expt <- exclude_genes_expt(zy_norm, ids = my_genes, method = "keep")
## Before removal, there were 8551 genes, now there are 6.
## There are 59 samples which kept less than 90 percent counts.
## TMRC20001 TMRC20065 TMRC20005 TMRC20066 TMRC20039 TMRC20037 TMRC20038 TMRC20067 
##    0.1313    0.1250    0.1321    0.1061    0.1301    0.1102    0.1131    0.1165 
## TMRC20068 TMRC20041 TMRC20015 TMRC20009 TMRC20010 TMRC20016 TMRC20011 TMRC20012 
##    0.1157    0.1182    0.1150    0.1137    0.1101    0.1062    0.1104    0.1209 
## TMRC20013 TMRC20017 TMRC20014 TMRC20018 TMRC20019 TMRC20070 TMRC20020 TMRC20021 
##    0.1208    0.1066    0.1092    0.1147    0.1225    0.1128    0.1103    0.1063 
## TMRC20022 TMRC20024 TMRC20036 TMRC20069 TMRC20033 TMRC20026 TMRC20031 TMRC20076 
##    0.1307    0.1126    0.1203    0.1164    0.1128    0.1386    0.1005    0.1203 
## TMRC20073 TMRC20055 TMRC20079 TMRC20071 TMRC20078 TMRC20042 TMRC20058 TMRC20072 
##    0.1228    0.1347    0.1269    0.1234    0.1341    0.1315    0.1182    0.1431 
## TMRC20059 TMRC20048 TMRC20060 TMRC20077 TMRC20074 TMRC20063 TMRC20053 TMRC20052 
##    0.1104    0.1033    0.1087    0.1221    0.1209    0.1168    0.1183    0.1106 
## TMRC20064 TMRC20075 TMRC20051 TMRC20050 TMRC20049 TMRC20062 TMRC20080 TMRC20043 
##    0.1140    0.1113    0.1283    0.1154    0.1394    0.1286    0.1155    0.1138 
## TMRC20054 TMRC20046 TMRC20044 
##    0.1278    0.1368    0.1338
zymo_heatmap <- plot_sample_heatmap(zymo_expt, row_label = my_names)
zymo_heatmap

6.7 Empirically observed Zymodeme genes from differential expression analysis

In contrast, the following plots take the set of genes which are shared among all differential expression methods (|lfc| >= 1.0 and adjp <= 0.05) and use them to make categories of genes which are increased in 2.3 or 2.2.

shared_zymo <- intersect_significant(zy_table)
## Deleting the file excel/intersect_significant.xlsx before writing the tables.
up_shared <- shared_zymo[["ups"]][[1]][["data"]][["all"]]
rownames(up_shared)
##  [1] "LPAL13_000033300" "LPAL13_000012000" "LPAL13_000038500" "LPAL13_000012100"
##  [5] "LPAL13_310031300" "LPAL13_000038400" "LPAL13_340039600" "LPAL13_050005000"
##  [9] "LPAL13_310031000" "LPAL13_310039200" "LPAL13_350063000" "LPAL13_180013900"
## [13] "LPAL13_140019300" "LPAL13_210015500" "LPAL13_340039700" "LPAL13_350013200"
## [17] "LPAL13_170015400" "LPAL13_250006300" "LPAL13_330024000" "LPAL13_350073400"
## [21] "LPAL13_140019100" "LPAL13_320038700" "LPAL13_210005000" "LPAL13_140019200"
## [25] "LPAL13_240009700" "LPAL13_000052700" "LPAL13_330021800" "LPAL13_160014500"
## [29] "LPAL13_230011200" "LPAL13_350073200" "LPAL13_330021900" "LPAL13_050009600"
## [33] "LPAL13_250025700" "LPAL13_230011500" "LPAL13_160014100" "LPAL13_040007800"
## [37] "LPAL13_230011400" "LPAL13_020006700" "LPAL13_310028500" "LPAL13_230011300"
## [41] "LPAL13_310032500" "LPAL13_140015200"
upshared_expt <- exclude_genes_expt(zy_norm, ids = rownames(up_shared), method = "keep")
## Before removal, there were 8551 genes, now there are 42.
## There are 59 samples which kept less than 90 percent counts.
## TMRC20001 TMRC20065 TMRC20005 TMRC20066 TMRC20039 TMRC20037 TMRC20038 TMRC20067 
##    0.3732    0.4452    0.1227    0.4003    0.1662    0.4253    0.5511    0.3428 
## TMRC20068 TMRC20041 TMRC20015 TMRC20009 TMRC20010 TMRC20016 TMRC20011 TMRC20012 
##    0.4034    0.1735    0.4406    0.1521    0.4268    0.3287    0.1662    0.1494 
## TMRC20013 TMRC20017 TMRC20014 TMRC20018 TMRC20019 TMRC20070 TMRC20020 TMRC20021 
##    0.3753    0.1677    0.1777    0.3503    0.1409    0.4173    0.1316    0.3789 
## TMRC20022 TMRC20024 TMRC20036 TMRC20069 TMRC20033 TMRC20026 TMRC20031 TMRC20076 
##    0.1415    0.1522    0.1971    0.1728    0.1518    0.1365    0.1279    0.1494 
## TMRC20073 TMRC20055 TMRC20079 TMRC20071 TMRC20078 TMRC20042 TMRC20058 TMRC20072 
##    0.4929    0.1837    0.5541    0.5077    0.1928    0.1522    0.6068    0.1794 
## TMRC20059 TMRC20048 TMRC20060 TMRC20077 TMRC20074 TMRC20063 TMRC20053 TMRC20052 
##    0.2983    0.3285    0.1431    0.1394    0.1635    0.1558    0.1790    0.4487 
## TMRC20064 TMRC20075 TMRC20051 TMRC20050 TMRC20049 TMRC20062 TMRC20080 TMRC20043 
##    0.4130    0.3533    0.6096    0.1523    0.1707    0.6237    0.4556    0.4232 
## TMRC20054 TMRC20046 TMRC20044 
##    0.5415    0.1726    0.1745

We can plot a quick heatmap to get a sense of the differences observed between the genes which are different between the two zymodemes.

6.7.1 Heatmap of zymodeme gene expression increased in 2.3 vs. 2.2

high_23_heatmap <- plot_sample_heatmap(upshared_expt, row_label = rownames(up_shared))
high_23_heatmap

6.7.2 Heatmap of zymodeme gene expression increased in 2.2 vs. 2.3

down_shared <- shared_zymo[["downs"]][[1]][["data"]][["all"]]
downshared_expt <- exclude_genes_expt(zy_norm, ids = rownames(down_shared), method = "keep")
## Before removal, there were 8551 genes, now there are 68.
## There are 59 samples which kept less than 90 percent counts.
## TMRC20001 TMRC20065 TMRC20005 TMRC20066 TMRC20039 TMRC20037 TMRC20038 TMRC20067 
##    0.2203    0.1876    0.6711    0.2263    0.6573    0.2104    0.2013    0.2387 
## TMRC20068 TMRC20041 TMRC20015 TMRC20009 TMRC20010 TMRC20016 TMRC20011 TMRC20012 
##    0.2022    0.6894    0.1876    0.6365    0.1681    0.2075    0.5677    0.5508 
## TMRC20013 TMRC20017 TMRC20014 TMRC20018 TMRC20019 TMRC20070 TMRC20020 TMRC20021 
##    0.1654    0.6549    0.6494    0.1595    0.6528    0.1856    0.6853    0.1536 
## TMRC20022 TMRC20024 TMRC20036 TMRC20069 TMRC20033 TMRC20026 TMRC20031 TMRC20076 
##    0.6763    0.7153    0.6761    0.6973    0.7206    0.6958    0.6111    0.6079 
## TMRC20073 TMRC20055 TMRC20079 TMRC20071 TMRC20078 TMRC20042 TMRC20058 TMRC20072 
##    0.1910    0.7106    0.1867    0.1696    0.5314    0.5479    0.2195    0.5354 
## TMRC20059 TMRC20048 TMRC20060 TMRC20077 TMRC20074 TMRC20063 TMRC20053 TMRC20052 
##    0.1378    0.1541    0.7602    0.5707    0.6624    0.6333    0.5717    0.1783 
## TMRC20064 TMRC20075 TMRC20051 TMRC20050 TMRC20049 TMRC20062 TMRC20080 TMRC20043 
##    0.1945    0.1802    0.1902    0.6151    0.6973    0.1885    0.1545    0.1727 
## TMRC20054 TMRC20046 TMRC20044 
##    0.2025    0.6338    0.6221
high_22_heatmap <- plot_sample_heatmap(downshared_expt, row_label = rownames(down_shared))
high_22_heatmap

7 SNP profiles

Now I will combine our previous samples and our new samples in the hopes of finding variant positions which help elucidate currently unknown aspects of either group via their clustering to known samples from the other group. In other words, we do not know the zymodeme annotations for the old samples nor the strain identities (or the shortcut ‘chronic vs. self-healing’) for the new samples. I hope to make educated guesses given the variant profiles. There are some differences in how the previous and current data sets were analyzed (though I have since redone the old samples so it should be trivial to remove those differences now).

I added our 2016 data to a specific TMRC2 sample sheet, dated 20191203. Thus I will load the data here. That previous data was mapped using tophat, so I will also need to make some changes to the gene names to accomodate the two mappings.

old_expt <- sm(create_expt("sample_sheets/tmrc2_samples_20191203.xlsx",
                           file_column = "tophat2file"))

tt <- lp_expt[["expressionset"]]
rownames(tt) <- gsub(pattern = "^exon_", replacement = "", x = rownames(tt))
rownames(tt) <- gsub(pattern = "\\.E1$", replacement = "", x = rownames(tt))
lp_expt$expressionset <- tt

tt <- old_expt$expressionset
rownames(tt) <- gsub(pattern = "^exon_", replacement = "", x = rownames(tt))
rownames(tt) <- gsub(pattern = "\\.1$", replacement = "", x = rownames(tt))
old_expt$expressionset <- tt
rm(tt)

7.1 Create the SNP expressionset

One other important caveat, we have a group of new samples which have not yet run through the variant search pipeline, so I need to remove them from consideration. Though it looks like they finished overnight…

## The next line drops the samples which are missing the SNP pipeline.
lp_snp <- subset_expt(lp_expt, subset="!is.na(pData(lp_expt)[['bcftable']])")
## subset_expt(): There were 71, now there are 66 samples.
new_snps <- sm(count_expt_snps(lp_snp, annot_column = "bcftable"))
old_snps <- sm(count_expt_snps(old_expt, annot_column = "bcftable", snp_column = 2))

both_snps <- combine_expts(new_snps, old_snps)
both_norm <- sm(normalize_expt(both_snps, transform = "log2", convert = "cpm", filter = TRUE))

## strains <- both_norm[["design"]][["strain"]]
both_strain <- set_expt_conditions(both_norm, fact = "strain")

The data structure ‘both_norm’ now contains our 2016 data along with the newer data collected since 2019.

7.2 Plot of SNP profiles for zymodemes

The following plot shows the SNP profiles of all samples (old and new) where the colors at the top show either the 2.2 strains (orange), 2.3 strains (green), the previous samples (purple), or the various lab strains (pink etc).

old_new_variant_heatmap <- plot_disheat(both_norm)
pp(file = "images/raw_snp_disheat.png", image = old_new_variant_heatmap,
   height = 12, width = 12)

The function get_snp_sets() takes the provided metadata factor (in this case ‘condition’) and looks for variants which are exclusive to each element in it. In this case, this is looking for differences between 2.2 and 2.3, as well as the set shared among them.

snp_sets <- get_snp_sets(both_snps, factor = "condition")
## The factor z2.3 has 26 rows.
## The factor z2.2 has 28 rows.
## The factor unknown has 7 rows.
## The factor z2.1 has 3 rows.
## The factor z2.4 has only 1 row.
## The factor null has only 1 row.
## The factor sh has 13 rows.
## The factor chr has 14 rows.
## The factor inf has 6 rows.
## Iterating over 727 elements.
both_expt <- combine_expts(lp_expt, old_expt)

snp_genes <- sm(snps_vs_genes(both_expt, snp_sets, expt_name_col = "chromosome"))
## I think we have some metrics here we can plot...
snp_subset <- sm(snp_subset_genes(
  both_expt, both_snps,
  genes = c("LPAL13_120010900", "LPAL13_340013000", "LPAL13_000054100",
            "LPAL13_140006100", "LPAL13_180018500", "LPAL13_320022300")))
zymo_heat <- plot_sample_heatmap(snp_subset, row_label = rownames(exprs(snp_subset)))
zymo_heat

Didn’t I create a set of densities by chromosome? Oh I think they come in from get_snp_sets()

7.3 SNPS associated with clinical response in the TMRC samples

clinical_sets <- get_snp_sets(new_snps, factor = "clinicalresponse")
## The factor cure has 26 rows.
## The factor failure has 21 rows.
## The factor laboratory line has only 1 row.
## The factor laboratory line miltefosine resistant has only 1 row.
## The factor nd has 13 rows.
## The factor reference strain has 4 rows.
## Iterating over 695 elements.
density_vec <- clinical_sets[["density"]]
chromosome_idx <- grep(pattern = "LpaL", x = names(density_vec))
density_df <- as.data.frame(density_vec[chromosome_idx])
density_df[["chr"]] <- rownames(density_df)
colnames(density_df) <- c("density_vec", "chr")
ggplot(density_df, aes_string(x = "chr", y = "density_vec")) +
  ggplot2::geom_col() +
  ggplot2::theme(axis.text = ggplot2::element_text(size = 10, colour = "black"),
                 axis.text.x = ggplot2::element_text(angle = 90, vjust = 0.5))

## clinical_written <- write_variants(new_snps)

7.3.1 Cross reference these variants by gene

clinical_genes <- sm(snps_vs_genes(lp_expt, clinical_sets, expt_name_col = "chromosome"))

snp_density <- merge(as.data.frame(clinical_genes[["summary_by_gene"]]),
                     as.data.frame(fData(lp_expt)),
                     by = "row.names")
snp_density <- snp_density[, c(1, 2, 4, 15)]
colnames(snp_density) <- c("name", "snps", "product", "length")
snp_density[["product"]] <- tolower(snp_density[["product"]])
snp_density[["length"]] <- as.numeric(snp_density[["length"]])
snp_density[["density"]] <- snp_density[["snps"]] / snp_density[["length"]]
snp_idx <- order(snp_density[["density"]], decreasing = TRUE)
snp_density <- snp_density[snp_idx, ]

removers <- c("amastin", "gp63", "leishmanolysin")
for (r in removers) {
  drop_idx <- grepl(pattern = r, x = snp_density[["product"]])
  snp_density <- snp_density[!drop_idx, ]
}
## Filter these for [A|a]mastin gp63 Leishmanolysin
clinical_snps <- snps_intersections(lp_expt, clinical_sets, chr_column = "chromosome")

fail_ref_snps <- as.data.frame(clinical_snps[["inters"]][["failure, reference strain"]])
cure_snps <- as.data.frame(clinical_snps[["inters"]][["cure"]])

head(fail_ref_snps)
##                                       seqnames  start    end width strand
## chr_LpaL13-02_pos_83184_ref_A_alt_G  LpaL13-02  83184  83185     2      +
## chr_LpaL13-05_pos_368787_ref_C_alt_T LpaL13-05 368787 368788     2      +
## chr_LpaL13-06_pos_2975_ref_T_alt_C   LpaL13-06   2975   2976     2      +
## chr_LpaL13-06_pos_66067_ref_T_alt_C  LpaL13-06  66067  66068     2      +
## chr_LpaL13-07_pos_461925_ref_A_alt_G LpaL13-07 461925 461926     2      +
## chr_LpaL13-10_pos_327353_ref_T_alt_C LpaL13-10 327353 327354     2      +
head(cure_snps)
##                                           seqnames  start    end width strand
## chr_LpaL13-03_pos_5444_ref_T_alt_C       LpaL13-03   5444   5445     2      +
## chr_LpaL13-08_pos_184791_ref_T_alt_A     LpaL13-08 184791 184792     2      +
## chr_LpaL13-12_pos_290500_ref_A_alt_G     LpaL13-12 290500 290501     2      +
## chr_LpaL13-20.1_pos_369935_ref_C_alt_T LpaL13-20.1 369935 369936     2      +
## chr_LpaL13-20.1_pos_370282_ref_C_alt_T LpaL13-20.1 370282 370283     2      +
## chr_LpaL13-20.1_pos_371356_ref_T_alt_C LpaL13-20.1 371356 371357     2      +
annot <- fData(lp_expt)
clinical_interest <- as.data.frame(clinical_snps[["gene_summaries"]][["cure"]])
clinical_interest <- merge(clinical_interest,
                           as.data.frame(clinical_snps[["gene_summaries"]][["failure, reference strain"]]),
                           by = "row.names")
rownames(clinical_interest) <- clinical_interest[["Row.names"]]
clinical_interest[["Row.names"]] <- NULL
colnames(clinical_interest) <- c("cure_snps","fail_snps")
annot <- merge(annot, clinical_interest, by = "row.names")
rownames(annot) <- annot[["Row.names"]]
annot[["Row.names"]] <- NULL
fData(lp_expt$expressionset) <- annot

8 Zymodeme for new samples

The heatmap produced here should show the variants only for the zymodeme genes.

8.1 Hunt for snp clusters

I am thinking that if we find clusters of locations which are variant, that might provide some PCR testing possibilities.

## Drop the 2.1, 2.4, unknown, and null
pruned_snps <- subset_expt(new_snps, subset="condition=='z2.2'|condition=='z2.3'")
## subset_expt(): There were 66, now there are 54 samples.
new_sets <- get_snp_sets(pruned_snps, factor = "zymodemecategorical")
## The factor z22 has 28 rows.
## The factor z23 has 26 rows.
## Iterating over 695 elements.
summary(new_sets)
##               Length Class      Mode     
## medians         3    data.frame list     
## possibilities   2    -none-     character
## intersections   3    -none-     list     
## chr_data      695    -none-     list     
## set_names       4    -none-     list     
## invert_names    4    -none-     list     
## density       695    -none-     numeric
## 1000000: 2.2
## 0100000: 2.3

summary(new_sets[["intersections"]][["10"]])
##    Length     Class      Mode 
##       888 character character
summary(new_sets[["intersections"]][["01"]])
##    Length     Class      Mode 
##     76419 character character

Thus we see that there are 511 variants associated with 2.2 and 49,790 associated with 2.3.

8.1.1 A small function for searching for potential PCR primers

The following function uses the positional data to look for sequential mismatches associated with zymodeme in the hopes that there will be some regions which would provide good potential targets for a PCR-based assay.

sequential_variants <- function(snp_sets, conditions = NULL, minimum = 3, maximum_separation = 3) {
  if (is.null(conditions)) {
    conditions <- 1
  }
  intersection_sets <- snp_sets[["intersections"]]
  intersection_names <- snp_sets[["set_names"]]
  chosen_intersection <- 1
  if (is.numeric(conditions)) {
    chosen_intersection <- conditions
  } else {
    intersection_idx <- intersection_names == conditions
    chosen_intersection <- names(intersection_names)[intersection_idx]
  }

  possible_positions <- intersection_sets[[chosen_intersection]]
  position_table <- data.frame(row.names = possible_positions)
  pat <- "^chr_(.+)_pos_(.+)_ref_.*$"
  position_table[["chr"]] <- gsub(pattern = pat, replacement = "\\1", x = rownames(position_table))
  position_table[["pos"]] <- as.numeric(gsub(pattern = pat, replacement = "\\2", x = rownames(position_table)))
  position_idx <- order(position_table[, "chr"], position_table[, "pos"])
  position_table <- position_table[position_idx, ]
  position_table[["dist"]] <- 0

  last_chr <- ""
  for (r in 1:nrow(position_table)) {
    this_chr <- position_table[r, "chr"]
    if (r == 1) {
      position_table[r, "dist"] <- position_table[r, "pos"]
      last_chr <- this_chr
      next
    }
    if (this_chr == last_chr) {
      position_table[r, "dist"] <- position_table[r, "pos"] - position_table[r - 1, "pos"]
    } else {
      position_table[r, "dist"] <- position_table[r, "pos"]
    }
    last_chr <- this_chr
  }

  ## Working interactively here.

  doubles <- position_table[["dist"]] == 1
  doubles <- position_table[doubles, ]
  write.csv(doubles, "doubles.csv")

  one_away <- position_table[["dist"]] == 2
  one_away <- position_table[one_away, ]
  write.csv(one_away, "one_away.csv")

  two_away <- position_table[["dist"]] == 3
  two_away <- position_table[two_away, ]
  write.csv(two_away, "two_away.csv")

  combined <- rbind(doubles, one_away)
  combined <- rbind(combined, two_away)
  position_idx <- order(combined[, "chr"], combined[, "pos"])
  combined <- combined[position_idx, ]

  this_chr <- ""
  for (r in 1:nrow(combined)) {
    this_chr <- combined[r, "chr"]
    if (r == 1) {
      combined[r, "dist_pair"] <- combined[r, "pos"]
      last_chr <- this_chr
      next
    }
    if (this_chr == last_chr) {
      combined[r, "dist_pair"] <- combined[r, "pos"] - combined[r - 1, "pos"]
    } else {
      combined[r, "dist_pair"] <- combined[r, "pos"]
    }
    last_chr <- this_chr
  }

  dist_pair_maximum <- 1000
  dist_pair_minimum <- 200
  dist_pair_idx <- combined[["dist_pair"]] <= dist_pair_maximum &
    combined[["dist_pair"]] >= dist_pair_minimum
  remaining <- combined[dist_pair_idx, ]
  no_weak_idx <- grepl(pattern="ref_(G|C)", x=rownames(remaining))
  remaining <- remaining[no_weak_idx, ]

  print(head(table(position_table[["dist"]])))
  sequentials <- position_table[["dist"]] <= maximum_separation
  message("There are ", sum(sequentials), " candidate regions.")

  ## The following can tell me how many runs of each length occurred, that is not quite what I want.
  ## Now use run length encoding to find the set of sequential sequentials!
  rle_result <- rle(sequentials)
  rle_values <- rle_result[["values"]]
  ## The following line is equivalent to just leaving values alone:
  ## true_values <- rle_result[["values"]] == TRUE
  rle_lengths <- rle_result[["lengths"]]
  true_sequentials <- rle_lengths[rle_values]
  rle_idx <- cumsum(rle_lengths)[which(rle_values)]

  position_table[["last_sequential"]] <- 0
  count <- 0
  for (r in rle_idx) {
    count <- count + 1
    position_table[r, "last_sequential"] <- true_sequentials[count]
  }
  message("The maximum sequential set is: ", max(position_table[["last_sequential"]]), ".")

  wanted_idx <- position_table[["last_sequential"]] >= minimum
  wanted <- position_table[wanted_idx, c("chr", "pos")]
  return(wanted)
}

zymo22_sequentials <- sequential_variants(new_sets, conditions = "z22", minimum=1, maximum_separation=2)
dim(zymo22_sequentials)
## 7 candidate regions for zymodeme 2.2 -- thus I am betting that the reference strain is a 2.2
zymo23_sequentials <- sequential_variants(new_sets, conditions = "z23",
                                          minimum = 2, maximum_separation = 2)
dim(zymo23_sequentials)
## In contrast, there are lots (587) of interesting regions for 2.3!

8.1.2 Extract a promising region from the genome

The first 4 candidate regions from my set of remaining: * Chr Pos. Distance * LpaL13-15 238433 448 * LpaL13-18 142844 613 * LpaL13-29 830342 252 * LpaL13-33 1331507 843

Lets define a couple of terms: * Third: Each of the 4 above positions. * Second: Third - Distance * End: Third + PrimerLen * Start: Second - Primerlen

In each instance, these are the last positions, so we want to grab three things:

  • The entire region from End -> Start, this way we can have a quick sanity check.
  • Start -> Second.
  • (Third -> End) <- Reverse complemented
## * LpaL13-15 238433 448
first_candidate_chr <- genome[["LpaL13_15"]]
primer_length <- 22
amplicon_length <- 448
first_candidate_third <- 238433
first_candidate_second <- first_candidate_third - amplicon_length
first_candidate_start <- first_candidate_second - primer_length
first_candidate_end <- first_candidate_third + primer_length
first_candidate_region <- subseq(first_candidate_chr, first_candidate_start, first_candidate_end)
first_candidate_region
first_candidate_5p <- subseq(first_candidate_chr, first_candidate_start, first_candidate_second)
as.character(first_candidate_5p)
first_candidate_3p <- spgs::reverseComplement(subseq(first_candidate_chr, first_candidate_third, first_candidate_end))
first_candidate_3p


## * LpaL13-18 142844 613
second_candidate_chr <- genome[["LpaL13_18"]]
primer_length <- 22
amplicon_length <- 613
second_candidate_third <- 142844
second_candidate_second <- second_candidate_third - amplicon_length
second_candidate_start <- second_candidate_second - primer_length
second_candidate_end <- second_candidate_third + primer_length
second_candidate_region <- subseq(second_candidate_chr, second_candidate_start, second_candidate_end)
second_candidate_region
second_candidate_5p <- subseq(second_candidate_chr, second_candidate_start, second_candidate_second)
as.character(second_candidate_5p)
second_candidate_3p <- spgs::reverseComplement(subseq(second_candidate_chr, second_candidate_third, second_candidate_end))
second_candidate_3p


## * LpaL13-29 830342 252
third_candidate_chr <- genome[["LpaL13_29"]]
primer_length <- 22
amplicon_length <- 252
third_candidate_third <- 830342
third_candidate_second <- third_candidate_third - amplicon_length
third_candidate_start <- third_candidate_second - primer_length
third_candidate_end <- third_candidate_third + primer_length
third_candidate_region <- subseq(third_candidate_chr, third_candidate_start, third_candidate_end)
third_candidate_region
third_candidate_5p <- subseq(third_candidate_chr, third_candidate_start, third_candidate_second)
as.character(third_candidate_5p)
third_candidate_3p <- spgs::reverseComplement(subseq(third_candidate_chr, third_candidate_third, third_candidate_end))
third_candidate_3p
## You are a garbage polypyrimidine tract.
## Which is actually interesting if the mutations mess it up.


## * LpaL13-33 1331507 843
fourth_candidate_chr <- genome[["LpaL13_33"]]
primer_length <- 22
amplicon_length <- 843
fourth_candidate_third <- 1331507
fourth_candidate_second <- fourth_candidate_third - amplicon_length
fourth_candidate_start <- fourth_candidate_second - primer_length
fourth_candidate_end <- fourth_candidate_third + primer_length
fourth_candidate_region <- subseq(fourth_candidate_chr, fourth_candidate_start, fourth_candidate_end)
fourth_candidate_region
fourth_candidate_5p <- subseq(fourth_candidate_chr, fourth_candidate_start, fourth_candidate_second)
as.character(fourth_candidate_5p)
fourth_candidate_3p <- spgs::reverseComplement(subseq(fourth_candidate_chr, fourth_candidate_third, fourth_candidate_end))
fourth_candidate_3p

8.2 Go hunting for Sanger sequencing regions

I made a fun little function which should find regions which have lots of variants associated with a given experimental factor.

pheno <- subset_expt(lp_expt, subset = "condition=='z2.2'|condition=='z2.3'")
pheno <- subset_expt(pheno, subset = "!is.na(pData(pheno)[['bcftable']])")
pheno_snps <- sm(count_expt_snps(pheno, annot_column = "bcftable"))

fun_stuff <- snp_density_primers(pheno_snps,
                                 bsgenome="BSGenome.Leishmania.panamensis.MHOMCOL81L13.v53",
                                 gff="reference/TriTrypDB-53_LpanamensisMHOMCOL81L13.gff")
drop_scaffolds <- grepl(x = rownames(fun_stuff$favorites), pattern = "SCAF")
favorite_primer_regions <- fun_stuff[["favorites"]][!drop_scaffolds, ]
favorite_primer_regions[["bin"]] <- rownames(favorite_primer_regions)
library(dplyr)
favorite_primer_regions <- favorite_primer_regions %>%
  relocate(bin)

8.3 Combine this table with 2.2/2.3 genes

Here is my note from our meeting:

Cross reference primers to DE genes of 2.2/2.3 and/or resistance/suscpetible, add a column to the primer spreadsheet with the DE genes (in retrospect I am guessing this actually means to put the logFC as a column.

One nice thing, I did a semantic removal on the lp_expt, so the set of logFC/pvalues should not have any of the offending types; thus I should be able to automagically get rid of them in the merge.

logfc <- zy_table[["data"]][["z23_vs_z22"]]
logfc_columns <- logfc[, c("deseq_logfc", "deseq_adjp")]
colnames(logfc_columns) <- c("z23_logfc", "z23_adjp")
new_table <- merge(favorite_primer_regions, logfc_columns,
                   by.x = "closest_gene_before_id", by.y = "row.names")
sus <- sus_table[["data"]][["sensitive_vs_resistant"]]
sus_columns <- sus[, c("deseq_logfc", "deseq_adjp")]
colnames(sus_columns) <- c("sus_logfc", "sus_adjp")
new_table <- merge(new_table, sus_columns,
                   by.x = "closest_gene_before_id", by.y = "row.names") %>%
  relocate(bin)
written <- write_xlsx(data=new_table,
                      excel="excel/favorite_primers_xref_zy_sus.
## Error: <text>:13:29: unexpected INCOMPLETE_STRING
## 12: written <- write_xlsx(data=new_table,
## 13:                       excel="excel/favorite_primers_xref_zy_sus.
##                                 ^

8.4 Make a heatmap describing the clustering of variants

We can cross reference the variants against the zymodeme status and plot a heatmap of the results and hopefully see how they separate.

## pruned_snps <- subset_expt(new_snps, subset="condition=='z2.2'|condition=='z2.3'")
snp_genes <- sm(snps_vs_genes(lp_expt, new_sets, expt_name_col = "chromosome"))

##new_zymo_norm <- normalize_expt(pruned_snps, filter = TRUE, convert = "cpm", norm = "quant", transform = TRUE)
##new_zymo_norm <- set_expt_conditions(new_zymo_norm, fact = "zymodemecategorical")
clinical_colors_v2 <- list(
    "z22" = "#0000cc",
    "z23" = "#cc0000")
new_zymo_norm <- normalize_expt(pruned_snps, filter = TRUE, convert = "cpm", norm = "quant", transform = TRUE) %>%
  set_expt_conditions(fact = "zymodemecategorical") %>%
  set_expt_colors(clinical_colors_v2)
## Removing 0 low-count genes (568413 remaining).
## transform_counts: Found 28454638 values equal to 0, adding 1 to the matrix.
zymo_heat <- plot_disheat(new_zymo_norm)
pp(file = "images/onlyz22_z23_snp_heatmap.png", image=zymo_heat[["plot"]])

zymo_heat[["plot"]]

8.4.1 Annotated heatmap of variants

Now let us try to make a heatmap which includes some of the annotation data.

des <- both_norm[["design"]]
undef_idx <- is.na(des[["strain"]])
des[undef_idx, "strain"] <- "unknown"

##hmcols <- colorRampPalette(c("yellow","black","darkblue"))(256)
correlations <- hpgl_cor(exprs(both_norm))

zymo_missing_idx <- is.na(des[["zymodemecategorical"]])
des[["zymodemecategorical"]] <- as.character(des[["cymodemecategorical"]])
## Error in `[[<-.data.frame`(`*tmp*`, "zymodemecategorical", value = character(0)): replacement has 0 rows, data has 99
des[["clinicalcategorical"]] <- as.character(des[["clinicalcategorical"]])
des[zymo_missing_idx, "zymodemecategorical"] <- "unknown"
mydendro <- list(
  "clustfun" = hclust,
  "lwd" = 2.0)
col_data <- as.data.frame(des[, c("zymodemecategorical", "clinicalcategorical")])

unknown_clinical <- is.na(col_data[["clinicalcategorical"]])
row_data <- as.data.frame(des[, c("strain")])
colnames(col_data) <- c("zymodeme", "outcome")
col_data[unknown_clinical, "outcome"] <- "undefined"

colnames(row_data) <- c("strain")
myannot <- list(
  "Col" = list("data" = col_data),
  "Row" = list("data" = row_data))
myclust <- list("cuth" = 1.0,
                "col" = BrewerClusterCol)
mylabs <- list(
  "Row" = list("nrow" = 4),
  "Col" = list("nrow" = 4))
hmcols <- colorRampPalette(c("darkblue", "beige"))(240)
map1 <- annHeatmap2(
  correlations,
  dendrogram = mydendro,
  annotation = myannot,
  cluster = myclust,
  labels = mylabs,
  ## The following controls if the picture is symmetric
  scale = "none",
  col = hmcols)
## Warning in breakColors(breaks, col): more colors than classes: ignoring 29 last
## colors
pp(file = "images/dendro_heatmap.png", image = map1, height = 20, width = 20)
## annotated Heatmap
## 
## Rows: 'dendrogram' with 2 branches and 99 members total, at height 6.095 
##   11  annotation variable(s)
## Cols: 'dendrogram' with 2 branches and 99 members total, at height 6.095 
##   10  annotation variable(s)

Print the larger heatmap so that all the labels appear. Keep in mind that as we get more samples, this image needs to continue getting bigger.

big heatmap

xref_prop <- table(pheno_snps[["conditions"]])
## Error in eval(quote(list(...)), env): object 'pheno_snps' not found
pheno_snps$conditions
## Error in eval(expr, envir, enclos): object 'pheno_snps' not found
idx_tbl <- exprs(pheno_snps) > 5
## Error in h(simpleError(msg, call)): error in evaluating the argument 'object' in selecting a method for function 'exprs': object 'pheno_snps' not found
new_tbl <- data.frame(row.names = rownames(exprs(pheno_snps)))
## Error in h(simpleError(msg, call)): error in evaluating the argument 'x' in selecting a method for function 'rownames': error in evaluating the argument 'object' in selecting a method for function 'exprs': object 'pheno_snps' not found
for (n in names(xref_prop)) {
  new_tbl[[n]] <- 0
  idx_cols <- which(pheno_snps[["conditions"]] == n)
  prop_col <- rowSums(idx_tbl[, idx_cols]) / xref_prop[n]
  new_tbl[n] <- prop_col
}
## Error in eval(expr, envir, enclos): object 'xref_prop' not found
keepers <- grepl(x = rownames(new_tbl), pattern = "LpaL13")
## Error in h(simpleError(msg, call)): error in evaluating the argument 'x' in selecting a method for function 'grepl': error in evaluating the argument 'x' in selecting a method for function 'rownames': object 'new_tbl' not found
new_tbl <- new_tbl[keepers, ]
## Error in eval(expr, envir, enclos): object 'new_tbl' not found
new_tbl[["strong22"]] <- 1.001 - new_tbl[["z2.2"]]
## Error in eval(expr, envir, enclos): object 'new_tbl' not found
new_tbl[["strong23"]] <- 1.001 - new_tbl[["z2.3"]]
## Error in eval(expr, envir, enclos): object 'new_tbl' not found
s22_na <- new_tbl[["strong22"]] > 1
## Error in eval(expr, envir, enclos): object 'new_tbl' not found
new_tbl[s22_na, "strong22"] <- 1
## Error in new_tbl[s22_na, "strong22"] <- 1: object 'new_tbl' not found
s23_na <- new_tbl[["strong23"]] > 1
## Error in eval(expr, envir, enclos): object 'new_tbl' not found
new_tbl[s23_na, "strong23"] <- 1
## Error in new_tbl[s23_na, "strong23"] <- 1: object 'new_tbl' not found
new_tbl[["SNP"]] <- rownames(new_tbl)
## Error in h(simpleError(msg, call)): error in evaluating the argument 'x' in selecting a method for function 'rownames': object 'new_tbl' not found
new_tbl[["Chromosome"]] <- gsub(x = new_tbl[["SNP"]], pattern = "chr_(.*)_pos_.*", replacement = "\\1")
## Error in h(simpleError(msg, call)): error in evaluating the argument 'x' in selecting a method for function 'gsub': object 'new_tbl' not found
new_tbl[["Position"]] <- gsub(x = new_tbl[["SNP"]], pattern = ".*_pos_(\\d+)_.*", replacement = "\\1")
## Error in h(simpleError(msg, call)): error in evaluating the argument 'x' in selecting a method for function 'gsub': object 'new_tbl' not found
new_tbl <- new_tbl[, c("SNP", "Chromosome", "Position", "strong22", "strong23")]
## Error in eval(expr, envir, enclos): object 'new_tbl' not found
library(CMplot)
## Much appreciate for using CMplot.
## Full description, Bug report, Suggestion and the latest codes:
## https://github.com/YinLiLin/CMplot
simplify <- new_tbl
## Error in eval(expr, envir, enclos): object 'new_tbl' not found
simplify[["strong22"]] <- NULL
## Error in simplify[["strong22"]] <- NULL: object 'simplify' not found
CMplot(simplify, bin.size = 100000)
## Error in is.data.frame(x): object 'simplify' not found
CMplot(new_tbl, plot.type="m", multracks=TRUE, threshold = c(0.01, 0.05),
       threshold.lwd=c(1,1), threshold.col=c("black","grey"),
       amplify=TRUE, bin.size=10000,
       chr.den.col=c("darkgreen", "yellow", "red"),
       signal.col=c("red", "green", "blue"),
       signal.cex=1, file="jpg", memo="", dpi=300, file.output=TRUE, verbose=TRUE)
## Error in is.data.frame(x): object 'new_tbl' not found

SNP Density Circular Manhattan Rectangular Manhattan QQ

8.5 Try out MatrixEQTL

This tool looks a little opaque, but provides sample data with things that make sense to me and should be pretty easy to recapitulate in our data.

  1. covariates.txt: Columns are samples, rows are things from pData – the most likely ones of interest for our data would be zymodeme, sensitivity
  2. geneloc.txt: columns are ‘geneid’, ‘chr’, ‘left’, ‘right’. I guess I can assume left and right are start/stop; in which case this is trivially acquirable from fData.
  3. ge.txt: This appears to be a log(rpkm/cpm) table with rows as genes and columns as samples
  4. snpsloc.txt: columns are ‘snpid’, ‘chr’, ‘pos’
  5. snps.txt: columns are samples, rows are the ids from snsploc, values a 0,1,2. I assume 0 is identical and 1..12 are the various A->TGC T->AGC C->AGT G->ACT
## For this, let us use the 'new_snps' data structure.
## Caveat here: these need to be coerced to numbers.
my_covariates <- pData(new_snps)[, c("zymodemecategorical", "clinicalcategorical")]
for (col in colnames(my_covariates)) {
  my_covariates[[col]] <- as.numeric(as.factor(my_covariates[[col]]))
}
my_covariates <- t(my_covariates)

my_geneloc <- fData(lp_expt)[, c("gid", "chromosome", "start", "end")]
colnames(my_geneloc) <- c("geneid", "chr", "left", "right")

my_ge <- exprs(normalize_expt(lp_expt, transform = "log2", filter = TRUE, convert = "cpm"))
used_samples <- tolower(colnames(my_ge)) %in% colnames(exprs(new_snps))
my_ge <- my_ge[, used_samples]

my_snpsloc <- data.frame(rownames = rownames(exprs(new_snps)))
## Oh, caveat here: Because of the way I stored the data,
## I could have duplicate rows which presumably will make matrixEQTL sad
my_snpsloc[["chr"]] <- gsub(pattern = "^chr_(.+)_pos(.+)_ref_.*$", replacement = "\\1",
                            x = rownames(my_snpsloc))
my_snpsloc[["pos"]] <- gsub(pattern = "^chr_(.+)_pos(.+)_ref_.*$", replacement = "\\2",
                            x = rownames(my_snpsloc))
test <- duplicated(my_snpsloc)
## Each duplicated row would be another variant at that position;
## so in theory we would do a rle to number them I am guessing
## However, I do not have different variants so I think I can ignore this for the moment
## but will need to make my matrix either 0 or 1.
if (sum(test) > 0) {
  message("There are: ", sum(duplicated), " duplicated entries.")
  keep_idx <- ! test
  my_snpsloc <- my_snpsloc[keep_idx, ]
}

my_snps <- exprs(new_snps)
one_idx <- my_snps > 0
my_snps[one_idx] <- 1

## Ok, at this point I think I have all the pieces which this method wants...
## Oh, no I guess not; it actually wants the data as a set of filenames...
library(MatrixEQTL)
write.table(my_snps, "eqtl/snps.tsv", na = "NA", col.names = TRUE, row.names = TRUE, sep = "\t", quote = TRUE)
## readr::write_tsv(my_snps, "eqtl/snps.tsv", )
write.table(my_snpsloc, "eqtl/snpsloc.tsv", na = "NA", col.names = TRUE, row.names = TRUE, sep = "\t", quote = TRUE)
## readr::write_tsv(my_snpsloc, "eqtl/snpsloc.tsv")
write.table(as.data.frame(my_ge), "eqtl/ge.tsv", na = "NA", col.names = TRUE, row.names = TRUE, sep = "\t", quote = TRUE)
## readr::write_tsv(as.data.frame(my_ge), "eqtl/ge.tsv")
write.table(as.data.frame(my_geneloc), "eqtl/geneloc.tsv", na = "NA", col.names = TRUE, row.names = TRUE, sep = "\t", quote = TRUE)
## readr::write_tsv(as.data.frame(my_geneloc), "eqtl/geneloc.tsv")
write.table(as.data.frame(my_covariates), "eqtl/covariates.tsv", na = "NA", col.names = TRUE, row.names = TRUE, sep = "\t", quote = TRUE)
## readr::write_tsv(as.data.frame(my_covariates), "eqtl/covariates.tsv")

useModel = modelLINEAR # modelANOVA, modelLINEAR, or modelLINEAR_CROSS

# Genotype file name
SNP_file_name = "eqtl/snps.tsv"
snps_location_file_name = "eqtl/snpsloc.tsv"
expression_file_name = "eqtl/ge.tsv"
gene_location_file_name = "eqtl/geneloc.tsv"
covariates_file_name = "eqtl/covariates.tsv"
# Output file name
output_file_name_cis = tempfile()
output_file_name_tra = tempfile()
# Only associations significant at this level will be saved
pvOutputThreshold_cis = 0.1
pvOutputThreshold_tra = 0.1
# Error covariance matrix
# Set to numeric() for identity.
errorCovariance = numeric()
# errorCovariance = read.table("Sample_Data/errorCovariance.txt");
# Distance for local gene-SNP pairs
cisDist = 1e6
## Load genotype data
snps = SlicedData$new()
snps$fileDelimiter = "\t"      # the TAB character
snps$fileOmitCharacters = "NA" # denote missing values;
snps$fileSkipRows = 1          # one row of column labels
snps$fileSkipColumns = 1       # one column of row labels
snps$fileSliceSize = 2000      # read file in slices of 2,000 rows
snps$LoadFile(SNP_file_name)
## Load gene expression data
gene = SlicedData$new()
gene$fileDelimiter = "\t"      # the TAB character
gene$fileOmitCharacters = "NA" # denote missing values;
gene$fileSkipRows = 1          # one row of column labels
gene$fileSkipColumns = 1       # one column of row labels
gene$fileSliceSize = 2000      # read file in slices of 2,000 rows
gene$LoadFile(expression_file_name)
## Load covariates
cvrt = SlicedData$new()
cvrt$fileDelimiter = "\t"      # the TAB character
cvrt$fileOmitCharacters = "NA" # denote missing values;
cvrt$fileSkipRows = 1          # one row of column labels
cvrt$fileSkipColumns = 1       # one column of row labels
if(length(covariates_file_name) > 0) {
  cvrt$LoadFile(covariates_file_name)
}
## Run the analysis
snpspos = read.table(snps_location_file_name, header = TRUE, stringsAsFactors = FALSE)
genepos = read.table(gene_location_file_name, header = TRUE, stringsAsFactors = FALSE)

me = Matrix_eQTL_main(
    snps = snps,
    gene = gene,
    cvrt = cvrt,
    output_file_name = output_file_name_tra,
    pvOutputThreshold = pvOutputThreshold_tra,
    useModel = useModel,
    errorCovariance = errorCovariance,
    verbose = TRUE,
    output_file_name.cis = output_file_name_cis,
    pvOutputThreshold.cis = pvOutputThreshold_cis,
    snpspos = snpspos,
    genepos = genepos,
    cisDist = cisDist,
    pvalue.hist = "qqplot",
    min.pv.by.genesnp = FALSE,
    noFDRsaveMemory = FALSE);
if (!isTRUE(get0("skip_load"))) {
  pander::pander(sessionInfo())
  message(paste0("This is hpgltools commit: ", get_git_commit()))
  message(paste0("Saving to ", savefile))
  tmp <- sm(saveme(filename = savefile))
}
## If you wish to reproduce this exact build of hpgltools, invoke the following:
## > git clone http://github.com/abelew/hpgltools.git
## > git reset 28ece304dcf5e4b5adbfb86dbece61884ad2d3b9
## This is hpgltools commit: Mon Jan 10 14:19:38 2022 -0500: 28ece304dcf5e4b5adbfb86dbece61884ad2d3b9
## Saving to tmrc2_02sample_estimation_v202201.rda.xz
tmp <- loadme(filename = savefile)
LS0tCnRpdGxlOiAiVE1SQzIgQ29tcHJlaGVuc2l2ZSBEYXRhIEFuYWx5c2lzOiAyMDIyMDEiCmF1dGhvcjogImF0YiBhYmVsZXdAZ21haWwuY29tIgpkYXRlOiAiYHIgU3lzLkRhdGUoKWAiCm91dHB1dDoKIGh0bWxfZG9jdW1lbnQ6CiAgY29kZV9kb3dubG9hZDogdHJ1ZQogIGNvZGVfZm9sZGluZzogc2hvdwogIGZpZ19jYXB0aW9uOiB0cnVlCiAgZmlnX2hlaWdodDogNwogIGZpZ193aWR0aDogNwogIGhpZ2hsaWdodDogZGVmYXVsdAogIGtlZXBfbWQ6IGZhbHNlCiAgbW9kZTogc2VsZmNvbnRhaW5lZAogIG51bWJlcl9zZWN0aW9uczogdHJ1ZQogIHNlbGZfY29udGFpbmVkOiB0cnVlCiAgdGhlbWU6IHJlYWRhYmxlCiAgdG9jOiB0cnVlCiAgdG9jX2Zsb2F0OgogICBjb2xsYXBzZWQ6IGZhbHNlCiAgIHNtb290aF9zY3JvbGw6IGZhbHNlCi0tLQoKPHN0eWxlPgogIGJvZHkgLm1haW4tY29udGFpbmVyIHsKICAgIG1heC13aWR0aDogMTYwMHB4OwogIH0KPC9zdHlsZT4KCmBgYHtyIG9wdGlvbnMsIGluY2x1ZGUgPSBGQUxTRX0KbGlicmFyeShocGdsdG9vbHMpCnR0IDwtIHNtKGRldnRvb2xzOjpsb2FkX2FsbCgifi9ocGdsdG9vbHMiKSkKa25pdHI6Om9wdHNfa25pdCRzZXQocHJvZ3Jlc3MgPSBUUlVFLAogICAgICAgICAgICAgICAgICAgICB2ZXJib3NlID0gVFJVRSwKICAgICAgICAgICAgICAgICAgICAgd2lkdGggPSA5MCwKICAgICAgICAgICAgICAgICAgICAgZWNobyA9IFRSVUUpCmtuaXRyOjpvcHRzX2NodW5rJHNldChlcnJvciA9IFRSVUUsCiAgICAgICAgICAgICAgICAgICAgICBmaWcud2lkdGggPSA4LAogICAgICAgICAgICAgICAgICAgICAgZmlnLmhlaWdodCA9IDgsCiAgICAgICAgICAgICAgICAgICAgICBkcGkgPSA5NikKb2xkX29wdGlvbnMgPC0gb3B0aW9ucyhkaWdpdHMgPSA0LAogICAgICAgICAgICAgICAgICAgICAgIHN0cmluZ3NBc0ZhY3RvcnMgPSBGQUxTRSwKICAgICAgICAgICAgICAgICAgICAgICBrbml0ci5kdXBsaWNhdGUubGFiZWwgPSAiYWxsb3ciKQpnZ3Bsb3QyOjp0aGVtZV9zZXQoZ2dwbG90Mjo6dGhlbWVfYncoYmFzZV9zaXplID0gMTIpKQp2ZXIgPC0gIjIwMjIwMSIKcnVuZGF0ZSA8LSBmb3JtYXQoU3lzLkRhdGUoKSwgZm9ybWF0ID0gIiVZJW0lZCIpCgojIyB0bXAgPC0gdHJ5KHNtKGxvYWRtZShmaWxlbmFtZSA9IGdzdWIocGF0dGVybiA9ICJcXC5SbWQiLCByZXBsYWNlID0gIlxcLnJkYVxcLnh6IiwgeCA9IHByZXZpb3VzX2ZpbGUpKSkpCnJtZF9maWxlIDwtIGdsdWU6OmdsdWUoInRtcmMyXzAyc2FtcGxlX2VzdGltYXRpb25fdnt2ZXJ9LlJtZCIpCnNhdmVmaWxlIDwtIGdzdWIocGF0dGVybiA9ICJcXC5SbWQiLCByZXBsYWNlID0gIlxcLnJkYVxcLnh6IiwgeCA9IHJtZF9maWxlKQoKbGlicmFyeShIZWF0cGx1cykKYGBgCgpgYGB7ciBjdXJyZW50X3NhbXBsZXNoZWV0fQpzYW1wbGVfc2hlZXQgPC0gZ2x1ZTo6Z2x1ZSgic2FtcGxlX3NoZWV0cy90bXJjMl9zYW1wbGVzXzIwMjIwMTA2Lnhsc3giKQpgYGAKCiMgSW50cm9kdWN0aW9uCgpUaGlzIGlzIG1vc3RseSBqdXN0IGEgcnVuIG9mIHRoaXMgd29ya3NoZWV0IHRvIHJlYWNxdWFpbnQgbXlzZWxmIHdpdGggaXQuCgpUaGlzIGRvY3VtZW50IGlzIGludGVuZGVkIHRvIHByb3ZpZGUgYSBnZW5lcmFsIG92ZXJ2aWV3IG9mIHRoZSBUTVJDMiBzYW1wbGVzCndoaWNoIGhhdmUgdGh1cyBmYXIgYmVlbiBzZXF1ZW5jZWQuICBJbiBzb21lIGNhc2VzLCB0aGlzIGluY2x1ZGVzIG9ubHkgdGhvc2UKc2FtcGxlcyBzdGFydGluZyBpbiAyMDE5OyBpbiBvdGhlciBpbnN0YW5jZXMgSSBhbSBpbmNsdWRpbmcgb3VyIHByZXZpb3VzCigyMDE1LTIwMTYpIHNhbXBsZXMuCgpJbiBhbGwgY2FzZXMgdGhlIHByb2Nlc3NpbmcgcGVyZm9ybWVkIHdhczoKCjEuICBEZWZhdWx0IHRyaW1taW5nIHdhcyBwZXJmb3JtZWQuCjIuICBIaXNhdDIgd2FzIHVzZWQgdG8gbWFwIHRoZSByZW1haW5pbmcgcmVhZHMgYWdhaW5zdCB0aGUgTGVpc2htYW5pYQogICAgcGFuYW1lbnNpcyBnZW5vbWUgcmV2aXNpb24gMzYuCjMuICBUaGUgYWxpZ25tZW50cyBmcm9tIGhpc2F0MiB3ZXJlIHVzZWQgdG8gY291bnQgcmVhZHMvZ2VuZSBhZ2FpbnN0IHRoZQogICAgcmV2aXNpb24gMzYgYW5ub3RhdGlvbnMgd2l0aCBodHNlcS4KNC4gIFRoZXNlIGFsaWdubWVudHMgd2VyZSBhbHNvIHBhc3NlZCB0byB0aGUgcGlsZXVwIGZ1bmN0aW9uYWxpdHkgb2Ygc2FtdG9vbHMKICAgIGFuZCB0aGUgdmNmL2JjZiB1dGlsaXRpZXMgaW4gb3JkZXIgdG8gbWFrZSBhIG1hdHJpeCBvZiBhbGwgb2JzZXJ2ZWQKICAgIGRpZmZlcmVuY2VzIGJldHdlZW4gZWFjaCBzYW1wbGUgd2l0aCByZXNwZWN0IHRvIHRoZSByZWZlcmVuY2UuCgpUaGUgYW5hbHlzZXMgaW4gdGhpcyBkb2N1bWVudCB1c2UgdGhlIG1hdHJpY2VzIG9mIGNvdW50cy9nZW5lIGZyb20gIzMgYW5kCnZhcmlhbnRzL3Bvc2l0aW9uIGZyb20gIzQgaW4gb3JkZXIgdG8gcHJvdmlkZSBzb21lIGltYWdlcyBhbmQgbWV0cmljcyBkZXNjcmliaW5nCnRoZSBzYW1wbGVzIHdlIGhhdmUgc2VxdWVuY2VkIHNvIGZhci4KCiMgQW5ub3RhdGlvbnMKCkV2ZXJ5dGhpbmcgd2hpY2ggZm9sbG93cyBkZXBlbmRzIG9uIHRoZSBFeGlzdGluZyBUcmlUcnlwREIgYW5ub3RhdGlvbnMgcmV2aXNpb24KNDYsIGNpcmNhIDIwMTkuICBUaGUgZm9sbG93aW5nIGJsb2NrIGxvYWRzIGEgZGF0YWJhc2Ugb2YgdGhlc2UgYW5ub3RhdGlvbnMgYW5kCnR1cm5zIGl0IGludG8gYSBtYXRyaXggd2hlcmUgdGhlIHJvd3MgYXJlIGdlbmVzIGFuZCBjb2x1bW5zIGFyZSBhbGwgdGhlCmFubm90YXRpb24gdHlwZXMgcHJvdmlkZWQgYnkgVHJpVHJ5cERCLgoKVGhlIHNhbWUgZGF0YWJhc2Ugd2FzIHVzZWQgdG8gY3JlYXRlIGEgbWF0cml4IG9mIG9ydGhvbG9nb3VzIGdlbmVzIGJldHdlZW4KTC5wYW5hbWVuc2lzIGFuZCBhbGwgb2YgdGhlIG90aGVyIHNwZWNpZXMgaW4gdGhlIFRyaVRyeXBEQi4KCmBgYHtyIGFubm90fQp0dCA8LSBzbShsaWJyYXJ5KEV1UGF0aERCKSkKdHQgPC0gc20obGlicmFyeShvcmcuTHBhbmFtZW5zaXMuTUhPTUNPTDgxTDEzLnY0Ni5lZy5kYikpCnBhbl9kYiA8LSBvcmcuTHBhbmFtZW5zaXMuTUhPTUNPTDgxTDEzLnY0Ni5lZy5kYgphbGxfZmllbGRzIDwtIGNvbHVtbnMocGFuX2RiKQoKYWxsX2xwX2Fubm90IDwtIHNtKGxvYWRfb3JnZGJfYW5ub3RhdGlvbnMoCiAgICBwYW5fZGIsCiAgICBrZXl0eXBlID0gImdpZCIsCiAgICBmaWVsZHMgPSBjKCJhbm5vdF9nZW5lX2VudHJlel9pZCIsICJhbm5vdF9nZW5lX25hbWUiLAogICAgICAgICAgICAgICAiYW5ub3Rfc3RyYW5kIiwgImFubm90X2Nocm9tb3NvbWUiLCAiYW5ub3RfY2RzX2xlbmd0aCIsCiAgICAgICAgICAgICAgICJhbm5vdF9nZW5lX3Byb2R1Y3QiKSkpJGdlbmVzCgpscF9nbyA8LSBzbShsb2FkX29yZ2RiX2dvKHBhbl9kYikpCmxwX2xlbmd0aHMgPC0gYWxsX2xwX2Fubm90WywgYygiZ2lkIiwgImFubm90X2Nkc19sZW5ndGgiKV0KY29sbmFtZXMobHBfbGVuZ3RocykgIDwtIGMoIklEIiwgImxlbmd0aCIpCmFsbF9scF9hbm5vdFtbImFubm90X2dlbmVfcHJvZHVjdCJdXSA8LSB0b2xvd2VyKGFsbF9scF9hbm5vdFtbImFubm90X2dlbmVfcHJvZHVjdCJdXSkKb3J0aG9zIDwtIHNtKEV1UGF0aERCOjpleHRyYWN0X2V1cGF0aF9vcnRob2xvZ3MoZGIgPSBwYW5fZGIpKQoKaGlzYXRfYW5ub3QgPC0gYWxsX2xwX2Fubm90CiMjIHJvd25hbWVzKGhpc2F0X2Fubm90KSA8LSBwYXN0ZTAoImV4b25fIiwgcm93bmFtZXMoaGlzYXRfYW5ub3QpLCAiLkUxIikKYGBgCgojIExvYWQgYSBnZW5vbWUKCmBgYHtyIGdlbm9tZX0KbWV0YSA8LSBFdVBhdGhEQjo6ZG93bmxvYWRfZXVwYXRoX21ldGFkYXRhKHdlYnNlcnZpY2U9InRyaXRyeXBkYiIpCmxwX2VudHJ5IDwtIEV1UGF0aERCOjpnZXRfZXVwYXRoX2VudHJ5KHNwZWNpZXM9IkxlaXNobWFuaWEgcGFuYW1lbnNpcyIsIG1ldGFkYXRhPW1ldGEpCmNvbG5hbWVzKGxwX2VudHJ5KQp0ZXN0aW5nX3BhbmFtZW5zaXMgPC0gIkJTR2Vub21lLkxlaXNobWFuaWEucGFuYW1lbnNpcy5NSE9NQ09MODFMMTMudjUzIgojIyB0ZXN0aW5nX3BhbmFtZW5zaXMgPC0gRXVQYXRoREI6Om1ha2VfZXVwYXRoX2JzZ2Vub21lKGVudHJ5PWxwX2VudHJ5LCBldV92ZXJzaW9uPSJ2NDYiKQpsaWJyYXJ5KGFzLmNoYXJhY3Rlcih0ZXN0aW5nX3BhbmFtZW5zaXMpLCBjaGFyYWN0ZXIub25seT1UUlVFKQpnZW5vbWUgPC0gZ2V0MChhcy5jaGFyYWN0ZXIodGVzdGluZ19wYW5hbWVuc2lzKSkKYGBgCgojIFRPRE86CgpSZXNlcXVlbmNlIHNhbXBsZXM6IFRNUkMyMDAwMiwgVE1SQzIwMDA2LCBUTVJDMjAwMDQgKG1heWJlIFRNUkMyMDAwOCBhbmQgVE1SQzIwMDI5KQoKIyBHZW5lcmF0ZSBFeHByZXNzaW9uc2V0cyBhbmQgU2FtcGxlIEVzdGltYXRpb24KClRoZSBwcm9jZXNzIG9mIHNhbXBsZSBlc3RpbWF0aW9uIHRha2VzIHR3byBwcmltYXJ5IGlucHV0czoKCjEuICBUaGUgc2FtcGxlIHNoZWV0LCB3aGljaCBjb250YWlucyBhbGwgdGhlIG1ldGFkYXRhIHdlIGN1cnJlbnRseSBoYXZlIG9uIGhhbmQsCiAgICBpbmNsdWRpbmcgZmlsZW5hbWVzIGZvciB0aGUgb3V0cHV0cyBvZiAjMyBhbmQgIzQgYWJvdmUuCjIuICBUaGUgZ2VuZSBhbm5vdGF0aW9ucy4KCkFuIGV4cHJlc3Npb25zZXQgaXMgYSBkYXRhIHN0cnVjdHVyZSB1c2VkIGluIFIgdG8gZXhhbWluZSBSTkFTZXEgZGF0YS4gIEl0CmlzIGNvbXByaXNlZCBvZiBhbm5vdGF0aW9ucywgbWV0YWRhdGEsIGFuZCBleHByZXNzaW9uIGRhdGEuICBJbiB0aGUgY2FzZSBvZiBvdXIKcHJvY2Vzc2luZyBwaXBlbGluZSwgdGhlIGxvY2F0aW9uIG9mIHRoZSBleHByZXNzaW9uIGRhdGEgaXMgcHJvdmlkZWQgYnkgdGhlCmZpbGVuYW1lcyBpbiB0aGUgbWV0YWRhdGEuCgpUaGUgZmlyc3QgbGluZXMgb2YgdGhlIGZvbGxvd2luZyBibG9jayBjcmVhdGUgdGhlIEV4cHJlc3Npb25zZXQuICBBbGwgb2YgdGhlCmZvbGxvd2luZyBsaW5lcyBwZXJmb3JtIHZhcmlvdXMgbm9ybWFsaXphdGlvbnMgYW5kIGdlbmVyYXRlIHBsb3RzIGZyb20gaXQuCgojIyBOb3RlcwoKVGhlIGZvbGxvd2luZyBzYW1wbGVzIGFyZSBtdWNoIGxvd2VyIGNvdmVyYWdlOgoKKiBUTVJDMjAwMDIKKiBUTVJDMjAwMDYKKiBUTVJDMjAwMDcKKiBUTVJDMjAwMDgKCjIwMjEwNjEwOiBJIG1hZGUgc29tZSBtYW51YWwgY2hhbmdlcyB0byB0aGUgc2FtcGxlIHNoZWV0IHdoaWNoIEkKZG93bmxvYWRlZCwgZmlsbGluZyBpbiBzb21lIHp5bW9kZW1lIHdpdGggJ3Vua25vd24nCgojIyBUT0RPOgoKMS4gIERvIHRoZSBtdWx0aS1nZW5lIGZhbWlseSByZW1vdmFsIHJpZ2h0IGhlcmUgaW5zdGVhZCBvZiB3YXkgZG93biBhdCB0aGUgYm90dG9tCjIuICBBZGQgenltb2RlbWUgc25wcyB0byB0aGUgYW5ub3RhdGlvbiBsYXRlci4KMy4gIFN0YXJ0IHBoeWxvZ2VuZXRpYyBhbmFseXNpcyBvZiB2YXJpYW50IHRhYmxlLgoKYGBge3IgbmV3X3NhbXBsZXNfaGlzYXR9CnNhbml0aXplX2NvbHVtbnMgPC0gYygicGFzc2FnZW51bWJlciIsICJjbGluaWNhbHJlc3BvbnNlIiwgImNsaW5pY2FsY2F0ZWdvcmljYWwiLAogICAgICAgICAgICAgICAgICAgICAgInp5bW9kZW1lY2F0ZWdvcmljYWwiLCAienltb2RlbWVjYXRlZ29yaWNhbCIpCmxwX2V4cHQgPC0gc20oY3JlYXRlX2V4cHQoc2FtcGxlX3NoZWV0LAogICAgICAgICAgICAgICAgICAgICAgICAgIGdlbmVfaW5mbyA9IGhpc2F0X2Fubm90LAogICAgICAgICAgICAgICAgICAgICAgICAgIGlkX2NvbHVtbiA9ICJocGdsaWRlbnRpZmllciIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgZmlsZV9jb2x1bW4gPSAibHBhbmFtZW5zaXN2MzZoaXNhdGZpbGUiKSkgJT4lCiAgc2V0X2V4cHRfY29uZGl0aW9ucyhmYWN0ID0gInp5bW9kZW1lY2F0ZWdvcmljYWwiKSAlPiUKICBzdWJzZXRfZXhwdChub256ZXJvID0gODU1MCkgJT4lCiAgc3Vic2V0X2V4cHQoY292ZXJhZ2UgPSA1MDAwMDAwKSAlPiUKICBzZW1hbnRpY19leHB0X2ZpbHRlcihzZW1hbnRpYyA9IGMoImFtYXN0aW4iLCAiZ3A2MyIsICJsZWlzaG1hbm9seXNpbiIpLAogICAgICAgICAgICAgICAgICAgICAgIHNlbWFudGljX2NvbHVtbiA9ICJhbm5vdF9nZW5lX3Byb2R1Y3QiKSAlPiUKICBzYW5pdGl6ZV9leHB0X21ldGFkYXRhKGNvbHVtbnMgPSBzYW5pdGl6ZV9jb2x1bW5zKSAlPiUKICBzZXRfZXhwdF9mYWN0b3JzKGNvbHVtbnMgPSBzYW5pdGl6ZV9jb2x1bW5zLCBjbGFzcyA9ICJmYWN0b3IiKQoKbGlic2l6ZXMgPC0gcGxvdF9saWJzaXplKGxwX2V4cHQpCnBwKGZpbGUgPSAiaW1hZ2VzL2xwX2V4cHRfbGlic2l6ZXMucG5nIiwgaW1hZ2UgPSBsaWJzaXplcyRwbG90LCB3aWR0aCA9IDE0LCBoZWlnaHQgPSA5KQojIyBJIHRoaW5rIHNhbXBsZXMgNywxMCBzaG91bGQgYmUgcmVtb3ZlZCBhdCBtaW5pbXVtLCBwcm9iYWJseSBhbHNvIDksMTEKbm9uemVybyA8LSBwbG90X25vbnplcm8obHBfZXhwdCkKcHAoZmlsZSA9ICJpbWFnZXMvbHBfbm9uemVyby5wbmciLCBpbWFnZSA9IG5vbnplcm8kcGxvdCwgd2lkdGggPSA5LCBoZWlnaHQgPSA5KQoKbHBfYm94IDwtIHBsb3RfYm94cGxvdChscF9leHB0KQpwcChmaWxlID0gImltYWdlcy9scF9leHB0X2JveHBsb3QucG5nIiwgaW1hZ2UgPSBscF9ib3gsIHdpZHRoID0gMTIsIGhlaWdodCA9IDkpCgpmaWx0ZXJfcGxvdCA8LSBwbG90X2xpYnNpemVfcHJlcG9zdChscF9leHB0KQpmaWx0ZXJfcGxvdCRsb3dnZW5lX3Bsb3QKZmlsdGVyX3Bsb3QkY291bnRfcGxvdApgYGAKCiMjIERpc3RyaWJ1dGlvbiBWaXN1YWxpemF0aW9uCgpOYWppYidzIGZhdm9yaXRlIHBsb3RzIGFyZSBvZiBjb3Vyc2UgdGhlIFBDQS9UTlNFLiAgVGhlc2UgYXJlIG5pY2UgdG8gbG9vayBhdCBpbgpvcmRlciB0byBnZXQgYSBzZW5zZSBvZiB0aGUgcmVsYXRpb25zaGlwcyBiZXR3ZWVuIHNhbXBsZXMuICBUaGV5IGFsc28gcHJvdmlkZSBhCmdvb2Qgb3Bwb3J0dW5pdHkgdG8gc2VlIHdoYXQgaGFwcGVucyB3aGVuIG9uZSBhcHBsaWVzIGRpZmZlcmVudCBub3JtYWxpemF0aW9ucywKc3Vycm9nYXRlIGFuYWx5c2VzLCBmaWx0ZXJzLCBldGMuICBJbiBhZGRpdGlvbiwgb25lIG1heSBzZXQgZGlmZmVyZW50CmV4cGVyaW1lbnRhbCBmYWN0b3JzIGFzIHRoZSBwcmltYXJ5ICdjb25kaXRpb24nICh1c3VhbGx5IHRoZSBjb2xvciBvZiBwbG90cykgYW5kCnN1cnJvZ2F0ZSAnYmF0Y2hlcycuCgojIyBCeSBTdXNjZXB0aWxpYml0eQoKQ29sdW1uICdRJyBpbiB0aGUgc2FtcGxlIHNoZWV0LCBtYWtlIGEgY2F0ZWdvcmljYWwgdmVyc2lvbiBvZiBpdCB3aXRoIHRoZXNlIHBhcmFtZXRlcnM6CgoqIDAgPD0geCA8PSAzNSBpcyByZXNpc3RhbnQKKiAzNiA8PSB4IDw9IDQ4IGlzIGFtYmlndW91cwoqIDQ5IDw9IHggaXMgc2Vuc2l0aXZlCgpgYGB7ciBzdXNjZXB0aWJpbGl0eX0Kc3RhcnRpbmcgPC0gYXMubnVtZXJpYyhwRGF0YShscF9leHB0KVtbInN1c2NlcHRpYmlsaXR5aW5mZWN0aW9ucmVkdWN0aW9uMzJ1Z21sc2J2aGlzdG9yaWNhbGRhdGEiXV0pCnN1c19jYXRlZ29yaWNhbCA8LSBzdGFydGluZwpuYV9pZHggPC0gaXMubmEoc3RhcnRpbmcpCnN1c19jYXRlZ29yaWNhbFtuYV9pZHhdIDwtICJ1bmtub3duIgoKcmVzaXN0X2lkeCA8LSBzdGFydGluZyA8PSAwLjM1CnN1c19jYXRlZ29yaWNhbFtyZXNpc3RfaWR4XSA8LSAicmVzaXN0YW50IgppbmRldGVybWluYW50X2lkeCA8LSBzdGFydGluZyA+PSAwLjM2ICYgc3RhcnRpbmcgPD0gMC40OApzdXNfY2F0ZWdvcmljYWxbaW5kZXRlcm1pbmFudF9pZHhdIDwtICJhbWJpZ3VvdXMiCnN1c2NlcHRpYmxlX2lkeCA8LSBzdGFydGluZyA+PSAwLjQ5CnN1c19jYXRlZ29yaWNhbFtzdXNjZXB0aWJsZV9pZHhdIDwtICJzZW5zaXRpdmUiCgpwRGF0YShscF9leHB0JGV4cHJlc3Npb25zZXQpW1sic3VzX2NhdGVnb3J5Il1dIDwtIHN1c19jYXRlZ29yaWNhbApgYGAKCmBgYHtyIHByZV9xdWVzdGlvbnN9CmNsaW5pY2FsX2NvbG9ycyA8LSBsaXN0KAogICAgIyMgICAgInoyLjEiID0gIiMwMDAwY2MiLAogICAgIyMgICAgInoyLjMiID0gIiM4NzQ0MDAiLAogICAgIyMgICAgInoyLjIiID0gIiNkZjcwMDAiLAogICAgIyMgICAgInoyLjQiID0gIiNjYzAwMDAiLAogICAgInoyLjEiID0gIiM4NzQ0MDAiLAogICAgInoyLjIiID0gIiMwMDAwY2MiLAogICAgInoyLjMiID0gIiNjYzAwMDAiLAogICAgInoyLjQiID0gIiNkZjcwMDAiLAogICAgInVua25vd24iID0gIiNjYmNiY2IiLAogICAgIm51bGwiID0gIiMwMDAwMDAiKQpjbGluaWNhbF9zYW1wbGVzIDwtIGxwX2V4cHQgJT4lCiAgc2V0X2V4cHRfYmF0Y2hlcyhmYWN0ID0gc3VzX2NhdGVnb3JpY2FsKSAlPiUKICBzZXRfZXhwdF9jb2xvcnMoY2xpbmljYWxfY29sb3JzKQoKY2xpbmljYWxfbm9ybSA8LSBzbShub3JtYWxpemVfZXhwdChjbGluaWNhbF9zYW1wbGVzLCBub3JtID0gInF1YW50IiwgdHJhbnNmb3JtID0gImxvZzIiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGNvbnZlcnQgPSAiY3BtIiwgYmF0Y2ggPSBGQUxTRSwgZmlsdGVyID0gVFJVRSkpCnp5bW9fcGNhIDwtIHBsb3RfcGNhKGNsaW5pY2FsX25vcm0sIHBsb3RfdGl0bGUgPSAiUENBIG9mIHBhcmFzaXRlIGV4cHJlc3Npb24gdmFsdWVzIiwKICAgICAgICAgICAgICAgICAgICAgcGxvdF9sYWJlbHMgPSBGQUxTRSkKcHAoZmlsZSA9ICJpbWFnZXMvenltb19wY2Ffc3VzX3NoYXBlLnBuZyIsIGltYWdlID0genltb19wY2EkcGxvdCkKCm9ubHlfdHdvX3R5cGVzIDwtIHN1YnNldF9leHB0KGNsaW5pY2FsX3NhbXBsZXMsIHN1YnNldCA9ICJjb25kaXRpb249PSd6Mi4zJ3xjb25kaXRpb249PSd6Mi4yJyIpCm9ubHlfdHdvX25vcm0gPC0gc20obm9ybWFsaXplX2V4cHQob25seV90d29fdHlwZXMsIG5vcm0gPSAicXVhbnQiLCB0cmFuc2Zvcm0gPSAibG9nMiIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgY29udmVydCA9ICJjcG0iLCBiYXRjaCA9IEZBTFNFLCBmaWx0ZXIgPSBUUlVFKSkKb25seXR3b19wY2EgPC0gcGxvdF9wY2Eob25seV90d29fbm9ybSwgcGxvdF90aXRsZSA9ICJQQ0Egb2YgejIuMiBhbmQgejIuMyBwYXJhc2l0ZSBleHByZXNzaW9uIHZhbHVlcyIsCiAgICAgICAgICAgICAgICAgICAgIHBsb3RfbGFiZWxzID0gRkFMU0UpCnBwKGZpbGUgPSAiaW1hZ2VzL3p5bW9fejIuMl96Mi4zX3BjYV9zdXNfc2hhcGUucG5nIiwgaW1hZ2UgPSBvbmx5dHdvX3BjYSRwbG90KQoKenltb18zZHBjYSA8LSBwbG90XzNkX3BjYSh6eW1vX3BjYSkKenltb18zZHBjYSRwbG90CgpjbGluaWNhbF9uIDwtIHNtKG5vcm1hbGl6ZV9leHB0KGNsaW5pY2FsX3NhbXBsZXMsIHRyYW5zZm9ybSA9ICJsb2cyIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBjb252ZXJ0ID0gImNwbSIsIGJhdGNoID0gRkFMU0UsIGZpbHRlciA9IFRSVUUpKQp6eW1vX3RzbmUgPC0gcGxvdF90c25lKGNsaW5pY2FsX24sIHBsb3RfdGl0bGUgPSAiVFNORSBvZiBwYXJhc2l0ZSBleHByZXNzaW9uIHZhbHVlcyIpCnp5bW9fdHNuZSRwbG90CgpjbGluaWNhbF9uYiA8LSBub3JtYWxpemVfZXhwdChjbGluaWNhbF9zYW1wbGVzLCBjb252ZXJ0ID0gImNwbSIsIHRyYW5zZm9ybSA9ICJsb2cyIiwKICAgICAgICAgICAgICAgICAgICAgICAgIGZpbHRlciA9IFRSVUUsIGJhdGNoID0gInN2YXNlcSIpCmNsaW5pY2FsX25iX3BjYSA8LSBwbG90X3BjYShjbGluaWNhbF9uYiwgcGxvdF90aXRsZSA9ICJQQ0Egb2YgcGFyYXNpdGUgZXhwcmVzc2lvbiB2YWx1ZXMiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgcGxvdF9sYWJlbHMgPSBGQUxTRSkKcHAoZmlsZSA9ICJpbWFnZXMvY2xpbmljYWxfbmJfcGNhX3N1c19zaGFwZS5wbmciLCBpbWFnZSA9IGNsaW5pY2FsX25iX3BjYSRwbG90KQoKY2xpbmljYWxfbmJfdHNuZSA8LSBwbG90X3RzbmUoY2xpbmljYWxfbmIsIHBsb3RfdGl0bGUgPSAiVFNORSBvZiBwYXJhc2l0ZSBleHByZXNzaW9uIHZhbHVlcyIpCmNsaW5pY2FsX25iX3RzbmUkcGxvdAoKY29yaGVhdCA8LSBwbG90X2NvcmhlYXQoY2xpbmljYWxfbm9ybSwgcGxvdF90aXRsZSA9ICJDb3JyZWxhdGlvbiBoZWF0bWFwIG9mIHBhcmFzaXRlCiAgICAgICAgICAgICAgICAgZXhwcmVzc2lvbiB2YWx1ZXMKIikKY29yaGVhdCRwbG90CgpwbG90X3NtKGNsaW5pY2FsX25vcm0pJHBsb3QKYGBgCgojIyBCeSBDdXJlL0ZhaWwgc3RhdHVzCgpgYGB7ciBjZl9zdGF0dXN9CmNmX2NvbG9ycyA8LSBsaXN0KAogICAgImN1cmUiID0gIiMwMDZmMDAiLAogICAgImZhaWwiID0gIiM5ZGZmYTAiLAogICAgInVua25vd24iID0gIiNjYmNiY2IiLAogICAgIm5vdGFwcGxpY2FibGUiID0gIiMwMDAwMDAiKQpjZl9leHB0IDwtIHNldF9leHB0X2NvbmRpdGlvbnMobHBfZXhwdCwgZmFjdCA9ICJjbGluaWNhbGNhdGVnb3JpY2FsIikgJT4lCiAgc2V0X2V4cHRfYmF0Y2hlcyhmYWN0ID0gc3VzX2NhdGVnb3JpY2FsKSAlPiUKICBzZXRfZXhwdF9jb2xvcnMoY2ZfY29sb3JzKQoKY2Zfbm9ybSA8LSBub3JtYWxpemVfZXhwdChjZl9leHB0LCBjb252ZXJ0ID0gImNwbSIsIHRyYW5zZm9ybSA9ICJsb2cyIiwKICAgICAgICAgICAgICAgICAgICAgICAgICBub3JtID0gInF1YW50IiwgZmlsdGVyID0gVFJVRSkKc3RhcnRfY2YgPC0gcGxvdF9wY2EoY2Zfbm9ybSwgcGxvdF90aXRsZSA9ICJQQ0Egb2YgcGFyYXNpdGUgZXhwcmVzc2lvbiB2YWx1ZXMiLAogICAgICAgICAgICAgICAgICAgICBwbG90X2xhYmVscyA9IEZBTFNFKQpwcChmaWxlID0gImltYWdlcy9jZl9zdXNfc2hhcGUucG5nIiwgaW1hZ2UgPSBzdGFydF9jZiRwbG90KQoKY2ZfbmIgPC0gbm9ybWFsaXplX2V4cHQoY2ZfZXhwdCwgY29udmVydCA9ICJjcG0iLCB0cmFuc2Zvcm0gPSAibG9nMiIsCiAgICAgICAgICAgICAgICAgICAgICAgIG5vcm0gPSAicXVhbnQiLCBmaWx0ZXIgPSBUUlVFLCBiYXRjaCA9ICJzdmFzZXEiKQpjZl9uYl9wY2EgPC0gcGxvdF9wY2EoY2ZfbmIsIHBsb3RfdGl0bGUgPSAiUENBIG9mIHBhcmFzaXRlIGV4cHJlc3Npb24gdmFsdWVzIiwKICAgICAgICAgICAgICAgICAgICAgIHBsb3RfbGFiZWxzID0gRkFMU0UpCnBwKGZpbGUgPSAiaW1hZ2VzL2NmX3N1c19zaGFyZV9uYi5wbmciLCBpbWFnZSA9IGNmX25iX3BjYSRwbG90KQoKY2Zfbm9ybSA8LSBub3JtYWxpemVfZXhwdChjZl9leHB0LCB0cmFuc2Zvcm0gPSAibG9nMiIsIGNvbnZlcnQgPSAiY3BtIiwKICAgICAgICAgICAgICAgICAgICAgICAgICBmaWx0ZXIgPSBUUlVFLCBub3JtID0gInF1YW50IikKCnRlc3QgPC0gcGNhX2luZm9ybWF0aW9uKGNmX25vcm0sCiAgICAgICAgICAgICAgICAgICAgICAgIGV4cHRfZmFjdG9ycyA9IGMoImNsaW5pY2FsY2F0ZWdvcmljYWwiLCAienltb2RlbWVjYXRlZ29yaWNhbCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgInBhdGhvZ2Vuc3RyYWluIiwgInBhc3NhZ2VudW1iZXIiKSwKICAgICAgICAgICAgICAgICAgICAgICAgbnVtX2NvbXBvbmVudHMgPSA2LCBwbG90X3BjYXMgPSBUUlVFKQp0ZXN0JGFub3ZhX3AKdGVzdCRjb3JfaGVhdG1hcApgYGAKCmBgYHtyIHN1c2NlcHRpYmlsaXR5X3BjYX0Kc3VzX2NvbG9ycyA8LSBsaXN0KAogICAgInJlc2lzdGFudCIgPSAiIzg1NjNhNyIsCiAgICAic2Vuc2l0aXZlIiA9ICIjOGQwMDAwIiwKICAgICJhbWJpZ3VvdXMiID0gIiNjYmNiY2IiLAogICAgInVua25vd24iID0gIiMwMDAwMDAiKQpzdXNfZXhwdCA8LSBzZXRfZXhwdF9jb25kaXRpb25zKGxwX2V4cHQsIGZhY3QgPSAic3VzX2NhdGVnb3J5IikgJT4lCiAgc2V0X2V4cHRfYmF0Y2hlcyhmYWN0ID0gInp5bW9kZW1lY2F0ZWdvcmljYWwiKSAlPiUKICBzZXRfZXhwdF9jb2xvcnMoY29sb3JzID0gc3VzX2NvbG9ycykgJT4lCiAgc3Vic2V0X2V4cHQoc3Vic2V0ID0gImJhdGNoIT0nejI0JyIpICU+JQogIHN1YnNldF9leHB0KHN1YnNldCA9ICJiYXRjaCE9J3oyMSciKQoKc3VzX25vcm0gPC0gbm9ybWFsaXplX2V4cHQoc3VzX2V4cHQsIHRyYW5zZm9ybSA9ICJsb2cyIiwgY29udmVydCA9ICJjcG0iLAogICAgICAgICAgICAgICAgICAgICAgICAgICBub3JtID0gInF1YW50IiwgZmlsdGVyID0gVFJVRSkKc3VzX3BjYSA8LSBwbG90X3BjYShzdXNfbm9ybSwgcGxvdF90aXRsZSA9ICJQQ0Egb2YgcGFyYXNpdGUgZXhwcmVzc2lvbiB2YWx1ZXMiLAogICAgICAgICAgICAgICAgICAgIHBsb3RfbGFiZWxzID0gRkFMU0UpCnBwKGZpbGUgPSAiaW1hZ2VzL3N1c19ub3JtX3BjYS5wbmciLCBpbWFnZSA9IHN1c19wY2FbWyJwbG90Il1dKQoKc3VzX25iIDwtIG5vcm1hbGl6ZV9leHB0KHN1c19leHB0LCB0cmFuc2Zvcm0gPSAibG9nMiIsIGNvbnZlcnQgPSAiY3BtIiwKICAgICAgICAgICAgICAgICAgICAgICAgIGJhdGNoID0gInN2YXNlcSIsIGZpbHRlciA9IFRSVUUpCnN1c19uYl9wY2EgPC0gcGxvdF9wY2Eoc3VzX25iLCBwbG90X3RpdGxlID0gIlBDQSBvZiBwYXJhc2l0ZSBleHByZXNzaW9uIHZhbHVlcyIsCiAgICAgICAgICAgICAgICAgICAgICAgcGxvdF9sYWJlbHMgPSBGQUxTRSkKcHAoZmlsZSA9ICJpbWFnZXMvc3VzX25iX3BjYS5wbmciLCBpbWFnZSA9IHN1c19uYl9wY2FbWyJwbG90Il1dKQpgYGAKCiMgWnltb2RlbWUgYW5hbHlzZXMKClRoZSBmb2xsb3dpbmcgc2VjdGlvbnMgcGVyZm9ybSBhIHNlcmllcyBvZiBhbmFseXNlcyB3aGljaCBzZWVrIHRvIGVsdWNpZGF0ZQpkaWZmZXJlbmNlcyBiZXR3ZWVuIHRoZSB6eW1vZGVtZXMgMi4yIGFuZCAyLjMgZWl0aGVyIHRocm91Z2ggZGlmZmVyZW50aWFsCmV4cHJlc3Npb24gb3IgdmFyaWFudCBwcm9maWxlcy4KCiMjIERpZmZlcmVudGlhbCBleHByZXNzaW9uCgojIyMgV2l0aCByZXNwZWN0IHRvIHp5bW9kZW1lIGF0dHJpYnV0aW9uCgpUT0RPOiBEbyB0aGlzIHdpdGggYW5kIHdpdGhvdXQgc3ZhIGFuZCBjb21wYXJlIHRoZSByZXN1bHRzLgoKYGBge3Igenltb19kZSwgZmlnLnNob3cgPSAiaGlkZSJ9Cnp5X2V4cHQgPC0gc3Vic2V0X2V4cHQobHBfZXhwdCwgc3Vic2V0ID0gImNvbmRpdGlvbj09J3oyLjInfGNvbmRpdGlvbj09J3oyLjMnIikKenlfbm9ybSA8LSBub3JtYWxpemVfZXhwdCh6eV9leHB0LCBmaWx0ZXIgPSBUUlVFLCBjb252ZXJ0ID0gImNwbSIsIG5vcm0gPSAicXVhbnQiKQp6eV9kZV9ub2JhdGNoIDwtIHNtKGFsbF9wYWlyd2lzZSh6eV9leHB0LCBmaWx0ZXIgPSBUUlVFLCBtb2RlbF9iYXRjaCA9ICJzdmFzZXEiKSkKenlfZGUgPC0gc20oYWxsX3BhaXJ3aXNlKHp5X2V4cHQsIGZpbHRlciA9IFRSVUUsIG1vZGVsX2JhdGNoID0gInN2YXNlcSIpKQp6eV90YWJsZSA8LSBzbShjb21iaW5lX2RlX3RhYmxlcyh6eV9kZSwgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC96eV90YWJsZXMtdnt2ZXJ9Lnhsc3giKSkpCnp5X3NpZyA8LSBzbShleHRyYWN0X3NpZ25pZmljYW50X2dlbmVzKHp5X3RhYmxlLCBleGNlbCA9IGdsdWU6OmdsdWUoImV4Y2VsL3p5X3NpZy12e3Zlcn0ueGxzeCIpKSkKYGBgCgojIyMgSW1hZ2VzIG9mIHp5bW9kZW1lIERFCgpgYGB7ciB6eW1vZF9kZV9waWN0dXJlc30KcHAoZmlsZSA9ICJpbWFnZXMvenltb19tYS5wbmciLCBpbWFnZSA9IHp5X3RhYmxlW1sicGxvdHMiXV1bWyJ6MjNfdnNfejIyIl1dW1siZGVzZXFfbWFfcGxvdHMiXV1bWyJwbG90Il1dKQpgYGAKCiMjIFdpdGggcmVzcGVjdCB0byBjdXJlL2ZhaWx1cmUKCkluIGNvbnRyYXN0LCB3ZSBjYW4gc2VhcmNoIGZvciBnZW5lcyB3aGljaCBhcmUgZGlmZmVyZW50aWFsbHkKZXhwcmVzc2VkIHdpdGggcmVzcGVjdCB0byBjdXJlL2ZhaWx1cmUgc3RhdHVzLgoKYGBge3IgY3VyZWZhaWxfZGUsIGZpZy5zaG93ID0gImhpZGUifQpjZl9kZSA8LSBzbShhbGxfcGFpcndpc2UoY2ZfZXhwdCwgZmlsdGVyID0gVFJVRSwgbW9kZWxfYmF0Y2ggPSAic3Zhc2VxIikpCmNmX3RhYmxlIDwtIHNtKGNvbWJpbmVfZGVfdGFibGVzKGNmX2RlLCBleGNlbCA9IGdsdWU6OmdsdWUoImV4Y2VsL2NmX3RhYmxlcy12e3Zlcn0ueGxzeCIpKSkKY2Zfc2lnIDwtIHNtKGV4dHJhY3Rfc2lnbmlmaWNhbnRfZ2VuZXMoY2ZfdGFibGUsIGV4Y2VsID0gZ2x1ZTo6Z2x1ZSgiZXhjZWwvY2Zfc2lnLXZ7dmVyfS54bHN4IikpKQpgYGAKCiMjIFdpdGggcmVzcGVjdCB0byBzdXNjZXB0aWJpbGl0eQoKRmluYWxseSwgd2UgY2FuIHVzZSBvdXIgY2F0ZWdvcnkgb2Ygc3VzY2VwdGliaWxpdHkgYW5kIGxvb2sgZm9yIGdlbmVzCndoaWNoIGNoYW5nZSBmcm9tIHNlbnNpdGl2ZSB0byByZXNpc3RhbnQuICBLZWVwIGluIG1pbmQsIHRob3VnaCwgdGhhdApmb3IgdGhlIG1vbWVudCB3ZSBoYXZlIGEgbG90IG9mIGFtYmlndW91cyBhbmQgdW5rbm93biBzdHJhaW5zLgoKYGBge3IgY3VyZWZhaWxfZGUwMiwgZmlnLnNob3cgPSAiaGlkZSJ9CnN1c19kZSA8LSBzbShhbGxfcGFpcndpc2Uoc3VzX2V4cHQsIGZpbHRlciA9IFRSVUUsIG1vZGVsX2JhdGNoID0gInN2YXNlcSIpKQpzdXNfdGFibGUgPC0gc20oY29tYmluZV9kZV90YWJsZXMoc3VzX2RlLCBleGNlbCA9IGdsdWU6OmdsdWUoImV4Y2VsL3N1c190YWJsZXMtdnt2ZXJ9Lnhsc3giKSkpCnN1c19zaWcgPC0gc20oZXh0cmFjdF9zaWduaWZpY2FudF9nZW5lcyhzdXNfdGFibGUsIGV4Y2VsID0gZ2x1ZTo6Z2x1ZSgiZXhjZWwvc3VzX3NpZy12e3Zlcn0ueGxzeCIpKSkKYGBgCgpgYGB7ciB6eW1vZF9kZV9waWN0dXJlczAxfQprbml0cjo6a2FibGUoaGVhZChzdXNfc2lnJGRlc2VxJHVwcyRzZW5zaXRpdmVfdnNfcmVzaXN0YW50LCBuID0gMjApKQoKa25pdHI6OmthYmxlKGhlYWQoc3VzX3NpZyRkZXNlcSRkb3ducyRzZW5zaXRpdmVfdnNfcmVzaXN0YW50LCBuID0gMjApKQoKc3VzX21hIDwtIHN1c190YWJsZVtbInBsb3RzIl1dW1sic2Vuc2l0aXZlX3ZzX3Jlc2lzdGFudCJdXVtbImRlc2VxX21hX3Bsb3RzIl1dW1sicGxvdCJdXQpwcChmaWxlID0gImltYWdlcy9zdXNfbWEucG5nIiwgaW1hZ2UgPSBzdXNfbWEpCgojIyB0ZXN0IDwtIGdncGx0KHN1c19tYSkKYGBgCgojIyBPbnRvbG9neSBzZWFyY2hlcwoKTm93IGxldCB1cyBsb29rIGZvciBvbnRvbG9neSBjYXRlZ29yaWVzIHdoaWNoIGFyZSBpbmNyZWFzZWQgaW4gdGhlIDIuMwpzYW1wbGVzIGZvbGxvd2VkIGJ5IHRoZSAyLjIgc2FtcGxlcy4KCmBgYHtyIGdvLCBzaWcuc2hvdyA9ICJoaWRlIn0KIyMgR2VuZSBjYXRlZ29yaWVzIG1vcmUgcmVwcmVzZW50ZWQgaW4gdGhlIDIuMyBncm91cC4KenlfZ29fdXAgPC0gc20oc2ltcGxlX2dvc2VxKHNpZ19nZW5lcyA9IHp5X3NpZ1tbImRlc2VxIl1dW1sidXBzIl1dW1sxXV0sCiAgICAgICAgICAgICAgICAgICAgICAgICAgICBnb19kYiA9IGxwX2dvLCBsZW5ndGhfZGIgPSBscF9sZW5ndGhzKSkKCiMjIEdlbmUgY2F0ZWdvcmllcyBtb3JlIHJlcHJlc2VudGVkIGluIHRoZSAyLjIgZ3JvdXAuCnp5X2dvX2Rvd24gPC0gc20oc2ltcGxlX2dvc2VxKHNpZ19nZW5lcyA9IHp5X3NpZ1tbImRlc2VxIl1dW1siZG93bnMiXV1bWzFdXSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgZ29fZGIgPSBscF9nbywgbGVuZ3RoX2RiID0gbHBfbGVuZ3RocykpCmBgYAoKIyMjIEEgY291cGxlIHBsb3RzIGZyb20gdGhlIGRpZmZlcmVudGlhbCBleHByZXNzaW9uCgojIyMjIE51bWJlciBvZiBnZW5lcyBpbiBhZ3JlZW1lbnQgYW1vbmcgREUgbWV0aG9kcywgMi4zIG1vcmUgdGhhbiAyLjIKCkluIHRoZSBmdW5jdGlvbiAnY29tYmluZWRfZGVfdGFibGVzKCknIGFib3ZlLCBvbmUgb2YgdGhlIHRhc2tzCnBlcmZvcm1lZCBpcyB0byBsb29rIGF0IHRoZSBhZ3JlZW1lbnQgYW1vbmcgREVTZXEyLCBsaW1tYSwgYW5kIGVkZ2VSLgpUaGUgZm9sbG93aW5nIHNob3cgYSBjb3VwbGUgb2YgdGhlc2UgZm9yIHRoZSBzZXQgb2YgZ2VuZXMgb2JzZXJ2ZWQKd2l0aCBhIGZvbGQtY2hhbmdlID49IHwyfCBhbmQgYWRqdXN0ZWQgcC12YWx1ZSA8PSAwLjA1LgoKYGBge3IgZGVfcGxvdHN9Cnp5X3RhYmxlW1sidmVubnMiXV1bWzFdXVtbInBfbGZjMSJdXVtbInVwX25vd2VpZ2h0Il1dCmBgYAoKIyMjIyBOdW1iZXIgb2YgZ2VuZXMgaW4gYWdyZWVtZW50IGFtb25nIERFIG1ldGhvZHMsIDIuMiBtb3JlIHRoYW4gMi4zCgpgYGB7ciBkZV9wbG90czAxfQp6eV90YWJsZVtbInZlbm5zIl1dW1sxXV1bWyJwX2xmYzEiXV1bWyJkb3duX25vd2VpZ2h0Il1dCmBgYAoKIyMjIyBnb3NlcSBvbnRvbG9neSBwbG90cyBvZiBncm91cHMgb2YgZ2VuZXMsIDIuMyBtb3JlIHRoYW4gMi4yCgpgYGB7ciBnb3NlcV91cH0KenlfZ29fdXAkcHZhbHVlX3Bsb3RzJGJwcF9wbG90X292ZXIKYGBgCgojIyMjIGdvc2VxIG9udG9sb2d5IHBsb3RzIG9mIGdyb3VwcyBvZiBnZW5lcywgMi4yIG1vcmUgdGhhbiAyLjMKCmBgYHtyIGdvc2VxX2Rvd259Cnp5X2dvX2Rvd24kcHZhbHVlX3Bsb3RzJGJwcF9wbG90X292ZXIKYGBgCgojIyBMb29rIGZvciBhZ3JlZW1lbnQgYmV0d2VlbiBzZW5zaXRpdml0eSBhbmQgenltb2RlbWVzCgpSZW1pbmQgbXlzZWxmLCB0aGUgZGF0YSBzdHJ1Y3R1cmVzIGFyZSAoenl8c3VzKV8oZGV8dGFibGV8c2lnKS4KCmBgYHtyIHNlbnNpdGl2ZV92c196eW1vfQp6eV9kZiA8LSB6eV90YWJsZVtbImRhdGEiXV1bWyJ6MjNfdnNfejIyIl1dCnN1c19kZiA8LSBzdXNfdGFibGVbWyJkYXRhIl1dW1sic2Vuc2l0aXZlX3ZzX3Jlc2lzdGFudCJdXQoKYm90aF9kZiA8LSBtZXJnZSh6eV9kZiwgc3VzX2RmLCBieSA9ICJyb3cubmFtZXMiKQpwbG90X2RmIDwtIGJvdGhfZGZbLCBjKCJkZXNlcV9sb2dmYy54IiwgImRlc2VxX2xvZ2ZjLnkiKV0Kcm93bmFtZXMocGxvdF9kZikgPC0gYm90aF9kZltbIlJvdy5uYW1lcyJdXQpjb2xuYW1lcyhwbG90X2RmKSA8LSBjKCJ6MjNfdnNfejIyIiwgInNlbnNpdGl2ZV92c19yZXNpc3RhbnQiKQoKY29tcGFyZSA8LSBwbG90X2xpbmVhcl9zY2F0dGVyKHBsb3RfZGYpCnBwKGZpbGUgPSAiaW1hZ2VzL2NvbXBhcmVfc3VzX3p5LnBuZyIsIGltYWdlID0gY29tcGFyZSRzY2F0dGVyKQpjb21wYXJlJGNvcgpgYGAKCiMjIFp5bW9kZW1lIGVuenltZSBnZW5lIElEcwoKTmFqaWIgcmVhZCBtZSBhbiBlbWFpbCBsaXN0aW5nIG9mZiB0aGUgZ2VuZSBuYW1lcyBhc3NvY2lhdGVkIHdpdGggdGhlIHp5bW9kZW1lCmNsYXNzaWZpY2F0aW9uLiAgSSB0b29rIHRob3NlIG5hbWVzIGFuZCBjcm9zcyByZWZlcmVuY2VkIHRoZW0gYWdhaW5zdCB0aGUKTGVpc2htYW5pYSBwYW5hbWVuc2lzIGdlbmUgYW5ub3RhdGlvbnMgYW5kIGZvdW5kIHRoZSBmb2xsb3dpbmc6CgpUaGV5IGFyZToKCjEuIEFMQVQ6IExQQUwxM18xMjAwMTA5MDAgLS0gYWxhbmluZSBhbWlub3RyYW5zZmVyYXNlCjIuIEFTQVQ6IExQQUwxM18zNDAwMTMwMDAgLS0gYXNwYXJ0YXRlIGFtaW5vdHJhbnNmZXJhc2UKMy4gRzZQRDogTFBBTDEzXzAwMDA1NDEwMCAtLSBnbHVjYXNlLTYtcGhvc3BoYXRlIDEtZGVoeWRyb2dlbmFzZQo0LiBOSDogTFBBTDEzXzE0MDA2MTAwLCBMUEFMMTNfMTgwMDE4NTAwIC0tIGlub3NpbmUtZ3VhbmluZSBudWNsZW9zaWRlIGh5ZHJvbGFzZQo1LiBNUEk6IExQQUwxM18zMjAwMjIzMDAgKG1heWJlKSAtLSBtYW5ub3NlIHBob3NwaGF0ZSBpc29tZXJhc2UgKEkgY2hvc2UgcGhvc3Bob21hbm5vc2UgaXNvbWVyYXNlKQoKR2l2ZW4gdGhlc2UgNiBnZW5lIElEcyAoTkggaGFzIHR3byBnZW5lIElEcyBhc3NvY2lhdGVkIHdpdGggaXQpLCBJIGNhbiBkbyBzb21lCmxvb2tpbmcgZm9yIHNwZWNpZmljIGRpZmZlcmVuY2VzIGFtb25nIHRoZSB2YXJpb3VzIHNhbXBsZXMuCgojIyMgRXhwcmVzc2lvbiBsZXZlbHMgb2Ygenltb2RlbWUgZ2VuZXMKClRoZSBmb2xsb3dpbmcgY3JlYXRlcyBhIGNvbG9yc3BhY2UgKHJlZCB0byBncmVlbikgaGVhdG1hcCBzaG93aW5nIHRoZSBvYnNlcnZlZApleHByZXNzaW9uIG9mIHRoZXNlIGdlbmVzIGluIGV2ZXJ5IHNhbXBsZS4KCmBgYHtyIHp5bW9kZW1lc30KbXlfZ2VuZXMgPC0gYygiTFBBTDEzXzEyMDAxMDkwMCIsICJMUEFMMTNfMzQwMDEzMDAwIiwgIkxQQUwxM18wMDAwNTQxMDAiLAogICAgICAgICAgICAgICJMUEFMMTNfMTQwMDA2MTAwIiwgIkxQQUwxM18xODAwMTg1MDAiLCAiTFBBTDEzXzMyMDAyMjMwMCIsCiAgICAgICAgICAgICAgIm90aGVyIikKbXlfbmFtZXMgPC0gYygiQUxBVCIsICJBU0FUIiwgIkc2UEQiLCAiTkh2MSIsICJOSHYyIiwgIk1QSSIsICJvdGhlciIpCgp6eW1vX2V4cHQgPC0gZXhjbHVkZV9nZW5lc19leHB0KHp5X25vcm0sIGlkcyA9IG15X2dlbmVzLCBtZXRob2QgPSAia2VlcCIpCnp5bW9faGVhdG1hcCA8LSBwbG90X3NhbXBsZV9oZWF0bWFwKHp5bW9fZXhwdCwgcm93X2xhYmVsID0gbXlfbmFtZXMpCnp5bW9faGVhdG1hcApgYGAKCiMjIEVtcGlyaWNhbGx5IG9ic2VydmVkIFp5bW9kZW1lIGdlbmVzIGZyb20gZGlmZmVyZW50aWFsIGV4cHJlc3Npb24gYW5hbHlzaXMKCkluIGNvbnRyYXN0LCB0aGUgZm9sbG93aW5nIHBsb3RzIHRha2UgdGhlIHNldCBvZiBnZW5lcyB3aGljaCBhcmUgc2hhcmVkIGFtb25nCmFsbCBkaWZmZXJlbnRpYWwgZXhwcmVzc2lvbiBtZXRob2RzICh8bGZjfCA+PSAxLjAgYW5kIGFkanAgPD0gMC4wNSkgYW5kIHVzZSB0aGVtCnRvIG1ha2UgY2F0ZWdvcmllcyBvZiBnZW5lcyB3aGljaCBhcmUgaW5jcmVhc2VkIGluIDIuMyBvciAyLjIuCgpgYGB7ciB6eW1vZGVtZV9nZW5lc19lbXBpcmljYWx9CnNoYXJlZF96eW1vIDwtIGludGVyc2VjdF9zaWduaWZpY2FudCh6eV90YWJsZSkKdXBfc2hhcmVkIDwtIHNoYXJlZF96eW1vW1sidXBzIl1dW1sxXV1bWyJkYXRhIl1dW1siYWxsIl1dCnJvd25hbWVzKHVwX3NoYXJlZCkKdXBzaGFyZWRfZXhwdCA8LSBleGNsdWRlX2dlbmVzX2V4cHQoenlfbm9ybSwgaWRzID0gcm93bmFtZXModXBfc2hhcmVkKSwgbWV0aG9kID0gImtlZXAiKQpgYGAKCldlIGNhbiBwbG90IGEgcXVpY2sgaGVhdG1hcCB0byBnZXQgYSBzZW5zZSBvZiB0aGUgZGlmZmVyZW5jZXMgb2JzZXJ2ZWQKYmV0d2VlbiB0aGUgZ2VuZXMgd2hpY2ggYXJlIGRpZmZlcmVudCBiZXR3ZWVuIHRoZSB0d28genltb2RlbWVzLgoKIyMjIEhlYXRtYXAgb2Ygenltb2RlbWUgZ2VuZSBleHByZXNzaW9uIGluY3JlYXNlZCBpbiAyLjMgdnMuIDIuMgoKYGBge3Igenltb2VtcHVwfQpoaWdoXzIzX2hlYXRtYXAgPC0gcGxvdF9zYW1wbGVfaGVhdG1hcCh1cHNoYXJlZF9leHB0LCByb3dfbGFiZWwgPSByb3duYW1lcyh1cF9zaGFyZWQpKQpoaWdoXzIzX2hlYXRtYXAKYGBgCgojIyMgSGVhdG1hcCBvZiB6eW1vZGVtZSBnZW5lIGV4cHJlc3Npb24gaW5jcmVhc2VkIGluIDIuMiB2cy4gMi4zCgpgYGB7ciB6eW1vZW1kb3dufQpkb3duX3NoYXJlZCA8LSBzaGFyZWRfenltb1tbImRvd25zIl1dW1sxXV1bWyJkYXRhIl1dW1siYWxsIl1dCmRvd25zaGFyZWRfZXhwdCA8LSBleGNsdWRlX2dlbmVzX2V4cHQoenlfbm9ybSwgaWRzID0gcm93bmFtZXMoZG93bl9zaGFyZWQpLCBtZXRob2QgPSAia2VlcCIpCmhpZ2hfMjJfaGVhdG1hcCA8LSBwbG90X3NhbXBsZV9oZWF0bWFwKGRvd25zaGFyZWRfZXhwdCwgcm93X2xhYmVsID0gcm93bmFtZXMoZG93bl9zaGFyZWQpKQpoaWdoXzIyX2hlYXRtYXAKYGBgCgojIFNOUCBwcm9maWxlcwoKTm93IEkgd2lsbCBjb21iaW5lIG91ciBwcmV2aW91cyBzYW1wbGVzIGFuZCBvdXIgbmV3IHNhbXBsZXMgaW4gdGhlCmhvcGVzIG9mIGZpbmRpbmcgdmFyaWFudCBwb3NpdGlvbnMgd2hpY2ggaGVscCBlbHVjaWRhdGUgY3VycmVudGx5CnVua25vd24gYXNwZWN0cyBvZiBlaXRoZXIgZ3JvdXAgdmlhIHRoZWlyIGNsdXN0ZXJpbmcgdG8ga25vd24gc2FtcGxlcwpmcm9tIHRoZSBvdGhlciBncm91cC4gSW4gb3RoZXIgd29yZHMsIHdlIGRvIG5vdCBrbm93IHRoZSB6eW1vZGVtZQphbm5vdGF0aW9ucyBmb3IgdGhlIG9sZCBzYW1wbGVzIG5vciB0aGUgc3RyYWluIGlkZW50aXRpZXMgKG9yIHRoZQpzaG9ydGN1dCAnY2hyb25pYyB2cy4gc2VsZi1oZWFsaW5nJykgZm9yIHRoZSBuZXcgc2FtcGxlcy4gSSBob3BlIHRvCm1ha2UgZWR1Y2F0ZWQgZ3Vlc3NlcyBnaXZlbiB0aGUgdmFyaWFudCBwcm9maWxlcy4gVGhlcmUgYXJlIHNvbWUKZGlmZmVyZW5jZXMgaW4gaG93IHRoZSBwcmV2aW91cyBhbmQgY3VycmVudCBkYXRhIHNldHMgd2VyZSBhbmFseXplZAoodGhvdWdoIEkgaGF2ZSBzaW5jZSByZWRvbmUgdGhlIG9sZCBzYW1wbGVzIHNvIGl0IHNob3VsZCBiZSB0cml2aWFsIHRvCnJlbW92ZSB0aG9zZSBkaWZmZXJlbmNlcyBub3cpLgoKSSBhZGRlZCBvdXIgMjAxNiBkYXRhIHRvIGEgc3BlY2lmaWMgVE1SQzIgc2FtcGxlIHNoZWV0LApkYXRlZCAyMDE5MTIwMy4gIFRodXMgSSB3aWxsIGxvYWQgdGhlIGRhdGEgaGVyZS4gIFRoYXQgcHJldmlvdXMgZGF0YQp3YXMgbWFwcGVkIHVzaW5nIHRvcGhhdCwgc28gSSB3aWxsIGFsc28gbmVlZCB0byBtYWtlIHNvbWUgY2hhbmdlcyB0bwp0aGUgZ2VuZSBuYW1lcyB0byBhY2NvbW9kYXRlIHRoZSB0d28gbWFwcGluZ3MuCgpgYGB7ciBvbGRuZXdfdmFyaWFudHN9Cm9sZF9leHB0IDwtIHNtKGNyZWF0ZV9leHB0KCJzYW1wbGVfc2hlZXRzL3RtcmMyX3NhbXBsZXNfMjAxOTEyMDMueGxzeCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgIGZpbGVfY29sdW1uID0gInRvcGhhdDJmaWxlIikpCgp0dCA8LSBscF9leHB0W1siZXhwcmVzc2lvbnNldCJdXQpyb3duYW1lcyh0dCkgPC0gZ3N1YihwYXR0ZXJuID0gIl5leG9uXyIsIHJlcGxhY2VtZW50ID0gIiIsIHggPSByb3duYW1lcyh0dCkpCnJvd25hbWVzKHR0KSA8LSBnc3ViKHBhdHRlcm4gPSAiXFwuRTEkIiwgcmVwbGFjZW1lbnQgPSAiIiwgeCA9IHJvd25hbWVzKHR0KSkKbHBfZXhwdCRleHByZXNzaW9uc2V0IDwtIHR0Cgp0dCA8LSBvbGRfZXhwdCRleHByZXNzaW9uc2V0CnJvd25hbWVzKHR0KSA8LSBnc3ViKHBhdHRlcm4gPSAiXmV4b25fIiwgcmVwbGFjZW1lbnQgPSAiIiwgeCA9IHJvd25hbWVzKHR0KSkKcm93bmFtZXModHQpIDwtIGdzdWIocGF0dGVybiA9ICJcXC4xJCIsIHJlcGxhY2VtZW50ID0gIiIsIHggPSByb3duYW1lcyh0dCkpCm9sZF9leHB0JGV4cHJlc3Npb25zZXQgPC0gdHQKcm0odHQpCmBgYAoKIyMgQ3JlYXRlIHRoZSBTTlAgZXhwcmVzc2lvbnNldAoKT25lIG90aGVyIGltcG9ydGFudCBjYXZlYXQsIHdlIGhhdmUgYSBncm91cCBvZiBuZXcgc2FtcGxlcyB3aGljaCBoYXZlCm5vdCB5ZXQgcnVuIHRocm91Z2ggdGhlIHZhcmlhbnQgc2VhcmNoIHBpcGVsaW5lLCBzbyBJIG5lZWQgdG8gcmVtb3ZlCnRoZW0gZnJvbSBjb25zaWRlcmF0aW9uLiAgVGhvdWdoIGl0IGxvb2tzIGxpa2UgdGhleSBmaW5pc2hlZCBvdmVybmlnaHQuLi4KCmBgYHtyIGNvdW50X2V4cHRfb2xkX25ld30KIyMgVGhlIG5leHQgbGluZSBkcm9wcyB0aGUgc2FtcGxlcyB3aGljaCBhcmUgbWlzc2luZyB0aGUgU05QIHBpcGVsaW5lLgpscF9zbnAgPC0gc3Vic2V0X2V4cHQobHBfZXhwdCwgc3Vic2V0PSIhaXMubmEocERhdGEobHBfZXhwdClbWydiY2Z0YWJsZSddXSkiKQpuZXdfc25wcyA8LSBzbShjb3VudF9leHB0X3NucHMobHBfc25wLCBhbm5vdF9jb2x1bW4gPSAiYmNmdGFibGUiKSkKb2xkX3NucHMgPC0gc20oY291bnRfZXhwdF9zbnBzKG9sZF9leHB0LCBhbm5vdF9jb2x1bW4gPSAiYmNmdGFibGUiLCBzbnBfY29sdW1uID0gMikpCgpib3RoX3NucHMgPC0gY29tYmluZV9leHB0cyhuZXdfc25wcywgb2xkX3NucHMpCmJvdGhfbm9ybSA8LSBzbShub3JtYWxpemVfZXhwdChib3RoX3NucHMsIHRyYW5zZm9ybSA9ICJsb2cyIiwgY29udmVydCA9ICJjcG0iLCBmaWx0ZXIgPSBUUlVFKSkKCiMjIHN0cmFpbnMgPC0gYm90aF9ub3JtW1siZGVzaWduIl1dW1sic3RyYWluIl1dCmJvdGhfc3RyYWluIDwtIHNldF9leHB0X2NvbmRpdGlvbnMoYm90aF9ub3JtLCBmYWN0ID0gInN0cmFpbiIpCmBgYAoKVGhlIGRhdGEgc3RydWN0dXJlICdib3RoX25vcm0nIG5vdyBjb250YWlucyBvdXIgMjAxNiBkYXRhIGFsb25nIHdpdGgKdGhlIG5ld2VyIGRhdGEgY29sbGVjdGVkIHNpbmNlIDIwMTkuCgojIyBQbG90IG9mIFNOUCBwcm9maWxlcyBmb3Igenltb2RlbWVzCgpUaGUgZm9sbG93aW5nIHBsb3Qgc2hvd3MgdGhlIFNOUCBwcm9maWxlcyBvZiBhbGwgc2FtcGxlcyAob2xkIGFuZCBuZXcpIHdoZXJlIHRoZQpjb2xvcnMgYXQgdGhlIHRvcCBzaG93IGVpdGhlciB0aGUgMi4yIHN0cmFpbnMgKG9yYW5nZSksIDIuMyBzdHJhaW5zIChncmVlbiksIHRoZQpwcmV2aW91cyBzYW1wbGVzIChwdXJwbGUpLCBvciB0aGUgdmFyaW91cyBsYWIgc3RyYWlucyAocGluayBldGMpLgoKYGBge3IgcGxvdHRpbmdfdmFyaWFudHN9Cm9sZF9uZXdfdmFyaWFudF9oZWF0bWFwIDwtIHBsb3RfZGlzaGVhdChib3RoX25vcm0pCnBwKGZpbGUgPSAiaW1hZ2VzL3Jhd19zbnBfZGlzaGVhdC5wbmciLCBpbWFnZSA9IG9sZF9uZXdfdmFyaWFudF9oZWF0bWFwLAogICBoZWlnaHQgPSAxMiwgd2lkdGggPSAxMikKYGBgCgpUaGUgZnVuY3Rpb24gZ2V0X3NucF9zZXRzKCkgdGFrZXMgdGhlIHByb3ZpZGVkIG1ldGFkYXRhIGZhY3RvciAoaW4KdGhpcyBjYXNlICdjb25kaXRpb24nKSBhbmQgbG9va3MgZm9yIHZhcmlhbnRzIHdoaWNoIGFyZSBleGNsdXNpdmUgdG8KZWFjaCBlbGVtZW50IGluIGl0LiAgSW4gdGhpcyBjYXNlLCB0aGlzIGlzIGxvb2tpbmcgZm9yIGRpZmZlcmVuY2VzCmJldHdlZW4gMi4yIGFuZCAyLjMsIGFzIHdlbGwgYXMgdGhlIHNldCBzaGFyZWQgYW1vbmcgdGhlbS4KCmBgYHtyIGdldF9zbnBfc2V0czF9CnNucF9zZXRzIDwtIGdldF9zbnBfc2V0cyhib3RoX3NucHMsIGZhY3RvciA9ICJjb25kaXRpb24iKQpib3RoX2V4cHQgPC0gY29tYmluZV9leHB0cyhscF9leHB0LCBvbGRfZXhwdCkKCnNucF9nZW5lcyA8LSBzbShzbnBzX3ZzX2dlbmVzKGJvdGhfZXhwdCwgc25wX3NldHMsIGV4cHRfbmFtZV9jb2wgPSAiY2hyb21vc29tZSIpKQojIyBJIHRoaW5rIHdlIGhhdmUgc29tZSBtZXRyaWNzIGhlcmUgd2UgY2FuIHBsb3QuLi4Kc25wX3N1YnNldCA8LSBzbShzbnBfc3Vic2V0X2dlbmVzKAogIGJvdGhfZXhwdCwgYm90aF9zbnBzLAogIGdlbmVzID0gYygiTFBBTDEzXzEyMDAxMDkwMCIsICJMUEFMMTNfMzQwMDEzMDAwIiwgIkxQQUwxM18wMDAwNTQxMDAiLAogICAgICAgICAgICAiTFBBTDEzXzE0MDAwNjEwMCIsICJMUEFMMTNfMTgwMDE4NTAwIiwgIkxQQUwxM18zMjAwMjIzMDAiKSkpCnp5bW9faGVhdCA8LSBwbG90X3NhbXBsZV9oZWF0bWFwKHNucF9zdWJzZXQsIHJvd19sYWJlbCA9IHJvd25hbWVzKGV4cHJzKHNucF9zdWJzZXQpKSkKenltb19oZWF0CmBgYAoKRGlkbid0IEkgY3JlYXRlIGEgc2V0IG9mIGRlbnNpdGllcyBieSBjaHJvbW9zb21lPwpPaCBJIHRoaW5rIHRoZXkgY29tZSBpbiBmcm9tIGdldF9zbnBfc2V0cygpCgojIyBTTlBTIGFzc29jaWF0ZWQgd2l0aCBjbGluaWNhbCByZXNwb25zZSBpbiB0aGUgVE1SQyBzYW1wbGVzCgpgYGB7ciBzbnBfY2xpbmljYWx9CmNsaW5pY2FsX3NldHMgPC0gZ2V0X3NucF9zZXRzKG5ld19zbnBzLCBmYWN0b3IgPSAiY2xpbmljYWxyZXNwb25zZSIpCgpkZW5zaXR5X3ZlYyA8LSBjbGluaWNhbF9zZXRzW1siZGVuc2l0eSJdXQpjaHJvbW9zb21lX2lkeCA8LSBncmVwKHBhdHRlcm4gPSAiTHBhTCIsIHggPSBuYW1lcyhkZW5zaXR5X3ZlYykpCmRlbnNpdHlfZGYgPC0gYXMuZGF0YS5mcmFtZShkZW5zaXR5X3ZlY1tjaHJvbW9zb21lX2lkeF0pCmRlbnNpdHlfZGZbWyJjaHIiXV0gPC0gcm93bmFtZXMoZGVuc2l0eV9kZikKY29sbmFtZXMoZGVuc2l0eV9kZikgPC0gYygiZGVuc2l0eV92ZWMiLCAiY2hyIikKZ2dwbG90KGRlbnNpdHlfZGYsIGFlc19zdHJpbmcoeCA9ICJjaHIiLCB5ID0gImRlbnNpdHlfdmVjIikpICsKICBnZ3Bsb3QyOjpnZW9tX2NvbCgpICsKICBnZ3Bsb3QyOjp0aGVtZShheGlzLnRleHQgPSBnZ3Bsb3QyOjplbGVtZW50X3RleHQoc2l6ZSA9IDEwLCBjb2xvdXIgPSAiYmxhY2siKSwKICAgICAgICAgICAgICAgICBheGlzLnRleHQueCA9IGdncGxvdDI6OmVsZW1lbnRfdGV4dChhbmdsZSA9IDkwLCB2anVzdCA9IDAuNSkpCgojIyBjbGluaWNhbF93cml0dGVuIDwtIHdyaXRlX3ZhcmlhbnRzKG5ld19zbnBzKQpgYGAKCiMjIyBDcm9zcyByZWZlcmVuY2UgdGhlc2UgdmFyaWFudHMgYnkgZ2VuZQoKYGBge3Igc25wX2NsYXNzaWZpY2F0aW9uc30KY2xpbmljYWxfZ2VuZXMgPC0gc20oc25wc192c19nZW5lcyhscF9leHB0LCBjbGluaWNhbF9zZXRzLCBleHB0X25hbWVfY29sID0gImNocm9tb3NvbWUiKSkKCnNucF9kZW5zaXR5IDwtIG1lcmdlKGFzLmRhdGEuZnJhbWUoY2xpbmljYWxfZ2VuZXNbWyJzdW1tYXJ5X2J5X2dlbmUiXV0pLAogICAgICAgICAgICAgICAgICAgICBhcy5kYXRhLmZyYW1lKGZEYXRhKGxwX2V4cHQpKSwKICAgICAgICAgICAgICAgICAgICAgYnkgPSAicm93Lm5hbWVzIikKc25wX2RlbnNpdHkgPC0gc25wX2RlbnNpdHlbLCBjKDEsIDIsIDQsIDE1KV0KY29sbmFtZXMoc25wX2RlbnNpdHkpIDwtIGMoIm5hbWUiLCAic25wcyIsICJwcm9kdWN0IiwgImxlbmd0aCIpCnNucF9kZW5zaXR5W1sicHJvZHVjdCJdXSA8LSB0b2xvd2VyKHNucF9kZW5zaXR5W1sicHJvZHVjdCJdXSkKc25wX2RlbnNpdHlbWyJsZW5ndGgiXV0gPC0gYXMubnVtZXJpYyhzbnBfZGVuc2l0eVtbImxlbmd0aCJdXSkKc25wX2RlbnNpdHlbWyJkZW5zaXR5Il1dIDwtIHNucF9kZW5zaXR5W1sic25wcyJdXSAvIHNucF9kZW5zaXR5W1sibGVuZ3RoIl1dCnNucF9pZHggPC0gb3JkZXIoc25wX2RlbnNpdHlbWyJkZW5zaXR5Il1dLCBkZWNyZWFzaW5nID0gVFJVRSkKc25wX2RlbnNpdHkgPC0gc25wX2RlbnNpdHlbc25wX2lkeCwgXQoKcmVtb3ZlcnMgPC0gYygiYW1hc3RpbiIsICJncDYzIiwgImxlaXNobWFub2x5c2luIikKZm9yIChyIGluIHJlbW92ZXJzKSB7CiAgZHJvcF9pZHggPC0gZ3JlcGwocGF0dGVybiA9IHIsIHggPSBzbnBfZGVuc2l0eVtbInByb2R1Y3QiXV0pCiAgc25wX2RlbnNpdHkgPC0gc25wX2RlbnNpdHlbIWRyb3BfaWR4LCBdCn0KIyMgRmlsdGVyIHRoZXNlIGZvciBbQXxhXW1hc3RpbiBncDYzIExlaXNobWFub2x5c2luCmBgYAoKCmBgYHtyIHNucF9pbnRlcnNlY3Rpb25zfQpjbGluaWNhbF9zbnBzIDwtIHNucHNfaW50ZXJzZWN0aW9ucyhscF9leHB0LCBjbGluaWNhbF9zZXRzLCBjaHJfY29sdW1uID0gImNocm9tb3NvbWUiKQoKZmFpbF9yZWZfc25wcyA8LSBhcy5kYXRhLmZyYW1lKGNsaW5pY2FsX3NucHNbWyJpbnRlcnMiXV1bWyJmYWlsdXJlLCByZWZlcmVuY2Ugc3RyYWluIl1dKQpjdXJlX3NucHMgPC0gYXMuZGF0YS5mcmFtZShjbGluaWNhbF9zbnBzW1siaW50ZXJzIl1dW1siY3VyZSJdXSkKCmhlYWQoZmFpbF9yZWZfc25wcykKaGVhZChjdXJlX3NucHMpCgphbm5vdCA8LSBmRGF0YShscF9leHB0KQpjbGluaWNhbF9pbnRlcmVzdCA8LSBhcy5kYXRhLmZyYW1lKGNsaW5pY2FsX3NucHNbWyJnZW5lX3N1bW1hcmllcyJdXVtbImN1cmUiXV0pCmNsaW5pY2FsX2ludGVyZXN0IDwtIG1lcmdlKGNsaW5pY2FsX2ludGVyZXN0LAogICAgICAgICAgICAgICAgICAgICAgICAgICBhcy5kYXRhLmZyYW1lKGNsaW5pY2FsX3NucHNbWyJnZW5lX3N1bW1hcmllcyJdXVtbImZhaWx1cmUsIHJlZmVyZW5jZSBzdHJhaW4iXV0pLAogICAgICAgICAgICAgICAgICAgICAgICAgICBieSA9ICJyb3cubmFtZXMiKQpyb3duYW1lcyhjbGluaWNhbF9pbnRlcmVzdCkgPC0gY2xpbmljYWxfaW50ZXJlc3RbWyJSb3cubmFtZXMiXV0KY2xpbmljYWxfaW50ZXJlc3RbWyJSb3cubmFtZXMiXV0gPC0gTlVMTApjb2xuYW1lcyhjbGluaWNhbF9pbnRlcmVzdCkgPC0gYygiY3VyZV9zbnBzIiwiZmFpbF9zbnBzIikKYW5ub3QgPC0gbWVyZ2UoYW5ub3QsIGNsaW5pY2FsX2ludGVyZXN0LCBieSA9ICJyb3cubmFtZXMiKQpyb3duYW1lcyhhbm5vdCkgPC0gYW5ub3RbWyJSb3cubmFtZXMiXV0KYW5ub3RbWyJSb3cubmFtZXMiXV0gPC0gTlVMTApmRGF0YShscF9leHB0JGV4cHJlc3Npb25zZXQpIDwtIGFubm90CmBgYAoKIyBaeW1vZGVtZSBmb3IgbmV3IHNhbXBsZXMKClRoZSBoZWF0bWFwIHByb2R1Y2VkIGhlcmUgc2hvdWxkIHNob3cgdGhlIHZhcmlhbnRzIG9ubHkgZm9yIHRoZSB6eW1vZGVtZSBnZW5lcy4KCiMjIEh1bnQgZm9yIHNucCBjbHVzdGVycwoKSSBhbSB0aGlua2luZyB0aGF0IGlmIHdlIGZpbmQgY2x1c3RlcnMgb2YgbG9jYXRpb25zIHdoaWNoIGFyZSB2YXJpYW50LCB0aGF0Cm1pZ2h0IHByb3ZpZGUgc29tZSBQQ1IgdGVzdGluZyBwb3NzaWJpbGl0aWVzLgoKYGBge3IgbmV3X3p5bW99CiMjIERyb3AgdGhlIDIuMSwgMi40LCB1bmtub3duLCBhbmQgbnVsbApwcnVuZWRfc25wcyA8LSBzdWJzZXRfZXhwdChuZXdfc25wcywgc3Vic2V0PSJjb25kaXRpb249PSd6Mi4yJ3xjb25kaXRpb249PSd6Mi4zJyIpCm5ld19zZXRzIDwtIGdldF9zbnBfc2V0cyhwcnVuZWRfc25wcywgZmFjdG9yID0gInp5bW9kZW1lY2F0ZWdvcmljYWwiKQpzdW1tYXJ5KG5ld19zZXRzKQojIyAxMDAwMDAwOiAyLjIKIyMgMDEwMDAwMDogMi4zCgpzdW1tYXJ5KG5ld19zZXRzW1siaW50ZXJzZWN0aW9ucyJdXVtbIjEwIl1dKQpzdW1tYXJ5KG5ld19zZXRzW1siaW50ZXJzZWN0aW9ucyJdXVtbIjAxIl1dKQpgYGAKClRodXMgd2Ugc2VlIHRoYXQgdGhlcmUgYXJlIDUxMSB2YXJpYW50cyBhc3NvY2lhdGVkIHdpdGggMi4yIGFuZCA0OSw3OTAgYXNzb2NpYXRlZCB3aXRoIDIuMy4KCiMjIyBBIHNtYWxsIGZ1bmN0aW9uIGZvciBzZWFyY2hpbmcgZm9yIHBvdGVudGlhbCBQQ1IgcHJpbWVycwoKVGhlIGZvbGxvd2luZyBmdW5jdGlvbiB1c2VzIHRoZSBwb3NpdGlvbmFsIGRhdGEgdG8gbG9vayBmb3Igc2VxdWVudGlhbAptaXNtYXRjaGVzIGFzc29jaWF0ZWQgd2l0aCB6eW1vZGVtZSBpbiB0aGUgaG9wZXMgdGhhdCB0aGVyZSB3aWxsIGJlCnNvbWUgcmVnaW9ucyB3aGljaCB3b3VsZCBwcm92aWRlIGdvb2QgcG90ZW50aWFsIHRhcmdldHMgZm9yIGEKUENSLWJhc2VkIGFzc2F5LgoKYGBge3Igc2VxdWVudGlhbF9zZWFyY2gsIGV2YWw9RkFMU0V9CnNlcXVlbnRpYWxfdmFyaWFudHMgPC0gZnVuY3Rpb24oc25wX3NldHMsIGNvbmRpdGlvbnMgPSBOVUxMLCBtaW5pbXVtID0gMywgbWF4aW11bV9zZXBhcmF0aW9uID0gMykgewogIGlmIChpcy5udWxsKGNvbmRpdGlvbnMpKSB7CiAgICBjb25kaXRpb25zIDwtIDEKICB9CiAgaW50ZXJzZWN0aW9uX3NldHMgPC0gc25wX3NldHNbWyJpbnRlcnNlY3Rpb25zIl1dCiAgaW50ZXJzZWN0aW9uX25hbWVzIDwtIHNucF9zZXRzW1sic2V0X25hbWVzIl1dCiAgY2hvc2VuX2ludGVyc2VjdGlvbiA8LSAxCiAgaWYgKGlzLm51bWVyaWMoY29uZGl0aW9ucykpIHsKICAgIGNob3Nlbl9pbnRlcnNlY3Rpb24gPC0gY29uZGl0aW9ucwogIH0gZWxzZSB7CiAgICBpbnRlcnNlY3Rpb25faWR4IDwtIGludGVyc2VjdGlvbl9uYW1lcyA9PSBjb25kaXRpb25zCiAgICBjaG9zZW5faW50ZXJzZWN0aW9uIDwtIG5hbWVzKGludGVyc2VjdGlvbl9uYW1lcylbaW50ZXJzZWN0aW9uX2lkeF0KICB9CgogIHBvc3NpYmxlX3Bvc2l0aW9ucyA8LSBpbnRlcnNlY3Rpb25fc2V0c1tbY2hvc2VuX2ludGVyc2VjdGlvbl1dCiAgcG9zaXRpb25fdGFibGUgPC0gZGF0YS5mcmFtZShyb3cubmFtZXMgPSBwb3NzaWJsZV9wb3NpdGlvbnMpCiAgcGF0IDwtICJeY2hyXyguKylfcG9zXyguKylfcmVmXy4qJCIKICBwb3NpdGlvbl90YWJsZVtbImNociJdXSA8LSBnc3ViKHBhdHRlcm4gPSBwYXQsIHJlcGxhY2VtZW50ID0gIlxcMSIsIHggPSByb3duYW1lcyhwb3NpdGlvbl90YWJsZSkpCiAgcG9zaXRpb25fdGFibGVbWyJwb3MiXV0gPC0gYXMubnVtZXJpYyhnc3ViKHBhdHRlcm4gPSBwYXQsIHJlcGxhY2VtZW50ID0gIlxcMiIsIHggPSByb3duYW1lcyhwb3NpdGlvbl90YWJsZSkpKQogIHBvc2l0aW9uX2lkeCA8LSBvcmRlcihwb3NpdGlvbl90YWJsZVssICJjaHIiXSwgcG9zaXRpb25fdGFibGVbLCAicG9zIl0pCiAgcG9zaXRpb25fdGFibGUgPC0gcG9zaXRpb25fdGFibGVbcG9zaXRpb25faWR4LCBdCiAgcG9zaXRpb25fdGFibGVbWyJkaXN0Il1dIDwtIDAKCiAgbGFzdF9jaHIgPC0gIiIKICBmb3IgKHIgaW4gMTpucm93KHBvc2l0aW9uX3RhYmxlKSkgewogICAgdGhpc19jaHIgPC0gcG9zaXRpb25fdGFibGVbciwgImNociJdCiAgICBpZiAociA9PSAxKSB7CiAgICAgIHBvc2l0aW9uX3RhYmxlW3IsICJkaXN0Il0gPC0gcG9zaXRpb25fdGFibGVbciwgInBvcyJdCiAgICAgIGxhc3RfY2hyIDwtIHRoaXNfY2hyCiAgICAgIG5leHQKICAgIH0KICAgIGlmICh0aGlzX2NociA9PSBsYXN0X2NocikgewogICAgICBwb3NpdGlvbl90YWJsZVtyLCAiZGlzdCJdIDwtIHBvc2l0aW9uX3RhYmxlW3IsICJwb3MiXSAtIHBvc2l0aW9uX3RhYmxlW3IgLSAxLCAicG9zIl0KICAgIH0gZWxzZSB7CiAgICAgIHBvc2l0aW9uX3RhYmxlW3IsICJkaXN0Il0gPC0gcG9zaXRpb25fdGFibGVbciwgInBvcyJdCiAgICB9CiAgICBsYXN0X2NociA8LSB0aGlzX2NocgogIH0KCiAgIyMgV29ya2luZyBpbnRlcmFjdGl2ZWx5IGhlcmUuCgogIGRvdWJsZXMgPC0gcG9zaXRpb25fdGFibGVbWyJkaXN0Il1dID09IDEKICBkb3VibGVzIDwtIHBvc2l0aW9uX3RhYmxlW2RvdWJsZXMsIF0KICB3cml0ZS5jc3YoZG91YmxlcywgImRvdWJsZXMuY3N2IikKCiAgb25lX2F3YXkgPC0gcG9zaXRpb25fdGFibGVbWyJkaXN0Il1dID09IDIKICBvbmVfYXdheSA8LSBwb3NpdGlvbl90YWJsZVtvbmVfYXdheSwgXQogIHdyaXRlLmNzdihvbmVfYXdheSwgIm9uZV9hd2F5LmNzdiIpCgogIHR3b19hd2F5IDwtIHBvc2l0aW9uX3RhYmxlW1siZGlzdCJdXSA9PSAzCiAgdHdvX2F3YXkgPC0gcG9zaXRpb25fdGFibGVbdHdvX2F3YXksIF0KICB3cml0ZS5jc3YodHdvX2F3YXksICJ0d29fYXdheS5jc3YiKQoKICBjb21iaW5lZCA8LSByYmluZChkb3VibGVzLCBvbmVfYXdheSkKICBjb21iaW5lZCA8LSByYmluZChjb21iaW5lZCwgdHdvX2F3YXkpCiAgcG9zaXRpb25faWR4IDwtIG9yZGVyKGNvbWJpbmVkWywgImNociJdLCBjb21iaW5lZFssICJwb3MiXSkKICBjb21iaW5lZCA8LSBjb21iaW5lZFtwb3NpdGlvbl9pZHgsIF0KCiAgdGhpc19jaHIgPC0gIiIKICBmb3IgKHIgaW4gMTpucm93KGNvbWJpbmVkKSkgewogICAgdGhpc19jaHIgPC0gY29tYmluZWRbciwgImNociJdCiAgICBpZiAociA9PSAxKSB7CiAgICAgIGNvbWJpbmVkW3IsICJkaXN0X3BhaXIiXSA8LSBjb21iaW5lZFtyLCAicG9zIl0KICAgICAgbGFzdF9jaHIgPC0gdGhpc19jaHIKICAgICAgbmV4dAogICAgfQogICAgaWYgKHRoaXNfY2hyID09IGxhc3RfY2hyKSB7CiAgICAgIGNvbWJpbmVkW3IsICJkaXN0X3BhaXIiXSA8LSBjb21iaW5lZFtyLCAicG9zIl0gLSBjb21iaW5lZFtyIC0gMSwgInBvcyJdCiAgICB9IGVsc2UgewogICAgICBjb21iaW5lZFtyLCAiZGlzdF9wYWlyIl0gPC0gY29tYmluZWRbciwgInBvcyJdCiAgICB9CiAgICBsYXN0X2NociA8LSB0aGlzX2NocgogIH0KCiAgZGlzdF9wYWlyX21heGltdW0gPC0gMTAwMAogIGRpc3RfcGFpcl9taW5pbXVtIDwtIDIwMAogIGRpc3RfcGFpcl9pZHggPC0gY29tYmluZWRbWyJkaXN0X3BhaXIiXV0gPD0gZGlzdF9wYWlyX21heGltdW0gJgogICAgY29tYmluZWRbWyJkaXN0X3BhaXIiXV0gPj0gZGlzdF9wYWlyX21pbmltdW0KICByZW1haW5pbmcgPC0gY29tYmluZWRbZGlzdF9wYWlyX2lkeCwgXQogIG5vX3dlYWtfaWR4IDwtIGdyZXBsKHBhdHRlcm49InJlZl8oR3xDKSIsIHg9cm93bmFtZXMocmVtYWluaW5nKSkKICByZW1haW5pbmcgPC0gcmVtYWluaW5nW25vX3dlYWtfaWR4LCBdCgogIHByaW50KGhlYWQodGFibGUocG9zaXRpb25fdGFibGVbWyJkaXN0Il1dKSkpCiAgc2VxdWVudGlhbHMgPC0gcG9zaXRpb25fdGFibGVbWyJkaXN0Il1dIDw9IG1heGltdW1fc2VwYXJhdGlvbgogIG1lc3NhZ2UoIlRoZXJlIGFyZSAiLCBzdW0oc2VxdWVudGlhbHMpLCAiIGNhbmRpZGF0ZSByZWdpb25zLiIpCgogICMjIFRoZSBmb2xsb3dpbmcgY2FuIHRlbGwgbWUgaG93IG1hbnkgcnVucyBvZiBlYWNoIGxlbmd0aCBvY2N1cnJlZCwgdGhhdCBpcyBub3QgcXVpdGUgd2hhdCBJIHdhbnQuCiAgIyMgTm93IHVzZSBydW4gbGVuZ3RoIGVuY29kaW5nIHRvIGZpbmQgdGhlIHNldCBvZiBzZXF1ZW50aWFsIHNlcXVlbnRpYWxzIQogIHJsZV9yZXN1bHQgPC0gcmxlKHNlcXVlbnRpYWxzKQogIHJsZV92YWx1ZXMgPC0gcmxlX3Jlc3VsdFtbInZhbHVlcyJdXQogICMjIFRoZSBmb2xsb3dpbmcgbGluZSBpcyBlcXVpdmFsZW50IHRvIGp1c3QgbGVhdmluZyB2YWx1ZXMgYWxvbmU6CiAgIyMgdHJ1ZV92YWx1ZXMgPC0gcmxlX3Jlc3VsdFtbInZhbHVlcyJdXSA9PSBUUlVFCiAgcmxlX2xlbmd0aHMgPC0gcmxlX3Jlc3VsdFtbImxlbmd0aHMiXV0KICB0cnVlX3NlcXVlbnRpYWxzIDwtIHJsZV9sZW5ndGhzW3JsZV92YWx1ZXNdCiAgcmxlX2lkeCA8LSBjdW1zdW0ocmxlX2xlbmd0aHMpW3doaWNoKHJsZV92YWx1ZXMpXQoKICBwb3NpdGlvbl90YWJsZVtbImxhc3Rfc2VxdWVudGlhbCJdXSA8LSAwCiAgY291bnQgPC0gMAogIGZvciAociBpbiBybGVfaWR4KSB7CiAgICBjb3VudCA8LSBjb3VudCArIDEKICAgIHBvc2l0aW9uX3RhYmxlW3IsICJsYXN0X3NlcXVlbnRpYWwiXSA8LSB0cnVlX3NlcXVlbnRpYWxzW2NvdW50XQogIH0KICBtZXNzYWdlKCJUaGUgbWF4aW11bSBzZXF1ZW50aWFsIHNldCBpczogIiwgbWF4KHBvc2l0aW9uX3RhYmxlW1sibGFzdF9zZXF1ZW50aWFsIl1dKSwgIi4iKQoKICB3YW50ZWRfaWR4IDwtIHBvc2l0aW9uX3RhYmxlW1sibGFzdF9zZXF1ZW50aWFsIl1dID49IG1pbmltdW0KICB3YW50ZWQgPC0gcG9zaXRpb25fdGFibGVbd2FudGVkX2lkeCwgYygiY2hyIiwgInBvcyIpXQogIHJldHVybih3YW50ZWQpCn0KCnp5bW8yMl9zZXF1ZW50aWFscyA8LSBzZXF1ZW50aWFsX3ZhcmlhbnRzKG5ld19zZXRzLCBjb25kaXRpb25zID0gInoyMiIsIG1pbmltdW09MSwgbWF4aW11bV9zZXBhcmF0aW9uPTIpCmRpbSh6eW1vMjJfc2VxdWVudGlhbHMpCiMjIDcgY2FuZGlkYXRlIHJlZ2lvbnMgZm9yIHp5bW9kZW1lIDIuMiAtLSB0aHVzIEkgYW0gYmV0dGluZyB0aGF0IHRoZSByZWZlcmVuY2Ugc3RyYWluIGlzIGEgMi4yCnp5bW8yM19zZXF1ZW50aWFscyA8LSBzZXF1ZW50aWFsX3ZhcmlhbnRzKG5ld19zZXRzLCBjb25kaXRpb25zID0gInoyMyIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG1pbmltdW0gPSAyLCBtYXhpbXVtX3NlcGFyYXRpb24gPSAyKQpkaW0oenltbzIzX3NlcXVlbnRpYWxzKQojIyBJbiBjb250cmFzdCwgdGhlcmUgYXJlIGxvdHMgKDU4Nykgb2YgaW50ZXJlc3RpbmcgcmVnaW9ucyBmb3IgMi4zIQpgYGAKCiMjIyBFeHRyYWN0IGEgcHJvbWlzaW5nIHJlZ2lvbiBmcm9tIHRoZSBnZW5vbWUKClRoZSBmaXJzdCA0IGNhbmRpZGF0ZSByZWdpb25zIGZyb20gbXkgc2V0IG9mIHJlbWFpbmluZzoKKiBDaHIgICAgICAgUG9zLiAgIERpc3RhbmNlCiogTHBhTDEzLTE1IDIzODQzMyA0NDgKKiBMcGFMMTMtMTggMTQyODQ0IDYxMwoqIExwYUwxMy0yOSA4MzAzNDIgMjUyCiogTHBhTDEzLTMzIDEzMzE1MDcgODQzCgpMZXRzIGRlZmluZSBhIGNvdXBsZSBvZiB0ZXJtczoKKiBUaGlyZDogRWFjaCBvZiB0aGUgNCBhYm92ZSBwb3NpdGlvbnMuCiogU2Vjb25kOiBUaGlyZCAtIERpc3RhbmNlCiogRW5kOiBUaGlyZCArIFByaW1lckxlbgoqIFN0YXJ0OiBTZWNvbmQgLSBQcmltZXJsZW4KCkluIGVhY2ggaW5zdGFuY2UsIHRoZXNlIGFyZSB0aGUgbGFzdCBwb3NpdGlvbnMsIHNvIHdlIHdhbnQgdG8gZ3JhYiB0aHJlZSB0aGluZ3M6CgoqIFRoZSBlbnRpcmUgcmVnaW9uIGZyb20gRW5kIC0+IFN0YXJ0LCB0aGlzIHdheSB3ZSBjYW4gaGF2ZSBhIHF1aWNrIHNhbml0eSBjaGVjay4KKiBTdGFydCAtPiBTZWNvbmQuCiogKFRoaXJkIC0+IEVuZCkgPC0gUmV2ZXJzZSBjb21wbGVtZW50ZWQKCmBgYHtyIGV4dHJhY3RfYnNnZW5vbWUsIGV2YWw9RkFMU0V9CiMjICogTHBhTDEzLTE1IDIzODQzMyA0NDgKZmlyc3RfY2FuZGlkYXRlX2NociA8LSBnZW5vbWVbWyJMcGFMMTNfMTUiXV0KcHJpbWVyX2xlbmd0aCA8LSAyMgphbXBsaWNvbl9sZW5ndGggPC0gNDQ4CmZpcnN0X2NhbmRpZGF0ZV90aGlyZCA8LSAyMzg0MzMKZmlyc3RfY2FuZGlkYXRlX3NlY29uZCA8LSBmaXJzdF9jYW5kaWRhdGVfdGhpcmQgLSBhbXBsaWNvbl9sZW5ndGgKZmlyc3RfY2FuZGlkYXRlX3N0YXJ0IDwtIGZpcnN0X2NhbmRpZGF0ZV9zZWNvbmQgLSBwcmltZXJfbGVuZ3RoCmZpcnN0X2NhbmRpZGF0ZV9lbmQgPC0gZmlyc3RfY2FuZGlkYXRlX3RoaXJkICsgcHJpbWVyX2xlbmd0aApmaXJzdF9jYW5kaWRhdGVfcmVnaW9uIDwtIHN1YnNlcShmaXJzdF9jYW5kaWRhdGVfY2hyLCBmaXJzdF9jYW5kaWRhdGVfc3RhcnQsIGZpcnN0X2NhbmRpZGF0ZV9lbmQpCmZpcnN0X2NhbmRpZGF0ZV9yZWdpb24KZmlyc3RfY2FuZGlkYXRlXzVwIDwtIHN1YnNlcShmaXJzdF9jYW5kaWRhdGVfY2hyLCBmaXJzdF9jYW5kaWRhdGVfc3RhcnQsIGZpcnN0X2NhbmRpZGF0ZV9zZWNvbmQpCmFzLmNoYXJhY3RlcihmaXJzdF9jYW5kaWRhdGVfNXApCmZpcnN0X2NhbmRpZGF0ZV8zcCA8LSBzcGdzOjpyZXZlcnNlQ29tcGxlbWVudChzdWJzZXEoZmlyc3RfY2FuZGlkYXRlX2NociwgZmlyc3RfY2FuZGlkYXRlX3RoaXJkLCBmaXJzdF9jYW5kaWRhdGVfZW5kKSkKZmlyc3RfY2FuZGlkYXRlXzNwCgoKIyMgKiBMcGFMMTMtMTggMTQyODQ0IDYxMwpzZWNvbmRfY2FuZGlkYXRlX2NociA8LSBnZW5vbWVbWyJMcGFMMTNfMTgiXV0KcHJpbWVyX2xlbmd0aCA8LSAyMgphbXBsaWNvbl9sZW5ndGggPC0gNjEzCnNlY29uZF9jYW5kaWRhdGVfdGhpcmQgPC0gMTQyODQ0CnNlY29uZF9jYW5kaWRhdGVfc2Vjb25kIDwtIHNlY29uZF9jYW5kaWRhdGVfdGhpcmQgLSBhbXBsaWNvbl9sZW5ndGgKc2Vjb25kX2NhbmRpZGF0ZV9zdGFydCA8LSBzZWNvbmRfY2FuZGlkYXRlX3NlY29uZCAtIHByaW1lcl9sZW5ndGgKc2Vjb25kX2NhbmRpZGF0ZV9lbmQgPC0gc2Vjb25kX2NhbmRpZGF0ZV90aGlyZCArIHByaW1lcl9sZW5ndGgKc2Vjb25kX2NhbmRpZGF0ZV9yZWdpb24gPC0gc3Vic2VxKHNlY29uZF9jYW5kaWRhdGVfY2hyLCBzZWNvbmRfY2FuZGlkYXRlX3N0YXJ0LCBzZWNvbmRfY2FuZGlkYXRlX2VuZCkKc2Vjb25kX2NhbmRpZGF0ZV9yZWdpb24Kc2Vjb25kX2NhbmRpZGF0ZV81cCA8LSBzdWJzZXEoc2Vjb25kX2NhbmRpZGF0ZV9jaHIsIHNlY29uZF9jYW5kaWRhdGVfc3RhcnQsIHNlY29uZF9jYW5kaWRhdGVfc2Vjb25kKQphcy5jaGFyYWN0ZXIoc2Vjb25kX2NhbmRpZGF0ZV81cCkKc2Vjb25kX2NhbmRpZGF0ZV8zcCA8LSBzcGdzOjpyZXZlcnNlQ29tcGxlbWVudChzdWJzZXEoc2Vjb25kX2NhbmRpZGF0ZV9jaHIsIHNlY29uZF9jYW5kaWRhdGVfdGhpcmQsIHNlY29uZF9jYW5kaWRhdGVfZW5kKSkKc2Vjb25kX2NhbmRpZGF0ZV8zcAoKCiMjICogTHBhTDEzLTI5IDgzMDM0MiAyNTIKdGhpcmRfY2FuZGlkYXRlX2NociA8LSBnZW5vbWVbWyJMcGFMMTNfMjkiXV0KcHJpbWVyX2xlbmd0aCA8LSAyMgphbXBsaWNvbl9sZW5ndGggPC0gMjUyCnRoaXJkX2NhbmRpZGF0ZV90aGlyZCA8LSA4MzAzNDIKdGhpcmRfY2FuZGlkYXRlX3NlY29uZCA8LSB0aGlyZF9jYW5kaWRhdGVfdGhpcmQgLSBhbXBsaWNvbl9sZW5ndGgKdGhpcmRfY2FuZGlkYXRlX3N0YXJ0IDwtIHRoaXJkX2NhbmRpZGF0ZV9zZWNvbmQgLSBwcmltZXJfbGVuZ3RoCnRoaXJkX2NhbmRpZGF0ZV9lbmQgPC0gdGhpcmRfY2FuZGlkYXRlX3RoaXJkICsgcHJpbWVyX2xlbmd0aAp0aGlyZF9jYW5kaWRhdGVfcmVnaW9uIDwtIHN1YnNlcSh0aGlyZF9jYW5kaWRhdGVfY2hyLCB0aGlyZF9jYW5kaWRhdGVfc3RhcnQsIHRoaXJkX2NhbmRpZGF0ZV9lbmQpCnRoaXJkX2NhbmRpZGF0ZV9yZWdpb24KdGhpcmRfY2FuZGlkYXRlXzVwIDwtIHN1YnNlcSh0aGlyZF9jYW5kaWRhdGVfY2hyLCB0aGlyZF9jYW5kaWRhdGVfc3RhcnQsIHRoaXJkX2NhbmRpZGF0ZV9zZWNvbmQpCmFzLmNoYXJhY3Rlcih0aGlyZF9jYW5kaWRhdGVfNXApCnRoaXJkX2NhbmRpZGF0ZV8zcCA8LSBzcGdzOjpyZXZlcnNlQ29tcGxlbWVudChzdWJzZXEodGhpcmRfY2FuZGlkYXRlX2NociwgdGhpcmRfY2FuZGlkYXRlX3RoaXJkLCB0aGlyZF9jYW5kaWRhdGVfZW5kKSkKdGhpcmRfY2FuZGlkYXRlXzNwCiMjIFlvdSBhcmUgYSBnYXJiYWdlIHBvbHlweXJpbWlkaW5lIHRyYWN0LgojIyBXaGljaCBpcyBhY3R1YWxseSBpbnRlcmVzdGluZyBpZiB0aGUgbXV0YXRpb25zIG1lc3MgaXQgdXAuCgoKIyMgKiBMcGFMMTMtMzMgMTMzMTUwNyA4NDMKZm91cnRoX2NhbmRpZGF0ZV9jaHIgPC0gZ2Vub21lW1siTHBhTDEzXzMzIl1dCnByaW1lcl9sZW5ndGggPC0gMjIKYW1wbGljb25fbGVuZ3RoIDwtIDg0Mwpmb3VydGhfY2FuZGlkYXRlX3RoaXJkIDwtIDEzMzE1MDcKZm91cnRoX2NhbmRpZGF0ZV9zZWNvbmQgPC0gZm91cnRoX2NhbmRpZGF0ZV90aGlyZCAtIGFtcGxpY29uX2xlbmd0aApmb3VydGhfY2FuZGlkYXRlX3N0YXJ0IDwtIGZvdXJ0aF9jYW5kaWRhdGVfc2Vjb25kIC0gcHJpbWVyX2xlbmd0aApmb3VydGhfY2FuZGlkYXRlX2VuZCA8LSBmb3VydGhfY2FuZGlkYXRlX3RoaXJkICsgcHJpbWVyX2xlbmd0aApmb3VydGhfY2FuZGlkYXRlX3JlZ2lvbiA8LSBzdWJzZXEoZm91cnRoX2NhbmRpZGF0ZV9jaHIsIGZvdXJ0aF9jYW5kaWRhdGVfc3RhcnQsIGZvdXJ0aF9jYW5kaWRhdGVfZW5kKQpmb3VydGhfY2FuZGlkYXRlX3JlZ2lvbgpmb3VydGhfY2FuZGlkYXRlXzVwIDwtIHN1YnNlcShmb3VydGhfY2FuZGlkYXRlX2NociwgZm91cnRoX2NhbmRpZGF0ZV9zdGFydCwgZm91cnRoX2NhbmRpZGF0ZV9zZWNvbmQpCmFzLmNoYXJhY3Rlcihmb3VydGhfY2FuZGlkYXRlXzVwKQpmb3VydGhfY2FuZGlkYXRlXzNwIDwtIHNwZ3M6OnJldmVyc2VDb21wbGVtZW50KHN1YnNlcShmb3VydGhfY2FuZGlkYXRlX2NociwgZm91cnRoX2NhbmRpZGF0ZV90aGlyZCwgZm91cnRoX2NhbmRpZGF0ZV9lbmQpKQpmb3VydGhfY2FuZGlkYXRlXzNwCmBgYAoKIyMgR28gaHVudGluZyBmb3IgU2FuZ2VyIHNlcXVlbmNpbmcgcmVnaW9ucwoKSSBtYWRlIGEgZnVuIGxpdHRsZSBmdW5jdGlvbiB3aGljaCBzaG91bGQgZmluZCByZWdpb25zIHdoaWNoIGhhdmUgbG90cyBvZiB2YXJpYW50cwphc3NvY2lhdGVkIHdpdGggYSBnaXZlbiBleHBlcmltZW50YWwgZmFjdG9yLgoKYGBge3Igc2FuZ2VyX2Z1biwgZXZhbD1GQUxTRX0KcGhlbm8gPC0gc3Vic2V0X2V4cHQobHBfZXhwdCwgc3Vic2V0ID0gImNvbmRpdGlvbj09J3oyLjInfGNvbmRpdGlvbj09J3oyLjMnIikKcGhlbm8gPC0gc3Vic2V0X2V4cHQocGhlbm8sIHN1YnNldCA9ICIhaXMubmEocERhdGEocGhlbm8pW1snYmNmdGFibGUnXV0pIikKcGhlbm9fc25wcyA8LSBzbShjb3VudF9leHB0X3NucHMocGhlbm8sIGFubm90X2NvbHVtbiA9ICJiY2Z0YWJsZSIpKQoKZnVuX3N0dWZmIDwtIHNucF9kZW5zaXR5X3ByaW1lcnMocGhlbm9fc25wcywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgYnNnZW5vbWU9IkJTR2Vub21lLkxlaXNobWFuaWEucGFuYW1lbnNpcy5NSE9NQ09MODFMMTMudjUzIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgZ2ZmPSJyZWZlcmVuY2UvVHJpVHJ5cERCLTUzX0xwYW5hbWVuc2lzTUhPTUNPTDgxTDEzLmdmZiIpCmRyb3Bfc2NhZmZvbGRzIDwtIGdyZXBsKHggPSByb3duYW1lcyhmdW5fc3R1ZmYkZmF2b3JpdGVzKSwgcGF0dGVybiA9ICJTQ0FGIikKZmF2b3JpdGVfcHJpbWVyX3JlZ2lvbnMgPC0gZnVuX3N0dWZmW1siZmF2b3JpdGVzIl1dWyFkcm9wX3NjYWZmb2xkcywgXQpmYXZvcml0ZV9wcmltZXJfcmVnaW9uc1tbImJpbiJdXSA8LSByb3duYW1lcyhmYXZvcml0ZV9wcmltZXJfcmVnaW9ucykKbGlicmFyeShkcGx5cikKZmF2b3JpdGVfcHJpbWVyX3JlZ2lvbnMgPC0gZmF2b3JpdGVfcHJpbWVyX3JlZ2lvbnMgJT4lCiAgcmVsb2NhdGUoYmluKQpgYGAKCiMjIENvbWJpbmUgdGhpcyB0YWJsZSB3aXRoIDIuMi8yLjMgZ2VuZXMKCkhlcmUgaXMgbXkgbm90ZSBmcm9tIG91ciBtZWV0aW5nOgoKQ3Jvc3MgcmVmZXJlbmNlIHByaW1lcnMgdG8gREUgZ2VuZXMgb2YgMi4yLzIuMyBhbmQvb3IgcmVzaXN0YW5jZS9zdXNjcGV0aWJsZSwKYWRkIGEgY29sdW1uIHRvIHRoZSBwcmltZXIgc3ByZWFkc2hlZXQgd2l0aCB0aGUgREUgZ2VuZXMgKGluIHJldHJvc3BlY3QgSSBhbSBndWVzc2luZwp0aGlzIGFjdHVhbGx5IG1lYW5zIHRvIHB1dCB0aGUgbG9nRkMgYXMgYSBjb2x1bW4uCgpPbmUgbmljZSB0aGluZywgSSBkaWQgYSBzZW1hbnRpYyByZW1vdmFsIG9uIHRoZSBscF9leHB0LCBzbyB0aGUgc2V0IG9mIGxvZ0ZDL3B2YWx1ZXMKc2hvdWxkIG5vdCBoYXZlIGFueSBvZiB0aGUgb2ZmZW5kaW5nIHR5cGVzOyB0aHVzIEkgc2hvdWxkIGJlIGFibGUgdG8gYXV0b21hZ2ljYWxseQpnZXQgcmlkIG9mIHRoZW0gaW4gdGhlIG1lcmdlLgoKYGBge3IgeHJlZl9wcmltZXJzX2RlZ30KbG9nZmMgPC0genlfdGFibGVbWyJkYXRhIl1dW1siejIzX3ZzX3oyMiJdXQpsb2dmY19jb2x1bW5zIDwtIGxvZ2ZjWywgYygiZGVzZXFfbG9nZmMiLCAiZGVzZXFfYWRqcCIpXQpjb2xuYW1lcyhsb2dmY19jb2x1bW5zKSA8LSBjKCJ6MjNfbG9nZmMiLCAiejIzX2FkanAiKQpuZXdfdGFibGUgPC0gbWVyZ2UoZmF2b3JpdGVfcHJpbWVyX3JlZ2lvbnMsIGxvZ2ZjX2NvbHVtbnMsCiAgICAgICAgICAgICAgICAgICBieS54ID0gImNsb3Nlc3RfZ2VuZV9iZWZvcmVfaWQiLCBieS55ID0gInJvdy5uYW1lcyIpCnN1cyA8LSBzdXNfdGFibGVbWyJkYXRhIl1dW1sic2Vuc2l0aXZlX3ZzX3Jlc2lzdGFudCJdXQpzdXNfY29sdW1ucyA8LSBzdXNbLCBjKCJkZXNlcV9sb2dmYyIsICJkZXNlcV9hZGpwIildCmNvbG5hbWVzKHN1c19jb2x1bW5zKSA8LSBjKCJzdXNfbG9nZmMiLCAic3VzX2FkanAiKQpuZXdfdGFibGUgPC0gbWVyZ2UobmV3X3RhYmxlLCBzdXNfY29sdW1ucywKICAgICAgICAgICAgICAgICAgIGJ5LnggPSAiY2xvc2VzdF9nZW5lX2JlZm9yZV9pZCIsIGJ5LnkgPSAicm93Lm5hbWVzIikgJT4lCiAgcmVsb2NhdGUoYmluKQp3cml0dGVuIDwtIHdyaXRlX3hsc3goZGF0YT1uZXdfdGFibGUsCiAgICAgICAgICAgICAgICAgICAgICBleGNlbD0iZXhjZWwvZmF2b3JpdGVfcHJpbWVyc194cmVmX3p5X3N1cy4KYGBgCgoKIyMgTWFrZSBhIGhlYXRtYXAgZGVzY3JpYmluZyB0aGUgY2x1c3RlcmluZyBvZiB2YXJpYW50cwoKV2UgY2FuIGNyb3NzIHJlZmVyZW5jZSB0aGUgdmFyaWFudHMgYWdhaW5zdCB0aGUgenltb2RlbWUgc3RhdHVzIGFuZApwbG90IGEgaGVhdG1hcCBvZiB0aGUgcmVzdWx0cyBhbmQgaG9wZWZ1bGx5IHNlZSBob3cgdGhleSBzZXBhcmF0ZS4KCmBgYHtyIHp5bW9faGVhdG1hcHN9CiMjIHBydW5lZF9zbnBzIDwtIHN1YnNldF9leHB0KG5ld19zbnBzLCBzdWJzZXQ9ImNvbmRpdGlvbj09J3oyLjInfGNvbmRpdGlvbj09J3oyLjMnIikKc25wX2dlbmVzIDwtIHNtKHNucHNfdnNfZ2VuZXMobHBfZXhwdCwgbmV3X3NldHMsIGV4cHRfbmFtZV9jb2wgPSAiY2hyb21vc29tZSIpKQoKIyNuZXdfenltb19ub3JtIDwtIG5vcm1hbGl6ZV9leHB0KHBydW5lZF9zbnBzLCBmaWx0ZXIgPSBUUlVFLCBjb252ZXJ0ID0gImNwbSIsIG5vcm0gPSAicXVhbnQiLCB0cmFuc2Zvcm0gPSBUUlVFKQojI25ld196eW1vX25vcm0gPC0gc2V0X2V4cHRfY29uZGl0aW9ucyhuZXdfenltb19ub3JtLCBmYWN0ID0gInp5bW9kZW1lY2F0ZWdvcmljYWwiKQpjbGluaWNhbF9jb2xvcnNfdjIgPC0gbGlzdCgKICAgICJ6MjIiID0gIiMwMDAwY2MiLAogICAgInoyMyIgPSAiI2NjMDAwMCIpCm5ld196eW1vX25vcm0gPC0gbm9ybWFsaXplX2V4cHQocHJ1bmVkX3NucHMsIGZpbHRlciA9IFRSVUUsIGNvbnZlcnQgPSAiY3BtIiwgbm9ybSA9ICJxdWFudCIsIHRyYW5zZm9ybSA9IFRSVUUpICU+JQogIHNldF9leHB0X2NvbmRpdGlvbnMoZmFjdCA9ICJ6eW1vZGVtZWNhdGVnb3JpY2FsIikgJT4lCiAgc2V0X2V4cHRfY29sb3JzKGNsaW5pY2FsX2NvbG9yc192MikKCnp5bW9faGVhdCA8LSBwbG90X2Rpc2hlYXQobmV3X3p5bW9fbm9ybSkKcHAoZmlsZSA9ICJpbWFnZXMvb25seXoyMl96MjNfc25wX2hlYXRtYXAucG5nIiwgaW1hZ2U9enltb19oZWF0W1sicGxvdCJdXSkKenltb19oZWF0W1sicGxvdCJdXQpgYGAKCiMjIyBBbm5vdGF0ZWQgaGVhdG1hcCBvZiB2YXJpYW50cwoKTm93IGxldCB1cyB0cnkgdG8gbWFrZSBhIGhlYXRtYXAgd2hpY2ggaW5jbHVkZXMgc29tZSBvZiB0aGUgYW5ub3RhdGlvbiBkYXRhLgoKYGBge3Igenltb19oZWF0X3BhbmVsX2dlbmVzfQpkZXMgPC0gYm90aF9ub3JtW1siZGVzaWduIl1dCnVuZGVmX2lkeCA8LSBpcy5uYShkZXNbWyJzdHJhaW4iXV0pCmRlc1t1bmRlZl9pZHgsICJzdHJhaW4iXSA8LSAidW5rbm93biIKCiMjaG1jb2xzIDwtIGNvbG9yUmFtcFBhbGV0dGUoYygieWVsbG93IiwiYmxhY2siLCJkYXJrYmx1ZSIpKSgyNTYpCmNvcnJlbGF0aW9ucyA8LSBocGdsX2NvcihleHBycyhib3RoX25vcm0pKQoKenltb19taXNzaW5nX2lkeCA8LSBpcy5uYShkZXNbWyJ6eW1vZGVtZWNhdGVnb3JpY2FsIl1dKQpkZXNbWyJ6eW1vZGVtZWNhdGVnb3JpY2FsIl1dIDwtIGFzLmNoYXJhY3RlcihkZXNbWyJjeW1vZGVtZWNhdGVnb3JpY2FsIl1dKQpkZXNbWyJjbGluaWNhbGNhdGVnb3JpY2FsIl1dIDwtIGFzLmNoYXJhY3RlcihkZXNbWyJjbGluaWNhbGNhdGVnb3JpY2FsIl1dKQpkZXNbenltb19taXNzaW5nX2lkeCwgInp5bW9kZW1lY2F0ZWdvcmljYWwiXSA8LSAidW5rbm93biIKbXlkZW5kcm8gPC0gbGlzdCgKICAiY2x1c3RmdW4iID0gaGNsdXN0LAogICJsd2QiID0gMi4wKQpjb2xfZGF0YSA8LSBhcy5kYXRhLmZyYW1lKGRlc1ssIGMoInp5bW9kZW1lY2F0ZWdvcmljYWwiLCAiY2xpbmljYWxjYXRlZ29yaWNhbCIpXSkKCnVua25vd25fY2xpbmljYWwgPC0gaXMubmEoY29sX2RhdGFbWyJjbGluaWNhbGNhdGVnb3JpY2FsIl1dKQpyb3dfZGF0YSA8LSBhcy5kYXRhLmZyYW1lKGRlc1ssIGMoInN0cmFpbiIpXSkKY29sbmFtZXMoY29sX2RhdGEpIDwtIGMoInp5bW9kZW1lIiwgIm91dGNvbWUiKQpjb2xfZGF0YVt1bmtub3duX2NsaW5pY2FsLCAib3V0Y29tZSJdIDwtICJ1bmRlZmluZWQiCgpjb2xuYW1lcyhyb3dfZGF0YSkgPC0gYygic3RyYWluIikKbXlhbm5vdCA8LSBsaXN0KAogICJDb2wiID0gbGlzdCgiZGF0YSIgPSBjb2xfZGF0YSksCiAgIlJvdyIgPSBsaXN0KCJkYXRhIiA9IHJvd19kYXRhKSkKbXljbHVzdCA8LSBsaXN0KCJjdXRoIiA9IDEuMCwKICAgICAgICAgICAgICAgICJjb2wiID0gQnJld2VyQ2x1c3RlckNvbCkKbXlsYWJzIDwtIGxpc3QoCiAgIlJvdyIgPSBsaXN0KCJucm93IiA9IDQpLAogICJDb2wiID0gbGlzdCgibnJvdyIgPSA0KSkKaG1jb2xzIDwtIGNvbG9yUmFtcFBhbGV0dGUoYygiZGFya2JsdWUiLCAiYmVpZ2UiKSkoMjQwKQptYXAxIDwtIGFubkhlYXRtYXAyKAogIGNvcnJlbGF0aW9ucywKICBkZW5kcm9ncmFtID0gbXlkZW5kcm8sCiAgYW5ub3RhdGlvbiA9IG15YW5ub3QsCiAgY2x1c3RlciA9IG15Y2x1c3QsCiAgbGFiZWxzID0gbXlsYWJzLAogICMjIFRoZSBmb2xsb3dpbmcgY29udHJvbHMgaWYgdGhlIHBpY3R1cmUgaXMgc3ltbWV0cmljCiAgc2NhbGUgPSAibm9uZSIsCiAgY29sID0gaG1jb2xzKQpwcChmaWxlID0gImltYWdlcy9kZW5kcm9faGVhdG1hcC5wbmciLCBpbWFnZSA9IG1hcDEsIGhlaWdodCA9IDIwLCB3aWR0aCA9IDIwKQpgYGAKClByaW50IHRoZSBsYXJnZXIgaGVhdG1hcCBzbyB0aGF0IGFsbCB0aGUgbGFiZWxzIGFwcGVhci4gIEtlZXAgaW4gbWluZAp0aGF0IGFzIHdlIGdldCBtb3JlIHNhbXBsZXMsIHRoaXMgaW1hZ2UgbmVlZHMgdG8gY29udGludWUgZ2V0dGluZwpiaWdnZXIuCgohW2JpZyBoZWF0bWFwXShpbWFnZXMvZGVuZHJvX2hlYXRtYXAucG5nKQoKCmBgYHtyIHRoZXJlc2FfaWRlYX0KeHJlZl9wcm9wIDwtIHRhYmxlKHBoZW5vX3NucHNbWyJjb25kaXRpb25zIl1dKQpwaGVub19zbnBzJGNvbmRpdGlvbnMKaWR4X3RibCA8LSBleHBycyhwaGVub19zbnBzKSA+IDUKbmV3X3RibCA8LSBkYXRhLmZyYW1lKHJvdy5uYW1lcyA9IHJvd25hbWVzKGV4cHJzKHBoZW5vX3NucHMpKSkKZm9yIChuIGluIG5hbWVzKHhyZWZfcHJvcCkpIHsKICBuZXdfdGJsW1tuXV0gPC0gMAogIGlkeF9jb2xzIDwtIHdoaWNoKHBoZW5vX3NucHNbWyJjb25kaXRpb25zIl1dID09IG4pCiAgcHJvcF9jb2wgPC0gcm93U3VtcyhpZHhfdGJsWywgaWR4X2NvbHNdKSAvIHhyZWZfcHJvcFtuXQogIG5ld190Ymxbbl0gPC0gcHJvcF9jb2wKfQprZWVwZXJzIDwtIGdyZXBsKHggPSByb3duYW1lcyhuZXdfdGJsKSwgcGF0dGVybiA9ICJMcGFMMTMiKQpuZXdfdGJsIDwtIG5ld190Ymxba2VlcGVycywgXQpuZXdfdGJsW1sic3Ryb25nMjIiXV0gPC0gMS4wMDEgLSBuZXdfdGJsW1siejIuMiJdXQpuZXdfdGJsW1sic3Ryb25nMjMiXV0gPC0gMS4wMDEgLSBuZXdfdGJsW1siejIuMyJdXQpzMjJfbmEgPC0gbmV3X3RibFtbInN0cm9uZzIyIl1dID4gMQpuZXdfdGJsW3MyMl9uYSwgInN0cm9uZzIyIl0gPC0gMQpzMjNfbmEgPC0gbmV3X3RibFtbInN0cm9uZzIzIl1dID4gMQpuZXdfdGJsW3MyM19uYSwgInN0cm9uZzIzIl0gPC0gMQoKbmV3X3RibFtbIlNOUCJdXSA8LSByb3duYW1lcyhuZXdfdGJsKQpuZXdfdGJsW1siQ2hyb21vc29tZSJdXSA8LSBnc3ViKHggPSBuZXdfdGJsW1siU05QIl1dLCBwYXR0ZXJuID0gImNocl8oLiopX3Bvc18uKiIsIHJlcGxhY2VtZW50ID0gIlxcMSIpCm5ld190YmxbWyJQb3NpdGlvbiJdXSA8LSBnc3ViKHggPSBuZXdfdGJsW1siU05QIl1dLCBwYXR0ZXJuID0gIi4qX3Bvc18oXFxkKylfLioiLCByZXBsYWNlbWVudCA9ICJcXDEiKQpuZXdfdGJsIDwtIG5ld190YmxbLCBjKCJTTlAiLCAiQ2hyb21vc29tZSIsICJQb3NpdGlvbiIsICJzdHJvbmcyMiIsICJzdHJvbmcyMyIpXQoKCmxpYnJhcnkoQ01wbG90KQpzaW1wbGlmeSA8LSBuZXdfdGJsCnNpbXBsaWZ5W1sic3Ryb25nMjIiXV0gPC0gTlVMTAoKCkNNcGxvdChzaW1wbGlmeSwgYmluLnNpemUgPSAxMDAwMDApCgpDTXBsb3QobmV3X3RibCwgcGxvdC50eXBlPSJtIiwgbXVsdHJhY2tzPVRSVUUsIHRocmVzaG9sZCA9IGMoMC4wMSwgMC4wNSksCiAgICAgICB0aHJlc2hvbGQubHdkPWMoMSwxKSwgdGhyZXNob2xkLmNvbD1jKCJibGFjayIsImdyZXkiKSwKICAgICAgIGFtcGxpZnk9VFJVRSwgYmluLnNpemU9MTAwMDAsCiAgICAgICBjaHIuZGVuLmNvbD1jKCJkYXJrZ3JlZW4iLCAieWVsbG93IiwgInJlZCIpLAogICAgICAgc2lnbmFsLmNvbD1jKCJyZWQiLCAiZ3JlZW4iLCAiYmx1ZSIpLAogICAgICAgc2lnbmFsLmNleD0xLCBmaWxlPSJqcGciLCBtZW1vPSIiLCBkcGk9MzAwLCBmaWxlLm91dHB1dD1UUlVFLCB2ZXJib3NlPVRSVUUpCmBgYAoKIVtTTlAgRGVuc2l0eV0oU05QLURlbnNpdHkucmF0aW8uanBnKQohW0NpcmN1bGFyIE1hbmhhdHRhbl0oQ2lyY3VsYXItTWFuaGF0dGFuLnJhdGlvLmpwZykKIVtSZWN0YW5ndWxhciBNYW5oYXR0YW5dKFJlY3Rhbmd1bGFyLU1hbmhhdHRhbi5yYXRpby5qcGcpCiFbUVFdKFFRcGxvdC5yYXRpby5qcGcpCgojIyBUcnkgb3V0IE1hdHJpeEVRVEwKClRoaXMgdG9vbCBsb29rcyBhIGxpdHRsZSBvcGFxdWUsIGJ1dCBwcm92aWRlcyBzYW1wbGUgZGF0YSB3aXRoIHRoaW5ncwp0aGF0IG1ha2Ugc2Vuc2UgdG8gbWUgYW5kIHNob3VsZCBiZSBwcmV0dHkgZWFzeSB0byByZWNhcGl0dWxhdGUgaW4gb3VyCmRhdGEuCgoxLiAgY292YXJpYXRlcy50eHQ6IENvbHVtbnMgYXJlIHNhbXBsZXMsIHJvd3MgYXJlIHRoaW5ncyBmcm9tIHBEYXRhIC0tIHRoZQogICAgbW9zdCBsaWtlbHkgb25lcyBvZiBpbnRlcmVzdCBmb3Igb3VyIGRhdGEgd291bGQgYmUgenltb2RlbWUsCiAgICBzZW5zaXRpdml0eQoyLiAgZ2VuZWxvYy50eHQ6IGNvbHVtbnMgYXJlICdnZW5laWQnLCAnY2hyJywgJ2xlZnQnLCAncmlnaHQnLiAgSQogICAgZ3Vlc3MgSSBjYW4gYXNzdW1lIGxlZnQgYW5kIHJpZ2h0IGFyZSBzdGFydC9zdG9wOyBpbiB3aGljaCBjYXNlCiAgICB0aGlzIGlzIHRyaXZpYWxseSBhY3F1aXJhYmxlIGZyb20gZkRhdGEuCjMuICBnZS50eHQ6IFRoaXMgYXBwZWFycyB0byBiZSBhIGxvZyhycGttL2NwbSkgdGFibGUgd2l0aCByb3dzIGFzIGdlbmVzIGFuZAogICAgY29sdW1ucyBhcyBzYW1wbGVzCjQuICBzbnBzbG9jLnR4dDogY29sdW1ucyBhcmUgJ3NucGlkJywgJ2NocicsICdwb3MnCjUuICBzbnBzLnR4dDogY29sdW1ucyBhcmUgc2FtcGxlcywgcm93cyBhcmUgdGhlIGlkcyBmcm9tIHNuc3Bsb2MsCiAgICB2YWx1ZXMgYSAwLDEsMi4gIEkgYXNzdW1lIDAgaXMgaWRlbnRpY2FsIGFuZCAxLi4xMiBhcmUgdGhlIHZhcmlvdXMKICAgIEEtPlRHQyBULT5BR0MgQy0+QUdUIEctPkFDVAoKYGBge3IgbWF0cml4ZXF0bCwgZXZhbD1GQUxTRX0KIyMgRm9yIHRoaXMsIGxldCB1cyB1c2UgdGhlICduZXdfc25wcycgZGF0YSBzdHJ1Y3R1cmUuCiMjIENhdmVhdCBoZXJlOiB0aGVzZSBuZWVkIHRvIGJlIGNvZXJjZWQgdG8gbnVtYmVycy4KbXlfY292YXJpYXRlcyA8LSBwRGF0YShuZXdfc25wcylbLCBjKCJ6eW1vZGVtZWNhdGVnb3JpY2FsIiwgImNsaW5pY2FsY2F0ZWdvcmljYWwiKV0KZm9yIChjb2wgaW4gY29sbmFtZXMobXlfY292YXJpYXRlcykpIHsKICBteV9jb3ZhcmlhdGVzW1tjb2xdXSA8LSBhcy5udW1lcmljKGFzLmZhY3RvcihteV9jb3ZhcmlhdGVzW1tjb2xdXSkpCn0KbXlfY292YXJpYXRlcyA8LSB0KG15X2NvdmFyaWF0ZXMpCgpteV9nZW5lbG9jIDwtIGZEYXRhKGxwX2V4cHQpWywgYygiZ2lkIiwgImNocm9tb3NvbWUiLCAic3RhcnQiLCAiZW5kIildCmNvbG5hbWVzKG15X2dlbmVsb2MpIDwtIGMoImdlbmVpZCIsICJjaHIiLCAibGVmdCIsICJyaWdodCIpCgpteV9nZSA8LSBleHBycyhub3JtYWxpemVfZXhwdChscF9leHB0LCB0cmFuc2Zvcm0gPSAibG9nMiIsIGZpbHRlciA9IFRSVUUsIGNvbnZlcnQgPSAiY3BtIikpCnVzZWRfc2FtcGxlcyA8LSB0b2xvd2VyKGNvbG5hbWVzKG15X2dlKSkgJWluJSBjb2xuYW1lcyhleHBycyhuZXdfc25wcykpCm15X2dlIDwtIG15X2dlWywgdXNlZF9zYW1wbGVzXQoKbXlfc25wc2xvYyA8LSBkYXRhLmZyYW1lKHJvd25hbWVzID0gcm93bmFtZXMoZXhwcnMobmV3X3NucHMpKSkKIyMgT2gsIGNhdmVhdCBoZXJlOiBCZWNhdXNlIG9mIHRoZSB3YXkgSSBzdG9yZWQgdGhlIGRhdGEsCiMjIEkgY291bGQgaGF2ZSBkdXBsaWNhdGUgcm93cyB3aGljaCBwcmVzdW1hYmx5IHdpbGwgbWFrZSBtYXRyaXhFUVRMIHNhZApteV9zbnBzbG9jW1siY2hyIl1dIDwtIGdzdWIocGF0dGVybiA9ICJeY2hyXyguKylfcG9zKC4rKV9yZWZfLiokIiwgcmVwbGFjZW1lbnQgPSAiXFwxIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHggPSByb3duYW1lcyhteV9zbnBzbG9jKSkKbXlfc25wc2xvY1tbInBvcyJdXSA8LSBnc3ViKHBhdHRlcm4gPSAiXmNocl8oLispX3BvcyguKylfcmVmXy4qJCIsIHJlcGxhY2VtZW50ID0gIlxcMiIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICB4ID0gcm93bmFtZXMobXlfc25wc2xvYykpCnRlc3QgPC0gZHVwbGljYXRlZChteV9zbnBzbG9jKQojIyBFYWNoIGR1cGxpY2F0ZWQgcm93IHdvdWxkIGJlIGFub3RoZXIgdmFyaWFudCBhdCB0aGF0IHBvc2l0aW9uOwojIyBzbyBpbiB0aGVvcnkgd2Ugd291bGQgZG8gYSBybGUgdG8gbnVtYmVyIHRoZW0gSSBhbSBndWVzc2luZwojIyBIb3dldmVyLCBJIGRvIG5vdCBoYXZlIGRpZmZlcmVudCB2YXJpYW50cyBzbyBJIHRoaW5rIEkgY2FuIGlnbm9yZSB0aGlzIGZvciB0aGUgbW9tZW50CiMjIGJ1dCB3aWxsIG5lZWQgdG8gbWFrZSBteSBtYXRyaXggZWl0aGVyIDAgb3IgMS4KaWYgKHN1bSh0ZXN0KSA+IDApIHsKICBtZXNzYWdlKCJUaGVyZSBhcmU6ICIsIHN1bShkdXBsaWNhdGVkKSwgIiBkdXBsaWNhdGVkIGVudHJpZXMuIikKICBrZWVwX2lkeCA8LSAhIHRlc3QKICBteV9zbnBzbG9jIDwtIG15X3NucHNsb2Nba2VlcF9pZHgsIF0KfQoKbXlfc25wcyA8LSBleHBycyhuZXdfc25wcykKb25lX2lkeCA8LSBteV9zbnBzID4gMApteV9zbnBzW29uZV9pZHhdIDwtIDEKCiMjIE9rLCBhdCB0aGlzIHBvaW50IEkgdGhpbmsgSSBoYXZlIGFsbCB0aGUgcGllY2VzIHdoaWNoIHRoaXMgbWV0aG9kIHdhbnRzLi4uCiMjIE9oLCBubyBJIGd1ZXNzIG5vdDsgaXQgYWN0dWFsbHkgd2FudHMgdGhlIGRhdGEgYXMgYSBzZXQgb2YgZmlsZW5hbWVzLi4uCmxpYnJhcnkoTWF0cml4RVFUTCkKd3JpdGUudGFibGUobXlfc25wcywgImVxdGwvc25wcy50c3YiLCBuYSA9ICJOQSIsIGNvbC5uYW1lcyA9IFRSVUUsIHJvdy5uYW1lcyA9IFRSVUUsIHNlcCA9ICJcdCIsIHF1b3RlID0gVFJVRSkKIyMgcmVhZHI6OndyaXRlX3RzdihteV9zbnBzLCAiZXF0bC9zbnBzLnRzdiIsICkKd3JpdGUudGFibGUobXlfc25wc2xvYywgImVxdGwvc25wc2xvYy50c3YiLCBuYSA9ICJOQSIsIGNvbC5uYW1lcyA9IFRSVUUsIHJvdy5uYW1lcyA9IFRSVUUsIHNlcCA9ICJcdCIsIHF1b3RlID0gVFJVRSkKIyMgcmVhZHI6OndyaXRlX3RzdihteV9zbnBzbG9jLCAiZXF0bC9zbnBzbG9jLnRzdiIpCndyaXRlLnRhYmxlKGFzLmRhdGEuZnJhbWUobXlfZ2UpLCAiZXF0bC9nZS50c3YiLCBuYSA9ICJOQSIsIGNvbC5uYW1lcyA9IFRSVUUsIHJvdy5uYW1lcyA9IFRSVUUsIHNlcCA9ICJcdCIsIHF1b3RlID0gVFJVRSkKIyMgcmVhZHI6OndyaXRlX3Rzdihhcy5kYXRhLmZyYW1lKG15X2dlKSwgImVxdGwvZ2UudHN2IikKd3JpdGUudGFibGUoYXMuZGF0YS5mcmFtZShteV9nZW5lbG9jKSwgImVxdGwvZ2VuZWxvYy50c3YiLCBuYSA9ICJOQSIsIGNvbC5uYW1lcyA9IFRSVUUsIHJvdy5uYW1lcyA9IFRSVUUsIHNlcCA9ICJcdCIsIHF1b3RlID0gVFJVRSkKIyMgcmVhZHI6OndyaXRlX3Rzdihhcy5kYXRhLmZyYW1lKG15X2dlbmVsb2MpLCAiZXF0bC9nZW5lbG9jLnRzdiIpCndyaXRlLnRhYmxlKGFzLmRhdGEuZnJhbWUobXlfY292YXJpYXRlcyksICJlcXRsL2NvdmFyaWF0ZXMudHN2IiwgbmEgPSAiTkEiLCBjb2wubmFtZXMgPSBUUlVFLCByb3cubmFtZXMgPSBUUlVFLCBzZXAgPSAiXHQiLCBxdW90ZSA9IFRSVUUpCiMjIHJlYWRyOjp3cml0ZV90c3YoYXMuZGF0YS5mcmFtZShteV9jb3ZhcmlhdGVzKSwgImVxdGwvY292YXJpYXRlcy50c3YiKQoKdXNlTW9kZWwgPSBtb2RlbExJTkVBUiAjIG1vZGVsQU5PVkEsIG1vZGVsTElORUFSLCBvciBtb2RlbExJTkVBUl9DUk9TUwoKIyBHZW5vdHlwZSBmaWxlIG5hbWUKU05QX2ZpbGVfbmFtZSA9ICJlcXRsL3NucHMudHN2IgpzbnBzX2xvY2F0aW9uX2ZpbGVfbmFtZSA9ICJlcXRsL3NucHNsb2MudHN2IgpleHByZXNzaW9uX2ZpbGVfbmFtZSA9ICJlcXRsL2dlLnRzdiIKZ2VuZV9sb2NhdGlvbl9maWxlX25hbWUgPSAiZXF0bC9nZW5lbG9jLnRzdiIKY292YXJpYXRlc19maWxlX25hbWUgPSAiZXF0bC9jb3ZhcmlhdGVzLnRzdiIKIyBPdXRwdXQgZmlsZSBuYW1lCm91dHB1dF9maWxlX25hbWVfY2lzID0gdGVtcGZpbGUoKQpvdXRwdXRfZmlsZV9uYW1lX3RyYSA9IHRlbXBmaWxlKCkKIyBPbmx5IGFzc29jaWF0aW9ucyBzaWduaWZpY2FudCBhdCB0aGlzIGxldmVsIHdpbGwgYmUgc2F2ZWQKcHZPdXRwdXRUaHJlc2hvbGRfY2lzID0gMC4xCnB2T3V0cHV0VGhyZXNob2xkX3RyYSA9IDAuMQojIEVycm9yIGNvdmFyaWFuY2UgbWF0cml4CiMgU2V0IHRvIG51bWVyaWMoKSBmb3IgaWRlbnRpdHkuCmVycm9yQ292YXJpYW5jZSA9IG51bWVyaWMoKQojIGVycm9yQ292YXJpYW5jZSA9IHJlYWQudGFibGUoIlNhbXBsZV9EYXRhL2Vycm9yQ292YXJpYW5jZS50eHQiKTsKIyBEaXN0YW5jZSBmb3IgbG9jYWwgZ2VuZS1TTlAgcGFpcnMKY2lzRGlzdCA9IDFlNgojIyBMb2FkIGdlbm90eXBlIGRhdGEKc25wcyA9IFNsaWNlZERhdGEkbmV3KCkKc25wcyRmaWxlRGVsaW1pdGVyID0gIlx0IiAgICAgICMgdGhlIFRBQiBjaGFyYWN0ZXIKc25wcyRmaWxlT21pdENoYXJhY3RlcnMgPSAiTkEiICMgZGVub3RlIG1pc3NpbmcgdmFsdWVzOwpzbnBzJGZpbGVTa2lwUm93cyA9IDEgICAgICAgICAgIyBvbmUgcm93IG9mIGNvbHVtbiBsYWJlbHMKc25wcyRmaWxlU2tpcENvbHVtbnMgPSAxICAgICAgICMgb25lIGNvbHVtbiBvZiByb3cgbGFiZWxzCnNucHMkZmlsZVNsaWNlU2l6ZSA9IDIwMDAgICAgICAjIHJlYWQgZmlsZSBpbiBzbGljZXMgb2YgMiwwMDAgcm93cwpzbnBzJExvYWRGaWxlKFNOUF9maWxlX25hbWUpCiMjIExvYWQgZ2VuZSBleHByZXNzaW9uIGRhdGEKZ2VuZSA9IFNsaWNlZERhdGEkbmV3KCkKZ2VuZSRmaWxlRGVsaW1pdGVyID0gIlx0IiAgICAgICMgdGhlIFRBQiBjaGFyYWN0ZXIKZ2VuZSRmaWxlT21pdENoYXJhY3RlcnMgPSAiTkEiICMgZGVub3RlIG1pc3NpbmcgdmFsdWVzOwpnZW5lJGZpbGVTa2lwUm93cyA9IDEgICAgICAgICAgIyBvbmUgcm93IG9mIGNvbHVtbiBsYWJlbHMKZ2VuZSRmaWxlU2tpcENvbHVtbnMgPSAxICAgICAgICMgb25lIGNvbHVtbiBvZiByb3cgbGFiZWxzCmdlbmUkZmlsZVNsaWNlU2l6ZSA9IDIwMDAgICAgICAjIHJlYWQgZmlsZSBpbiBzbGljZXMgb2YgMiwwMDAgcm93cwpnZW5lJExvYWRGaWxlKGV4cHJlc3Npb25fZmlsZV9uYW1lKQojIyBMb2FkIGNvdmFyaWF0ZXMKY3ZydCA9IFNsaWNlZERhdGEkbmV3KCkKY3ZydCRmaWxlRGVsaW1pdGVyID0gIlx0IiAgICAgICMgdGhlIFRBQiBjaGFyYWN0ZXIKY3ZydCRmaWxlT21pdENoYXJhY3RlcnMgPSAiTkEiICMgZGVub3RlIG1pc3NpbmcgdmFsdWVzOwpjdnJ0JGZpbGVTa2lwUm93cyA9IDEgICAgICAgICAgIyBvbmUgcm93IG9mIGNvbHVtbiBsYWJlbHMKY3ZydCRmaWxlU2tpcENvbHVtbnMgPSAxICAgICAgICMgb25lIGNvbHVtbiBvZiByb3cgbGFiZWxzCmlmKGxlbmd0aChjb3ZhcmlhdGVzX2ZpbGVfbmFtZSkgPiAwKSB7CiAgY3ZydCRMb2FkRmlsZShjb3ZhcmlhdGVzX2ZpbGVfbmFtZSkKfQojIyBSdW4gdGhlIGFuYWx5c2lzCnNucHNwb3MgPSByZWFkLnRhYmxlKHNucHNfbG9jYXRpb25fZmlsZV9uYW1lLCBoZWFkZXIgPSBUUlVFLCBzdHJpbmdzQXNGYWN0b3JzID0gRkFMU0UpCmdlbmVwb3MgPSByZWFkLnRhYmxlKGdlbmVfbG9jYXRpb25fZmlsZV9uYW1lLCBoZWFkZXIgPSBUUlVFLCBzdHJpbmdzQXNGYWN0b3JzID0gRkFMU0UpCgptZSA9IE1hdHJpeF9lUVRMX21haW4oCiAgICBzbnBzID0gc25wcywKICAgIGdlbmUgPSBnZW5lLAogICAgY3ZydCA9IGN2cnQsCiAgICBvdXRwdXRfZmlsZV9uYW1lID0gb3V0cHV0X2ZpbGVfbmFtZV90cmEsCiAgICBwdk91dHB1dFRocmVzaG9sZCA9IHB2T3V0cHV0VGhyZXNob2xkX3RyYSwKICAgIHVzZU1vZGVsID0gdXNlTW9kZWwsCiAgICBlcnJvckNvdmFyaWFuY2UgPSBlcnJvckNvdmFyaWFuY2UsCiAgICB2ZXJib3NlID0gVFJVRSwKICAgIG91dHB1dF9maWxlX25hbWUuY2lzID0gb3V0cHV0X2ZpbGVfbmFtZV9jaXMsCiAgICBwdk91dHB1dFRocmVzaG9sZC5jaXMgPSBwdk91dHB1dFRocmVzaG9sZF9jaXMsCiAgICBzbnBzcG9zID0gc25wc3BvcywKICAgIGdlbmVwb3MgPSBnZW5lcG9zLAogICAgY2lzRGlzdCA9IGNpc0Rpc3QsCiAgICBwdmFsdWUuaGlzdCA9ICJxcXBsb3QiLAogICAgbWluLnB2LmJ5LmdlbmVzbnAgPSBGQUxTRSwKICAgIG5vRkRSc2F2ZU1lbW9yeSA9IEZBTFNFKTsKYGBgCgoKCmBgYHtyIHNhdmVtZX0KaWYgKCFpc1RSVUUoZ2V0MCgic2tpcF9sb2FkIikpKSB7CiAgcGFuZGVyOjpwYW5kZXIoc2Vzc2lvbkluZm8oKSkKICBtZXNzYWdlKHBhc3RlMCgiVGhpcyBpcyBocGdsdG9vbHMgY29tbWl0OiAiLCBnZXRfZ2l0X2NvbW1pdCgpKSkKICBtZXNzYWdlKHBhc3RlMCgiU2F2aW5nIHRvICIsIHNhdmVmaWxlKSkKICB0bXAgPC0gc20oc2F2ZW1lKGZpbGVuYW1lID0gc2F2ZWZpbGUpKQp9CmBgYAoKYGBge3IgbG9hZG1lX2FmdGVyLCBldmFsID0gRkFMU0V9CnRtcCA8LSBsb2FkbWUoZmlsZW5hbWUgPSBzYXZlZmlsZSkKYGBgCg==