+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
++Stat: Import data (phyloseq_import_data.py +v5.1.0) +
+## Import packages
+library(phyloseq)
+library(ggplot2)
+library(ape)
+
+## Settin variables
+ ## The ASV abundance matrix with taxonomy annotation file (biom format)
+ # biomfile <- ""
+
+ ## The sample metadata file(TSV format)
+ # samplefile <- ""
+
+ ## (optional) the ASV tree file (nwk format). Write "None" if you do not have any tree
+ # treefile <- ""
+
+ ## The ordered taxonomic levels stored in BIOM. Each level is separated by one space.
+ ## default : "Kingdom Phylum Class Order Family Genus Species"
+ # ranks <- ""
+
+ ## Do you want to normalise your data ? "True" or "False"
+ # normalisation <- ""
+
+## Create input and parameters dataframe
+ # params <- data.frame( "biomfile" = biomfile, "samplefile" = samplefile, "tree" = tree, "ranks" = ranks, "normalisation" = normalisation)
+
+
+## Import data
+biomfile <- params$biomfile
+data <- import_biom(biomfile)
+sampledata <- read.csv(params$samplefile, sep = "\t", row.names = 1, check.names = FALSE)
+
+# if taxonomy starts with k__ it means that its Greengenes like format
+# import need to be done using parse_taxonomy_greengenes function
+# in this case user taxonomic rank names are ignored
+tax <- tax_table(data)[[1]]
+if ((gregexpr('k__', tax))[[1]][1]>0) {
+ cat("Warning : Taxonomic affiliations come from Greengenes database, user specified ranks names are ignored.")
+ data <- import_biom(biomfile, parseFunction = parse_taxonomy_greengenes)
+} else {
+## else, custumize rank name with the user specified ranks variable
+ new_rank <- as.list(strsplit(params$ranks, " ")[[1]])
+ colnames(tax_table(data)) <- new_rank
+}Warning : Taxonomic affiliations come from Greengenes database, user specified ranks names are ignored.## add sample name to metadata, as SampleID variable
+sampledata$SampleID <- rownames(sampledata)
+sample_data(data) <- sampledata
+
+## add tree metadata if available
+if (params$treefile != "None"){
+ treefile <- read.tree(params$treefile)
+ phy_tree(data) <- treefile
+}
+
+## change de sample metadata order as in input samplefile
+for ( variable in sample_variables(data)){
+ variable.order = as.vector(unique(sampledata[,variable]))
+ sample_data(data)[,variable] <- factor(get_variable(data, variable),levels=variable.order)
+}
+
+## remove empty samples
+empty_samples <- sample_names(data)[which(sample_sums(data)==0)]
+sample_to_keep <- sample_names(data)[which(sample_sums(data)>0)]
+data <- prune_samples(sample_to_keep, data)
+
+empty_taxa <- taxa_names(data)[which(taxa_sums(data)==0)]
+taxa_to_keep <- taxa_names(data)[which(taxa_sums(data)>0)]
+data <- prune_taxa(taxa_to_keep, data)
+
+## abundance normalisation
+if(params$normalisation){ data <- rarefy_even_depth(data, rngseed = 1121983)}
+
+## save phyloseq object in Rdata file
+save(data, file=params$outputRdata)
+
+
+
+
+
+
+
+
+
+
+Summary
+ + + +dataphyloseq-class experiment-level object
+otu_table() OTU Table: [ 507 taxa and 64 samples ]
+sample_data() Sample Data: [ 64 samples by 4 sample variables ]
+tax_table() Taxonomy Table: [ 507 taxa by 7 taxonomic ranks ]
+phy_tree() Phylogenetic Tree: [ 507 tips and 506 internal nodes ]if(length(empty_samples) > 0) {cat(paste('Remove empty samples: ', paste(empty_samples, collapse=",")))}
+if(length(empty_taxa) > 0) {cat(paste('Remove empty taxa: ', paste(empty_taxa, collapse=",")))}Remove empty taxa: otu_01781if(params$normalisation){cat(paste('Number of sequences in each sample after normalisation: ', head(sample_sums(data))[[1]]))}
+
+Ranks Names
+ + + +cat(paste('Rank names : ',paste(rank_names(data),collapse=', ')))Rank names : Kingdom, Phylum, Class, Order, Family, Genus, Species
+
+Sample metadata
+ + + +variables <- sample_variables(data)
+cat(paste('Sample variables: ', paste(variables, collapse=', ' )))Sample variables: EnvType, Description, FoodType, SampleIDfor (var in variables){
+ cat(paste(var,': ',paste(levels(factor(get_variable(data, varName = var))),collapse=', '), '\n\n'))
+}EnvType : BoeufHache, VeauHache, DesLardons, MerguezVolaille, SaumonFume, FiletSaumon, FiletCabillaud, Crevette
+
+Description : LOT1, LOT3, LOT4, LOT5, LOT6, LOT7, LOT8, LOT10, LOT2, LOT9
+
+FoodType : Meat, Seafood
+
+SampleID : BHT0.LOT01, BHT0.LOT03, BHT0.LOT04, BHT0.LOT05, BHT0.LOT06, BHT0.LOT07, BHT0.LOT08, BHT0.LOT10, VHT0.LOT01, VHT0.LOT02, VHT0.LOT03, VHT0.LOT04, VHT0.LOT06, VHT0.LOT07, VHT0.LOT08, VHT0.LOT10, DLT0.LOT01, DLT0.LOT03, DLT0.LOT04, DLT0.LOT05, DLT0.LOT06, DLT0.LOT07, DLT0.LOT08, DLT0.LOT10, MVT0.LOT01, MVT0.LOT03, MVT0.LOT05, MVT0.LOT06, MVT0.LOT07, MVT0.LOT08, MVT0.LOT09, MVT0.LOT10, SFT0.LOT01, SFT0.LOT02, SFT0.LOT03, SFT0.LOT04, SFT0.LOT05, SFT0.LOT06, SFT0.LOT07, SFT0.LOT08, FST0.LOT01, FST0.LOT02, FST0.LOT03, FST0.LOT05, FST0.LOT06, FST0.LOT07, FST0.LOT08, FST0.LOT10, FCT0.LOT01, FCT0.LOT02, FCT0.LOT03, FCT0.LOT05, FCT0.LOT06, FCT0.LOT07, FCT0.LOT08, FCT0.LOT10, CDT0.LOT02, CDT0.LOT04, CDT0.LOT05, CDT0.LOT06, CDT0.LOT07, CDT0.LOT08, CDT0.LOT09, CDT0.LOT10
+
+Plot tree
+ + + +if(params$treefile!="None"){
+ p <- plot_tree(data, color=rank_names(data)[2]) +
+ ggtitle(paste("Phylogenetic tree colored by", rank_names(data)[2])) +
+ theme(plot.title = element_text(hjust = 0.5))
+ plot(p)
+}
if(params$treefile=="None"){
+ cat("There is no phylogenetic tree in the object you have provided.")
+}
+
+Reproducibility token
+ + + +sessioninfo::session_info()─ Session info ───────────────────────────────────────────────────────────────
+ setting value
+ version R version 4.1.2 (2021-11-01)
+ os Ubuntu 24.04.2 LTS
+ system x86_64, linux-gnu
+ ui X11
+ language fr_FR:en
+ collate en_US.utf8
+ ctype en_US.utf8
+ tz Europe/Paris
+ date 2026-01-14
+ pandoc 2.19.2 @ /home/maria/miniforge3/envs/frogs@5.1.0/bin/ (via rmarkdown)
+
+─ Packages ───────────────────────────────────────────────────────────────────
+ package * version date (UTC) lib source
+ ade4 1.7-22 2023-02-06 [1] CRAN (R 4.1.3)
+ ape * 5.7-1 2023-03-13 [1] CRAN (R 4.1.3)
+ Biobase 2.54.0 2021-10-26 [1] Bioconductor
+ BiocGenerics 0.40.0 2021-10-26 [1] Bioconductor
+ biomformat 1.22.0 2021-10-26 [1] Bioconductor
+ Biostrings 2.62.0 2021-10-26 [1] Bioconductor
+ bitops 1.0-7 2021-04-24 [1] CRAN (R 4.1.3)
+ bslib 0.5.0 2023-06-09 [1] CRAN (R 4.1.3)
+ cachem 1.0.8 2023-05-01 [1] CRAN (R 4.1.3)
+ cli 3.6.1 2023-03-23 [1] CRAN (R 4.1.3)
+ cluster 2.1.4 2022-08-22 [1] CRAN (R 4.1.3)
+ codetools 0.2-19 2023-02-01 [1] CRAN (R 4.1.3)
+ colorspace 2.1-0 2023-01-23 [1] CRAN (R 4.1.3)
+ crayon 1.5.2 2022-09-29 [1] CRAN (R 4.1.3)
+ data.table 1.14.8 2023-02-17 [1] CRAN (R 4.1.3)
+ digest 0.6.31 2022-12-11 [1] CRAN (R 4.1.3)
+ dplyr 1.1.2 2023-04-20 [1] CRAN (R 4.1.3)
+ evaluate 0.21 2023-05-05 [1] CRAN (R 4.1.3)
+ fansi 1.0.4 2023-01-22 [1] CRAN (R 4.1.3)
+ farver 2.1.1 2022-07-06 [1] CRAN (R 4.1.3)
+ fastmap 1.1.1 2023-02-24 [1] CRAN (R 4.1.3)
+ foreach 1.5.2 2022-02-02 [1] CRAN (R 4.1.3)
+ generics 0.1.3 2022-07-05 [1] CRAN (R 4.1.3)
+ GenomeInfoDb 1.30.1 2022-01-30 [1] Bioconductor
+ GenomeInfoDbData 1.2.7 2026-01-14 [1] Bioconductor
+ ggplot2 * 3.4.2 2023-04-03 [1] CRAN (R 4.1.3)
+ glue 1.6.2 2022-02-24 [1] CRAN (R 4.1.3)
+ gtable 0.3.3 2023-03-21 [1] CRAN (R 4.1.3)
+ highr 0.10 2022-12-22 [1] CRAN (R 4.1.3)
+ htmltools 0.5.5 2023-03-23 [1] CRAN (R 4.1.3)
+ igraph 1.3.5 2022-09-22 [1] CRAN (R 4.1.3)
+ IRanges 2.28.0 2021-10-26 [1] Bioconductor
+ iterators 1.0.14 2022-02-05 [1] CRAN (R 4.1.3)
+ jquerylib 0.1.4 2021-04-26 [1] CRAN (R 4.1.3)
+ jsonlite 1.8.5 2023-06-05 [1] CRAN (R 4.1.3)
+ knitr 1.43 2023-05-25 [1] CRAN (R 4.1.3)
+ labeling 0.4.2 2020-10-20 [1] CRAN (R 4.1.3)
+ lattice 0.21-8 2023-04-05 [1] CRAN (R 4.1.3)
+ lifecycle 1.0.3 2022-10-07 [1] CRAN (R 4.1.3)
+ magrittr 2.0.3 2022-03-30 [1] CRAN (R 4.1.3)
+ MASS 7.3-58.3 2023-03-07 [1] CRAN (R 4.1.3)
+ Matrix 1.5-4.1 2023-05-18 [1] CRAN (R 4.1.3)
+ mgcv 1.8-42 2023-03-02 [1] CRAN (R 4.1.3)
+ multtest 2.50.0 2021-10-26 [1] Bioconductor
+ munsell 0.5.0 2018-06-12 [1] CRAN (R 4.1.3)
+ nlme 3.1-162 2023-01-31 [1] CRAN (R 4.1.3)
+ permute 0.9-7 2022-01-27 [1] CRAN (R 4.1.3)
+ phyloseq * 1.38.0 2021-10-26 [1] Bioconductor
+ pillar 1.9.0 2023-03-22 [1] CRAN (R 4.1.3)
+ pkgconfig 2.0.3 2019-09-22 [1] CRAN (R 4.1.3)
+ plyr 1.8.8 2022-11-11 [1] CRAN (R 4.1.3)
+ R6 2.5.1 2021-08-19 [1] CRAN (R 4.1.3)
+ Rcpp 1.0.10 2023-01-22 [1] CRAN (R 4.1.3)
+ RCurl 1.98-1.12 2023-03-27 [1] CRAN (R 4.1.3)
+ reshape2 1.4.4 2020-04-09 [1] CRAN (R 4.1.3)
+ rhdf5 2.38.1 2022-03-10 [1] Bioconductor
+ rhdf5filters 1.6.0 2021-10-26 [1] Bioconductor
+ Rhdf5lib 1.16.0 2021-10-26 [1] Bioconductor
+ rlang 1.1.1 2023-04-28 [1] CRAN (R 4.1.3)
+ rmarkdown 2.22 2023-06-01 [1] CRAN (R 4.1.3)
+ rstudioapi 0.14 2022-08-22 [1] CRAN (R 4.1.3)
+ S4Vectors 0.32.4 2022-03-24 [1] Bioconductor
+ sass 0.4.6 2023-05-03 [1] CRAN (R 4.1.3)
+ scales 1.2.1 2022-08-20 [1] CRAN (R 4.1.3)
+ sessioninfo 1.2.2 2021-12-06 [1] CRAN (R 4.1.3)
+ stringi 1.7.6 2021-11-29 [1] CRAN (R 4.1.1)
+ stringr 1.5.0 2022-12-02 [1] CRAN (R 4.1.3)
+ survival 3.5-5 2023-03-12 [1] CRAN (R 4.1.3)
+ tibble 3.2.1 2023-03-20 [1] CRAN (R 4.1.3)
+ tidyselect 1.2.0 2022-10-10 [1] CRAN (R 4.1.3)
+ utf8 1.2.3 2023-01-31 [1] CRAN (R 4.1.3)
+ vctrs 0.6.2 2023-04-19 [1] CRAN (R 4.1.3)
+ vegan 2.6-4 2022-10-11 [1] CRAN (R 4.1.3)
+ withr 2.5.0 2022-03-03 [1] CRAN (R 4.1.3)
+ xfun 0.39 2023-04-20 [1] CRAN (R 4.1.3)
+ XVector 0.34.0 2021-10-26 [1] Bioconductor
+ yaml 2.3.7 2023-01-23 [1] CRAN (R 4.1.3)
+ zlibbioc 1.40.0 2021-10-26 [1] Bioconductor
+
+ [1] /home/maria/miniforge3/envs/frogs@5.1.0/lib/R/library
+
+──────────────────────────────────────────────────────────────────────────────LS0tCmNvcHlyaWdodDogIkNvcHlyaWdodCAoQykgMjAyNSBJTlJBRSIKbGljZW5zZTogIkdOVSBHZW5lcmFsIFB1YmxpYyBMaWNlbnNlIgpvdXRwdXQ6IAogIGh0bWxfbm90ZWJvb2s6CiAgICBjb2RlX2ZvbGRpbmc6IGhpZGUKcGFyYW1zOgogICMgUGF0aCBBU1YgYWJ1bmRhbmNlIGJpb20gZmlsZQogIGJpb21maWxlOiAKICAgIHZhbHVlOiB4CiAgIyBQYXRoIHRvIHNhbXBsZSBtZXRhZGF0YSBUU1YgZmlsZQogIHNhbXBsZWZpbGU6CiAgICB2YWx1ZTogeAogICMgKG9wdGlvbmFsKSBQYXRoIHRvIEFTVnMgdHJlZSBOV0sgZmlsZQogIHRyZWVmaWxlOgogICAgdmFsdWU6IHgKICAjIChvcHRpb25hbCkgbm9ybWFsaXNhdGlvbiBzdHJpbmcgb3B0aW9uIDogVFJVRSBvciBGQUxTRQogIG5vcm1hbGlzYXRpb246CiAgICB2YWx1ZTogeAogICMgUGF0aCB0byBwaHlsb3NlcSBvYmplY3Qgc3RvcmUgaW4gUmRhdGEgZmlsZQogIG91dHB1dFJkYXRhOgogICAgdmFsdWU6IHgKICAjIExpc3Qgb2YgdGF4b25vbWljIHJhbmtzIG5hbWVzCiAgcmFua3MgOiAKICAgIHZhbHVlOiB4CiAgIyBQYXRoIHRvIHRoZSBwaHlsb3NlcS1leHRlbmQgUiBmdW5jdGlvbnMgOiBodHRwczovL2dpdGh1Yi5jb20vbWFoZW5kcmEtbWFyaWFkYXNzb3UvcGh5bG9zZXEtZXh0ZW5kZWQsIHR5cGljYWxseSBzdG9yZWQgaW4gRlJPR1MvbGliL2V4dGVybmFsLWxpYgogIGxpYmRpcjoKICAgIHZhbHVlOiB4CiAgdmVyc2lvbjoKICAgIHZhbHVlOiB4CiAgdG9vbDoKICAgIHZhbHVlOiB4Ci0tLQoKPG1ldGEgaHR0cC1lcXVpdj0iQ29udGVudC1UeXBlIiBjb250ZW50PSJ0ZXh0L2h0bWw7IGNoYXJzZXQ9VVRGLTgiIC8+CjxtZXRhIG5hbWU9ImF1dGhvciIgY29udGVudD0iVGEgVGhpIE5nYW4gLSBTSUdFTkFFL0dBQkkgJiBNYXJpYSBCZXJuYXJkIC0gU0lHRU5BRS9HQUJJIiAvPgo8bWV0YSBuYW1lPSJ2ZXJzaW9uIiBjb250ZW50PSI1LjEuMCIgLz4KPG1ldGEgbmFtZT0iY29weXJpZ2h0IiBjb250ZW50PSJDb3B5cmlnaHQgKEMpIDIwMjUgSU5SQUUiIC8+CgpgYGB7ciwgZWNobz1GQUxTRSwgcmVzdWx0cz0iYXNpcyJ9CkNVUlJFTlRfRElSIDwtIGtuaXRyOjpvcHRzX2tuaXQkZ2V0KCJvdXRwdXQuZGlyIikKaWYgKGlzLm51bGwoQ1VSUkVOVF9ESVIpKSB7CiAgQ1VSUkVOVF9ESVIgPC0gZ2V0d2QoKQp9ClRIRU1FX0RJUiA8LSBub3JtYWxpemVQYXRoKAogIGZpbGUucGF0aChkaXJuYW1lKENVUlJFTlRfRElSKSwgInN0YXRpYyIpLAogIHdpbnNsYXNoID0gIi8iLAogIG11c3RXb3JrID0gRkFMU0UKKQoKanNfdGhlbWVfZmlsZSA8LSBub3JtYWxpemVQYXRoKAogIGZpbGUucGF0aChUSEVNRV9ESVIsICJqcyIsICJ0aGVtZS5qcyIpLAogIHdpbnNsYXNoID0gIi8iLAogIG11c3RXb3JrID0gRkFMU0UKKQoKbGluZXMgPC0gcmVhZExpbmVzKGpzX3RoZW1lX2ZpbGUsIHdhcm4gPSBGQUxTRSkKCnNlcCA8LSB3aGljaChncmVwbCgiXi8vIyMiLCB0cmltd3MobGluZXMpKSkKCmlmIChsZW5ndGgoc2VwKSA8IDIpIHsKICBzdG9wKCJUYWcgbm90IGZvdW5kIGluIHRoZW1lLmpzOiAnLy8jIyciKQp9Cgpqc19jb2RlIDwtIGxpbmVzWyhzZXBbMV0gKyAxKTooc2VwWzJdIC0gMSldCgpjYXQoIjxzY3JpcHQ+XG4iKQpjYXQoanNfY29kZSwgc2VwID0gIlxuIikKY2F0KCJcbjwvc2NyaXB0PiIpCmBgYAoKPGRpdiBzdHlsZT0iZGlzcGxheTpmbGV4OyBhbGlnbi1pdGVtczpjZW50ZXI7IGdhcDoxNXB4OyBtYXJnaW4tYm90dG9tOjIwcHg7Ij4KICA8aW1nIGlkPSJsb2dvIiBzcmM9ImRhdGE6aW1hZ2UvcG5nO2Jhc2U2NCwiIGFsdD0iRlJPR1MgbG9nbyIgc3R5bGU9ImhlaWdodDoxODBweDsiPgogIDxoMSBzdHlsZT0ibWFyZ2luOjA7Ij5TdGF0OiBJbXBvcnQgZGF0YSA8aT48c3BhbiBjbGFzcz0idGV4dC1hY2NlbnQiPihgciBwYXJhbXMkdG9vbGAgdmByIHBhcmFtcyR2ZXJzaW9uYCk8L3NwYW4+PC9pPjwvaDE+CjwvZGl2PgoKPHNjcmlwdD4KCiQoZnVuY3Rpb24gKCkgewogIHVwZGF0ZV9sb2dvKENVUlJFTlRfVEhFTUUpOwp9KTsKCgo8L3NjcmlwdD4KCjxzdHlsZSB0eXBlPSJ0ZXh0L2NzcyI+Cjpyb290IHsKCQkJLS1mcm9nc0NvbG9yOiAjMDNBNUE4OyAKCQkJLS1mcm9nc0NvbG9ySG92ZXI6IGhzbChmcm9tIHZhcigtLWZyb2dzQ29sb3IpIGggY2FsYyhzICsgNCkgY2FsYyhsIC0gMykgLyAxKTsKCQkJLS1mcm9nc1ByZUNvbG9yOiBoc2woZnJvbSB2YXIoLS1mcm9nc0NvbG9yKSBoIHMgbCAvIDAuMSk7CgkJfQoudGV4dC1hY2NlbnQgewogIGNvbG9yOiAjNTk1YzVmOwp9Ci5tYWluLWNvbnRhaW5lcnsKCW1heC13aWR0aDogNzUlICFpbXBvcnRhbnQ7Cn0KaDF7Cgljb2xvcjogdmFyKC0tZnJvZ3NDb2xvcik7Cn0KaDJ7Cgljb2xvcjogdmFyKC0tZnJvZ3NDb2xvcik7Cn0KYSB7Cgljb2xvcjogdmFyKC0tZnJvZ3NDb2xvcik7Cn0KYTpob3ZlcnsKCWNvbG9yOiB2YXIoLS1mcm9nc0NvbG9ySG92ZXIpOwp9Ci5uYXYtcGlsbHMgPiBsaS5hY3RpdmUgPiBhLCAubmF2LXBpbGxzID4gbGkuYWN0aXZlID4gYTpmb2N1c3sKCWNvbG9yOiAjZmZmOwoJYmFja2dyb3VuZC1jb2xvcjogdmFyKC0tZnJvZ3NDb2xvcik7Cglib3JkZXItY29sb3I6ICNkZWUyZTYgI2RlZTJlNiAjZmZmOwp9Ci5uYXYtcGlsbHMgPiBsaS5hY3RpdmUgPiBhOmhvdmVyIHsKCWJhY2tncm91bmQtY29sb3I6IHZhcigtLWZyb2dzQ29sb3JIb3Zlcik7Cn0KLmJ1dHRvbiB7CiAgICBiYWNrZ3JvdW5kLWNvbG9yOiB2YXIoLS1mcm9nc0NvbG9yKSA7CiAgICBib3JkZXIgICAgICAgICAgOiBub25lOwogICAgY29sb3IgICAgICAgICAgIDogd2hpdGU7CiAgICBwYWRkaW5nICAgICAgICAgOiA1cHggMTBweDsKICAgIHRleHQtYWxpZ24gICAgICA6IGNlbnRlcjsKICAgIHRleHQtZGVjb3JhdGlvbiA6IG5vbmU7CiAgICBkaXNwbGF5ICAgICAgICAgOiBpbmxpbmUtYmxvY2s7CiAgICBmb250LXNpemUgICAgICAgOiAxMnB4OwogICAgbWFyZ2luICAgICAgICAgIDogNHB4IDJweDsKICAgIGN1cnNvciAgICAgICAgICA6IHBvaW50ZXI7CiAgICBib3JkZXItcmFkaXVzICAgOiA4cHg7Cn0KaDQgeyAKICAgIGRpc3BsYXkgICAgICA6IGJsb2NrOwogICAgZm9udC1zaXplICAgIDogMWVtOwogICAgbWFyZ2luLXRvcCAgIDogMS4zM2VtOwogICAgbWFyZ2luLWJvdHRvbTogMS4zM2VtOwogICAgbWFyZ2luLWxlZnQgIDogMDsKICAgIG1hcmdpbi1yaWdodCA6IDA7CiAgICBmb250LXdlaWdodCAgOiBib2xkOwogICAgY29sb3IgICAgICAgIDogdmFyKC0tZnJvZ3NDb2xvcik7Cn0KY29kZS5yeyAvKiBDb2RlIGJsb2NrICovCiAgZm9udC1zaXplOiAxMXB4Owp9CnByZXsKICBmb250LXNpemU6IDExcHggOwogIGJhY2tncm91bmQtY29sb3I6IHZhcigtLWZyb2dzUHJlQ29sb3IpICFpbXBvcnRhbnQ7Cn0KI3RoZW1lY2hvaWNlewogIGJvcmRlci1yYWRpdXM6IDAuMjVlbTsKICBiYWNrZ3JvdW5kLWNvbG9yOiB3aGl0ZTsKICBib3JkZXI6IDFweCBzb2xpZCAjZGVlMmU2OwogIGZvbnQtd2VpZ2h0OiA0MDA7CiAgbGluZS1oZWlnaHQ6IDEuNTsKICBwYWRkaW5nOiAuMzc1cmVtIDIuMjVyZW0gLjM3NXJlbSAuNzVyZW07Cn0KPC9zdHlsZT4KIApgYGB7ciBwYWNrYWdlX2ltcG9ydCwgbWVzc2FnZT1GQUxTRSwgd2FybmluZz1GQUxTRX0KCiMjIEltcG9ydCBwYWNrYWdlcwpsaWJyYXJ5KHBoeWxvc2VxKQpsaWJyYXJ5KGdncGxvdDIpCmxpYnJhcnkoYXBlKQoKIyMgU2V0dGluIHZhcmlhYmxlcwogICMjIFRoZSBBU1YgYWJ1bmRhbmNlIG1hdHJpeCB3aXRoIHRheG9ub215IGFubm90YXRpb24gZmlsZSAoYmlvbSBmb3JtYXQpCiAgIyBiaW9tZmlsZSA8LSAiIgoKICAjIyBUaGUgc2FtcGxlIG1ldGFkYXRhIGZpbGUoVFNWIGZvcm1hdCkKICAjIHNhbXBsZWZpbGUgPC0gIiIKCiAgIyMgKG9wdGlvbmFsKSB0aGUgQVNWIHRyZWUgZmlsZSAobndrIGZvcm1hdCkuIFdyaXRlICJOb25lIiBpZiB5b3UgZG8gbm90IGhhdmUgYW55IHRyZWUKICAjIHRyZWVmaWxlIDwtICIiCgogICMjIFRoZSBvcmRlcmVkIHRheG9ub21pYyBsZXZlbHMgc3RvcmVkIGluIEJJT00uIEVhY2ggbGV2ZWwgaXMgc2VwYXJhdGVkIGJ5IG9uZSBzcGFjZS4KICAjIyBkZWZhdWx0IDogIktpbmdkb20gUGh5bHVtIENsYXNzIE9yZGVyIEZhbWlseSBHZW51cyBTcGVjaWVzIgogICMgcmFua3MgPC0gIiIKCiAgIyMgRG8geW91IHdhbnQgdG8gbm9ybWFsaXNlIHlvdXIgZGF0YSA/ICJUcnVlIiBvciAiRmFsc2UiCiAgIyBub3JtYWxpc2F0aW9uIDwtICIiCgojIyBDcmVhdGUgaW5wdXQgYW5kIHBhcmFtZXRlcnMgZGF0YWZyYW1lCiAgIyBwYXJhbXMgPC0gZGF0YS5mcmFtZSggImJpb21maWxlIiA9IGJpb21maWxlLCAic2FtcGxlZmlsZSIgPSBzYW1wbGVmaWxlLCAidHJlZSIgPSB0cmVlLCAicmFua3MiID0gcmFua3MsICJub3JtYWxpc2F0aW9uIiA9IG5vcm1hbGlzYXRpb24pCgoKIyMgSW1wb3J0IGRhdGEKYmlvbWZpbGUgPC0gcGFyYW1zJGJpb21maWxlCmRhdGEgICAgIDwtIGltcG9ydF9iaW9tKGJpb21maWxlKQpzYW1wbGVkYXRhIDwtIHJlYWQuY3N2KHBhcmFtcyRzYW1wbGVmaWxlLCBzZXAgPSAiXHQiLCByb3cubmFtZXMgPSAxLCBjaGVjay5uYW1lcyA9IEZBTFNFKQoKIyBpZiB0YXhvbm9teSBzdGFydHMgd2l0aCBrX18gaXQgbWVhbnMgdGhhdCBpdHMgR3JlZW5nZW5lcyBsaWtlIGZvcm1hdAojIGltcG9ydCBuZWVkIHRvIGJlIGRvbmUgdXNpbmcgcGFyc2VfdGF4b25vbXlfZ3JlZW5nZW5lcyBmdW5jdGlvbgojIGluIHRoaXMgY2FzZSB1c2VyIHRheG9ub21pYyByYW5rIG5hbWVzIGFyZSBpZ25vcmVkCnRheCAgICAgIDwtIHRheF90YWJsZShkYXRhKVtbMV1dCmlmICgoZ3JlZ2V4cHIoJ2tfXycsIHRheCkpW1sxXV1bMV0+MCkgeyAKICBjYXQoIldhcm5pbmcgOiBUYXhvbm9taWMgYWZmaWxpYXRpb25zIGNvbWUgZnJvbSBHcmVlbmdlbmVzIGRhdGFiYXNlLCB1c2VyIHNwZWNpZmllZCByYW5rcyBuYW1lcyBhcmUgaWdub3JlZC4iKQogIGRhdGEgPC0gaW1wb3J0X2Jpb20oYmlvbWZpbGUsIHBhcnNlRnVuY3Rpb24gPSBwYXJzZV90YXhvbm9teV9ncmVlbmdlbmVzKQp9IGVsc2UgewojIyBlbHNlLCBjdXN0dW1pemUgcmFuayBuYW1lIHdpdGggdGhlIHVzZXIgc3BlY2lmaWVkIHJhbmtzIHZhcmlhYmxlCiAgbmV3X3JhbmsgPC0gYXMubGlzdChzdHJzcGxpdChwYXJhbXMkcmFua3MsICIgIilbWzFdXSkKICBjb2xuYW1lcyh0YXhfdGFibGUoZGF0YSkpIDwtIG5ld19yYW5rCn0KCiMjIGFkZCBzYW1wbGUgbmFtZSB0byBtZXRhZGF0YSwgYXMgU2FtcGxlSUQgdmFyaWFibGUKc2FtcGxlZGF0YSRTYW1wbGVJRCA8LSByb3duYW1lcyhzYW1wbGVkYXRhKQpzYW1wbGVfZGF0YShkYXRhKSA8LSBzYW1wbGVkYXRhCgojIyBhZGQgdHJlZSBtZXRhZGF0YSBpZiBhdmFpbGFibGUKaWYgKHBhcmFtcyR0cmVlZmlsZSAhPSAiTm9uZSIpewogIHRyZWVmaWxlIDwtIHJlYWQudHJlZShwYXJhbXMkdHJlZWZpbGUpCiAgcGh5X3RyZWUoZGF0YSkgPC0gdHJlZWZpbGUKfQoKIyMgY2hhbmdlIGRlIHNhbXBsZSBtZXRhZGF0YSBvcmRlciBhcyBpbiBpbnB1dCBzYW1wbGVmaWxlCmZvciAoIHZhcmlhYmxlIGluIHNhbXBsZV92YXJpYWJsZXMoZGF0YSkpewogIHZhcmlhYmxlLm9yZGVyID0gYXMudmVjdG9yKHVuaXF1ZShzYW1wbGVkYXRhWyx2YXJpYWJsZV0pKQogIHNhbXBsZV9kYXRhKGRhdGEpWyx2YXJpYWJsZV0gPC0gZmFjdG9yKGdldF92YXJpYWJsZShkYXRhLCB2YXJpYWJsZSksbGV2ZWxzPXZhcmlhYmxlLm9yZGVyKQp9CgojIyByZW1vdmUgZW1wdHkgc2FtcGxlcwplbXB0eV9zYW1wbGVzIDwtIHNhbXBsZV9uYW1lcyhkYXRhKVt3aGljaChzYW1wbGVfc3VtcyhkYXRhKT09MCldCnNhbXBsZV90b19rZWVwIDwtIHNhbXBsZV9uYW1lcyhkYXRhKVt3aGljaChzYW1wbGVfc3VtcyhkYXRhKT4wKV0KZGF0YSA8LSBwcnVuZV9zYW1wbGVzKHNhbXBsZV90b19rZWVwLCBkYXRhKQoKZW1wdHlfdGF4YSA8LSB0YXhhX25hbWVzKGRhdGEpW3doaWNoKHRheGFfc3VtcyhkYXRhKT09MCldCnRheGFfdG9fa2VlcCA8LSB0YXhhX25hbWVzKGRhdGEpW3doaWNoKHRheGFfc3VtcyhkYXRhKT4wKV0KZGF0YSA8LSBwcnVuZV90YXhhKHRheGFfdG9fa2VlcCwgZGF0YSkKCiMjIGFidW5kYW5jZSBub3JtYWxpc2F0aW9uCmlmKHBhcmFtcyRub3JtYWxpc2F0aW9uKXsgZGF0YSA8LSByYXJlZnlfZXZlbl9kZXB0aChkYXRhLCBybmdzZWVkID0gMTEyMTk4Myl9CgojIyBzYXZlIHBoeWxvc2VxIG9iamVjdCBpbiBSZGF0YSBmaWxlCnNhdmUoZGF0YSwgZmlsZT1wYXJhbXMkb3V0cHV0UmRhdGEpCmBgYAoKIyB7LnRhYnNldCAudGFic2V0LWZhZGUgLnRhYnNldC1waWxsc30KCmBgYHs9aHRtbH0KPGRpdiBjbGFzcz0icm93Ij4KICA8ZGl2IHN0eWxlPSJmbG9hdDpyaWdodCI+CiAgPHNlbGVjdCBpZD0idGhlbWVjaG9pY2UiCiAgICAgICAgICBjbGFzcz0iZm9ybS1zZWxlY3QgZm9ybS1zZWxlY3Qtc20iCiAgICAgICAgICBzdHlsZT0id2lkdGg6IGF1dG87IgogICAgICAgICAgb25jaGFuZ2U9InVwZGF0ZV90aGVtZV9SbWQodGhpcy52YWx1ZSkiCiAgICAgICAgICBhcmlhLWxhYmVsPSJTd2l0Y2ggdGhlbWUiPgogICAgPG9wdGlvbiBzZWxlY3RlZCBkaXNhYmxlZCB2YWx1ZT0iIj5Td2l0Y2ggdGhlbWU8L29wdGlvbj4KICAgIDxvcHRpb24gZGlzYWJsZWQgdmFsdWU9IkRlZmF1bHRUaGVtZSI+RGVmYXVsdDwvb3B0aW9uPgogICAgPG9wdGlvbiB2YWx1ZT0iQ29yYWxUaGVtZSI+Q29yYWw8L29wdGlvbj4KICAgIDxvcHRpb24gdmFsdWU9IkdvbGRUaGVtZSI+R29sZDwvb3B0aW9uPgogICAgPG9wdGlvbiB2YWx1ZT0iU3RlZWxUaGVtZSI+U3RlZWw8L29wdGlvbj4KICA8L3NlbGVjdD4KICA8L2Rpdj4KPC9kaXY+CmBgYAoKIyMgU3VtbWFyeQpgYGB7ciBzdW1tYXJ5fQpkYXRhCmlmKGxlbmd0aChlbXB0eV9zYW1wbGVzKSA+IDApIHtjYXQocGFzdGUoJ1JlbW92ZSBlbXB0eSBzYW1wbGVzOiAnLCBwYXN0ZShlbXB0eV9zYW1wbGVzLCBjb2xsYXBzZT0iLCIpKSl9CmlmKGxlbmd0aChlbXB0eV90YXhhKSA+IDApIHtjYXQocGFzdGUoJ1JlbW92ZSBlbXB0eSB0YXhhOiAnLCBwYXN0ZShlbXB0eV90YXhhLCBjb2xsYXBzZT0iLCIpKSl9CmlmKHBhcmFtcyRub3JtYWxpc2F0aW9uKXtjYXQocGFzdGUoJ051bWJlciBvZiBzZXF1ZW5jZXMgaW4gZWFjaCBzYW1wbGUgYWZ0ZXIgbm9ybWFsaXNhdGlvbjogJywgaGVhZChzYW1wbGVfc3VtcyhkYXRhKSlbWzFdXSkpfQpgYGAKCiMjIFJhbmtzIE5hbWVzCmBgYHtyIHJhbmtzfQpjYXQocGFzdGUoJ1JhbmsgbmFtZXMgOiAnLHBhc3RlKHJhbmtfbmFtZXMoZGF0YSksY29sbGFwc2U9JywgJykpKQpgYGAKCiMjIFNhbXBsZSBtZXRhZGF0YQpgYGB7ciBzYW1wbGV9CnZhcmlhYmxlcyA8LSBzYW1wbGVfdmFyaWFibGVzKGRhdGEpCmNhdChwYXN0ZSgnU2FtcGxlIHZhcmlhYmxlczogJywgcGFzdGUodmFyaWFibGVzLCBjb2xsYXBzZT0nLCAnICkpKQoKZm9yICh2YXIgaW4gdmFyaWFibGVzKXsKICBjYXQocGFzdGUodmFyLCc6ICcscGFzdGUobGV2ZWxzKGZhY3RvcihnZXRfdmFyaWFibGUoZGF0YSwgdmFyTmFtZSA9IHZhcikpKSxjb2xsYXBzZT0nLCAnKSwgJ1xuXG4nKSkKfQpgYGAKCiMjIFBsb3QgdHJlZQpgYGB7ciB0cmVlLCBmaWcud2lkdGg9MTAsIGZpZy5oZWlnaHQ9OH0KCmlmKHBhcmFtcyR0cmVlZmlsZSE9Ik5vbmUiKXsKICBwIDwtIHBsb3RfdHJlZShkYXRhLCBjb2xvcj1yYW5rX25hbWVzKGRhdGEpWzJdKSArIAogICAgICAgICAgZ2d0aXRsZShwYXN0ZSgiUGh5bG9nZW5ldGljIHRyZWUgY29sb3JlZCBieSIsIHJhbmtfbmFtZXMoZGF0YSlbMl0pKSArIAogICAgICAgICAgdGhlbWUocGxvdC50aXRsZSA9IGVsZW1lbnRfdGV4dChoanVzdCA9IDAuNSkpCiAgcGxvdChwKQp9CgppZihwYXJhbXMkdHJlZWZpbGU9PSJOb25lIil7CiAgY2F0KCJUaGVyZSBpcyBubyBwaHlsb2dlbmV0aWMgdHJlZSBpbiB0aGUgb2JqZWN0IHlvdSBoYXZlIHByb3ZpZGVkLiIpCn0KYGBgCgojIyBSZXByb2R1Y2liaWxpdHkgdG9rZW4KCmBgYHtyIHNlc3Npb24sIGVjaG89VFJVRSwgZXZhbD1UUlVFfQpzZXNzaW9uaW5mbzo6c2Vzc2lvbl9pbmZvKCkKYGBg
+
+
+
+
