1 Rework DIA-Umpire to feed OpenSWATH.

In the following blocks I want to use DIA Umpire to create transition libraries for openswath, then I want to run openswath and score the runs.

1.1 Invoke DIA Umpire

cd ~/scratch/proteomics/mycobacterium_tuberculosis_2018
module add openms

type="mzXML"
export VERSION="20190228"
basedir="${HOME}/scratch/proteomics/mycobacterium_tuberculosis_2018"
base_input="${basedir}/results/01${type}/dia/${VERSION}/"
umpire_inputs=$(/usr/bin/find "${base_input}" -name "*.${type}" | sort)
echo "Checking in: ${umpire_inputs}"
for input in ${umpire_inputs};
do
    in_name=$(basename $input ".${type}")
    out_name="${in_name}_Q1.mgf"
    if [[ ! -f "${base_input}/${out_name}" ]]; then
        echo "The output file: ${out_name} already exists."
    else
        java -jar DIA_Umpire_SE.jar ${input} diaumpire_se.params
    fi
done

1.2 Convert DIA Umpire results

msconvert --mzXML results/02mgf/*.mgf
mv *.mzXML results/03_dia_umpire_mzxml

1.3 Search the Umpire results

comet \
    -Pparameters/comet_dia_umpire_params.txt \
    results/03_dia_umpire_mzxml/*.mzXML

1.4 Merge them

xinteract \
    -dDECOY_ \
    -OARPd \
    -Ninteract.comet.pep.xml \
    results/03_dia_umpire_mzxml/*.pep.xml

mv interact.comet.pep.xml results/04_dia_umpire_xinteract

1.5 Combine the statistics

ProteinProphet \
    results/04_dia_umpire_xinteract/interact.comet.pep.xml \
    results/05_dia_umpire_prophet/combined.prot.xml

InterProphetParser \
    DECOY=DECOY \
    results/04_dia_umpire_xinteract/interact.comet.pep.xml \
    results/05_dia_umpire_prophet/iProphet.pep.xml

Mayu.pl \
    -A results/05_dia_umpire_prophet/iProphet.pep.xml \
    -C reference/mtb_irt.fasta \
    -E DECOY

mayu_output <- "../2019-05-12_12.37.03_main_1.07.csv"
number <- hpgltools::extract_mayu_pps_fdr(mayu_output)
message("The number is: ", number)
## 0.43291

## Rerunning because writing the file failed.
spectrast \
    -cNSpecLib -cICID-QTOF \
    -cf "Protein! ~ DECOY_" \
    -cP0.4237 \
    -c_IRTreference/irt.txt \
    -c_IRR results/05_dia_umpire_prophet/iProphet.pep.xml

spectrast \
    -cNSpecLib_cons \
    -cICID-QTOF \
    -cAC SpecLib.splib

spectrast2tsv.py \
    -l 350,2000 \
    -s b,y \
    -x 1,2 \
    -o 6 \
    -n 6 \
    -p 0.05 \
    -d -e \
    -k openswath \
    -w windows/2018_0817BrikenTrypsinDIA19.txt \
    -a SpecLib_cons_openswath.tsv \
    SpecLib_cons.sptxt

TargetedFileConverter \
    -in SpecLib_cons_openswath.tsv \
    -in_type tsv \
    -out SpecLib_cons_openswath.TraML \
    -out_type TraML

OpenSwathDecoyGenerator \
    -in SpecLib_cons_openswath.TraML \
    -out SpecLib_cons_openswath_decoy.TraML \
    -method shuffle
##    -exclude_similar \
##    -similarity_threshold 0.05 \
##    -identity_threshold 0.7

TargetedFileConverter \
    -in SpecLib_cons_openswath_decoy.TraML \
    -in_type TraML \
    -out SpecLib_cons_openswath_decoy.tsv \
    -out_type tsv

TargetedFileConverter \
    -in SpecLib_cons_openswath_decoy.TraML \
    -in_type TraML \
    -out SpecLib_cons_openswath_decoy.pqp \
    -out_type pqp

export VERSION=${VERSION:-20190327}
echo "Loading environment modules and parameters for version: ${VERSION}."
source "parameters/${VERSION}_settings.sh"

echo "Invoking the OpenSwathWorkflow using local comet-derived transitions."
type="diaumpire"
input_type="mzXML"
export TRANSITION_PREFIX="SpecLib_cons_openswath_decoy"
echo "Checking in, the transition library is: ${TRANSITION_PREFIX}.pqp"
base_mzxmldir="results/01${input_type}/dia/${VERSION}"
swath_inputs=$(/usr/bin/find "${base_mzxmldir}" -name *.${input_type} -print | sort)
echo "Checking in, the inputs are: ${swath_inputs}"
mkdir -p "${SWATH_OUTDIR}_${type}"
pypdir="${PYPROPHET_OUTDIR}_${type}"
mkdir -p "${pypdir}"
for input in ${swath_inputs}
do
    name=$(basename "${input}" ".${input_type}")
    echo "Starting openswath run, library type ${type} for ${name} using ${MZ_WINDOWS} windows at $(date)."
    swath_output_prefix="${SWATH_OUTDIR}_${type}/${name}_${DDA_METHOD}"
    pyprophet_output_prefix="${PYPROPHET_OUTDIR}_${type}/${name}_${DDA_METHOD}"
    echo "Deleting previous swath output file: ${swath_output_prefix}.osw"
    rm -f "${swath_output_prefix}.osw"
    rm -f "${swath_output_prefix}.tsv"
    OpenSwathWorkflow \
        -in "${input}" \
        -force \
        -sort_swath_maps \
        -min_upper_edge_dist 1 \
        -mz_correction_function "quadratic_regression_delta_ppm" \
        -Scoring:TransitionGroupPicker:background_subtraction "original" \
        -Scoring:stop_report_after_feature "5" \
        -swath_windows_file "windows/openswath_${name}.txt" \
        -tr "${TRANSITION_PREFIX}.pqp" \
        -out_tsv "${swath_output_prefix}.tsv"
    OpenSwathWorkflow \
        -in "${input}" \
        -force \
        -sort_swath_maps \
        -min_upper_edge_dist 1 \
        -mz_correction_function "quadratic_regression_delta_ppm" \
        -Scoring:TransitionGroupPicker:background_subtraction "original" \
        -Scoring:stop_report_after_feature "5" \
        -swath_windows_file "windows/openswath_${name}.txt" \
        -tr "${TRANSITION_PREFIX}.pqp" \
        -out_osw "${swath_output_prefix}.osw"
    ##2>"${swath_output_prefix}_osw.log" 1>&2
done
swath_out=$(dirname ${swath_output_prefix})
pyprophet_out="$(dirname "${pyprophet_output_prefix}")/openswath_merged.osw"
echo "Merging osw files to ${pyprophet_out}"
pyprophet merge \
          --template "${TRANSITION_PREFIX}.pqp" \
          --out="${pyprophet_out}" \
          ${swath_out}/*.osw
pyprophet score --in="${pyprophet_out}"
pyprophet export --in="${pyprophet_out}" --out "test.tsv"
## pyprophet always exports to the current working directory.
final_name="$(dirname ${pyprophet_out})/$(basename ${pyprophet_out} ".osw").tsv"
echo $final_name
mv "test.tsv"
ls -ld "${pyprophet_out}"

tric_tb="${TRIC_OUTDIR}_tuberculist"
mkdir -p "${tric_tb}"
feature_alignment.py \
    --force \
    --in "./${pypdir}/"*.tsv \
    --out "${tric_tb}/${SEARCH_METHOD}_${DDA_METHOD}.tsv" \
    --out_matrix "${tric_tb}/${DDA_METHOD}_outmatrix.tsv" \
    --out_meta "${tric_tb}/${DDA_METHOD}_meta.tsv"
2>"${tric_tb}/feature_alignment.err" \
 1>"${tric_tb}/feature_alignment.out"
echo "Wrote final output to ${tric_tb}/${SEARCH_METHOD}_${DDA_METHOD}.tsv"

2 DEP usage

Thanks to Vivek, I now am aware of DEP, which does everything I wish MSstats did. The matrix given to me by tric’s feature_alignment.py I think gives me what DEP requires, along with my annotations and sample sheet.

Let us see if this is true.

2.1 Protein annotations

mtb_gff <- "reference/mycobacterium_tuberculosis_h37rv_2.gff.gz"
mtb_genome <- "reference/mtuberculosis_h37rv_genbank.fasta"
mtb_cds <- "reference/mtb_cds.fasta"

mtb_annotations <- sm(load_gff_annotations(mtb_gff, type="gene"))
colnames(mtb_annotations) <- gsub(pattern="\\.", replacement="", x=colnames(mtb_annotations))
mtb_annotations[["description"]] <- gsub(pattern="\\+", replacement=" ",
                                         x=mtb_annotations[["description"]])
mtb_annotations[["function"]] <- gsub(pattern="\\+", replacement=" ",
                                      x=mtb_annotations[["function"]])
rownames(mtb_annotations) <- mtb_annotations[["ID"]]

mtb_microbes <- load_microbesonline_annotations(id=83332)

## The species being downloaded is: Mycobacterium tuberculosis H37Rv

2.2 Preprocess intensities in preparation for DEP

tric_data <- read.csv(
  paste0("results/tric/", ver, "/whole_8mz_dia_umpire/comet_HCD.tsv"), sep="\t")
tric_data[["ProteinName"]] <- gsub(pattern="^(.*)_.*$", replacement="\\1",
                                   x=tric_data[["ProteinName"]])
sample_annot <- extract_metadata(paste0("sample_sheets/Mtb_dia_samples_", ver, ".xlsx"))
kept <- ! grepl(x=rownames(sample_annot), pattern="^s\\.\\.")
sample_annot <- sample_annot[kept, ]
devtools::load_all("~/scratch/git/SWATH2stats_myforked")

## Loading SWATH2stats

s2s_exp <- SWATH2stats::sample_annotation(data=tric_data, check_files=FALSE,
                             sample_annotation=sample_annot,
                             fullpeptidename_column="fullpeptidename")

## Not checking that the files are identical between the annotation and data.

decoy_lists <- assess_decoy_rate(s2s_exp)

## Number of non-decoy peptides: 17081
## Number of decoy peptides: 1801
## Decoy rate: 0.1054

## This seems a bit high to me, yesno?
fdr_overall <- assess_fdr_overall(s2s_exp, output="Rconsole", plot=TRUE)

byrun_fdr <- assess_fdr_byrun(s2s_exp, FFT=0.7, plot=TRUE, output="Rconsole")

## The average FDR by run on assay level is 0.014

## The average FDR by run on peptide level is 0.016

## The average FDR by run on protein level is 0.001

chosen_mscore <- mscore4assayfdr(s2s_exp, FFT=0.7, fdr_target=0.02)

## Target assay FDR: 0.02

## Required overall m-score cutoff: 0.0031623
## achieving assay FDR: 0.0194

prot_score <- mscore4protfdr(s2s_exp, FFT=0.7, fdr_target=0.02)

## Target protein FDR: 0.02

## Required overall m-score cutoff: 0.01
## achieving protein FDR: 0.00115

filtered_ms <- filter_mscore(s2s_exp, chosen_mscore)

## Original dimension: 221952, new dimension: 211415, difference: 10537.

filtered_fq <- filter_mscore_freqobs(s2s_exp, 0.01, 0.8, rm.decoy=FALSE)

## Peptides need to have been quantified in more conditions than: 48 in order to pass this percentage-based threshold.

## Fraction of peptides selected: 0.00058

## Original dimension: 224796, new dimension: 680, difference: 224116.

filtered_ms_fdr <- filter_mscore_fdr(filtered_ms, FFT=0.7,
                                     overall_protein_fdr_target=prot_score,
                                     upper_overall_peptide_fdr_limit=0.05)

## Target protein FDR: 0.01

## Required overall m-score cutoff: 0.01
## achieving protein FDR: 0

## filter_mscore_fdr is filtering the data...

## finding m-score cutoff to achieve desired protein FDR in protein master list..

## finding m-score cutoff to achieve desired global peptide FDR..

## Target peptide FDR: 0.05

## Required overall m-score cutoff: 0.01

## Achieving peptide FDR: 0

## Proteins selected: 
## Total proteins selected: 2412
## Final target proteins: 2412
## Final decoy proteins: 0

## Peptides mapping to these protein entries selected:
## Total mapping peptides: 16868
## Final target peptides: 16868
## Final decoy peptides: 0

## Total peptides selected from:
## Total peptides: 16868
## Final target peptides: 16868
## Final decoy peptides: 0

## Individual run FDR quality of the peptides was not calculated
## as not every run contains a decoy.

## The decoys have been removed from the returned data.

filtered_ms_fdr_pr <- filter_proteotypic_peptides(filtered_ms_fdr)

## Number of proteins detected: 2363
## Protein identifiers: Rv0577, Rv0242c, Rv3012c, Rv2467, Rv3715c, Rv2220

## Number of proteins detected that are supported by a proteotypic peptide: 2337
## Number of proteotypic peptides detected: 16728

filtered_ms_fdr_pr_all <- filter_all_peptides(filtered_ms_fdr_pr)

## Number of proteins detected: 2337

## First 6 protein identifiers: Rv0577, Rv0242c, Rv3012c, Rv2467, Rv3715c, Rv2220

filtered_ms_fdr_pr_all_str <- filter_on_max_peptides(data=filtered_ms_fdr_pr_all, n_peptides=10)

## Before filtering:
##   Number of proteins: 2337
##   Number of peptides: 16728
## 
## Percentage of peptides removed: 25.94%
## 
## After filtering:
##   Number of proteins: 2331
##   Number of peptides: 12388

filtered_all_filters <- filter_on_min_peptides(data=filtered_ms_fdr_pr_all_str, n_peptides=3)

## Before filtering:
##   Number of proteins: 2331
##   Number of peptides: 12388
## 
## Percentage of peptides removed: 0%
## 
## After filtering:
##   Number of proteins: 2284
##   Number of peptides: 12388

matrix_prefix <- file.path("results", "swath2stats", ver)
if (!file.exists(matrix_prefix)) {
  dir.create(matrix_prefix)
}
protein_matrix_all <- write_matrix_proteins(
  s2s_exp, write.csv=TRUE,
  filename=file.path(matrix_prefix, "protein_all.csv"))

## Protein overview matrix results/swath2stats/20190327/protein_all.csv written to working folder.

dim(protein_matrix_all)

## [1] 2434   45

protein_matrix_mscore <- write_matrix_proteins(
  filtered_ms, write.csv=TRUE,
  filename=file.path(matrix_prefix, "protein_matrix_mscore.csv"))

## Protein overview matrix results/swath2stats/20190327/protein_matrix_mscore.csv written to working folder.

dim(protein_matrix_mscore)

## [1] 2412   45

peptide_matrix_mscore <- write_matrix_peptides(
  filtered_ms, write.csv=TRUE,
  filename=file.path(matrix_prefix, "peptide_matrix_mscore.csv"))

## Peptide overview matrix results/swath2stats/20190327/peptide_matrix_mscore.csv written to working folder.

dim(peptide_matrix_mscore)

## [1] 16868    45

protein_matrix_filtered <- write_matrix_proteins(
  filtered_all_filters, write.csv=TRUE,
  filename=file.path(matrix_prefix, "protein_matrix_filtered.csv"))

## Protein overview matrix results/swath2stats/20190327/protein_matrix_filtered.csv written to working folder.

dim(protein_matrix_filtered)

## [1] 2284   45

peptide_matrix_filtered <- write_matrix_peptides(
  filtered_all_filters, write.csv=TRUE,
  filename=file.path(matrix_prefix, "peptide_matrix_filtered.csv"))

## Peptide overview matrix results/swath2stats/20190327/peptide_matrix_filtered.csv written to working folder.

dim(peptide_matrix_filtered)

## [1] 144819     45

3 Pass the data to DEP and see what happens.

library(DEP)

## 
## Attaching package: 'DEP'

## The following objects are masked from 'package:hpgltools':
## 
##     plot_heatmap, plot_pca

intensities <- protein_matrix_filtered
cols <- gsub(x=colnames(intensities), pattern="^.*(2018.*$)", replacement="s\\1")
cols[[1]] <- "Protein"
colnames(intensities) <- cols

sample_annot[["label"]] <- rownames(sample_annot)
unique_intensities <- make_unique(intensities, "Protein", "Protein", delim=";")
intensity_columns <- grep(pattern="2018", x=cols)

sample_annot[["condition"]] <- as.character(sample_annot[["condition"]])
sample_annot[["label"]] <- as.character(sample_annot[["label"]])
sample_annot[["replicate"]] <- make.names(paste0(sample_annot[["condition"]], "_", sample_annot[["batch"]]), unique=TRUE)
design <- sample_annot[, c("label", "condition", "replicate")]

mtb_se <- DEP::make_se(unique_intensities, intensity_columns, design)

DEP::plot_frequency(mtb_se)

mtb_filt <- DEP::filter_missval(mtb_se, thr=1)
DEP::plot_numbers(mtb_filt)

DEP::plot_coverage(mtb_filt)

mtb_norm <- DEP::normalize_vsn(mtb_filt)

## vsn2: 1628 x 43 matrix (1 stratum).

## Please use 'meanSdPlot' to verify the fit.

DEP::plot_normalization(mtb_filt, mtb_norm)

DEP::plot_missval(mtb_filt)

DEP::plot_detect(mtb_filt)

mtb_imp <- DEP::impute(mtb_norm, fun="MinProb", q=0.01)

## Loading required package: imputeLCMD

## Loading required package: tmvtnorm

## Loading required package: mvtnorm

## Loading required package: Matrix

## Loading required package: stats4

## Loading required package: gmm

## Loading required package: sandwich

## Loading required package: norm

## Loading required package: pcaMethods

## 
## Attaching package: 'pcaMethods'

## The following object is masked from 'package:stats':
## 
##     loadings

## Loading required package: impute

## [1] 1.557

DEP::plot_imputation(mtb_norm, mtb_imp)

mtb_diff <- DEP::test_diff(mtb_imp, type="all")

## Tested contrasts: wt_filtrate_vs_wt_whole, wt_filtrate_vs_delta_filtrate, wt_filtrate_vs_comp_filtrate, wt_filtrate_vs_delta_whole, wt_filtrate_vs_comp_whole, wt_whole_vs_delta_filtrate, wt_whole_vs_comp_filtrate, wt_whole_vs_delta_whole, wt_whole_vs_comp_whole, delta_filtrate_vs_comp_filtrate, delta_filtrate_vs_delta_whole, delta_filtrate_vs_comp_whole, comp_filtrate_vs_delta_whole, comp_filtrate_vs_comp_whole, delta_whole_vs_comp_whole

mtb_diff <- DEP::test_diff(mtb_imp, type="manual",
                           test=c("wt_filtrate_vs_wt_whole",
                                  "delta_filtrate_vs_wt_filtrate",
                                  "comp_filtrate_vs_wt_filtrate",
                                  "wt_filtrate_vs_delta_filtrate",
                                  "wt_filtrate_vs_comp_filtrate",
                                  "wt_whole_vs_delta_whole",
                                  "wt_whole_vs_comp_whole"))

## Tested contrasts: wt_filtrate_vs_wt_whole, delta_filtrate_vs_wt_filtrate, comp_filtrate_vs_wt_filtrate, wt_filtrate_vs_delta_filtrate, wt_filtrate_vs_comp_filtrate, wt_whole_vs_delta_whole, wt_whole_vs_comp_whole

mtb_dep <- DEP::add_rejections(mtb_diff, alpha=0.05, lfc=0.6)
## mtb_pca <- DEP::plot_pca(mtb_dep)
## The PCA plotter provided by DEP has some problems.

DEP::plot_cor(mtb_dep)

DEP::plot_heatmap(mtb_dep, type="centered", kmeans=TRUE)

DEP::plot_heatmap(mtb_dep, type="contrast", kmeans=TRUE)

DEP::plot_volcano(mtb_dep, contrast="wt_whole_vs_delta_whole")

DEP::plot_volcano(mtb_dep, contrast="wt_filtrate_vs_delta_filtrate")

DEP::plot_single(mtb_dep, proteins=c("Rv0287"))

DEP::plot_single(mtb_dep, proteins=c("Rv0287"), type="centered")

DEP::plot_cond(mtb_dep)

mtb_result <- DEP::get_results(mtb_dep)

if (!isTRUE(get0("skip_load"))) {
  message(paste0("This is hpgltools commit: ", get_git_commit()))
  this_save <- paste0(gsub(pattern="\\.Rmd", replace="", x=rmd_file), "-v", ver, ".rda.xz")
  message(paste0("Saving to ", this_save))
  tmp <- sm(saveme(filename=this_save))
  pander::pander(sessionInfo())
}

## If you wish to reproduce this exact build of hpgltools, invoke the following:

## > git clone http://github.com/abelew/hpgltools.git

## > git reset 0abc58e173be7300595d30d407b7efd4e4a512d6

## This is hpgltools commit: Thu May 9 14:56:34 2019 -0400: 0abc58e173be7300595d30d407b7efd4e4a512d6

## Saving to dia_umpire_20190308-v20190327.rda.xz

R version 3.5.3 (2019-03-11)

Platform: x86_64-pc-linux-gnu (64-bit)

locale: LC_CTYPE=en_US.UTF-8, LC_NUMERIC=C, LC_TIME=en_US.UTF-8, LC_COLLATE=en_US.UTF-8, LC_MONETARY=en_US.UTF-8, LC_MESSAGES=en_US.UTF-8, LC_PAPER=en_US.UTF-8, LC_NAME=C, LC_ADDRESS=C, LC_TELEPHONE=C, LC_MEASUREMENT=en_US.UTF-8 and LC_IDENTIFICATION=C

attached base packages: stats4, parallel, stats, graphics, grDevices, utils, datasets, methods and base

other attached packages: imputeLCMD(v.2.0), impute(v.1.56.0), pcaMethods(v.1.74.0), norm(v.1.0-9.5), tmvtnorm(v.1.4-10), gmm(v.1.6-2), sandwich(v.2.5-1), Matrix(v.1.2-17), mvtnorm(v.1.0-10), DEP(v.1.4.1), SWATH2stats(v.1.13.5), testthat(v.2.1.1), hpgltools(v.1.0), Biobase(v.2.42.0) and BiocGenerics(v.0.28.0)

loaded via a namespace (and not attached): shinydashboard(v.0.7.1), tidyselect(v.0.2.5), lme4(v.1.1-21), htmlwidgets(v.1.3), RSQLite(v.2.1.1), AnnotationDbi(v.1.44.0), grid(v.3.5.3), BiocParallel(v.1.16.6), devtools(v.2.0.2), munsell(v.0.5.0), codetools(v.0.2-16), preprocessCore(v.1.44.0), DT(v.0.6), withr(v.2.1.2), colorspace(v.1.4-1), GOSemSim(v.2.8.0), knitr(v.1.23), rstudioapi(v.0.10), DOSE(v.3.8.2), mzID(v.1.20.1), labeling(v.0.3), urltools(v.1.7.3), GenomeInfoDbData(v.1.2.0), polyclip(v.1.10-0), bit64(v.0.9-7), farver(v.1.1.0), rprojroot(v.1.3-2), xfun(v.0.7), R6(v.2.4.0), doParallel(v.1.0.14), GenomeInfoDb(v.1.18.2), bitops(v.1.0-6), fgsea(v.1.8.0), gridGraphics(v.0.4-1), DelayedArray(v.0.8.0), assertthat(v.0.2.1), promises(v.1.0.1), scales(v.1.0.0), ggraph(v.1.0.2), enrichplot(v.1.2.0), gtable(v.0.3.0), affy(v.1.60.0), sva(v.3.30.1), processx(v.3.3.1), rlang(v.0.3.4), genefilter(v.1.64.0), mzR(v.2.16.2), GlobalOptions(v.0.1.0), splines(v.3.5.3), rtracklayer(v.1.42.2), lazyeval(v.0.2.2), selectr(v.0.4-1), europepmc(v.0.3), BiocManager(v.1.30.4), yaml(v.2.2.0), reshape2(v.1.4.3), GenomicFeatures(v.1.34.8), backports(v.1.1.4), httpuv(v.1.5.1), qvalue(v.2.14.1), clusterProfiler(v.3.10.1), tools(v.3.5.3), usethis(v.1.5.0), ggplotify(v.0.0.3), ggplot2(v.3.1.1), affyio(v.1.52.0), gplots(v.3.0.1.1), RColorBrewer(v.1.1-2), sessioninfo(v.1.1.1), MSnbase(v.2.8.3), ggridges(v.0.5.1), Rcpp(v.1.0.1), plyr(v.1.8.4), base64enc(v.0.1-3), progress(v.1.2.2), zlibbioc(v.1.28.0), purrr(v.0.3.2), RCurl(v.1.95-4.12), ps(v.1.3.0), prettyunits(v.1.0.2), GetoptLong(v.0.1.7), viridis(v.0.5.1), cowplot(v.0.9.4), zoo(v.1.8-5), S4Vectors(v.0.20.1), cluster(v.2.0.9), SummarizedExperiment(v.1.12.0), ggrepel(v.0.8.1), colorRamps(v.2.3), fs(v.1.3.1), variancePartition(v.1.12.3), magrittr(v.1.5), data.table(v.1.12.2), DO.db(v.2.9), openxlsx(v.4.1.0), circlize(v.0.4.6), triebeard(v.0.3.0), packrat(v.0.5.0), ProtGenerics(v.1.14.0), matrixStats(v.0.54.0), pkgload(v.1.0.2), mime(v.0.6), hms(v.0.4.2), evaluate(v.0.13), xtable(v.1.8-4), pbkrtest(v.0.4-7), XML(v.3.98-1.19), IRanges(v.2.16.0), gridExtra(v.2.3), shape(v.1.4.4), compiler(v.3.5.3), biomaRt(v.2.38.0), tibble(v.2.1.1), KernSmooth(v.2.23-15), ncdf4(v.1.16.1), crayon(v.1.3.4), minqa(v.1.2.4), htmltools(v.0.3.6), later(v.0.8.0), mgcv(v.1.8-28), tidyr(v.0.8.3), DBI(v.1.0.0), tweenr(v.1.0.1), ComplexHeatmap(v.1.20.0), MASS(v.7.3-51.4), boot(v.1.3-22), readr(v.1.3.1), cli(v.1.1.0), vsn(v.3.50.0), gdata(v.2.18.0), igraph(v.1.2.4.1), GenomicRanges(v.1.34.0), pkgconfig(v.2.0.2), rvcheck(v.0.1.3), GenomicAlignments(v.1.18.1), MALDIquant(v.1.19.3), xml2(v.1.2.0), foreach(v.1.4.4), annotate(v.1.60.1), XVector(v.0.22.0), rvest(v.0.3.4), stringr(v.1.4.0), callr(v.3.2.0), digest(v.0.6.19), Biostrings(v.2.50.2), rmarkdown(v.1.12), fastmatch(v.1.1-0), curl(v.3.3), shiny(v.1.3.2), Rsamtools(v.1.34.1), gtools(v.3.8.1), rjson(v.0.2.20), nloptr(v.1.2.1), nlme(v.3.1-140), jsonlite(v.1.6), desc(v.1.2.0), viridisLite(v.0.3.0), limma(v.3.38.3), pillar(v.1.4.0), lattice(v.0.20-38), httr(v.1.4.0), pkgbuild(v.1.0.3), survival(v.2.44-1.1), GO.db(v.3.7.0), glue(v.1.3.1), remotes(v.2.0.4), fdrtool(v.1.2.15), zip(v.2.0.2), UpSetR(v.1.3.3), iterators(v.1.0.10), pander(v.0.6.3), bit(v.1.1-14), ggforce(v.0.2.2), stringi(v.1.4.3), blob(v.1.1.1), caTools(v.1.17.1.2), memoise(v.1.1.0) and dplyr(v.0.8.1)

tmp <- loadme(filename=this_save)

LS0tCnRpdGxlOiAiTS4gdHViZXJjdWxvc2lzIDIwMTkwMzI3OiBESUEtVW1waXJlIGJhc2VkIE9wZW5TV0FUSCB3b3JrZmxvdy4iCmF1dGhvcjogImF0YiBhYmVsZXdAZ21haWwuY29tIgpkYXRlOiAiYHIgU3lzLkRhdGUoKWAiCm91dHB1dDoKICBodG1sX2RvY3VtZW50OgogICAgY29kZV9kb3dubG9hZDogdHJ1ZQogICAgY29kZV9mb2xkaW5nOiBzaG93CiAgICBmaWdfY2FwdGlvbjogdHJ1ZQogICAgZmlnX2hlaWdodDogNwogICAgZmlnX3dpZHRoOiA3CiAgICBoaWdobGlnaHQ6IHRhbmdvCiAgICBrZWVwX21kOiBmYWxzZQogICAgbW9kZTogc2VsZmNvbnRhaW5lZAogICAgbnVtYmVyX3NlY3Rpb25zOiB0cnVlCiAgICBzZWxmX2NvbnRhaW5lZDogdHJ1ZQogICAgdGhlbWU6IHJlYWRhYmxlCiAgICB0b2M6IHRydWUKICAgIHRvY19mbG9hdDoKICAgICAgY29sbGFwc2VkOiBmYWxzZQogICAgICBzbW9vdGhfc2Nyb2xsOiBmYWxzZQogIHJtZGZvcm1hdHM6OnJlYWR0aGVkb3duOgogICAgY29kZV9kb3dubG9hZDogdHJ1ZQogICAgY29kZV9mb2xkaW5nOiBzaG93CiAgICBkZl9wcmludDogcGFnZWQKICAgIGZpZ19jYXB0aW9uOiB0cnVlCiAgICBmaWdfaGVpZ2h0OiA3CiAgICBmaWdfd2lkdGg6IDcKICAgIGhpZ2hsaWdodDogdGFuZ28KICAgIHdpZHRoOiAzMDAKICAgIGtlZXBfbWQ6IGZhbHNlCiAgICBtb2RlOiBzZWxmY29udGFpbmVkCiAgICB0b2NfZmxvYXQ6IHRydWUKICBCaW9jU3R5bGU6Omh0bWxfZG9jdW1lbnQ6CiAgICBjb2RlX2Rvd25sb2FkOiB0cnVlCiAgICBjb2RlX2ZvbGRpbmc6IHNob3cKICAgIGZpZ19jYXB0aW9uOiB0cnVlCiAgICBmaWdfaGVpZ2h0OiA3CiAgICBmaWdfd2lkdGg6IDcKICAgIGhpZ2hsaWdodDogdGFuZ28KICAgIGtlZXBfbWQ6IGZhbHNlCiAgICBtb2RlOiBzZWxmY29udGFpbmVkCiAgICB0b2NfZmxvYXQ6IHRydWUKLS0tCgo8c3R5bGUgdHlwZT0idGV4dC9jc3MiPgpib2R5LCB0ZCB7CiAgZm9udC1zaXplOiAxNnB4Owp9CmNvZGUucnsKICBmb250LXNpemU6IDE2cHg7Cn0KcHJlIHsKIGZvbnQtc2l6ZTogMTZweAp9Cjwvc3R5bGU+CgpgYGB7ciBvcHRpb25zLCBpbmNsdWRlPUZBTFNFfQpsaWJyYXJ5KCJocGdsdG9vbHMiKQp0dCA8LSBkZXZ0b29sczo6bG9hZF9hbGwoIi9kYXRhL2hwZ2x0b29scyIpCmtuaXRyOjpvcHRzX2tuaXQkc2V0KHdpZHRoPTEyMCwKICAgICAgICAgICAgICAgICAgICAgcHJvZ3Jlc3M9VFJVRSwKICAgICAgICAgICAgICAgICAgICAgdmVyYm9zZT1UUlVFLAogICAgICAgICAgICAgICAgICAgICBlY2hvPVRSVUUpCmtuaXRyOjpvcHRzX2NodW5rJHNldChlcnJvcj1UUlVFLAogICAgICAgICAgICAgICAgICAgICAgZHBpPTk2KQpvbGRfb3B0aW9ucyA8LSBvcHRpb25zKGRpZ2l0cz00LAogICAgICAgICAgICAgICAgICAgICAgIHN0cmluZ3NBc0ZhY3RvcnM9RkFMU0UsCiAgICAgICAgICAgICAgICAgICAgICAga25pdHIuZHVwbGljYXRlLmxhYmVsPSJhbGxvdyIpCmdncGxvdDI6OnRoZW1lX3NldChnZ3Bsb3QyOjp0aGVtZV9idyhiYXNlX3NpemU9MTApKQpydW5kYXRlIDwtIGZvcm1hdChTeXMuRGF0ZSgpLCBmb3JtYXQ9IiVZJW0lZCIpCnByZXZpb3VzX2ZpbGUgPC0gIjAyX2VzdGltYXRpb25faW5mZWN0aW9uXzIwMTgwODIyLlJtZCIKdmVyIDwtICIyMDE5MDMyNyIKCiMjdG1wIDwtIHNtKGxvYWRtZShmaWxlbmFtZT1wYXN0ZTAoZ3N1YihwYXR0ZXJuPSJcXC5SbWQiLCByZXBsYWNlPSIiLCB4PXByZXZpb3VzX2ZpbGUpLCAiLXYiLCB2ZXIsICIucmRhLnh6IikpKQpybWRfZmlsZSA8LSAiZGlhX3VtcGlyZV8yMDE5MDMwOC5SbWQiCmBgYAoKIyBSZXdvcmsgRElBLVVtcGlyZSB0byBmZWVkIE9wZW5TV0FUSC4KCkluIHRoZSBmb2xsb3dpbmcgYmxvY2tzIEkgd2FudCB0byB1c2UgRElBIFVtcGlyZSB0byBjcmVhdGUgdHJhbnNpdGlvbiBsaWJyYXJpZXMgZm9yIG9wZW5zd2F0aCwKdGhlbiBJIHdhbnQgdG8gcnVuIG9wZW5zd2F0aCBhbmQgc2NvcmUgdGhlIHJ1bnMuCgojIyBJbnZva2UgRElBIFVtcGlyZQoKYGBge2Jhc2ggdW1waXJlLCBldmFsPUZBTFNFfQpjZCB+L3NjcmF0Y2gvcHJvdGVvbWljcy9teWNvYmFjdGVyaXVtX3R1YmVyY3Vsb3Npc18yMDE4Cm1vZHVsZSBhZGQgb3Blbm1zCgp0eXBlPSJtelhNTCIKZXhwb3J0IFZFUlNJT049IjIwMTkwMjI4IgpiYXNlZGlyPSIke0hPTUV9L3NjcmF0Y2gvcHJvdGVvbWljcy9teWNvYmFjdGVyaXVtX3R1YmVyY3Vsb3Npc18yMDE4IgpiYXNlX2lucHV0PSIke2Jhc2VkaXJ9L3Jlc3VsdHMvMDEke3R5cGV9L2RpYS8ke1ZFUlNJT059LyIKdW1waXJlX2lucHV0cz0kKC91c3IvYmluL2ZpbmQgIiR7YmFzZV9pbnB1dH0iIC1uYW1lICIqLiR7dHlwZX0iIHwgc29ydCkKZWNobyAiQ2hlY2tpbmcgaW46ICR7dW1waXJlX2lucHV0c30iCmZvciBpbnB1dCBpbiAke3VtcGlyZV9pbnB1dHN9OwpkbwogICAgaW5fbmFtZT0kKGJhc2VuYW1lICRpbnB1dCAiLiR7dHlwZX0iKQogICAgb3V0X25hbWU9IiR7aW5fbmFtZX1fUTEubWdmIgogICAgaWYgW1sgISAtZiAiJHtiYXNlX2lucHV0fS8ke291dF9uYW1lfSIgXV07IHRoZW4KICAgICAgICBlY2hvICJUaGUgb3V0cHV0IGZpbGU6ICR7b3V0X25hbWV9IGFscmVhZHkgZXhpc3RzLiIKICAgIGVsc2UKICAgICAgICBqYXZhIC1qYXIgRElBX1VtcGlyZV9TRS5qYXIgJHtpbnB1dH0gZGlhdW1waXJlX3NlLnBhcmFtcwogICAgZmkKZG9uZQpgYGAKCiMjIENvbnZlcnQgRElBIFVtcGlyZSByZXN1bHRzCgpgYGB7YmFzaCBjb252ZXJ0LCBldmFsPUZBTFNFfQptc2NvbnZlcnQgLS1telhNTCByZXN1bHRzLzAybWdmLyoubWdmCm12ICoubXpYTUwgcmVzdWx0cy8wM19kaWFfdW1waXJlX216eG1sCmBgYAoKIyMgU2VhcmNoIHRoZSBVbXBpcmUgcmVzdWx0cwoKYGBge2Jhc2ggc2VhcmNoX3VtcGlyZSwgZXZhbD1GQUxTRX0KY29tZXQgXAogICAgLVBwYXJhbWV0ZXJzL2NvbWV0X2RpYV91bXBpcmVfcGFyYW1zLnR4dCBcCiAgICByZXN1bHRzLzAzX2RpYV91bXBpcmVfbXp4bWwvKi5telhNTApgYGAKCiMjIE1lcmdlIHRoZW0KCmBgYHtiYXNoIG1lcmdlX3hpbnRlcmFjdCwgZXZhbD1GQUxTRX0KeGludGVyYWN0IFwKICAgIC1kREVDT1lfIFwKICAgIC1PQVJQZCBcCiAgICAtTmludGVyYWN0LmNvbWV0LnBlcC54bWwgXAogICAgcmVzdWx0cy8wM19kaWFfdW1waXJlX216eG1sLyoucGVwLnhtbAoKbXYgaW50ZXJhY3QuY29tZXQucGVwLnhtbCByZXN1bHRzLzA0X2RpYV91bXBpcmVfeGludGVyYWN0CmBgYAoKIyMgQ29tYmluZSB0aGUgc3RhdGlzdGljcwoKYGBge2Jhc2ggcHJvdGVpbl9wcm9waGV0LCBldmFsPUZBTFNFfQpQcm90ZWluUHJvcGhldCBcCiAgICByZXN1bHRzLzA0X2RpYV91bXBpcmVfeGludGVyYWN0L2ludGVyYWN0LmNvbWV0LnBlcC54bWwgXAogICAgcmVzdWx0cy8wNV9kaWFfdW1waXJlX3Byb3BoZXQvY29tYmluZWQucHJvdC54bWwKCkludGVyUHJvcGhldFBhcnNlciBcCiAgICBERUNPWT1ERUNPWSBcCiAgICByZXN1bHRzLzA0X2RpYV91bXBpcmVfeGludGVyYWN0L2ludGVyYWN0LmNvbWV0LnBlcC54bWwgXAogICAgcmVzdWx0cy8wNV9kaWFfdW1waXJlX3Byb3BoZXQvaVByb3BoZXQucGVwLnhtbAoKTWF5dS5wbCBcCiAgICAtQSByZXN1bHRzLzA1X2RpYV91bXBpcmVfcHJvcGhldC9pUHJvcGhldC5wZXAueG1sIFwKICAgIC1DIHJlZmVyZW5jZS9tdGJfaXJ0LmZhc3RhIFwKICAgIC1FIERFQ09ZCmBgYAoKYGBge3IgZXh0cmFjdF9wY3RfbWF5dSwgZXZhbD1GQUxTRX0KbWF5dV9vdXRwdXQgPC0gIi4uLzIwMTktMDUtMTJfMTIuMzcuMDNfbWFpbl8xLjA3LmNzdiIKbnVtYmVyIDwtIGhwZ2x0b29sczo6ZXh0cmFjdF9tYXl1X3Bwc19mZHIobWF5dV9vdXRwdXQpCm1lc3NhZ2UoIlRoZSBudW1iZXIgaXM6ICIsIG51bWJlcikKIyMgMC40MzI5MQpgYGAKCgpgYGB7YmFzaCB1bXBpcmVfY29udGQsIGV2YWw9RkFMU0V9CiMjIFJlcnVubmluZyBiZWNhdXNlIHdyaXRpbmcgdGhlIGZpbGUgZmFpbGVkLgpzcGVjdHJhc3QgXAogICAgLWNOU3BlY0xpYiAtY0lDSUQtUVRPRiBcCiAgICAtY2YgIlByb3RlaW4hIH4gREVDT1lfIiBcCiAgICAtY1AwLjQyMzcgXAogICAgLWNfSVJUcmVmZXJlbmNlL2lydC50eHQgXAogICAgLWNfSVJSIHJlc3VsdHMvMDVfZGlhX3VtcGlyZV9wcm9waGV0L2lQcm9waGV0LnBlcC54bWwKCnNwZWN0cmFzdCBcCiAgICAtY05TcGVjTGliX2NvbnMgXAogICAgLWNJQ0lELVFUT0YgXAogICAgLWNBQyBTcGVjTGliLnNwbGliCgpzcGVjdHJhc3QydHN2LnB5IFwKICAgIC1sIDM1MCwyMDAwIFwKICAgIC1zIGIseSBcCiAgICAteCAxLDIgXAogICAgLW8gNiBcCiAgICAtbiA2IFwKICAgIC1wIDAuMDUgXAogICAgLWQgLWUgXAogICAgLWsgb3BlbnN3YXRoIFwKICAgIC13IHdpbmRvd3MvMjAxOF8wODE3QnJpa2VuVHJ5cHNpbkRJQTE5LnR4dCBcCiAgICAtYSBTcGVjTGliX2NvbnNfb3BlbnN3YXRoLnRzdiBcCiAgICBTcGVjTGliX2NvbnMuc3B0eHQKClRhcmdldGVkRmlsZUNvbnZlcnRlciBcCiAgICAtaW4gU3BlY0xpYl9jb25zX29wZW5zd2F0aC50c3YgXAogICAgLWluX3R5cGUgdHN2IFwKICAgIC1vdXQgU3BlY0xpYl9jb25zX29wZW5zd2F0aC5UcmFNTCBcCiAgICAtb3V0X3R5cGUgVHJhTUwKCk9wZW5Td2F0aERlY295R2VuZXJhdG9yIFwKICAgIC1pbiBTcGVjTGliX2NvbnNfb3BlbnN3YXRoLlRyYU1MIFwKICAgIC1vdXQgU3BlY0xpYl9jb25zX29wZW5zd2F0aF9kZWNveS5UcmFNTCBcCiAgICAtbWV0aG9kIHNodWZmbGUKIyMgICAgLWV4Y2x1ZGVfc2ltaWxhciBcCiMjICAgIC1zaW1pbGFyaXR5X3RocmVzaG9sZCAwLjA1IFwKIyMgICAgLWlkZW50aXR5X3RocmVzaG9sZCAwLjcKClRhcmdldGVkRmlsZUNvbnZlcnRlciBcCiAgICAtaW4gU3BlY0xpYl9jb25zX29wZW5zd2F0aF9kZWNveS5UcmFNTCBcCiAgICAtaW5fdHlwZSBUcmFNTCBcCiAgICAtb3V0IFNwZWNMaWJfY29uc19vcGVuc3dhdGhfZGVjb3kudHN2IFwKICAgIC1vdXRfdHlwZSB0c3YKClRhcmdldGVkRmlsZUNvbnZlcnRlciBcCiAgICAtaW4gU3BlY0xpYl9jb25zX29wZW5zd2F0aF9kZWNveS5UcmFNTCBcCiAgICAtaW5fdHlwZSBUcmFNTCBcCiAgICAtb3V0IFNwZWNMaWJfY29uc19vcGVuc3dhdGhfZGVjb3kucHFwIFwKICAgIC1vdXRfdHlwZSBwcXAKCmV4cG9ydCBWRVJTSU9OPSR7VkVSU0lPTjotMjAxOTAzMjd9CmVjaG8gIkxvYWRpbmcgZW52aXJvbm1lbnQgbW9kdWxlcyBhbmQgcGFyYW1ldGVycyBmb3IgdmVyc2lvbjogJHtWRVJTSU9OfS4iCnNvdXJjZSAicGFyYW1ldGVycy8ke1ZFUlNJT059X3NldHRpbmdzLnNoIgoKZWNobyAiSW52b2tpbmcgdGhlIE9wZW5Td2F0aFdvcmtmbG93IHVzaW5nIGxvY2FsIGNvbWV0LWRlcml2ZWQgdHJhbnNpdGlvbnMuIgp0eXBlPSJkaWF1bXBpcmUiCmlucHV0X3R5cGU9Im16WE1MIgpleHBvcnQgVFJBTlNJVElPTl9QUkVGSVg9IlNwZWNMaWJfY29uc19vcGVuc3dhdGhfZGVjb3kiCmVjaG8gIkNoZWNraW5nIGluLCB0aGUgdHJhbnNpdGlvbiBsaWJyYXJ5IGlzOiAke1RSQU5TSVRJT05fUFJFRklYfS5wcXAiCmJhc2VfbXp4bWxkaXI9InJlc3VsdHMvMDEke2lucHV0X3R5cGV9L2RpYS8ke1ZFUlNJT059Igpzd2F0aF9pbnB1dHM9JCgvdXNyL2Jpbi9maW5kICIke2Jhc2VfbXp4bWxkaXJ9IiAtbmFtZSAqLiR7aW5wdXRfdHlwZX0gLXByaW50IHwgc29ydCkKZWNobyAiQ2hlY2tpbmcgaW4sIHRoZSBpbnB1dHMgYXJlOiAke3N3YXRoX2lucHV0c30iCm1rZGlyIC1wICIke1NXQVRIX09VVERJUn1fJHt0eXBlfSIKcHlwZGlyPSIke1BZUFJPUEhFVF9PVVRESVJ9XyR7dHlwZX0iCm1rZGlyIC1wICIke3B5cGRpcn0iCmZvciBpbnB1dCBpbiAke3N3YXRoX2lucHV0c30KZG8KICAgIG5hbWU9JChiYXNlbmFtZSAiJHtpbnB1dH0iICIuJHtpbnB1dF90eXBlfSIpCiAgICBlY2hvICJTdGFydGluZyBvcGVuc3dhdGggcnVuLCBsaWJyYXJ5IHR5cGUgJHt0eXBlfSBmb3IgJHtuYW1lfSB1c2luZyAke01aX1dJTkRPV1N9IHdpbmRvd3MgYXQgJChkYXRlKS4iCiAgICBzd2F0aF9vdXRwdXRfcHJlZml4PSIke1NXQVRIX09VVERJUn1fJHt0eXBlfS8ke25hbWV9XyR7RERBX01FVEhPRH0iCiAgICBweXByb3BoZXRfb3V0cHV0X3ByZWZpeD0iJHtQWVBST1BIRVRfT1VURElSfV8ke3R5cGV9LyR7bmFtZX1fJHtEREFfTUVUSE9EfSIKICAgIGVjaG8gIkRlbGV0aW5nIHByZXZpb3VzIHN3YXRoIG91dHB1dCBmaWxlOiAke3N3YXRoX291dHB1dF9wcmVmaXh9Lm9zdyIKICAgIHJtIC1mICIke3N3YXRoX291dHB1dF9wcmVmaXh9Lm9zdyIKICAgIHJtIC1mICIke3N3YXRoX291dHB1dF9wcmVmaXh9LnRzdiIKICAgIE9wZW5Td2F0aFdvcmtmbG93IFwKICAgICAgICAtaW4gIiR7aW5wdXR9IiBcCiAgICAgICAgLWZvcmNlIFwKICAgICAgICAtc29ydF9zd2F0aF9tYXBzIFwKICAgICAgICAtbWluX3VwcGVyX2VkZ2VfZGlzdCAxIFwKICAgICAgICAtbXpfY29ycmVjdGlvbl9mdW5jdGlvbiAicXVhZHJhdGljX3JlZ3Jlc3Npb25fZGVsdGFfcHBtIiBcCiAgICAgICAgLVNjb3Jpbmc6VHJhbnNpdGlvbkdyb3VwUGlja2VyOmJhY2tncm91bmRfc3VidHJhY3Rpb24gIm9yaWdpbmFsIiBcCiAgICAgICAgLVNjb3Jpbmc6c3RvcF9yZXBvcnRfYWZ0ZXJfZmVhdHVyZSAiNSIgXAogICAgICAgIC1zd2F0aF93aW5kb3dzX2ZpbGUgIndpbmRvd3Mvb3BlbnN3YXRoXyR7bmFtZX0udHh0IiBcCiAgICAgICAgLXRyICIke1RSQU5TSVRJT05fUFJFRklYfS5wcXAiIFwKICAgICAgICAtb3V0X3RzdiAiJHtzd2F0aF9vdXRwdXRfcHJlZml4fS50c3YiCiAgICBPcGVuU3dhdGhXb3JrZmxvdyBcCiAgICAgICAgLWluICIke2lucHV0fSIgXAogICAgICAgIC1mb3JjZSBcCiAgICAgICAgLXNvcnRfc3dhdGhfbWFwcyBcCiAgICAgICAgLW1pbl91cHBlcl9lZGdlX2Rpc3QgMSBcCiAgICAgICAgLW16X2NvcnJlY3Rpb25fZnVuY3Rpb24gInF1YWRyYXRpY19yZWdyZXNzaW9uX2RlbHRhX3BwbSIgXAogICAgICAgIC1TY29yaW5nOlRyYW5zaXRpb25Hcm91cFBpY2tlcjpiYWNrZ3JvdW5kX3N1YnRyYWN0aW9uICJvcmlnaW5hbCIgXAogICAgICAgIC1TY29yaW5nOnN0b3BfcmVwb3J0X2FmdGVyX2ZlYXR1cmUgIjUiIFwKICAgICAgICAtc3dhdGhfd2luZG93c19maWxlICJ3aW5kb3dzL29wZW5zd2F0aF8ke25hbWV9LnR4dCIgXAogICAgICAgIC10ciAiJHtUUkFOU0lUSU9OX1BSRUZJWH0ucHFwIiBcCiAgICAgICAgLW91dF9vc3cgIiR7c3dhdGhfb3V0cHV0X3ByZWZpeH0ub3N3IgogICAgIyMyPiIke3N3YXRoX291dHB1dF9wcmVmaXh9X29zdy5sb2ciIDE+JjIKZG9uZQpzd2F0aF9vdXQ9JChkaXJuYW1lICR7c3dhdGhfb3V0cHV0X3ByZWZpeH0pCnB5cHJvcGhldF9vdXQ9IiQoZGlybmFtZSAiJHtweXByb3BoZXRfb3V0cHV0X3ByZWZpeH0iKS9vcGVuc3dhdGhfbWVyZ2VkLm9zdyIKZWNobyAiTWVyZ2luZyBvc3cgZmlsZXMgdG8gJHtweXByb3BoZXRfb3V0fSIKcHlwcm9waGV0IG1lcmdlIFwKICAgICAgICAgIC0tdGVtcGxhdGUgIiR7VFJBTlNJVElPTl9QUkVGSVh9LnBxcCIgXAogICAgICAgICAgLS1vdXQ9IiR7cHlwcm9waGV0X291dH0iIFwKICAgICAgICAgICR7c3dhdGhfb3V0fS8qLm9zdwpweXByb3BoZXQgc2NvcmUgLS1pbj0iJHtweXByb3BoZXRfb3V0fSIKcHlwcm9waGV0IGV4cG9ydCAtLWluPSIke3B5cHJvcGhldF9vdXR9IiAtLW91dCAidGVzdC50c3YiCiMjIHB5cHJvcGhldCBhbHdheXMgZXhwb3J0cyB0byB0aGUgY3VycmVudCB3b3JraW5nIGRpcmVjdG9yeS4KZmluYWxfbmFtZT0iJChkaXJuYW1lICR7cHlwcm9waGV0X291dH0pLyQoYmFzZW5hbWUgJHtweXByb3BoZXRfb3V0fSAiLm9zdyIpLnRzdiIKZWNobyAkZmluYWxfbmFtZQptdiAidGVzdC50c3YiCmxzIC1sZCAiJHtweXByb3BoZXRfb3V0fSIKCnRyaWNfdGI9IiR7VFJJQ19PVVRESVJ9X3R1YmVyY3VsaXN0Igpta2RpciAtcCAiJHt0cmljX3RifSIKZmVhdHVyZV9hbGlnbm1lbnQucHkgXAogICAgLS1mb3JjZSBcCiAgICAtLWluICIuLyR7cHlwZGlyfS8iKi50c3YgXAogICAgLS1vdXQgIiR7dHJpY190Yn0vJHtTRUFSQ0hfTUVUSE9EfV8ke0REQV9NRVRIT0R9LnRzdiIgXAogICAgLS1vdXRfbWF0cml4ICIke3RyaWNfdGJ9LyR7RERBX01FVEhPRH1fb3V0bWF0cml4LnRzdiIgXAogICAgLS1vdXRfbWV0YSAiJHt0cmljX3RifS8ke0REQV9NRVRIT0R9X21ldGEudHN2IgoyPiIke3RyaWNfdGJ9L2ZlYXR1cmVfYWxpZ25tZW50LmVyciIgXAogMT4iJHt0cmljX3RifS9mZWF0dXJlX2FsaWdubWVudC5vdXQiCmVjaG8gIldyb3RlIGZpbmFsIG91dHB1dCB0byAke3RyaWNfdGJ9LyR7U0VBUkNIX01FVEhPRH1fJHtEREFfTUVUSE9EfS50c3YiCmBgYAoKIyBERVAgdXNhZ2UKClRoYW5rcyB0byBWaXZlaywgSSBub3cgYW0gYXdhcmUgb2YgREVQLCB3aGljaCBkb2VzIGV2ZXJ5dGhpbmcgSSB3aXNoIE1Tc3RhdHMgZGlkLgpUaGUgbWF0cml4IGdpdmVuIHRvIG1lIGJ5IHRyaWMncyBmZWF0dXJlX2FsaWdubWVudC5weSBJIHRoaW5rIGdpdmVzIG1lIHdoYXQgREVQCnJlcXVpcmVzLCBhbG9uZyB3aXRoIG15IGFubm90YXRpb25zIGFuZCBzYW1wbGUgc2hlZXQuCgpMZXQgdXMgc2VlIGlmIHRoaXMgaXMgdHJ1ZS4KCiMjIFByb3RlaW4gYW5ub3RhdGlvbnMKCmBgYHtyIHByb3RlaW5fYW5ub3RhdGlvbnN9Cm10Yl9nZmYgPC0gInJlZmVyZW5jZS9teWNvYmFjdGVyaXVtX3R1YmVyY3Vsb3Npc19oMzdydl8yLmdmZi5neiIKbXRiX2dlbm9tZSA8LSAicmVmZXJlbmNlL210dWJlcmN1bG9zaXNfaDM3cnZfZ2VuYmFuay5mYXN0YSIKbXRiX2NkcyA8LSAicmVmZXJlbmNlL210Yl9jZHMuZmFzdGEiCgptdGJfYW5ub3RhdGlvbnMgPC0gc20obG9hZF9nZmZfYW5ub3RhdGlvbnMobXRiX2dmZiwgdHlwZT0iZ2VuZSIpKQpjb2xuYW1lcyhtdGJfYW5ub3RhdGlvbnMpIDwtIGdzdWIocGF0dGVybj0iXFwuIiwgcmVwbGFjZW1lbnQ9IiIsIHg9Y29sbmFtZXMobXRiX2Fubm90YXRpb25zKSkKbXRiX2Fubm90YXRpb25zW1siZGVzY3JpcHRpb24iXV0gPC0gZ3N1YihwYXR0ZXJuPSJcXCsiLCByZXBsYWNlbWVudD0iICIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgeD1tdGJfYW5ub3RhdGlvbnNbWyJkZXNjcmlwdGlvbiJdXSkKbXRiX2Fubm90YXRpb25zW1siZnVuY3Rpb24iXV0gPC0gZ3N1YihwYXR0ZXJuPSJcXCsiLCByZXBsYWNlbWVudD0iICIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgeD1tdGJfYW5ub3RhdGlvbnNbWyJmdW5jdGlvbiJdXSkKcm93bmFtZXMobXRiX2Fubm90YXRpb25zKSA8LSBtdGJfYW5ub3RhdGlvbnNbWyJJRCJdXQoKbXRiX21pY3JvYmVzIDwtIGxvYWRfbWljcm9iZXNvbmxpbmVfYW5ub3RhdGlvbnMoaWQ9ODMzMzIpCmBgYAoKIyMgUHJlcHJvY2VzcyBpbnRlbnNpdGllcyBpbiBwcmVwYXJhdGlvbiBmb3IgREVQCgpgYGB7ciBzd2F0aDJzdGF0cywgZmlnLnNob3c9ImhpZGUifQp0cmljX2RhdGEgPC0gcmVhZC5jc3YoCiAgcGFzdGUwKCJyZXN1bHRzL3RyaWMvIiwgdmVyLCAiL3dob2xlXzhtel9kaWFfdW1waXJlL2NvbWV0X0hDRC50c3YiKSwgc2VwPSJcdCIpCnRyaWNfZGF0YVtbIlByb3RlaW5OYW1lIl1dIDwtIGdzdWIocGF0dGVybj0iXiguKilfLiokIiwgcmVwbGFjZW1lbnQ9IlxcMSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgeD10cmljX2RhdGFbWyJQcm90ZWluTmFtZSJdXSkKc2FtcGxlX2Fubm90IDwtIGV4dHJhY3RfbWV0YWRhdGEocGFzdGUwKCJzYW1wbGVfc2hlZXRzL010Yl9kaWFfc2FtcGxlc18iLCB2ZXIsICIueGxzeCIpKQprZXB0IDwtICEgZ3JlcGwoeD1yb3duYW1lcyhzYW1wbGVfYW5ub3QpLCBwYXR0ZXJuPSJec1xcLlxcLiIpCnNhbXBsZV9hbm5vdCA8LSBzYW1wbGVfYW5ub3Rba2VwdCwgXQpkZXZ0b29sczo6bG9hZF9hbGwoIn4vc2NyYXRjaC9naXQvU1dBVEgyc3RhdHNfbXlmb3JrZWQiKQpzMnNfZXhwIDwtIFNXQVRIMnN0YXRzOjpzYW1wbGVfYW5ub3RhdGlvbihkYXRhPXRyaWNfZGF0YSwgY2hlY2tfZmlsZXM9RkFMU0UsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgc2FtcGxlX2Fubm90YXRpb249c2FtcGxlX2Fubm90LAogICAgICAgICAgICAgICAgICAgICAgICAgICAgIGZ1bGxwZXB0aWRlbmFtZV9jb2x1bW49ImZ1bGxwZXB0aWRlbmFtZSIpCgpkZWNveV9saXN0cyA8LSBhc3Nlc3NfZGVjb3lfcmF0ZShzMnNfZXhwKQojIyBUaGlzIHNlZW1zIGEgYml0IGhpZ2ggdG8gbWUsIHllc25vPwpmZHJfb3ZlcmFsbCA8LSBhc3Nlc3NfZmRyX292ZXJhbGwoczJzX2V4cCwgb3V0cHV0PSJSY29uc29sZSIsIHBsb3Q9VFJVRSkKCmJ5cnVuX2ZkciA8LSBhc3Nlc3NfZmRyX2J5cnVuKHMyc19leHAsIEZGVD0wLjcsIHBsb3Q9VFJVRSwgb3V0cHV0PSJSY29uc29sZSIpCmNob3Nlbl9tc2NvcmUgPC0gbXNjb3JlNGFzc2F5ZmRyKHMyc19leHAsIEZGVD0wLjcsIGZkcl90YXJnZXQ9MC4wMikKcHJvdF9zY29yZSA8LSBtc2NvcmU0cHJvdGZkcihzMnNfZXhwLCBGRlQ9MC43LCBmZHJfdGFyZ2V0PTAuMDIpCgpmaWx0ZXJlZF9tcyA8LSBmaWx0ZXJfbXNjb3JlKHMyc19leHAsIGNob3Nlbl9tc2NvcmUpCmZpbHRlcmVkX2ZxIDwtIGZpbHRlcl9tc2NvcmVfZnJlcW9icyhzMnNfZXhwLCAwLjAxLCAwLjgsIHJtLmRlY295PUZBTFNFKQpmaWx0ZXJlZF9tc19mZHIgPC0gZmlsdGVyX21zY29yZV9mZHIoZmlsdGVyZWRfbXMsIEZGVD0wLjcsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvdmVyYWxsX3Byb3RlaW5fZmRyX3RhcmdldD1wcm90X3Njb3JlLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgdXBwZXJfb3ZlcmFsbF9wZXB0aWRlX2Zkcl9saW1pdD0wLjA1KQpmaWx0ZXJlZF9tc19mZHJfcHIgPC0gZmlsdGVyX3Byb3Rlb3R5cGljX3BlcHRpZGVzKGZpbHRlcmVkX21zX2ZkcikKZmlsdGVyZWRfbXNfZmRyX3ByX2FsbCA8LSBmaWx0ZXJfYWxsX3BlcHRpZGVzKGZpbHRlcmVkX21zX2Zkcl9wcikKZmlsdGVyZWRfbXNfZmRyX3ByX2FsbF9zdHIgPC0gZmlsdGVyX29uX21heF9wZXB0aWRlcyhkYXRhPWZpbHRlcmVkX21zX2Zkcl9wcl9hbGwsIG5fcGVwdGlkZXM9MTApCmZpbHRlcmVkX2FsbF9maWx0ZXJzIDwtIGZpbHRlcl9vbl9taW5fcGVwdGlkZXMoZGF0YT1maWx0ZXJlZF9tc19mZHJfcHJfYWxsX3N0ciwgbl9wZXB0aWRlcz0zKQoKbWF0cml4X3ByZWZpeCA8LSBmaWxlLnBhdGgoInJlc3VsdHMiLCAic3dhdGgyc3RhdHMiLCB2ZXIpCmlmICghZmlsZS5leGlzdHMobWF0cml4X3ByZWZpeCkpIHsKICBkaXIuY3JlYXRlKG1hdHJpeF9wcmVmaXgpCn0KcHJvdGVpbl9tYXRyaXhfYWxsIDwtIHdyaXRlX21hdHJpeF9wcm90ZWlucygKICBzMnNfZXhwLCB3cml0ZS5jc3Y9VFJVRSwKICBmaWxlbmFtZT1maWxlLnBhdGgobWF0cml4X3ByZWZpeCwgInByb3RlaW5fYWxsLmNzdiIpKQpkaW0ocHJvdGVpbl9tYXRyaXhfYWxsKQpwcm90ZWluX21hdHJpeF9tc2NvcmUgPC0gd3JpdGVfbWF0cml4X3Byb3RlaW5zKAogIGZpbHRlcmVkX21zLCB3cml0ZS5jc3Y9VFJVRSwKICBmaWxlbmFtZT1maWxlLnBhdGgobWF0cml4X3ByZWZpeCwgInByb3RlaW5fbWF0cml4X21zY29yZS5jc3YiKSkKZGltKHByb3RlaW5fbWF0cml4X21zY29yZSkKcGVwdGlkZV9tYXRyaXhfbXNjb3JlIDwtIHdyaXRlX21hdHJpeF9wZXB0aWRlcygKICBmaWx0ZXJlZF9tcywgd3JpdGUuY3N2PVRSVUUsCiAgZmlsZW5hbWU9ZmlsZS5wYXRoKG1hdHJpeF9wcmVmaXgsICJwZXB0aWRlX21hdHJpeF9tc2NvcmUuY3N2IikpCmRpbShwZXB0aWRlX21hdHJpeF9tc2NvcmUpCnByb3RlaW5fbWF0cml4X2ZpbHRlcmVkIDwtIHdyaXRlX21hdHJpeF9wcm90ZWlucygKICBmaWx0ZXJlZF9hbGxfZmlsdGVycywgd3JpdGUuY3N2PVRSVUUsCiAgZmlsZW5hbWU9ZmlsZS5wYXRoKG1hdHJpeF9wcmVmaXgsICJwcm90ZWluX21hdHJpeF9maWx0ZXJlZC5jc3YiKSkKZGltKHByb3RlaW5fbWF0cml4X2ZpbHRlcmVkKQpwZXB0aWRlX21hdHJpeF9maWx0ZXJlZCA8LSB3cml0ZV9tYXRyaXhfcGVwdGlkZXMoCiAgZmlsdGVyZWRfYWxsX2ZpbHRlcnMsIHdyaXRlLmNzdj1UUlVFLAogIGZpbGVuYW1lPWZpbGUucGF0aChtYXRyaXhfcHJlZml4LCAicGVwdGlkZV9tYXRyaXhfZmlsdGVyZWQuY3N2IikpCmRpbShwZXB0aWRlX21hdHJpeF9maWx0ZXJlZCkKYGBgCgojIFBhc3MgdGhlIGRhdGEgdG8gREVQIGFuZCBzZWUgd2hhdCBoYXBwZW5zLgoKYGBge3Igc3RhcnRfREVQfQpsaWJyYXJ5KERFUCkKaW50ZW5zaXRpZXMgPC0gcHJvdGVpbl9tYXRyaXhfZmlsdGVyZWQKY29scyA8LSBnc3ViKHg9Y29sbmFtZXMoaW50ZW5zaXRpZXMpLCBwYXR0ZXJuPSJeLiooMjAxOC4qJCkiLCByZXBsYWNlbWVudD0ic1xcMSIpCmNvbHNbWzFdXSA8LSAiUHJvdGVpbiIKY29sbmFtZXMoaW50ZW5zaXRpZXMpIDwtIGNvbHMKCnNhbXBsZV9hbm5vdFtbImxhYmVsIl1dIDwtIHJvd25hbWVzKHNhbXBsZV9hbm5vdCkKdW5pcXVlX2ludGVuc2l0aWVzIDwtIG1ha2VfdW5pcXVlKGludGVuc2l0aWVzLCAiUHJvdGVpbiIsICJQcm90ZWluIiwgZGVsaW09IjsiKQppbnRlbnNpdHlfY29sdW1ucyA8LSBncmVwKHBhdHRlcm49IjIwMTgiLCB4PWNvbHMpCgpzYW1wbGVfYW5ub3RbWyJjb25kaXRpb24iXV0gPC0gYXMuY2hhcmFjdGVyKHNhbXBsZV9hbm5vdFtbImNvbmRpdGlvbiJdXSkKc2FtcGxlX2Fubm90W1sibGFiZWwiXV0gPC0gYXMuY2hhcmFjdGVyKHNhbXBsZV9hbm5vdFtbImxhYmVsIl1dKQpzYW1wbGVfYW5ub3RbWyJyZXBsaWNhdGUiXV0gPC0gbWFrZS5uYW1lcyhwYXN0ZTAoc2FtcGxlX2Fubm90W1siY29uZGl0aW9uIl1dLCAiXyIsIHNhbXBsZV9hbm5vdFtbImJhdGNoIl1dKSwgdW5pcXVlPVRSVUUpCmRlc2lnbiA8LSBzYW1wbGVfYW5ub3RbLCBjKCJsYWJlbCIsICJjb25kaXRpb24iLCAicmVwbGljYXRlIildCgptdGJfc2UgPC0gREVQOjptYWtlX3NlKHVuaXF1ZV9pbnRlbnNpdGllcywgaW50ZW5zaXR5X2NvbHVtbnMsIGRlc2lnbikKCkRFUDo6cGxvdF9mcmVxdWVuY3kobXRiX3NlKQptdGJfZmlsdCA8LSBERVA6OmZpbHRlcl9taXNzdmFsKG10Yl9zZSwgdGhyPTEpCkRFUDo6cGxvdF9udW1iZXJzKG10Yl9maWx0KQpERVA6OnBsb3RfY292ZXJhZ2UobXRiX2ZpbHQpCgptdGJfbm9ybSA8LSBERVA6Om5vcm1hbGl6ZV92c24obXRiX2ZpbHQpCkRFUDo6cGxvdF9ub3JtYWxpemF0aW9uKG10Yl9maWx0LCBtdGJfbm9ybSkKCkRFUDo6cGxvdF9taXNzdmFsKG10Yl9maWx0KQpERVA6OnBsb3RfZGV0ZWN0KG10Yl9maWx0KQoKbXRiX2ltcCA8LSBERVA6OmltcHV0ZShtdGJfbm9ybSwgZnVuPSJNaW5Qcm9iIiwgcT0wLjAxKQpERVA6OnBsb3RfaW1wdXRhdGlvbihtdGJfbm9ybSwgbXRiX2ltcCkKCm10Yl9kaWZmIDwtIERFUDo6dGVzdF9kaWZmKG10Yl9pbXAsIHR5cGU9ImFsbCIpCm10Yl9kaWZmIDwtIERFUDo6dGVzdF9kaWZmKG10Yl9pbXAsIHR5cGU9Im1hbnVhbCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgIHRlc3Q9Yygid3RfZmlsdHJhdGVfdnNfd3Rfd2hvbGUiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgImRlbHRhX2ZpbHRyYXRlX3ZzX3d0X2ZpbHRyYXRlIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJjb21wX2ZpbHRyYXRlX3ZzX3d0X2ZpbHRyYXRlIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJ3dF9maWx0cmF0ZV92c19kZWx0YV9maWx0cmF0ZSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAid3RfZmlsdHJhdGVfdnNfY29tcF9maWx0cmF0ZSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAid3Rfd2hvbGVfdnNfZGVsdGFfd2hvbGUiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgInd0X3dob2xlX3ZzX2NvbXBfd2hvbGUiKSkKCm10Yl9kZXAgPC0gREVQOjphZGRfcmVqZWN0aW9ucyhtdGJfZGlmZiwgYWxwaGE9MC4wNSwgbGZjPTAuNikKIyMgbXRiX3BjYSA8LSBERVA6OnBsb3RfcGNhKG10Yl9kZXApCiMjIFRoZSBQQ0EgcGxvdHRlciBwcm92aWRlZCBieSBERVAgaGFzIHNvbWUgcHJvYmxlbXMuCgpERVA6OnBsb3RfY29yKG10Yl9kZXApCkRFUDo6cGxvdF9oZWF0bWFwKG10Yl9kZXAsIHR5cGU9ImNlbnRlcmVkIiwga21lYW5zPVRSVUUpCkRFUDo6cGxvdF9oZWF0bWFwKG10Yl9kZXAsIHR5cGU9ImNvbnRyYXN0Iiwga21lYW5zPVRSVUUpCgpERVA6OnBsb3Rfdm9sY2FubyhtdGJfZGVwLCBjb250cmFzdD0id3Rfd2hvbGVfdnNfZGVsdGFfd2hvbGUiKQpERVA6OnBsb3Rfdm9sY2FubyhtdGJfZGVwLCBjb250cmFzdD0id3RfZmlsdHJhdGVfdnNfZGVsdGFfZmlsdHJhdGUiKQoKREVQOjpwbG90X3NpbmdsZShtdGJfZGVwLCBwcm90ZWlucz1jKCJSdjAyODciKSkKREVQOjpwbG90X3NpbmdsZShtdGJfZGVwLCBwcm90ZWlucz1jKCJSdjAyODciKSwgdHlwZT0iY2VudGVyZWQiKQpERVA6OnBsb3RfY29uZChtdGJfZGVwKQoKbXRiX3Jlc3VsdCA8LSBERVA6OmdldF9yZXN1bHRzKG10Yl9kZXApCmBgYAoKCmBgYHtyIHNhdmVtZX0KaWYgKCFpc1RSVUUoZ2V0MCgic2tpcF9sb2FkIikpKSB7CiAgbWVzc2FnZShwYXN0ZTAoIlRoaXMgaXMgaHBnbHRvb2xzIGNvbW1pdDogIiwgZ2V0X2dpdF9jb21taXQoKSkpCiAgdGhpc19zYXZlIDwtIHBhc3RlMChnc3ViKHBhdHRlcm49IlxcLlJtZCIsIHJlcGxhY2U9IiIsIHg9cm1kX2ZpbGUpLCAiLXYiLCB2ZXIsICIucmRhLnh6IikKICBtZXNzYWdlKHBhc3RlMCgiU2F2aW5nIHRvICIsIHRoaXNfc2F2ZSkpCiAgdG1wIDwtIHNtKHNhdmVtZShmaWxlbmFtZT10aGlzX3NhdmUpKQogIHBhbmRlcjo6cGFuZGVyKHNlc3Npb25JbmZvKCkpCn0KYGBgCgpgYGB7ciBsb2FkbWUsIGV2YWw9RkFMU0V9CnRtcCA8LSBsb2FkbWUoZmlsZW5hbWU9dGhpc19zYXZlKQpgYGAK

M. tuberculosis 20190327: DIA-Umpire based OpenSWATH workflow.

atb abelew@gmail.com

2019-05-21