metagene: a package to produce metagene plots
19 May 2021
Contents
Package: metagene
Modified: 18 september, 2015
Compiled: Wed May 19 18:05:04 2021
License: Artistic-2.0 | file LICENSE
1 Introduction
This package produces metagene-like plots to compare the behavior of DNA-interacting proteins at selected groups of features. A typical analysis can be done in viscinity of transcription start sites (TSS) of genes or at any regions of interest (such as enhancers). Multiple combinations of group of features and/or group of bam files can be compared in a single analysis. Bootstraping analysis is used to compare the groups and locate regions with statistically different enrichment profiles. In order to increase the sensitivity of the analysis, alignment data is used instead of peaks produced with peak callers (i.e.: MACS2 or PICS). The metagene package uses bootstrap to obtain a better estimation of the mean enrichment and the confidence interval for every group of samples.
This vignette will introduce all the main features of the metagene package.
2 Loading the metagene package
library(metagene)3 Inputs
3.1 Alignment files (BAM files)
There is no hard limit in the number of BAM files that can be included in an
analysis (but with too many BAM files, memory may become an issue). BAM files
must be indexed. For instance, if you use a file names file.bam, a file
named file.bam.bai or file.baimust be present in the same directory.
The path (relative or absolute) to the BAM files must be in a vector:
bam_files <- get_demo_bam_files()
bam_files## [1] "/tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/align1_rep1.bam"
## [2] "/tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/align1_rep2.bam"
## [3] "/tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/align2_rep1.bam"
## [4] "/tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/align2_rep2.bam"
## [5] "/tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/ctrl.bam"For this demo, we have 2 samples (each with 2 replicates). It is also possible to use a named vector to add your own names to each BAM files:
named_bam_files <- bam_files
names(named_bam_files) <- letters[seq_along(bam_files)]
named_bam_files## a
## "/tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/align1_rep1.bam"
## b
## "/tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/align1_rep2.bam"
## c
## "/tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/align2_rep1.bam"
## d
## "/tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/align2_rep2.bam"
## e
## "/tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/ctrl.bam"Using named BAM files can simplify the use of the metagene helper functions and the creation of the design.
3.2 Genomic regions
3.2.1 BED files
To compare custom regions of interest, it is possible to use a list of one or more BED files.
regions <- get_demo_regions()
regions## [1] "/tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/list1.bed"
## [2] "/tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/list2.bed"The name of the files (without the extension) will be used to name each groups.
metagene also support the
narrowPeak
and the broadPeak.
3.2.2 GRanges or GRangesList objects - Regions
As an alternative to a list of BED files, GRanges or GRangesList objects
can be used.
3.2.3 Available datasets
Some common datasets are already available with the metagene package:
promoters_hg19promoters_hg18promoters_mm10promoters_mm9
data(promoters_hg19)
promoters_hg19## GRanges object with 23056 ranges and 1 metadata column:
## seqnames ranges strand | gene_id
## <Rle> <IRanges> <Rle> | <character>
## 1 chr19 58873215-58875214 - | 1
## 10 chr8 18247755-18249754 + | 10
## 100 chr20 43279377-43281376 - | 100
## 1000 chr18 25756446-25758445 - | 1000
## 10000 chr1 244005887-244007886 - | 10000
## ... ... ... ... . ...
## 9991 chr9 115094945-115096944 - | 9991
## 9992 chr21 35735323-35737322 + | 9992
## 9993 chr22 19108968-19110967 - | 9993
## 9994 chr6 90538619-90540618 + | 9994
## 9997 chr22 50963906-50965905 - | 9997
## -------
## seqinfo: 93 sequences (1 circular) from hg19 genomeFor more details about each datasets, please refer to their documentation
(i.e.:?promoters_hg19).
3.3 Design groups
A design group contains a set of BAM files that, when put together, represent a logical analysis. Furthermore, a design group contains the relationship between every BAM files present. Samples (with or without replicates) and controls can be assigned to a same design group. There can be as many groups as necessary. A BAM file can be assigned to more than one group.
To represent the relationship between every BAM files, design groups must have the following columns:
- The list of paths to every BAM files related to an analysis.
- One column per group of files (replicates and/or controls).
There is two possible way to create design groups, by reading a file or by directly creating a design object in R.
3.3.1 Design groups from a file
Design groups can be loaded into the metagene package by using a text file. As the relationship between BAM files as to be specified, the following columns are mandatory:
- First column: The list of paths (absolute or relative) to every BAM files for all the design groups, the BAM filenames or the BAM names (if a named BAM. file was used).
- Following columns: One column per design group (replicates and/or controls).
The column can take only three values:
- 0: ignore file
- 1: input
- 2: control
The file must also contain a header. It is recommanded to use Samples for the
name of the first column, but the value is not checked. The other columns in
the design file will be used for naming design groups, and must be unique.
fileDesign <- system.file("extdata/design.txt", package="metagene")
design <- read.table(fileDesign, header=TRUE, stringsAsFactors=FALSE)
design$Samples <- paste(system.file("extdata", package="metagene"),
design$Samples, sep="/")
kable(design)| Samples | align1 | align2 |
|---|---|---|
| /tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/align1_rep1.bam | 1 | 0 |
| /tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/align1_rep2.bam | 1 | 0 |
| /tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/align2_rep1.bam | 0 | 1 |
| /tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/align2_rep2.bam | 0 | 1 |
| /tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/ctrl.bam | 2 | 2 |
3.3.2 Design groups from R
It is not obligatory to use a design file, you can create the design
data.frame using your prefered method (as long as the restrictions on the
values mentioned previously are respected).
For instance, the previous design data.frame could have been create directly in R:
design <- data.frame(Samples = c("align1_rep1.bam", "align1_rep2.bam",
"align2_rep1.bam", "align2_rep2.bam", "ctrl.bam"),
align1 = c(1,1,0,0,2), align2 = c(0,0,1,1,2))
design$Samples <- paste0(system.file("extdata", package="metagene"), "/",
design$Samples)
kable(design)| Samples | align1 | align2 |
|---|---|---|
| /tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/align1_rep1.bam | 1 | 0 |
| /tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/align1_rep2.bam | 1 | 0 |
| /tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/align2_rep1.bam | 0 | 1 |
| /tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/align2_rep2.bam | 0 | 1 |
| /tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/ctrl.bam | 2 | 2 |
4 Analysis steps
A typical metagene analysis will consist steps:
- Extraction the read count of every BAM files in selected regions.
- Conversion in coverage.
- Noise removal
- Normalization of the coverage values.
- Table production.
- Data frame production.
- Generation of the metagene plot.
4.1 Minimal analysis
A minimal metagene analysis can be performed in 2 steps:
- Initialization (the
newfunction). plot
regions <- get_demo_regions()
bam_files <- get_demo_bam_files()
# Initialization
mg <- metagene$new(regions = regions, bam_files = bam_files)
# Plotting
mg$plot(title = "Demo metagene plot")## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory
## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory## produce data table : ChIP-Seq## produce data frame : ChIP-Seq## Plot : ChIP-Seq
As you can see, it is not mandatory to explicitly call each step of the
metagene analysis. For instance, in the previous example, the plot function
call the other steps automatically with default values (the next section will
describe the steps in more details).
In this specific case, the plot is messy since by default metagene will produce a curve for each possible combinations of BAM file and regions. Since we have 5 BAM files and 2 regions, this gives us 10 curves.
If we want more control on how every step of the analysis are performed, we have to call each functions directly.
4.2 Complete analysis
In order to fully control every step of a metagene analysis, it is important to understand how a complete analysis is performed. If we are satisfied with the default values, it is not mandatory to explicitly call every step (as was shown in the previous section).
4.2.1 Initialization
During this step, the coverages for every regions specified are extracted from
every BAM files. More specifically, a new GRanges is created by combining
all the regions specified with the regions param of the new function.
regions <- get_demo_regions()
bam_files <- get_demo_bam_files()
mg <- metagene$new(regions = regions, bam_files = bam_files)4.2.2 Producing the table
To produce the table, coverages (produced from Genomics regions (.BED),
Alignment Files (.BAM) and Design Sheet) was treated for noise removal
and normalized. Furthermore, to reduce the computation time during the
following steps, the positions are also binned. Regions, designs, bins,
associated values and orientation of strands are pulled into a data.table
called ‘table’ and accessible thanks to the getter get_table.
We can control the size of the bins with the bin_count argument. By
default, a bin_count of 100 will be used during this step.
mg$produce_table()## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory
## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory## produce data table : ChIP-SeqWe can also use the design we produced earlier to remove background signal and combine replicates:
mg$produce_table(design = design)## produce data table : ChIP-Seq4.2.3 Producing the data.frame
The metagene plot are produced using the ggplot2 package, which require a
data.frame as input. During this step, the values of the ribbon are
calculated. Metagene uses “bootstrap” to obtain a better estimation of the
mean of enrichment for every positions in each groups.
mg$produce_data_frame()## produce data frame : ChIP-Seq4.2.4 Plotting
During this step, metagene will use the data.frame to plot the calculated
values using ggplot2. We show a subset of the regions by using the
region_names and design_names parameter. The region_names correspond to
the names of the regions used during the initialization. The design_name
will vary depending if a design was added. If no design was added, this param
correspond to the BAM name or BAM filenames. Otherwise, we have to use the
names of the columns from the design.
mg$plot(region_names = "list1", title = "Demo plot subset")## Plot : ChIP-Seq
5 Manipulating the metagene objects
5.1 Getters
Multiple getters functions are available to access the data that is stored in a
metagene object.
5.1.1 get_table
To get the data.table containing regions, designs, bins, values at bins and orientation of strands.
mg <- get_demo_metagene()
mg$produce_table()## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory
## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory## produce data table : ChIP-Seqmg$get_table()## region design bin value strand
## 1: list1 align1_rep1 1 0.00 *
## 2: list1 align1_rep1 2 0.75 *
## 3: list1 align1_rep1 3 1.55 *
## 4: list1 align1_rep1 4 1.15 *
## 5: list1 align1_rep1 5 1.25 *
## ---
## 49996: list2 ctrl 96 0.00 *
## 49997: list2 ctrl 97 0.00 *
## 49998: list2 ctrl 98 0.00 *
## 49999: list2 ctrl 99 0.00 *
## 50000: list2 ctrl 100 0.00 *5.1.2 get_matrices
To get the data.table as matrices (the former data structure)
mg <- get_demo_metagene()
mg$produce_table()## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory
## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory## produce data table : ChIP-Seqm <- mg$get_matrices()
# m$list1$ctrl$input to access to region 'list1' and 'ctrl' design5.1.3 get_data_frame
get_data_frame = function(region_names = NULL, design_names = NULL) To get the data.frame containing regions and design
mg <- get_demo_metagene()
mg$produce_table()## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory
## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory## produce data table : ChIP-Seqmg$produce_data_frame()## produce data frame : ChIP-Seqmg$get_data_frame()## region design bin value strand qinf qsup
## 1 list1 align1_rep1 1 0.449 * 0.40626500 0.50566250
## 2 list1 align1_rep1 2 0.552 * 0.50269125 0.61146625
## 3 list1 align1_rep1 3 0.609 * 0.54186625 0.65801125
## 4 list1 align1_rep1 4 0.534 * 0.50469625 0.56690125
## 5 list1 align1_rep1 5 0.522 * 0.46821750 0.56119875
## 6 list1 align1_rep1 6 0.589 * 0.54775000 0.63215000
## 7 list1 align1_rep1 7 0.704 * 0.64611750 0.79647125
## 8 list1 align1_rep1 8 0.591 * 0.53193875 0.63332250
## 9 list1 align1_rep1 9 0.562 * 0.50202375 0.60981375
## 10 list1 align1_rep1 10 0.488 * 0.44708375 0.53143125
## 11 list1 align1_rep1 11 0.614 * 0.56079750 0.67288625
## 12 list1 align1_rep1 12 0.745 * 0.66685875 0.83056500
## 13 list1 align1_rep1 13 0.772 * 0.70260250 0.91077750
## 14 list1 align1_rep1 14 0.659 * 0.57649500 0.74469750
## 15 list1 align1_rep1 15 0.481 * 0.42810625 0.54512375
## 16 list1 align1_rep1 16 0.558 * 0.51856625 0.59794875
## 17 list1 align1_rep1 17 0.614 * 0.57704250 0.65053125
## 18 list1 align1_rep1 18 0.708 * 0.64868375 0.77067625
## 19 list1 align1_rep1 19 0.847 * 0.78266000 0.91433125
## 20 list1 align1_rep1 20 0.745 * 0.69052000 0.79611625
## 21 list1 align1_rep1 21 0.731 * 0.67233250 0.76582750
## 22 list1 align1_rep1 22 0.795 * 0.74179250 0.85845875
## 23 list1 align1_rep1 23 0.772 * 0.71703500 0.85761375
## 24 list1 align1_rep1 24 0.866 * 0.81089875 0.93111625
## 25 list1 align1_rep1 25 0.872 * 0.83640500 0.94473375
## 26 list1 align1_rep1 26 0.841 * 0.77644500 0.89402500
## 27 list1 align1_rep1 27 0.808 * 0.74618875 0.92421000
## 28 list1 align1_rep1 28 0.835 * 0.74232000 0.90883125
## 29 list1 align1_rep1 29 0.915 * 0.80338125 0.99074000
## 30 list1 align1_rep1 30 1.123 * 1.02327750 1.27105500
## 31 list1 align1_rep1 31 1.161 * 0.97320500 1.27363000
## 32 list1 align1_rep1 32 1.122 * 0.94990500 1.28195250
## 33 list1 align1_rep1 33 1.301 * 1.17477500 1.37233625
## 34 list1 align1_rep1 34 1.235 * 1.14676375 1.31868875
## 35 list1 align1_rep1 35 1.285 * 1.22805750 1.34115125
## 36 list1 align1_rep1 36 1.497 * 1.41330625 1.59192000
## 37 list1 align1_rep1 37 1.480 * 1.36268000 1.59290750
## 38 list1 align1_rep1 38 1.433 * 1.33113500 1.61859875
## 39 list1 align1_rep1 39 1.234 * 1.11964125 1.37299250
## 40 list1 align1_rep1 40 1.227 * 1.15369750 1.28723000
## 41 list1 align1_rep1 41 1.302 * 1.22183250 1.38807875
## 42 list1 align1_rep1 42 1.267 * 1.15743750 1.34071750
## 43 list1 align1_rep1 43 1.584 * 1.46939750 1.68998000
## 44 list1 align1_rep1 44 1.602 * 1.52015375 1.67807125
## 45 list1 align1_rep1 45 1.727 * 1.62902375 1.81561625
## 46 list1 align1_rep1 46 2.020 * 1.94740250 2.10723250
## 47 list1 align1_rep1 47 2.120 * 2.00463625 2.20142750
## 48 list1 align1_rep1 48 2.044 * 1.94335375 2.12461000
## 49 list1 align1_rep1 49 2.196 * 2.10252375 2.32661250
## 50 list1 align1_rep1 50 2.340 * 2.24115625 2.44270500
## 51 list1 align1_rep1 51 2.367 * 2.27018125 2.45785500
## 52 list1 align1_rep1 52 2.374 * 2.27672625 2.47803125
## 53 list1 align1_rep1 53 2.086 * 1.95027875 2.19515500
## 54 list1 align1_rep1 54 2.053 * 1.94416125 2.19298375
## 55 list1 align1_rep1 55 2.126 * 2.00572625 2.23785375
## 56 list1 align1_rep1 56 1.710 * 1.61504625 1.80769375
## 57 list1 align1_rep1 57 1.429 * 1.35034375 1.54629500
## 58 list1 align1_rep1 58 1.425 * 1.34431750 1.51256750
## 59 list1 align1_rep1 59 1.618 * 1.52934500 1.72112875
## 60 list1 align1_rep1 60 1.707 * 1.60167625 1.79401750
## 61 list1 align1_rep1 61 1.716 * 1.56517125 1.80916250
## 62 list1 align1_rep1 62 1.699 * 1.54709500 1.80391500
## 63 list1 align1_rep1 63 1.383 * 1.27742625 1.48617375
## 64 list1 align1_rep1 64 1.316 * 1.21961000 1.43805000
## 65 list1 align1_rep1 65 1.293 * 1.21101375 1.36691625
## 66 list1 align1_rep1 66 1.136 * 1.07482375 1.19040125
## 67 list1 align1_rep1 67 1.095 * 1.01282125 1.14271375
## 68 list1 align1_rep1 68 1.028 * 0.96442750 1.08180375
## 69 list1 align1_rep1 69 0.895 * 0.81484625 0.94638250
## 70 list1 align1_rep1 70 0.825 * 0.78043000 0.86521375
## 71 list1 align1_rep1 71 0.898 * 0.83389625 0.94443875
## 72 list1 align1_rep1 72 0.856 * 0.79721625 0.92007125
## 73 list1 align1_rep1 73 0.879 * 0.83165500 0.93093875
## 74 list1 align1_rep1 74 0.899 * 0.83891500 0.95838750
## 75 list1 align1_rep1 75 0.958 * 0.89915000 0.99954250
## 76 list1 align1_rep1 76 1.048 * 0.96514250 1.11053750
## 77 list1 align1_rep1 77 0.963 * 0.90559375 1.03857500
## 78 list1 align1_rep1 78 0.941 * 0.88659125 1.00012375
## 79 list1 align1_rep1 79 0.823 * 0.76883000 0.87416000
## 80 list1 align1_rep1 80 0.564 * 0.52155750 0.60329375
## 81 list1 align1_rep1 81 0.486 * 0.43378000 0.53366250
## 82 list1 align1_rep1 82 0.730 * 0.68888625 0.77919250
## 83 list1 align1_rep1 83 0.820 * 0.76926000 0.85933750
## 84 list1 align1_rep1 84 0.711 * 0.66022875 0.75578375
## 85 list1 align1_rep1 85 0.702 * 0.64510875 0.75465000
## 86 list1 align1_rep1 86 0.664 * 0.59841375 0.72975250
## 87 list1 align1_rep1 87 0.508 * 0.46430000 0.54992375
## 88 list1 align1_rep1 88 0.538 * 0.49657625 0.57120750
## 89 list1 align1_rep1 89 0.582 * 0.54791250 0.62072125
## 90 list1 align1_rep1 90 0.701 * 0.65202000 0.76675250
## 91 list1 align1_rep1 91 0.628 * 0.58238250 0.67972375
## 92 list1 align1_rep1 92 0.649 * 0.61871875 0.69037625
## 93 list1 align1_rep1 93 0.900 * 0.83684625 0.95703625
## 94 list1 align1_rep1 94 0.830 * 0.76429000 0.90837750
## 95 list1 align1_rep1 95 0.590 * 0.54175750 0.65276750
## 96 list1 align1_rep1 96 0.557 * 0.51185250 0.61802125
## 97 list1 align1_rep1 97 0.575 * 0.53569375 0.61488500
## 98 list1 align1_rep1 98 0.525 * 0.46755125 0.56373625
## 99 list1 align1_rep1 99 0.404 * 0.37062750 0.43852000
## 100 list1 align1_rep1 100 0.537 * 0.49404875 0.58099750
## 101 list1 align1_rep2 1 0.092 * 0.06066500 0.12381750
## 102 list1 align1_rep2 2 0.161 * 0.12860125 0.21020875
## 103 list1 align1_rep2 3 0.201 * 0.16586875 0.21888750
## 104 list1 align1_rep2 4 0.147 * 0.13056875 0.17684000
## 105 list1 align1_rep2 5 0.077 * 0.06555875 0.09168500
## 106 list1 align1_rep2 6 0.180 * 0.16535750 0.20291625
## 107 list1 align1_rep2 7 0.267 * 0.23428625 0.29845375
## 108 list1 align1_rep2 8 0.233 * 0.19738500 0.27969375
## 109 list1 align1_rep2 9 0.172 * 0.14705500 0.20016625
## 110 list1 align1_rep2 10 0.253 * 0.22434000 0.28784500
## 111 list1 align1_rep2 11 0.258 * 0.22599625 0.29876750
## 112 list1 align1_rep2 12 0.176 * 0.14371500 0.21012000
## 113 list1 align1_rep2 13 0.193 * 0.16845000 0.22846625
## 114 list1 align1_rep2 14 0.216 * 0.18634625 0.24677625
## 115 list1 align1_rep2 15 0.164 * 0.13840000 0.18906000
## 116 list1 align1_rep2 16 0.123 * 0.10361875 0.14612000
## 117 list1 align1_rep2 17 0.110 * 0.09026500 0.12162625
## 118 list1 align1_rep2 18 0.133 * 0.11339000 0.15024750
## 119 list1 align1_rep2 19 0.169 * 0.15113375 0.18761625
## 120 list1 align1_rep2 20 0.265 * 0.22179250 0.29075375
## 121 list1 align1_rep2 21 0.306 * 0.25964625 0.34779375
## 122 list1 align1_rep2 22 0.399 * 0.36616750 0.44510125
## 123 list1 align1_rep2 23 0.326 * 0.29859000 0.36631375
## 124 list1 align1_rep2 24 0.179 * 0.15772625 0.19920000
## 125 list1 align1_rep2 25 0.251 * 0.22316375 0.27799375
## 126 list1 align1_rep2 26 0.369 * 0.32436750 0.42996250
## 127 list1 align1_rep2 27 0.372 * 0.33316000 0.39962625
## 128 list1 align1_rep2 28 0.320 * 0.27926000 0.35598375
## 129 list1 align1_rep2 29 0.254 * 0.22558375 0.28147375
## 130 list1 align1_rep2 30 0.270 * 0.24206375 0.30244750
## 131 list1 align1_rep2 31 0.382 * 0.33360000 0.45553750
## 132 list1 align1_rep2 32 0.402 * 0.34845625 0.46082125
## 133 list1 align1_rep2 33 0.354 * 0.31916125 0.39139000
## 134 list1 align1_rep2 34 0.499 * 0.45145500 0.54142375
## 135 list1 align1_rep2 35 0.635 * 0.57902000 0.69465375
## 136 list1 align1_rep2 36 0.915 * 0.80854250 0.99433125
## 137 list1 align1_rep2 37 1.010 * 0.92665500 1.10215625
## 138 list1 align1_rep2 38 0.811 * 0.74227875 0.89114125
## 139 list1 align1_rep2 39 0.678 * 0.62120375 0.73277500
## 140 list1 align1_rep2 40 0.727 * 0.65789000 0.82245875
## 141 list1 align1_rep2 41 0.769 * 0.70013000 0.84041750
## 142 list1 align1_rep2 42 0.760 * 0.69615875 0.81119500
## 143 list1 align1_rep2 43 0.647 * 0.58273750 0.68909750
## 144 list1 align1_rep2 44 0.849 * 0.78805625 0.89336250
## 145 list1 align1_rep2 45 1.043 * 0.97912375 1.11073375
## 146 list1 align1_rep2 46 1.037 * 0.95079625 1.12439375
## 147 list1 align1_rep2 47 1.105 * 1.03618625 1.17172000
## 148 list1 align1_rep2 48 1.042 * 0.96295000 1.08441875
## 149 list1 align1_rep2 49 0.911 * 0.84245000 0.99967875
## 150 list1 align1_rep2 50 0.962 * 0.87428125 1.05952625
## 151 list1 align1_rep2 51 1.255 * 1.15751250 1.36179750
## 152 list1 align1_rep2 52 1.254 * 1.15844750 1.36251500
## 153 list1 align1_rep2 53 1.109 * 1.05152625 1.17627875
## 154 list1 align1_rep2 54 1.160 * 1.09980875 1.22621500
## 155 list1 align1_rep2 55 1.112 * 1.06576875 1.17275500
## 156 list1 align1_rep2 56 0.912 * 0.86658500 0.96516375
## 157 list1 align1_rep2 57 0.804 * 0.74820375 0.84356625
## 158 list1 align1_rep2 58 0.721 * 0.67627375 0.75664250
## 159 list1 align1_rep2 59 0.542 * 0.50537750 0.58680125
## 160 list1 align1_rep2 60 0.425 * 0.37689000 0.47602750
## 161 list1 align1_rep2 61 0.479 * 0.42070500 0.52678625
## 162 list1 align1_rep2 62 0.729 * 0.68390875 0.79691250
## 163 list1 align1_rep2 63 0.796 * 0.73034875 0.83698625
## 164 list1 align1_rep2 64 0.647 * 0.58292750 0.70523750
## 165 list1 align1_rep2 65 0.603 * 0.56746500 0.64669250
## 166 list1 align1_rep2 66 0.399 * 0.36832875 0.45282500
## 167 list1 align1_rep2 67 0.430 * 0.40206875 0.47267375
## 168 list1 align1_rep2 68 0.498 * 0.46551125 0.54062000
## 169 list1 align1_rep2 69 0.338 * 0.31734125 0.37603250
## 170 list1 align1_rep2 70 0.440 * 0.40454750 0.48167500
## 171 list1 align1_rep2 71 0.454 * 0.41744000 0.48862625
## 172 list1 align1_rep2 72 0.445 * 0.40629625 0.47024375
## 173 list1 align1_rep2 73 0.386 * 0.35586125 0.41496375
## 174 list1 align1_rep2 74 0.373 * 0.34159875 0.41908250
## 175 list1 align1_rep2 75 0.352 * 0.31343500 0.38641625
## 176 list1 align1_rep2 76 0.182 * 0.16209625 0.20415750
## 177 list1 align1_rep2 77 0.200 * 0.17660000 0.21784875
## 178 list1 align1_rep2 78 0.282 * 0.25702625 0.31989125
## 179 list1 align1_rep2 79 0.196 * 0.17925125 0.21788750
## 180 list1 align1_rep2 80 0.184 * 0.16254375 0.20544125
## 181 list1 align1_rep2 81 0.271 * 0.24877375 0.29537000
## 182 list1 align1_rep2 82 0.261 * 0.22915500 0.28155250
## 183 list1 align1_rep2 83 0.308 * 0.27197250 0.34231375
## 184 list1 align1_rep2 84 0.314 * 0.27862125 0.34594750
## 185 list1 align1_rep2 85 0.225 * 0.19245250 0.26841500
## 186 list1 align1_rep2 86 0.190 * 0.16499000 0.21088250
## 187 list1 align1_rep2 87 0.290 * 0.25949125 0.32778875
## 188 list1 align1_rep2 88 0.375 * 0.34247750 0.42105500
## 189 list1 align1_rep2 89 0.263 * 0.23486375 0.29906500
## 190 list1 align1_rep2 90 0.148 * 0.12806125 0.16756375
## 191 list1 align1_rep2 91 0.258 * 0.22778000 0.29141000
## 192 list1 align1_rep2 92 0.267 * 0.22542250 0.30156875
## 193 list1 align1_rep2 93 0.164 * 0.14053000 0.18277750
## 194 list1 align1_rep2 94 0.162 * 0.13536875 0.17878125
## 195 list1 align1_rep2 95 0.133 * 0.11285750 0.15223250
## 196 list1 align1_rep2 96 0.098 * 0.08138375 0.11741875
## 197 list1 align1_rep2 97 0.117 * 0.08487500 0.14326625
## 198 list1 align1_rep2 98 0.179 * 0.15646250 0.19550750
## 199 list1 align1_rep2 99 0.190 * 0.16591875 0.21058375
## 200 list1 align1_rep2 100 0.185 * 0.16112500 0.21171750
## 201 list1 align2_rep1 1 0.087 * 0.07045625 0.10396750
## 202 list1 align2_rep1 2 0.148 * 0.13012250 0.17394000
## 203 list1 align2_rep1 3 0.175 * 0.14701125 0.20330875
## 204 list1 align2_rep1 4 0.125 * 0.10622750 0.13736500
## 205 list1 align2_rep1 5 0.097 * 0.07974250 0.11416000
## 206 list1 align2_rep1 6 0.060 * 0.04799500 0.06968250
## 207 list1 align2_rep1 7 0.130 * 0.10872625 0.15307875
## 208 list1 align2_rep1 8 0.199 * 0.17607625 0.22675000
## 209 list1 align2_rep1 9 0.166 * 0.14643000 0.18418000
## 210 list1 align2_rep1 10 0.104 * 0.08810625 0.11604250
## 211 list1 align2_rep1 11 0.049 * 0.03926750 0.06443375
## 212 list1 align2_rep1 12 0.031 * 0.02733500 0.03838625
## 213 list1 align2_rep1 13 0.117 * 0.09802500 0.14015125
## 214 list1 align2_rep1 14 0.170 * 0.15518250 0.19083375
## 215 list1 align2_rep1 15 0.163 * 0.14670875 0.18169250
## 216 list1 align2_rep1 16 0.127 * 0.10415250 0.14874250
## 217 list1 align2_rep1 17 0.123 * 0.10064000 0.14282500
## 218 list1 align2_rep1 18 0.183 * 0.16102250 0.20268875
## 219 list1 align2_rep1 19 0.201 * 0.17624875 0.22037000
## 220 list1 align2_rep1 20 0.146 * 0.12726750 0.16540000
## 221 list1 align2_rep1 21 0.153 * 0.13481125 0.17393500
## 222 list1 align2_rep1 22 0.169 * 0.15121750 0.20514375
## 223 list1 align2_rep1 23 0.163 * 0.12821125 0.19085250
## 224 list1 align2_rep1 24 0.151 * 0.12787375 0.17268250
## 225 list1 align2_rep1 25 0.222 * 0.19558000 0.25517375
## 226 list1 align2_rep1 26 0.241 * 0.21971125 0.26067875
## 227 list1 align2_rep1 27 0.230 * 0.20945750 0.25145125
## 228 list1 align2_rep1 28 0.249 * 0.22665500 0.26358875
## 229 list1 align2_rep1 29 0.280 * 0.24457250 0.31573375
## 230 list1 align2_rep1 30 0.255 * 0.23469875 0.27939000
## 231 list1 align2_rep1 31 0.173 * 0.15298500 0.20311375
## 232 list1 align2_rep1 32 0.164 * 0.14268000 0.18799750
## 233 list1 align2_rep1 33 0.275 * 0.24629625 0.30770000
## 234 list1 align2_rep1 34 0.247 * 0.21837500 0.27467625
## 235 list1 align2_rep1 35 0.245 * 0.22449125 0.26742500
## 236 list1 align2_rep1 36 0.299 * 0.25873125 0.32317250
## 237 list1 align2_rep1 37 0.337 * 0.31554750 0.36502000
## 238 list1 align2_rep1 38 0.313 * 0.28352875 0.34597250
## 239 list1 align2_rep1 39 0.257 * 0.22169625 0.28025625
## 240 list1 align2_rep1 40 0.252 * 0.23200500 0.27913750
## 241 list1 align2_rep1 41 0.305 * 0.27454875 0.33345375
## 242 list1 align2_rep1 42 0.371 * 0.33341250 0.39699750
## 243 list1 align2_rep1 43 0.332 * 0.30203000 0.36662625
## 244 list1 align2_rep1 44 0.345 * 0.32086250 0.37590750
## 245 list1 align2_rep1 45 0.318 * 0.29204875 0.34352875
## 246 list1 align2_rep1 46 0.282 * 0.26099875 0.30824750
## 247 list1 align2_rep1 47 0.302 * 0.26943500 0.33703375
## 248 list1 align2_rep1 48 0.374 * 0.34922250 0.39637750
## 249 list1 align2_rep1 49 0.306 * 0.27729625 0.32678125
## 250 list1 align2_rep1 50 0.411 * 0.37497000 0.44416000
## 251 list1 align2_rep1 51 0.649 * 0.59912625 0.69599125
## 252 list1 align2_rep1 52 0.614 * 0.57282000 0.66676875
## 253 list1 align2_rep1 53 0.476 * 0.43517750 0.50669875
## 254 list1 align2_rep1 54 0.430 * 0.39725625 0.46868625
## 255 list1 align2_rep1 55 0.436 * 0.40853375 0.47040625
## 256 list1 align2_rep1 56 0.426 * 0.39286500 0.46471250
## 257 list1 align2_rep1 57 0.396 * 0.36386125 0.45035375
## 258 list1 align2_rep1 58 0.395 * 0.36806500 0.42155375
## 259 list1 align2_rep1 59 0.379 * 0.34643625 0.42558750
## 260 list1 align2_rep1 60 0.304 * 0.27442375 0.33178500
## 261 list1 align2_rep1 61 0.283 * 0.26126125 0.30987375
## 262 list1 align2_rep1 62 0.371 * 0.33776000 0.40363625
## 263 list1 align2_rep1 63 0.286 * 0.26096875 0.31623875
## 264 list1 align2_rep1 64 0.211 * 0.18578250 0.23014375
## 265 list1 align2_rep1 65 0.259 * 0.22963000 0.29162500
## 266 list1 align2_rep1 66 0.258 * 0.23146500 0.28695500
## 267 list1 align2_rep1 67 0.252 * 0.23201125 0.27949250
## 268 list1 align2_rep1 68 0.200 * 0.17785500 0.21792375
## 269 list1 align2_rep1 69 0.140 * 0.12468625 0.15650125
## 270 list1 align2_rep1 70 0.092 * 0.08215875 0.10360250
## 271 list1 align2_rep1 71 0.140 * 0.12077250 0.15423875
## 272 list1 align2_rep1 72 0.238 * 0.20625750 0.28482250
## 273 list1 align2_rep1 73 0.255 * 0.22625625 0.29186750
## 274 list1 align2_rep1 74 0.266 * 0.24052250 0.30746375
## 275 list1 align2_rep1 75 0.256 * 0.22984000 0.28844875
## 276 list1 align2_rep1 76 0.222 * 0.19371250 0.24643250
## 277 list1 align2_rep1 77 0.157 * 0.13541250 0.17866875
## 278 list1 align2_rep1 78 0.199 * 0.17570625 0.22252875
## 279 list1 align2_rep1 79 0.294 * 0.26453000 0.33033750
## 280 list1 align2_rep1 80 0.232 * 0.20588750 0.26138500
## 281 list1 align2_rep1 81 0.177 * 0.15544625 0.20024125
## 282 list1 align2_rep1 82 0.158 * 0.13372750 0.17681125
## 283 list1 align2_rep1 83 0.166 * 0.14947875 0.18125000
## 284 list1 align2_rep1 84 0.168 * 0.14504125 0.18484750
## 285 list1 align2_rep1 85 0.247 * 0.22545625 0.27760500
## 286 list1 align2_rep1 86 0.276 * 0.24284000 0.30062125
## 287 list1 align2_rep1 87 0.162 * 0.14561250 0.17516625
## 288 list1 align2_rep1 88 0.064 * 0.04914750 0.07638250
## 289 list1 align2_rep1 89 0.056 * 0.04593375 0.06606875
## 290 list1 align2_rep1 90 0.129 * 0.11533000 0.14581125
## 291 list1 align2_rep1 91 0.191 * 0.17741750 0.21411000
## 292 list1 align2_rep1 92 0.199 * 0.17367250 0.22230375
## 293 list1 align2_rep1 93 0.167 * 0.14668500 0.18469875
## 294 list1 align2_rep1 94 0.150 * 0.13134750 0.16443125
## 295 list1 align2_rep1 95 0.127 * 0.11015875 0.14275625
## 296 list1 align2_rep1 96 0.143 * 0.12130250 0.16305500
## 297 list1 align2_rep1 97 0.130 * 0.11699000 0.14629375
## 298 list1 align2_rep1 98 0.118 * 0.09952250 0.13681750
## 299 list1 align2_rep1 99 0.116 * 0.08379875 0.13955125
## 300 list1 align2_rep1 100 0.191 * 0.16772750 0.21750250
## 301 list1 align2_rep2 1 0.150 * 0.11433125 0.19528125
## 302 list1 align2_rep2 2 0.132 * 0.09672500 0.17464500
## 303 list1 align2_rep2 3 0.139 * 0.10862875 0.17069375
## 304 list1 align2_rep2 4 0.320 * 0.25945000 0.38009000
## 305 list1 align2_rep2 5 0.677 * 0.59009625 0.80637625
## 306 list1 align2_rep2 6 0.463 * 0.39040250 0.53544875
## 307 list1 align2_rep2 7 0.033 * 0.02697000 0.04136625
## 308 list1 align2_rep2 8 0.080 * 0.06000000 0.10310000
## 309 list1 align2_rep2 9 0.048 * 0.03260625 0.06429000
## 310 list1 align2_rep2 10 0.020 * 0.01222500 0.02677500
## 311 list1 align2_rep2 11 0.015 * 0.01065875 0.01852125
## 312 list1 align2_rep2 12 0.021 * 0.01431375 0.03084250
## 313 list1 align2_rep2 13 0.029 * 0.01995125 0.03937125
## 314 list1 align2_rep2 14 0.153 * 0.12488125 0.18688125
## 315 list1 align2_rep2 15 0.220 * 0.17460000 0.26875000
## 316 list1 align2_rep2 16 0.177 * 0.14507875 0.20933500
## 317 list1 align2_rep2 17 0.342 * 0.27835250 0.39620750
## 318 list1 align2_rep2 18 0.435 * 0.31359500 0.53542000
## 319 list1 align2_rep2 19 0.603 * 0.51795000 0.74629500
## 320 list1 align2_rep2 20 0.400 * 0.28917750 0.49642750
## 321 list1 align2_rep2 21 0.277 * 0.22591625 0.32771250
## 322 list1 align2_rep2 22 0.287 * 0.19003875 0.36615375
## 323 list1 align2_rep2 23 0.277 * 0.22152875 0.33765500
## 324 list1 align2_rep2 24 0.480 * 0.40620500 0.58203500
## 325 list1 align2_rep2 25 0.477 * 0.41500250 0.57655625
## 326 list1 align2_rep2 26 0.364 * 0.30015125 0.42646500
## 327 list1 align2_rep2 27 0.373 * 0.30584375 0.48737000
## 328 list1 align2_rep2 28 0.558 * 0.46808750 0.72727000
## 329 list1 align2_rep2 29 0.397 * 0.32106250 0.46886625
## 330 list1 align2_rep2 30 0.266 * 0.20389625 0.35225375
## 331 list1 align2_rep2 31 0.259 * 0.20261500 0.31036875
## 332 list1 align2_rep2 32 0.555 * 0.44248500 0.69347125
## 333 list1 align2_rep2 33 0.247 * 0.20654750 0.29003250
## 334 list1 align2_rep2 34 0.163 * 0.13613375 0.19466750
## 335 list1 align2_rep2 35 0.275 * 0.22612375 0.33637375
## 336 list1 align2_rep2 36 0.155 * 0.13893500 0.17960250
## 337 list1 align2_rep2 37 0.412 * 0.37301375 0.47652000
## 338 list1 align2_rep2 38 0.509 * 0.42199000 0.57333000
## 339 list1 align2_rep2 39 0.128 * 0.09993250 0.15865000
## 340 list1 align2_rep2 40 0.261 * 0.22739125 0.30338125
## 341 list1 align2_rep2 41 0.901 * 0.79134750 1.01688375
## 342 list1 align2_rep2 42 0.811 * 0.72302375 0.92077875
## 343 list1 align2_rep2 43 0.915 * 0.76224375 1.05698125
## 344 list1 align2_rep2 44 1.259 * 1.12192250 1.39448375
## 345 list1 align2_rep2 45 1.376 * 1.27109375 1.53683875
## 346 list1 align2_rep2 46 1.493 * 1.34748125 1.74295500
## 347 list1 align2_rep2 47 1.185 * 1.07406875 1.31928250
## 348 list1 align2_rep2 48 1.324 * 1.17738000 1.43885125
## 349 list1 align2_rep2 49 2.221 * 1.96893750 2.48386750
## 350 list1 align2_rep2 50 2.301 * 2.06371125 2.51392750
## 351 list1 align2_rep2 51 1.610 * 1.44682500 1.71474125
## 352 list1 align2_rep2 52 1.464 * 1.28046875 1.63153750
## 353 list1 align2_rep2 53 1.342 * 1.18369500 1.57314375
## 354 list1 align2_rep2 54 1.217 * 1.08532375 1.33817875
## 355 list1 align2_rep2 55 1.312 * 1.18478250 1.44607625
## 356 list1 align2_rep2 56 1.316 * 1.20200625 1.47637000
## 357 list1 align2_rep2 57 1.322 * 1.16473375 1.45444250
## 358 list1 align2_rep2 58 0.944 * 0.84119125 1.05375500
## 359 list1 align2_rep2 59 0.709 * 0.64555375 0.78716375
## 360 list1 align2_rep2 60 0.642 * 0.53608500 0.76408250
## 361 list1 align2_rep2 61 0.807 * 0.69071625 0.93295625
## 362 list1 align2_rep2 62 0.846 * 0.73733625 0.98883250
## 363 list1 align2_rep2 63 0.693 * 0.59754750 0.79304625
## 364 list1 align2_rep2 64 1.135 * 0.94503000 1.27938625
## 365 list1 align2_rep2 65 0.830 * 0.76659750 0.95336250
## 366 list1 align2_rep2 66 0.297 * 0.25552125 0.33865250
## 367 list1 align2_rep2 67 0.370 * 0.33190000 0.42343250
## 368 list1 align2_rep2 68 0.341 * 0.27885500 0.38800500
## 369 list1 align2_rep2 69 0.306 * 0.23781000 0.40508250
## 370 list1 align2_rep2 70 0.371 * 0.29895125 0.47000000
## 371 list1 align2_rep2 71 0.412 * 0.36633500 0.47923875
## 372 list1 align2_rep2 72 0.378 * 0.32710875 0.44063375
## 373 list1 align2_rep2 73 0.346 * 0.29063875 0.41915375
## 374 list1 align2_rep2 74 0.292 * 0.22937875 0.35148625
## 375 list1 align2_rep2 75 0.140 * 0.11257000 0.18454750
## 376 list1 align2_rep2 76 0.069 * 0.04140000 0.09582375
## 377 list1 align2_rep2 77 0.105 * 0.06742500 0.13674375
## 378 list1 align2_rep2 78 0.137 * 0.08819125 0.18878375
## 379 list1 align2_rep2 79 0.137 * 0.11442750 0.16554750
## 380 list1 align2_rep2 80 0.472 * 0.39744375 0.55906875
## 381 list1 align2_rep2 81 0.456 * 0.38113375 0.51682250
## 382 list1 align2_rep2 82 0.157 * 0.12145625 0.20216000
## 383 list1 align2_rep2 83 0.046 * 0.02890500 0.05841000
## 384 list1 align2_rep2 84 0.151 * 0.10818250 0.18373000
## 385 list1 align2_rep2 85 0.193 * 0.15918625 0.25068125
## 386 list1 align2_rep2 86 0.220 * 0.18524750 0.26657750
## 387 list1 align2_rep2 87 0.499 * 0.39846125 0.60886500
## 388 list1 align2_rep2 88 0.337 * 0.22834875 0.42002750
## 389 list1 align2_rep2 89 0.140 * 0.10786250 0.17910750
## 390 list1 align2_rep2 90 0.203 * 0.16452500 0.24477750
## 391 list1 align2_rep2 91 0.378 * 0.31453625 0.44423000
## 392 list1 align2_rep2 92 0.556 * 0.44023625 0.63550625
## 393 list1 align2_rep2 93 0.284 * 0.21087250 0.36032250
## 394 list1 align2_rep2 94 0.081 * 0.05261625 0.10875375
## 395 list1 align2_rep2 95 0.138 * 0.10254750 0.17178000
## 396 list1 align2_rep2 96 0.361 * 0.27461500 0.50440500
## 397 list1 align2_rep2 97 0.298 * 0.24387750 0.35793125
## 398 list1 align2_rep2 98 0.606 * 0.51645000 0.75081250
## 399 list1 align2_rep2 99 0.536 * 0.44003000 0.61472750
## 400 list1 align2_rep2 100 0.174 * 0.14021000 0.21396250
## 401 list1 ctrl 1 0.584 * 0.53829375 0.62061500
## 402 list1 ctrl 2 0.515 * 0.47031875 0.55940625
## 403 list1 ctrl 3 0.592 * 0.54872125 0.62753250
## 404 list1 ctrl 4 0.654 * 0.58128000 0.70613125
## 405 list1 ctrl 5 0.651 * 0.58922625 0.70443625
## 406 list1 ctrl 6 0.660 * 0.62938375 0.70468875
## 407 list1 ctrl 7 0.614 * 0.56810500 0.65812000
## 408 list1 ctrl 8 0.574 * 0.54176875 0.60875250
## 409 list1 ctrl 9 0.567 * 0.52167625 0.60489875
## 410 list1 ctrl 10 0.568 * 0.53407625 0.60563250
## 411 list1 ctrl 11 0.644 * 0.60215625 0.68008250
## 412 list1 ctrl 12 0.504 * 0.48022000 0.57097500
## 413 list1 ctrl 13 0.487 * 0.44505250 0.53394750
## 414 list1 ctrl 14 0.530 * 0.48206750 0.55769375
## 415 list1 ctrl 15 0.620 * 0.56022375 0.66398750
## 416 list1 ctrl 16 0.708 * 0.66097625 0.75346875
## 417 list1 ctrl 17 0.798 * 0.73790875 0.84820500
## 418 list1 ctrl 18 0.796 * 0.75146625 0.84118875
## 419 list1 ctrl 19 0.855 * 0.80060875 0.91496500
## 420 list1 ctrl 20 0.790 * 0.72994500 0.83202125
## 421 list1 ctrl 21 0.792 * 0.74160125 0.87766500
## 422 list1 ctrl 22 0.854 * 0.76016875 0.90338750
## 423 list1 ctrl 23 0.786 * 0.73270125 0.84613875
## 424 list1 ctrl 24 0.984 * 0.92391000 1.06104750
## 425 list1 ctrl 25 1.143 * 1.04531375 1.22022625
## 426 list1 ctrl 26 0.919 * 0.83372625 1.03549000
## 427 list1 ctrl 27 0.940 * 0.83896250 1.04449750
## 428 list1 ctrl 28 1.023 * 0.92742625 1.13120625
## 429 list1 ctrl 29 1.112 * 1.06362750 1.17368875
## 430 list1 ctrl 30 1.076 * 1.01844125 1.12696625
## 431 list1 ctrl 31 1.020 * 0.94588625 1.06766125
## 432 list1 ctrl 32 1.217 * 1.15739500 1.29982000
## 433 list1 ctrl 33 1.373 * 1.30355250 1.46312500
## 434 list1 ctrl 34 1.663 * 1.61445375 1.76480250
## 435 list1 ctrl 35 1.750 * 1.65438625 1.83240125
## 436 list1 ctrl 36 1.727 * 1.61514375 1.80554500
## 437 list1 ctrl 37 1.786 * 1.68765125 1.90781875
## 438 list1 ctrl 38 2.013 * 1.90332000 2.17578875
## 439 list1 ctrl 39 2.331 * 2.19805000 2.46585875
## 440 list1 ctrl 40 2.336 * 2.21065625 2.46169875
## 441 list1 ctrl 41 2.092 * 1.99340000 2.19618375
## 442 list1 ctrl 42 2.298 * 2.20229625 2.39054375
## 443 list1 ctrl 43 2.821 * 2.71422875 2.97064125
## 444 list1 ctrl 44 3.097 * 2.94106250 3.27878125
## 445 list1 ctrl 45 3.621 * 3.43600375 3.74016875
## 446 list1 ctrl 46 4.533 * 4.34706375 4.78179000
## 447 list1 ctrl 47 4.631 * 4.44241375 4.80593625
## 448 list1 ctrl 48 4.739 * 4.59530375 4.93530625
## 449 list1 ctrl 49 4.906 * 4.70220500 5.07722875
## 450 list1 ctrl 50 5.254 * 5.03593125 5.42192875
## 451 list1 ctrl 51 5.477 * 5.25016750 5.61831375
## 452 list1 ctrl 52 4.813 * 4.63543875 5.07026125
## 453 list1 ctrl 53 4.509 * 4.34007125 4.73219250
## 454 list1 ctrl 54 3.935 * 3.77377750 4.16446500
## 455 list1 ctrl 55 3.140 * 3.00196875 3.28657375
## 456 list1 ctrl 56 2.912 * 2.79107125 3.03701625
## 457 list1 ctrl 57 2.945 * 2.80649250 3.08159375
## 458 list1 ctrl 58 2.748 * 2.61087500 2.84916000
## 459 list1 ctrl 59 2.555 * 2.41355125 2.66061750
## 460 list1 ctrl 60 2.708 * 2.56737500 2.82127125
## 461 list1 ctrl 61 2.500 * 2.34344250 2.71276000
## 462 list1 ctrl 62 2.279 * 2.16693000 2.41514250
## 463 list1 ctrl 63 2.672 * 2.54231125 2.77967000
## 464 list1 ctrl 64 2.878 * 2.70043250 3.10231125
## 465 list1 ctrl 65 2.269 * 2.15159375 2.42702750
## 466 list1 ctrl 66 1.775 * 1.69339375 1.91705375
## 467 list1 ctrl 67 1.684 * 1.57463750 1.82965125
## 468 list1 ctrl 68 1.558 * 1.46377375 1.63401625
## 469 list1 ctrl 69 1.326 * 1.23619250 1.40341125
## 470 list1 ctrl 70 1.319 * 1.22089750 1.44143125
## 471 list1 ctrl 71 1.376 * 1.25692375 1.52875000
## 472 list1 ctrl 72 1.121 * 1.06089250 1.19779000
## 473 list1 ctrl 73 1.067 * 1.01048750 1.12360750
## 474 list1 ctrl 74 1.111 * 1.06102125 1.19186125
## 475 list1 ctrl 75 1.069 * 1.01594625 1.11104000
## 476 list1 ctrl 76 1.086 * 1.03909250 1.15187625
## 477 list1 ctrl 77 0.938 * 0.88129250 1.00403000
## 478 list1 ctrl 78 0.836 * 0.76304000 0.89531625
## 479 list1 ctrl 79 0.831 * 0.76975500 0.89074250
## 480 list1 ctrl 80 0.868 * 0.77987500 0.93736875
## 481 list1 ctrl 81 0.918 * 0.82936875 0.99007000
## 482 list1 ctrl 82 0.690 * 0.64276250 0.73399000
## 483 list1 ctrl 83 0.764 * 0.71651250 0.81329375
## 484 list1 ctrl 84 0.957 * 0.90002375 1.02990375
## 485 list1 ctrl 85 0.839 * 0.78220125 0.88268375
## 486 list1 ctrl 86 0.984 * 0.91186375 1.03555500
## 487 list1 ctrl 87 1.081 * 0.99129250 1.17212750
## 488 list1 ctrl 88 0.868 * 0.78152750 0.93587125
## 489 list1 ctrl 89 0.779 * 0.73266000 0.82474125
## 490 list1 ctrl 90 1.013 * 0.95754250 1.06675250
## 491 list1 ctrl 91 0.984 * 0.93447500 1.04849500
## 492 list1 ctrl 92 0.995 * 0.94161125 1.06167375
## 493 list1 ctrl 93 1.035 * 0.96106250 1.08389625
## 494 list1 ctrl 94 0.810 * 0.77199125 0.85698625
## 495 list1 ctrl 95 0.902 * 0.83808875 0.98693875
## 496 list1 ctrl 96 1.034 * 0.90993000 1.14143500
## 497 list1 ctrl 97 0.903 * 0.82777125 0.99422125
## 498 list1 ctrl 98 0.808 * 0.72004625 0.87678125
## 499 list1 ctrl 99 0.876 * 0.83715750 0.95494375
## 500 list1 ctrl 100 0.655 * 0.58989750 0.72428125
## 501 list2 align1_rep1 1 0.525 * 0.47166750 0.58705875
## 502 list2 align1_rep1 2 0.495 * 0.43081500 0.56817000
## 503 list2 align1_rep1 3 0.469 * 0.42715625 0.53439125
## 504 list2 align1_rep1 4 0.442 * 0.37196250 0.48740875
## 505 list2 align1_rep1 5 0.482 * 0.38132875 0.56643250
## 506 list2 align1_rep1 6 0.627 * 0.50191500 0.75118125
## 507 list2 align1_rep1 7 0.719 * 0.58677750 0.93772750
## 508 list2 align1_rep1 8 0.704 * 0.53917000 0.88457000
## 509 list2 align1_rep1 9 0.777 * 0.66983000 0.89226750
## 510 list2 align1_rep1 10 0.825 * 0.65219875 1.02668000
## 511 list2 align1_rep1 11 0.668 * 0.56909500 0.78207625
## 512 list2 align1_rep1 12 0.564 * 0.45899250 0.66300125
## 513 list2 align1_rep1 13 0.497 * 0.43784375 0.58682750
## 514 list2 align1_rep1 14 0.580 * 0.50738750 0.65913000
## 515 list2 align1_rep1 15 0.652 * 0.49333750 0.76075000
## 516 list2 align1_rep1 16 0.419 * 0.36255750 0.49452875
## 517 list2 align1_rep1 17 0.443 * 0.39113875 0.51258125
## 518 list2 align1_rep1 18 0.517 * 0.42692250 0.61399125
## 519 list2 align1_rep1 19 0.437 * 0.37201250 0.48415125
## 520 list2 align1_rep1 20 0.500 * 0.43892750 0.53101500
## 521 list2 align1_rep1 21 0.702 * 0.65743125 0.74819250
## 522 list2 align1_rep1 22 0.666 * 0.61844125 0.72955250
## 523 list2 align1_rep1 23 0.572 * 0.51292500 0.62291000
## 524 list2 align1_rep1 24 0.594 * 0.54501000 0.65687125
## 525 list2 align1_rep1 25 0.507 * 0.47194500 0.55213375
## 526 list2 align1_rep1 26 0.497 * 0.44197375 0.52815375
## 527 list2 align1_rep1 27 0.602 * 0.55197500 0.64805000
## 528 list2 align1_rep1 28 0.518 * 0.48604875 0.56176625
## 529 list2 align1_rep1 29 0.353 * 0.30716000 0.39237750
## 530 list2 align1_rep1 30 0.329 * 0.29115625 0.35567750
## 531 list2 align1_rep1 31 0.398 * 0.36628125 0.43587000
## 532 list2 align1_rep1 32 0.426 * 0.39166000 0.45921750
## 533 list2 align1_rep1 33 0.404 * 0.37844500 0.44220500
## 534 list2 align1_rep1 34 0.448 * 0.41244625 0.48698875
## 535 list2 align1_rep1 35 0.622 * 0.55475125 0.69338125
## 536 list2 align1_rep1 36 0.621 * 0.56410750 0.69825750
## 537 list2 align1_rep1 37 0.686 * 0.63125125 0.75304375
## 538 list2 align1_rep1 38 0.706 * 0.63460750 0.76682750
## 539 list2 align1_rep1 39 0.596 * 0.54085125 0.64824375
## 540 list2 align1_rep1 40 0.655 * 0.59007500 0.73609000
## 541 list2 align1_rep1 41 0.687 * 0.62314625 0.74552125
## 542 list2 align1_rep1 42 0.784 * 0.71714250 0.86630875
## 543 list2 align1_rep1 43 0.687 * 0.64538000 0.74190625
## 544 list2 align1_rep1 44 0.623 * 0.55472750 0.67438375
## 545 list2 align1_rep1 45 0.844 * 0.77121750 0.90939250
## 546 list2 align1_rep1 46 0.989 * 0.92240375 1.06856375
## 547 list2 align1_rep1 47 1.054 * 0.96962125 1.15390500
## 548 list2 align1_rep1 48 1.022 * 0.92110250 1.15518625
## 549 list2 align1_rep1 49 1.129 * 1.00834625 1.22264250
## 550 list2 align1_rep1 50 1.097 * 1.00870750 1.20008125
## 551 list2 align1_rep1 51 0.865 * 0.81379875 0.92615750
## 552 list2 align1_rep1 52 0.872 * 0.81252250 0.93062000
## 553 list2 align1_rep1 53 0.926 * 0.87589000 0.98377625
## 554 list2 align1_rep1 54 0.811 * 0.77340625 0.85836625
## 555 list2 align1_rep1 55 0.660 * 0.61915750 0.71331375
## 556 list2 align1_rep1 56 0.695 * 0.65182250 0.75098375
## 557 list2 align1_rep1 57 0.696 * 0.63469000 0.76492375
## 558 list2 align1_rep1 58 0.737 * 0.68021875 0.77833625
## 559 list2 align1_rep1 59 0.593 * 0.55197625 0.64057750
## 560 list2 align1_rep1 60 0.756 * 0.66787875 0.81139250
## 561 list2 align1_rep1 61 0.943 * 0.85072000 1.04444250
## 562 list2 align1_rep1 62 0.850 * 0.77832375 0.94776375
## 563 list2 align1_rep1 63 0.826 * 0.75851125 0.91511750
## 564 list2 align1_rep1 64 0.568 * 0.51124875 0.62760750
## 565 list2 align1_rep1 65 0.590 * 0.53782250 0.63394750
## 566 list2 align1_rep1 66 0.691 * 0.63396625 0.74152125
## 567 list2 align1_rep1 67 0.668 * 0.60675500 0.72693750
## 568 list2 align1_rep1 68 0.566 * 0.52277625 0.62171000
## 569 list2 align1_rep1 69 0.512 * 0.47621375 0.57281375
## 570 list2 align1_rep1 70 0.630 * 0.57289000 0.68094750
## 571 list2 align1_rep1 71 0.720 * 0.66524125 0.77428875
## 572 list2 align1_rep1 72 0.652 * 0.60678875 0.70142750
## 573 list2 align1_rep1 73 0.460 * 0.41177750 0.50180375
## 574 list2 align1_rep1 74 0.388 * 0.34934625 0.41775750
## 575 list2 align1_rep1 75 0.589 * 0.52860875 0.63920250
## 576 list2 align1_rep1 76 0.691 * 0.63554625 0.76461000
## 577 list2 align1_rep1 77 0.645 * 0.59896375 0.69044125
## 578 list2 align1_rep1 78 0.520 * 0.48052750 0.55302375
## 579 list2 align1_rep1 79 0.590 * 0.53718500 0.66891125
## 580 list2 align1_rep1 80 0.609 * 0.56581000 0.65537375
## 581 list2 align1_rep1 81 0.530 * 0.46075000 0.58869625
## 582 list2 align1_rep1 82 0.513 * 0.46423500 0.56294875
## 583 list2 align1_rep1 83 0.491 * 0.42970500 0.54639250
## 584 list2 align1_rep1 84 0.433 * 0.40050000 0.48225375
## 585 list2 align1_rep1 85 0.399 * 0.36564125 0.45764500
## 586 list2 align1_rep1 86 0.394 * 0.34447875 0.43500375
## 587 list2 align1_rep1 87 0.454 * 0.42612375 0.49684875
## 588 list2 align1_rep1 88 0.499 * 0.44855000 0.54237125
## 589 list2 align1_rep1 89 0.541 * 0.45715250 0.59160625
## 590 list2 align1_rep1 90 0.465 * 0.41963500 0.50197625
## 591 list2 align1_rep1 91 0.518 * 0.47375625 0.57753625
## 592 list2 align1_rep1 92 0.433 * 0.39186250 0.47727125
## 593 list2 align1_rep1 93 0.477 * 0.42407375 0.53260125
## 594 list2 align1_rep1 94 0.479 * 0.43481375 0.53200500
## 595 list2 align1_rep1 95 0.598 * 0.51857375 0.67622250
## 596 list2 align1_rep1 96 0.736 * 0.65947500 0.84909375
## 597 list2 align1_rep1 97 0.519 * 0.45595000 0.58547625
## 598 list2 align1_rep1 98 0.637 * 0.54774875 0.74869375
## 599 list2 align1_rep1 99 0.771 * 0.67636750 0.87599875
## 600 list2 align1_rep1 100 0.677 * 0.60452375 0.75997375
## 601 list2 align1_rep2 1 0.062 * 0.04881250 0.07510500
## 602 list2 align1_rep2 2 0.112 * 0.09604500 0.13964250
## 603 list2 align1_rep2 3 0.141 * 0.11874000 0.16255000
## 604 list2 align1_rep2 4 0.112 * 0.09271500 0.13707500
## 605 list2 align1_rep2 5 0.155 * 0.12085500 0.20616375
## 606 list2 align1_rep2 6 0.228 * 0.17337500 0.26383875
## 607 list2 align1_rep2 7 0.219 * 0.17928625 0.26590375
## 608 list2 align1_rep2 8 0.207 * 0.14845625 0.24807875
## 609 list2 align1_rep2 9 0.264 * 0.18074750 0.32731625
## 610 list2 align1_rep2 10 0.245 * 0.16887875 0.31996750
## 611 list2 align1_rep2 11 0.309 * 0.25292000 0.36141625
## 612 list2 align1_rep2 12 0.279 * 0.23233625 0.31295375
## 613 list2 align1_rep2 13 0.206 * 0.17092375 0.23372000
## 614 list2 align1_rep2 14 0.243 * 0.20393125 0.27633375
## 615 list2 align1_rep2 15 0.231 * 0.19991250 0.26114750
## 616 list2 align1_rep2 16 0.106 * 0.08830000 0.12034500
## 617 list2 align1_rep2 17 0.087 * 0.07745000 0.10087750
## 618 list2 align1_rep2 18 0.218 * 0.18181000 0.24442500
## 619 list2 align1_rep2 19 0.296 * 0.26365375 0.32770250
## 620 list2 align1_rep2 20 0.250 * 0.22132625 0.27572250
## 621 list2 align1_rep2 21 0.196 * 0.16960125 0.23271250
## 622 list2 align1_rep2 22 0.180 * 0.15228375 0.22253375
## 623 list2 align1_rep2 23 0.182 * 0.16087000 0.21320375
## 624 list2 align1_rep2 24 0.202 * 0.17569500 0.22939750
## 625 list2 align1_rep2 25 0.164 * 0.13358375 0.19047375
## 626 list2 align1_rep2 26 0.178 * 0.14680625 0.21184500
## 627 list2 align1_rep2 27 0.176 * 0.13289375 0.20720125
## 628 list2 align1_rep2 28 0.150 * 0.12448000 0.17257750
## 629 list2 align1_rep2 29 0.146 * 0.13300125 0.16442000
## 630 list2 align1_rep2 30 0.195 * 0.16679500 0.22285875
## 631 list2 align1_rep2 31 0.171 * 0.14388750 0.20166250
## 632 list2 align1_rep2 32 0.155 * 0.13930625 0.17706875
## 633 list2 align1_rep2 33 0.212 * 0.18836750 0.24044125
## 634 list2 align1_rep2 34 0.218 * 0.19255750 0.24067375
## 635 list2 align1_rep2 35 0.198 * 0.17101375 0.23139375
## 636 list2 align1_rep2 36 0.211 * 0.18420625 0.23721875
## 637 list2 align1_rep2 37 0.282 * 0.25246500 0.29852750
## 638 list2 align1_rep2 38 0.289 * 0.25287750 0.32035250
## 639 list2 align1_rep2 39 0.306 * 0.27578750 0.33394250
## 640 list2 align1_rep2 40 0.261 * 0.22585250 0.31412625
## 641 list2 align1_rep2 41 0.223 * 0.19125000 0.24746750
## 642 list2 align1_rep2 42 0.426 * 0.37494500 0.45346375
## 643 list2 align1_rep2 43 0.374 * 0.32688250 0.40193750
## 644 list2 align1_rep2 44 0.340 * 0.30444875 0.37312500
## 645 list2 align1_rep2 45 0.486 * 0.43706125 0.55433250
## 646 list2 align1_rep2 46 0.428 * 0.38942125 0.47631000
## 647 list2 align1_rep2 47 0.386 * 0.31891125 0.44401125
## 648 list2 align1_rep2 48 0.398 * 0.34455125 0.43758625
## 649 list2 align1_rep2 49 0.386 * 0.34213250 0.42725625
## 650 list2 align1_rep2 50 0.447 * 0.39926625 0.47936625
## 651 list2 align1_rep2 51 0.417 * 0.37382000 0.45629500
## 652 list2 align1_rep2 52 0.378 * 0.34695125 0.42203625
## 653 list2 align1_rep2 53 0.299 * 0.26300625 0.32943000
## 654 list2 align1_rep2 54 0.310 * 0.27058750 0.36252250
## 655 list2 align1_rep2 55 0.356 * 0.32696125 0.39278500
## 656 list2 align1_rep2 56 0.444 * 0.39494000 0.48189375
## 657 list2 align1_rep2 57 0.367 * 0.33493875 0.39422375
## 658 list2 align1_rep2 58 0.285 * 0.25768250 0.31276875
## 659 list2 align1_rep2 59 0.271 * 0.24921875 0.29702750
## 660 list2 align1_rep2 60 0.234 * 0.19890875 0.26405250
## 661 list2 align1_rep2 61 0.317 * 0.26325875 0.35668375
## 662 list2 align1_rep2 62 0.388 * 0.35069625 0.43226875
## 663 list2 align1_rep2 63 0.314 * 0.28415125 0.35273750
## 664 list2 align1_rep2 64 0.289 * 0.24578750 0.32158375
## 665 list2 align1_rep2 65 0.316 * 0.27981000 0.34357875
## 666 list2 align1_rep2 66 0.388 * 0.33941250 0.44430500
## 667 list2 align1_rep2 67 0.361 * 0.32625375 0.38851625
## 668 list2 align1_rep2 68 0.252 * 0.21847625 0.27332750
## 669 list2 align1_rep2 69 0.217 * 0.18840000 0.24028125
## 670 list2 align1_rep2 70 0.246 * 0.22485250 0.27590375
## 671 list2 align1_rep2 71 0.227 * 0.19898500 0.25409500
## 672 list2 align1_rep2 72 0.179 * 0.15645500 0.20069250
## 673 list2 align1_rep2 73 0.162 * 0.13139500 0.18347250
## 674 list2 align1_rep2 74 0.189 * 0.16124875 0.21920000
## 675 list2 align1_rep2 75 0.216 * 0.19113500 0.24718500
## 676 list2 align1_rep2 76 0.279 * 0.23088625 0.31336375
## 677 list2 align1_rep2 77 0.346 * 0.29115875 0.40223125
## 678 list2 align1_rep2 78 0.188 * 0.15814500 0.21514000
## 679 list2 align1_rep2 79 0.148 * 0.12738125 0.17235375
## 680 list2 align1_rep2 80 0.350 * 0.31991875 0.39661375
## 681 list2 align1_rep2 81 0.290 * 0.27086750 0.31863875
## 682 list2 align1_rep2 82 0.131 * 0.11784375 0.14242500
## 683 list2 align1_rep2 83 0.080 * 0.07223375 0.09918875
## 684 list2 align1_rep2 84 0.094 * 0.08040625 0.11066250
## 685 list2 align1_rep2 85 0.133 * 0.10517250 0.15773125
## 686 list2 align1_rep2 86 0.128 * 0.10277125 0.14866000
## 687 list2 align1_rep2 87 0.166 * 0.13254500 0.20321500
## 688 list2 align1_rep2 88 0.199 * 0.17181500 0.23047500
## 689 list2 align1_rep2 89 0.148 * 0.11870375 0.18399250
## 690 list2 align1_rep2 90 0.108 * 0.09558250 0.12401875
## 691 list2 align1_rep2 91 0.168 * 0.14216625 0.18615500
## 692 list2 align1_rep2 92 0.220 * 0.18895875 0.25327625
## 693 list2 align1_rep2 93 0.260 * 0.23071500 0.27858125
## 694 list2 align1_rep2 94 0.311 * 0.26427250 0.35888250
## 695 list2 align1_rep2 95 0.221 * 0.19160125 0.26007750
## 696 list2 align1_rep2 96 0.179 * 0.16171000 0.19761500
## 697 list2 align1_rep2 97 0.296 * 0.25286000 0.33862750
## 698 list2 align1_rep2 98 0.323 * 0.25714500 0.37834000
## 699 list2 align1_rep2 99 0.251 * 0.20565625 0.30164500
## 700 list2 align1_rep2 100 0.332 * 0.29375500 0.37544000
## 701 list2 align2_rep1 1 0.055 * 0.04407375 0.06749250
## 702 list2 align2_rep1 2 0.052 * 0.03961000 0.06640875
## 703 list2 align2_rep1 3 0.068 * 0.05269500 0.08255000
## 704 list2 align2_rep1 4 0.063 * 0.05066500 0.07657000
## 705 list2 align2_rep1 5 0.082 * 0.06832375 0.09622500
## 706 list2 align2_rep1 6 0.111 * 0.09518750 0.12605500
## 707 list2 align2_rep1 7 0.152 * 0.13574750 0.17116750
## 708 list2 align2_rep1 8 0.153 * 0.12960125 0.17930625
## 709 list2 align2_rep1 9 0.156 * 0.13069375 0.17687625
## 710 list2 align2_rep1 10 0.105 * 0.09105000 0.12295500
## 711 list2 align2_rep1 11 0.063 * 0.05430000 0.07751000
## 712 list2 align2_rep1 12 0.123 * 0.10807500 0.13745250
## 713 list2 align2_rep1 13 0.144 * 0.12382750 0.16199750
## 714 list2 align2_rep1 14 0.140 * 0.12261500 0.15403750
## 715 list2 align2_rep1 15 0.167 * 0.14832875 0.18993625
## 716 list2 align2_rep1 16 0.159 * 0.14620625 0.17720500
## 717 list2 align2_rep1 17 0.150 * 0.13297375 0.17267625
## 718 list2 align2_rep1 18 0.137 * 0.11418375 0.15443875
## 719 list2 align2_rep1 19 0.080 * 0.06749000 0.09238250
## 720 list2 align2_rep1 20 0.072 * 0.05853500 0.08432750
## 721 list2 align2_rep1 21 0.131 * 0.09576250 0.15375125
## 722 list2 align2_rep1 22 0.160 * 0.14002875 0.18546125
## 723 list2 align2_rep1 23 0.084 * 0.07214000 0.09973625
## 724 list2 align2_rep1 24 0.063 * 0.05354625 0.07539875
## 725 list2 align2_rep1 25 0.093 * 0.07738875 0.10749250
## 726 list2 align2_rep1 26 0.071 * 0.06174000 0.08637750
## 727 list2 align2_rep1 27 0.075 * 0.06053875 0.09590625
## 728 list2 align2_rep1 28 0.121 * 0.10749500 0.13547125
## 729 list2 align2_rep1 29 0.106 * 0.09236375 0.12461375
## 730 list2 align2_rep1 30 0.112 * 0.09373375 0.12932375
## 731 list2 align2_rep1 31 0.127 * 0.10772875 0.14828250
## 732 list2 align2_rep1 32 0.054 * 0.03850250 0.06588875
## 733 list2 align2_rep1 33 0.020 * 0.01326000 0.02591500
## 734 list2 align2_rep1 34 0.057 * 0.04578500 0.06994375
## 735 list2 align2_rep1 35 0.114 * 0.09658250 0.13478375
## 736 list2 align2_rep1 36 0.168 * 0.14986000 0.19592125
## 737 list2 align2_rep1 37 0.175 * 0.15626125 0.19553625
## 738 list2 align2_rep1 38 0.105 * 0.09149625 0.11921375
## 739 list2 align2_rep1 39 0.123 * 0.10663750 0.13670375
## 740 list2 align2_rep1 40 0.130 * 0.11102125 0.14317125
## 741 list2 align2_rep1 41 0.099 * 0.08445375 0.12231625
## 742 list2 align2_rep1 42 0.125 * 0.10976500 0.15261125
## 743 list2 align2_rep1 43 0.106 * 0.08630500 0.12275375
## 744 list2 align2_rep1 44 0.173 * 0.15108625 0.19548250
## 745 list2 align2_rep1 45 0.226 * 0.19961250 0.25190000
## 746 list2 align2_rep1 46 0.169 * 0.15373375 0.19220375
## 747 list2 align2_rep1 47 0.251 * 0.22887000 0.27117625
## 748 list2 align2_rep1 48 0.241 * 0.21144250 0.26948250
## 749 list2 align2_rep1 49 0.240 * 0.21985125 0.27398000
## 750 list2 align2_rep1 50 0.183 * 0.15911750 0.20648500
## 751 list2 align2_rep1 51 0.181 * 0.16662625 0.20506375
## 752 list2 align2_rep1 52 0.240 * 0.20700125 0.26312000
## 753 list2 align2_rep1 53 0.149 * 0.12790125 0.17046625
## 754 list2 align2_rep1 54 0.148 * 0.12163875 0.16376000
## 755 list2 align2_rep1 55 0.223 * 0.19865125 0.24670000
## 756 list2 align2_rep1 56 0.170 * 0.14739500 0.18083000
## 757 list2 align2_rep1 57 0.127 * 0.10712875 0.15813500
## 758 list2 align2_rep1 58 0.107 * 0.08604125 0.13231875
## 759 list2 align2_rep1 59 0.113 * 0.08737375 0.13844250
## 760 list2 align2_rep1 60 0.133 * 0.10994125 0.15695125
## 761 list2 align2_rep1 61 0.102 * 0.08475375 0.11498250
## 762 list2 align2_rep1 62 0.186 * 0.16908375 0.20595000
## 763 list2 align2_rep1 63 0.205 * 0.18399125 0.22477375
## 764 list2 align2_rep1 64 0.162 * 0.13615250 0.18884375
## 765 list2 align2_rep1 65 0.134 * 0.11794500 0.15667750
## 766 list2 align2_rep1 66 0.102 * 0.08527375 0.11462750
## 767 list2 align2_rep1 67 0.032 * 0.02428000 0.03819875
## 768 list2 align2_rep1 68 0.081 * 0.06645875 0.09435375
## 769 list2 align2_rep1 69 0.158 * 0.14369500 0.17282500
## 770 list2 align2_rep1 70 0.151 * 0.12902000 0.17401375
## 771 list2 align2_rep1 71 0.097 * 0.08522500 0.11235125
## 772 list2 align2_rep1 72 0.080 * 0.06316750 0.09904500
## 773 list2 align2_rep1 73 0.132 * 0.11544375 0.15931250
## 774 list2 align2_rep1 74 0.121 * 0.10524625 0.14646875
## 775 list2 align2_rep1 75 0.121 * 0.10307500 0.14214125
## 776 list2 align2_rep1 76 0.181 * 0.14995500 0.20905875
## 777 list2 align2_rep1 77 0.144 * 0.12443000 0.16092125
## 778 list2 align2_rep1 78 0.105 * 0.08738125 0.11707500
## 779 list2 align2_rep1 79 0.104 * 0.08941125 0.11503875
## 780 list2 align2_rep1 80 0.190 * 0.16077250 0.21812250
## 781 list2 align2_rep1 81 0.237 * 0.21344000 0.26812625
## 782 list2 align2_rep1 82 0.188 * 0.16717375 0.20875875
## 783 list2 align2_rep1 83 0.129 * 0.11344875 0.14160375
## 784 list2 align2_rep1 84 0.099 * 0.07863500 0.11336000
## 785 list2 align2_rep1 85 0.136 * 0.11604000 0.15661375
## 786 list2 align2_rep1 86 0.130 * 0.11096500 0.15516125
## 787 list2 align2_rep1 87 0.093 * 0.08016250 0.11159375
## 788 list2 align2_rep1 88 0.098 * 0.08134000 0.12287750
## 789 list2 align2_rep1 89 0.073 * 0.05889000 0.09063875
## 790 list2 align2_rep1 90 0.046 * 0.03803750 0.05511750
## 791 list2 align2_rep1 91 0.049 * 0.03687750 0.05958250
## 792 list2 align2_rep1 92 0.046 * 0.03441250 0.05563000
## 793 list2 align2_rep1 93 0.052 * 0.03747125 0.05984250
## 794 list2 align2_rep1 94 0.081 * 0.06555875 0.09495750
## 795 list2 align2_rep1 95 0.070 * 0.05671500 0.08811000
## 796 list2 align2_rep1 96 0.074 * 0.05728000 0.08801250
## 797 list2 align2_rep1 97 0.097 * 0.07844500 0.11253750
## 798 list2 align2_rep1 98 0.078 * 0.05905125 0.09622250
## 799 list2 align2_rep1 99 0.063 * 0.04763625 0.08282500
## 800 list2 align2_rep1 100 0.049 * 0.03868625 0.05756875
## 801 list2 align2_rep2 1 0.170 * 0.12947500 0.20166500
## 802 list2 align2_rep2 2 0.202 * 0.16120250 0.26117000
## 803 list2 align2_rep2 3 0.173 * 0.12978125 0.22113875
## 804 list2 align2_rep2 4 0.187 * 0.15419375 0.21866375
## 805 list2 align2_rep2 5 0.366 * 0.30362750 0.45787750
## 806 list2 align2_rep2 6 0.241 * 0.16941375 0.30541750
## 807 list2 align2_rep2 7 0.146 * 0.11179875 0.20317750
## 808 list2 align2_rep2 8 0.207 * 0.15634375 0.25532375
## 809 list2 align2_rep2 9 0.172 * 0.14591875 0.21249125
## 810 list2 align2_rep2 10 0.279 * 0.21269125 0.37419000
## 811 list2 align2_rep2 11 0.266 * 0.19148500 0.38085500
## 812 list2 align2_rep2 12 0.169 * 0.11221500 0.21663875
## 813 list2 align2_rep2 13 0.180 * 0.13122500 0.21977500
## 814 list2 align2_rep2 14 0.031 * 0.02414125 0.03933250
## 815 list2 align2_rep2 15 0.040 * 0.02445000 0.06110000
## 816 list2 align2_rep2 16 0.014 * 0.00715750 0.01960000
## 817 list2 align2_rep2 17 0.000 * 0.00000000 0.00000000
## 818 list2 align2_rep2 18 0.000 * 0.00000000 0.00000000
## 819 list2 align2_rep2 19 0.000 * 0.00000000 0.00000000
## 820 list2 align2_rep2 20 0.000 * 0.00000000 0.00000000
## 821 list2 align2_rep2 21 0.051 * 0.03232125 0.07025250
## 822 list2 align2_rep2 22 0.142 * 0.10855250 0.17818125
## 823 list2 align2_rep2 23 0.095 * 0.06020625 0.13822500
## 824 list2 align2_rep2 24 0.015 * 0.00991875 0.02308125
## 825 list2 align2_rep2 25 0.060 * 0.03067500 0.08332500
## 826 list2 align2_rep2 26 0.139 * 0.10831375 0.17998375
## 827 list2 align2_rep2 27 0.172 * 0.11853500 0.24853000
## 828 list2 align2_rep2 28 0.025 * 0.01334375 0.03443750
## 829 list2 align2_rep2 29 0.000 * 0.00000000 0.00000000
## 830 list2 align2_rep2 30 0.038 * 0.02555500 0.05467250
## 831 list2 align2_rep2 31 0.082 * 0.06300000 0.10858250
## 832 list2 align2_rep2 32 0.250 * 0.18201250 0.33137500
## 833 list2 align2_rep2 33 0.283 * 0.22130375 0.38654750
## 834 list2 align2_rep2 34 0.235 * 0.17036875 0.29596875
## 835 list2 align2_rep2 35 0.216 * 0.14384000 0.27759000
## 836 list2 align2_rep2 36 0.050 * 0.03596750 0.07539500
## 837 list2 align2_rep2 37 0.099 * 0.07068000 0.13176625
## 838 list2 align2_rep2 38 0.294 * 0.22683000 0.39772750
## 839 list2 align2_rep2 39 0.161 * 0.13026750 0.20753000
## 840 list2 align2_rep2 40 0.227 * 0.16600625 0.29163625
## 841 list2 align2_rep2 41 0.128 * 0.10213500 0.15815000
## 842 list2 align2_rep2 42 0.073 * 0.05327125 0.09521000
## 843 list2 align2_rep2 43 0.179 * 0.13973250 0.24285000
## 844 list2 align2_rep2 44 0.115 * 0.07883625 0.16932125
## 845 list2 align2_rep2 45 0.313 * 0.22728000 0.38399875
## 846 list2 align2_rep2 46 0.186 * 0.12108500 0.22822000
## 847 list2 align2_rep2 47 0.110 * 0.08681250 0.13616750
## 848 list2 align2_rep2 48 0.391 * 0.30909625 0.46390375
## 849 list2 align2_rep2 49 0.303 * 0.25731375 0.35984875
## 850 list2 align2_rep2 50 0.190 * 0.13568125 0.23524375
## 851 list2 align2_rep2 51 0.344 * 0.27938250 0.39124500
## 852 list2 align2_rep2 52 0.372 * 0.28593625 0.47257250
## 853 list2 align2_rep2 53 0.188 * 0.14541750 0.25668625
## 854 list2 align2_rep2 54 0.229 * 0.16317250 0.27287500
## 855 list2 align2_rep2 55 0.290 * 0.20855500 0.36727250
## 856 list2 align2_rep2 56 0.189 * 0.14917250 0.25969625
## 857 list2 align2_rep2 57 0.196 * 0.14152000 0.25565500
## 858 list2 align2_rep2 58 0.323 * 0.21247375 0.47892500
## 859 list2 align2_rep2 59 0.205 * 0.13111875 0.25898125
## 860 list2 align2_rep2 60 0.137 * 0.10352875 0.16259625
## 861 list2 align2_rep2 61 0.275 * 0.22016875 0.35498750
## 862 list2 align2_rep2 62 0.384 * 0.31047500 0.45279000
## 863 list2 align2_rep2 63 0.273 * 0.23297875 0.32616625
## 864 list2 align2_rep2 64 0.159 * 0.12138625 0.20147375
## 865 list2 align2_rep2 65 0.312 * 0.25774500 0.36894250
## 866 list2 align2_rep2 66 0.418 * 0.35736250 0.48395875
## 867 list2 align2_rep2 67 0.248 * 0.20012250 0.31602625
## 868 list2 align2_rep2 68 0.083 * 0.06902375 0.09888125
## 869 list2 align2_rep2 69 0.156 * 0.11420250 0.20027250
## 870 list2 align2_rep2 70 0.151 * 0.13185250 0.18761625
## 871 list2 align2_rep2 71 0.169 * 0.12242625 0.20400875
## 872 list2 align2_rep2 72 0.152 * 0.08903000 0.18853000
## 873 list2 align2_rep2 73 0.025 * 0.01690750 0.03619250
## 874 list2 align2_rep2 74 0.032 * 0.02307000 0.04186500
## 875 list2 align2_rep2 75 0.128 * 0.09745500 0.16602500
## 876 list2 align2_rep2 76 0.296 * 0.23658000 0.37885375
## 877 list2 align2_rep2 77 0.259 * 0.19798625 0.30556750
## 878 list2 align2_rep2 78 0.175 * 0.12544500 0.21694375
## 879 list2 align2_rep2 79 0.394 * 0.30637250 0.46873250
## 880 list2 align2_rep2 80 0.385 * 0.29780625 0.48957000
## 881 list2 align2_rep2 81 0.295 * 0.22012250 0.40564250
## 882 list2 align2_rep2 82 0.478 * 0.40222750 0.61692125
## 883 list2 align2_rep2 83 0.252 * 0.17862250 0.33691000
## 884 list2 align2_rep2 84 0.018 * 0.01260000 0.02520000
## 885 list2 align2_rep2 85 0.130 * 0.07912500 0.17038750
## 886 list2 align2_rep2 86 0.243 * 0.20752000 0.30158000
## 887 list2 align2_rep2 87 0.214 * 0.15279250 0.27708750
## 888 list2 align2_rep2 88 0.030 * 0.02257000 0.03711500
## 889 list2 align2_rep2 89 0.147 * 0.10707625 0.19274500
## 890 list2 align2_rep2 90 0.281 * 0.22939125 0.33552625
## 891 list2 align2_rep2 91 0.251 * 0.20710000 0.32966250
## 892 list2 align2_rep2 92 0.145 * 0.09966250 0.18058250
## 893 list2 align2_rep2 93 0.080 * 0.04890000 0.10310000
## 894 list2 align2_rep2 94 0.024 * 0.01649000 0.03022000
## 895 list2 align2_rep2 95 0.088 * 0.07196000 0.10640000
## 896 list2 align2_rep2 96 0.361 * 0.26617875 0.48757250
## 897 list2 align2_rep2 97 0.586 * 0.45113375 0.74314625
## 898 list2 align2_rep2 98 0.313 * 0.25311750 0.38314500
## 899 list2 align2_rep2 99 0.466 * 0.32066375 0.59205000
## 900 list2 align2_rep2 100 0.415 * 0.33138000 0.50820750
## 901 list2 ctrl 1 0.365 * 0.33655875 0.40806625
## 902 list2 ctrl 2 0.296 * 0.26663250 0.32104250
## 903 list2 ctrl 3 0.311 * 0.27189750 0.35307250
## 904 list2 ctrl 4 0.416 * 0.36138250 0.46161875
## 905 list2 ctrl 5 0.536 * 0.48076375 0.59965125
## 906 list2 ctrl 6 0.357 * 0.32319375 0.39623625
## 907 list2 ctrl 7 0.278 * 0.23213625 0.30516500
## 908 list2 ctrl 8 0.340 * 0.27710500 0.40646375
## 909 list2 ctrl 9 0.403 * 0.35632500 0.49927750
## 910 list2 ctrl 10 0.584 * 0.49471125 0.63407750
## 911 list2 ctrl 11 0.500 * 0.43374500 0.55396250
## 912 list2 ctrl 12 0.448 * 0.39611875 0.50072625
## 913 list2 ctrl 13 0.461 * 0.41823000 0.49972250
## 914 list2 ctrl 14 0.493 * 0.44908375 0.54571875
## 915 list2 ctrl 15 0.430 * 0.40356250 0.49190750
## 916 list2 ctrl 16 0.388 * 0.35707500 0.43735750
## 917 list2 ctrl 17 0.466 * 0.41200875 0.50938250
## 918 list2 ctrl 18 0.387 * 0.32255250 0.43760000
## 919 list2 ctrl 19 0.341 * 0.28290500 0.41054875
## 920 list2 ctrl 20 0.379 * 0.33646875 0.40846375
## 921 list2 ctrl 21 0.340 * 0.30740500 0.37326125
## 922 list2 ctrl 22 0.265 * 0.22953000 0.29654625
## 923 list2 ctrl 23 0.323 * 0.27906500 0.35821375
## 924 list2 ctrl 24 0.492 * 0.42567625 0.53518750
## 925 list2 ctrl 25 0.459 * 0.41037375 0.50349125
## 926 list2 ctrl 26 0.288 * 0.26727000 0.31859625
## 927 list2 ctrl 27 0.358 * 0.31941375 0.39092625
## 928 list2 ctrl 28 0.389 * 0.36708000 0.42657500
## 929 list2 ctrl 29 0.373 * 0.32861500 0.41168875
## 930 list2 ctrl 30 0.374 * 0.34958125 0.40368125
## 931 list2 ctrl 31 0.444 * 0.39225000 0.49747625
## 932 list2 ctrl 32 0.381 * 0.33447500 0.44210625
## 933 list2 ctrl 33 0.481 * 0.43601250 0.53397125
## 934 list2 ctrl 34 0.519 * 0.47233000 0.57550375
## 935 list2 ctrl 35 0.576 * 0.51014625 0.63862375
## 936 list2 ctrl 36 0.589 * 0.52072375 0.67489625
## 937 list2 ctrl 37 0.744 * 0.67554375 0.88336750
## 938 list2 ctrl 38 0.742 * 0.66523000 0.83539875
## 939 list2 ctrl 39 0.729 * 0.67287875 0.79096375
## 940 list2 ctrl 40 0.870 * 0.76487750 0.95000000
## 941 list2 ctrl 41 0.719 * 0.65916625 0.78219250
## 942 list2 ctrl 42 0.616 * 0.54213000 0.69888375
## 943 list2 ctrl 43 0.791 * 0.67088500 0.88395250
## 944 list2 ctrl 44 0.907 * 0.82698000 0.98935000
## 945 list2 ctrl 45 0.984 * 0.89462250 1.07022500
## 946 list2 ctrl 46 1.127 * 1.03708375 1.22816375
## 947 list2 ctrl 47 1.295 * 1.19999125 1.40812500
## 948 list2 ctrl 48 1.275 * 1.17031750 1.39350250
## 949 list2 ctrl 49 1.160 * 1.01661250 1.26384750
## 950 list2 ctrl 50 1.024 * 0.92919625 1.10214750
## 951 list2 ctrl 51 0.964 * 0.86169500 1.05108375
## 952 list2 ctrl 52 0.858 * 0.78890125 0.94288625
## 953 list2 ctrl 53 0.808 * 0.75355500 0.87126125
## 954 list2 ctrl 54 0.888 * 0.83680375 0.94939750
## 955 list2 ctrl 55 0.897 * 0.84335625 0.95057250
## 956 list2 ctrl 56 0.697 * 0.63118000 0.75161875
## 957 list2 ctrl 57 0.753 * 0.65949625 0.83672750
## 958 list2 ctrl 58 0.737 * 0.65431500 0.81789750
## 959 list2 ctrl 59 0.552 * 0.50228375 0.60522875
## 960 list2 ctrl 60 0.520 * 0.46430750 0.57451375
## 961 list2 ctrl 61 0.798 * 0.74610750 0.86052375
## 962 list2 ctrl 62 1.014 * 0.94678125 1.12346250
## 963 list2 ctrl 63 0.752 * 0.68453125 0.81110750
## 964 list2 ctrl 64 0.589 * 0.53658750 0.63468000
## 965 list2 ctrl 65 0.544 * 0.50771875 0.62701750
## 966 list2 ctrl 66 0.483 * 0.44558500 0.52777000
## 967 list2 ctrl 67 0.566 * 0.52498000 0.61449000
## 968 list2 ctrl 68 0.663 * 0.62279000 0.72062625
## 969 list2 ctrl 69 0.521 * 0.48032625 0.55781625
## 970 list2 ctrl 70 0.484 * 0.39682250 0.54094375
## 971 list2 ctrl 71 0.412 * 0.36366750 0.48286875
## 972 list2 ctrl 72 0.385 * 0.34053000 0.44250375
## 973 list2 ctrl 73 0.407 * 0.36316375 0.44614250
## 974 list2 ctrl 74 0.540 * 0.49153000 0.59080500
## 975 list2 ctrl 75 0.574 * 0.51531375 0.63220625
## 976 list2 ctrl 76 0.625 * 0.55535250 0.70138125
## 977 list2 ctrl 77 0.655 * 0.54349625 0.73879125
## 978 list2 ctrl 78 0.558 * 0.50804125 0.61436000
## 979 list2 ctrl 79 0.557 * 0.49919000 0.61815500
## 980 list2 ctrl 80 0.529 * 0.48136625 0.57790250
## 981 list2 ctrl 81 0.586 * 0.53538125 0.64846750
## 982 list2 ctrl 82 0.664 * 0.58009750 0.75707125
## 983 list2 ctrl 83 0.529 * 0.45893500 0.59803250
## 984 list2 ctrl 84 0.324 * 0.28662500 0.35028125
## 985 list2 ctrl 85 0.433 * 0.40857875 0.47152625
## 986 list2 ctrl 86 0.504 * 0.47443000 0.53364125
## 987 list2 ctrl 87 0.347 * 0.31964875 0.38402250
## 988 list2 ctrl 88 0.383 * 0.35306125 0.41493750
## 989 list2 ctrl 89 0.415 * 0.37648000 0.46087625
## 990 list2 ctrl 90 0.349 * 0.31078375 0.37814750
## 991 list2 ctrl 91 0.478 * 0.44627875 0.52290250
## 992 list2 ctrl 92 0.593 * 0.54600750 0.67730750
## 993 list2 ctrl 93 0.384 * 0.33653250 0.42527750
## 994 list2 ctrl 94 0.390 * 0.35802250 0.42197875
## 995 list2 ctrl 95 0.544 * 0.48059750 0.62634500
## 996 list2 ctrl 96 0.513 * 0.42739625 0.58603625
## 997 list2 ctrl 97 0.615 * 0.51542875 0.71969375
## 998 list2 ctrl 98 0.705 * 0.52907500 0.82677625
## 999 list2 ctrl 99 0.708 * 0.58902500 0.89336125
## 1000 list2 ctrl 100 0.652 * 0.55527250 0.72730750
## group
## 1 align1_rep1_list1
## 2 align1_rep1_list1
## 3 align1_rep1_list1
## 4 align1_rep1_list1
## 5 align1_rep1_list1
## 6 align1_rep1_list1
## 7 align1_rep1_list1
## 8 align1_rep1_list1
## 9 align1_rep1_list1
## 10 align1_rep1_list1
## 11 align1_rep1_list1
## 12 align1_rep1_list1
## 13 align1_rep1_list1
## 14 align1_rep1_list1
## 15 align1_rep1_list1
## 16 align1_rep1_list1
## 17 align1_rep1_list1
## 18 align1_rep1_list1
## 19 align1_rep1_list1
## 20 align1_rep1_list1
## 21 align1_rep1_list1
## 22 align1_rep1_list1
## 23 align1_rep1_list1
## 24 align1_rep1_list1
## 25 align1_rep1_list1
## 26 align1_rep1_list1
## 27 align1_rep1_list1
## 28 align1_rep1_list1
## 29 align1_rep1_list1
## 30 align1_rep1_list1
## 31 align1_rep1_list1
## 32 align1_rep1_list1
## 33 align1_rep1_list1
## 34 align1_rep1_list1
## 35 align1_rep1_list1
## 36 align1_rep1_list1
## 37 align1_rep1_list1
## 38 align1_rep1_list1
## 39 align1_rep1_list1
## 40 align1_rep1_list1
## 41 align1_rep1_list1
## 42 align1_rep1_list1
## 43 align1_rep1_list1
## 44 align1_rep1_list1
## 45 align1_rep1_list1
## 46 align1_rep1_list1
## 47 align1_rep1_list1
## 48 align1_rep1_list1
## 49 align1_rep1_list1
## 50 align1_rep1_list1
## 51 align1_rep1_list1
## 52 align1_rep1_list1
## 53 align1_rep1_list1
## 54 align1_rep1_list1
## 55 align1_rep1_list1
## 56 align1_rep1_list1
## 57 align1_rep1_list1
## 58 align1_rep1_list1
## 59 align1_rep1_list1
## 60 align1_rep1_list1
## 61 align1_rep1_list1
## 62 align1_rep1_list1
## 63 align1_rep1_list1
## 64 align1_rep1_list1
## 65 align1_rep1_list1
## 66 align1_rep1_list1
## 67 align1_rep1_list1
## 68 align1_rep1_list1
## 69 align1_rep1_list1
## 70 align1_rep1_list1
## 71 align1_rep1_list1
## 72 align1_rep1_list1
## 73 align1_rep1_list1
## 74 align1_rep1_list1
## 75 align1_rep1_list1
## 76 align1_rep1_list1
## 77 align1_rep1_list1
## 78 align1_rep1_list1
## 79 align1_rep1_list1
## 80 align1_rep1_list1
## 81 align1_rep1_list1
## 82 align1_rep1_list1
## 83 align1_rep1_list1
## 84 align1_rep1_list1
## 85 align1_rep1_list1
## 86 align1_rep1_list1
## 87 align1_rep1_list1
## 88 align1_rep1_list1
## 89 align1_rep1_list1
## 90 align1_rep1_list1
## 91 align1_rep1_list1
## 92 align1_rep1_list1
## 93 align1_rep1_list1
## 94 align1_rep1_list1
## 95 align1_rep1_list1
## 96 align1_rep1_list1
## 97 align1_rep1_list1
## 98 align1_rep1_list1
## 99 align1_rep1_list1
## 100 align1_rep1_list1
## 101 align1_rep2_list1
## 102 align1_rep2_list1
## 103 align1_rep2_list1
## 104 align1_rep2_list1
## 105 align1_rep2_list1
## 106 align1_rep2_list1
## 107 align1_rep2_list1
## 108 align1_rep2_list1
## 109 align1_rep2_list1
## 110 align1_rep2_list1
## 111 align1_rep2_list1
## 112 align1_rep2_list1
## 113 align1_rep2_list1
## 114 align1_rep2_list1
## 115 align1_rep2_list1
## 116 align1_rep2_list1
## 117 align1_rep2_list1
## 118 align1_rep2_list1
## 119 align1_rep2_list1
## 120 align1_rep2_list1
## 121 align1_rep2_list1
## 122 align1_rep2_list1
## 123 align1_rep2_list1
## 124 align1_rep2_list1
## 125 align1_rep2_list1
## 126 align1_rep2_list1
## 127 align1_rep2_list1
## 128 align1_rep2_list1
## 129 align1_rep2_list1
## 130 align1_rep2_list1
## 131 align1_rep2_list1
## 132 align1_rep2_list1
## 133 align1_rep2_list1
## 134 align1_rep2_list1
## 135 align1_rep2_list1
## 136 align1_rep2_list1
## 137 align1_rep2_list1
## 138 align1_rep2_list1
## 139 align1_rep2_list1
## 140 align1_rep2_list1
## 141 align1_rep2_list1
## 142 align1_rep2_list1
## 143 align1_rep2_list1
## 144 align1_rep2_list1
## 145 align1_rep2_list1
## 146 align1_rep2_list1
## 147 align1_rep2_list1
## 148 align1_rep2_list1
## 149 align1_rep2_list1
## 150 align1_rep2_list1
## 151 align1_rep2_list1
## 152 align1_rep2_list1
## 153 align1_rep2_list1
## 154 align1_rep2_list1
## 155 align1_rep2_list1
## 156 align1_rep2_list1
## 157 align1_rep2_list1
## 158 align1_rep2_list1
## 159 align1_rep2_list1
## 160 align1_rep2_list1
## 161 align1_rep2_list1
## 162 align1_rep2_list1
## 163 align1_rep2_list1
## 164 align1_rep2_list1
## 165 align1_rep2_list1
## 166 align1_rep2_list1
## 167 align1_rep2_list1
## 168 align1_rep2_list1
## 169 align1_rep2_list1
## 170 align1_rep2_list1
## 171 align1_rep2_list1
## 172 align1_rep2_list1
## 173 align1_rep2_list1
## 174 align1_rep2_list1
## 175 align1_rep2_list1
## 176 align1_rep2_list1
## 177 align1_rep2_list1
## 178 align1_rep2_list1
## 179 align1_rep2_list1
## 180 align1_rep2_list1
## 181 align1_rep2_list1
## 182 align1_rep2_list1
## 183 align1_rep2_list1
## 184 align1_rep2_list1
## 185 align1_rep2_list1
## 186 align1_rep2_list1
## 187 align1_rep2_list1
## 188 align1_rep2_list1
## 189 align1_rep2_list1
## 190 align1_rep2_list1
## 191 align1_rep2_list1
## 192 align1_rep2_list1
## 193 align1_rep2_list1
## 194 align1_rep2_list1
## 195 align1_rep2_list1
## 196 align1_rep2_list1
## 197 align1_rep2_list1
## 198 align1_rep2_list1
## 199 align1_rep2_list1
## 200 align1_rep2_list1
## 201 align2_rep1_list1
## 202 align2_rep1_list1
## 203 align2_rep1_list1
## 204 align2_rep1_list1
## 205 align2_rep1_list1
## 206 align2_rep1_list1
## 207 align2_rep1_list1
## 208 align2_rep1_list1
## 209 align2_rep1_list1
## 210 align2_rep1_list1
## 211 align2_rep1_list1
## 212 align2_rep1_list1
## 213 align2_rep1_list1
## 214 align2_rep1_list1
## 215 align2_rep1_list1
## 216 align2_rep1_list1
## 217 align2_rep1_list1
## 218 align2_rep1_list1
## 219 align2_rep1_list1
## 220 align2_rep1_list1
## 221 align2_rep1_list1
## 222 align2_rep1_list1
## 223 align2_rep1_list1
## 224 align2_rep1_list1
## 225 align2_rep1_list1
## 226 align2_rep1_list1
## 227 align2_rep1_list1
## 228 align2_rep1_list1
## 229 align2_rep1_list1
## 230 align2_rep1_list1
## 231 align2_rep1_list1
## 232 align2_rep1_list1
## 233 align2_rep1_list1
## 234 align2_rep1_list1
## 235 align2_rep1_list1
## 236 align2_rep1_list1
## 237 align2_rep1_list1
## 238 align2_rep1_list1
## 239 align2_rep1_list1
## 240 align2_rep1_list1
## 241 align2_rep1_list1
## 242 align2_rep1_list1
## 243 align2_rep1_list1
## 244 align2_rep1_list1
## 245 align2_rep1_list1
## 246 align2_rep1_list1
## 247 align2_rep1_list1
## 248 align2_rep1_list1
## 249 align2_rep1_list1
## 250 align2_rep1_list1
## 251 align2_rep1_list1
## 252 align2_rep1_list1
## 253 align2_rep1_list1
## 254 align2_rep1_list1
## 255 align2_rep1_list1
## 256 align2_rep1_list1
## 257 align2_rep1_list1
## 258 align2_rep1_list1
## 259 align2_rep1_list1
## 260 align2_rep1_list1
## 261 align2_rep1_list1
## 262 align2_rep1_list1
## 263 align2_rep1_list1
## 264 align2_rep1_list1
## 265 align2_rep1_list1
## 266 align2_rep1_list1
## 267 align2_rep1_list1
## 268 align2_rep1_list1
## 269 align2_rep1_list1
## 270 align2_rep1_list1
## 271 align2_rep1_list1
## 272 align2_rep1_list1
## 273 align2_rep1_list1
## 274 align2_rep1_list1
## 275 align2_rep1_list1
## 276 align2_rep1_list1
## 277 align2_rep1_list1
## 278 align2_rep1_list1
## 279 align2_rep1_list1
## 280 align2_rep1_list1
## 281 align2_rep1_list1
## 282 align2_rep1_list1
## 283 align2_rep1_list1
## 284 align2_rep1_list1
## 285 align2_rep1_list1
## 286 align2_rep1_list1
## 287 align2_rep1_list1
## 288 align2_rep1_list1
## 289 align2_rep1_list1
## 290 align2_rep1_list1
## 291 align2_rep1_list1
## 292 align2_rep1_list1
## 293 align2_rep1_list1
## 294 align2_rep1_list1
## 295 align2_rep1_list1
## 296 align2_rep1_list1
## 297 align2_rep1_list1
## 298 align2_rep1_list1
## 299 align2_rep1_list1
## 300 align2_rep1_list1
## 301 align2_rep2_list1
## 302 align2_rep2_list1
## 303 align2_rep2_list1
## 304 align2_rep2_list1
## 305 align2_rep2_list1
## 306 align2_rep2_list1
## 307 align2_rep2_list1
## 308 align2_rep2_list1
## 309 align2_rep2_list1
## 310 align2_rep2_list1
## 311 align2_rep2_list1
## 312 align2_rep2_list1
## 313 align2_rep2_list1
## 314 align2_rep2_list1
## 315 align2_rep2_list1
## 316 align2_rep2_list1
## 317 align2_rep2_list1
## 318 align2_rep2_list1
## 319 align2_rep2_list1
## 320 align2_rep2_list1
## 321 align2_rep2_list1
## 322 align2_rep2_list1
## 323 align2_rep2_list1
## 324 align2_rep2_list1
## 325 align2_rep2_list1
## 326 align2_rep2_list1
## 327 align2_rep2_list1
## 328 align2_rep2_list1
## 329 align2_rep2_list1
## 330 align2_rep2_list1
## 331 align2_rep2_list1
## 332 align2_rep2_list1
## 333 align2_rep2_list1
## 334 align2_rep2_list1
## 335 align2_rep2_list1
## 336 align2_rep2_list1
## 337 align2_rep2_list1
## 338 align2_rep2_list1
## 339 align2_rep2_list1
## 340 align2_rep2_list1
## 341 align2_rep2_list1
## 342 align2_rep2_list1
## 343 align2_rep2_list1
## 344 align2_rep2_list1
## 345 align2_rep2_list1
## 346 align2_rep2_list1
## 347 align2_rep2_list1
## 348 align2_rep2_list1
## 349 align2_rep2_list1
## 350 align2_rep2_list1
## 351 align2_rep2_list1
## 352 align2_rep2_list1
## 353 align2_rep2_list1
## 354 align2_rep2_list1
## 355 align2_rep2_list1
## 356 align2_rep2_list1
## 357 align2_rep2_list1
## 358 align2_rep2_list1
## 359 align2_rep2_list1
## 360 align2_rep2_list1
## 361 align2_rep2_list1
## 362 align2_rep2_list1
## 363 align2_rep2_list1
## 364 align2_rep2_list1
## 365 align2_rep2_list1
## 366 align2_rep2_list1
## 367 align2_rep2_list1
## 368 align2_rep2_list1
## 369 align2_rep2_list1
## 370 align2_rep2_list1
## 371 align2_rep2_list1
## 372 align2_rep2_list1
## 373 align2_rep2_list1
## 374 align2_rep2_list1
## 375 align2_rep2_list1
## 376 align2_rep2_list1
## 377 align2_rep2_list1
## 378 align2_rep2_list1
## 379 align2_rep2_list1
## 380 align2_rep2_list1
## 381 align2_rep2_list1
## 382 align2_rep2_list1
## 383 align2_rep2_list1
## 384 align2_rep2_list1
## 385 align2_rep2_list1
## 386 align2_rep2_list1
## 387 align2_rep2_list1
## 388 align2_rep2_list1
## 389 align2_rep2_list1
## 390 align2_rep2_list1
## 391 align2_rep2_list1
## 392 align2_rep2_list1
## 393 align2_rep2_list1
## 394 align2_rep2_list1
## 395 align2_rep2_list1
## 396 align2_rep2_list1
## 397 align2_rep2_list1
## 398 align2_rep2_list1
## 399 align2_rep2_list1
## 400 align2_rep2_list1
## 401 ctrl_list1
## 402 ctrl_list1
## 403 ctrl_list1
## 404 ctrl_list1
## 405 ctrl_list1
## 406 ctrl_list1
## 407 ctrl_list1
## 408 ctrl_list1
## 409 ctrl_list1
## 410 ctrl_list1
## 411 ctrl_list1
## 412 ctrl_list1
## 413 ctrl_list1
## 414 ctrl_list1
## 415 ctrl_list1
## 416 ctrl_list1
## 417 ctrl_list1
## 418 ctrl_list1
## 419 ctrl_list1
## 420 ctrl_list1
## 421 ctrl_list1
## 422 ctrl_list1
## 423 ctrl_list1
## 424 ctrl_list1
## 425 ctrl_list1
## 426 ctrl_list1
## 427 ctrl_list1
## 428 ctrl_list1
## 429 ctrl_list1
## 430 ctrl_list1
## 431 ctrl_list1
## 432 ctrl_list1
## 433 ctrl_list1
## 434 ctrl_list1
## 435 ctrl_list1
## 436 ctrl_list1
## 437 ctrl_list1
## 438 ctrl_list1
## 439 ctrl_list1
## 440 ctrl_list1
## 441 ctrl_list1
## 442 ctrl_list1
## 443 ctrl_list1
## 444 ctrl_list1
## 445 ctrl_list1
## 446 ctrl_list1
## 447 ctrl_list1
## 448 ctrl_list1
## 449 ctrl_list1
## 450 ctrl_list1
## 451 ctrl_list1
## 452 ctrl_list1
## 453 ctrl_list1
## 454 ctrl_list1
## 455 ctrl_list1
## 456 ctrl_list1
## 457 ctrl_list1
## 458 ctrl_list1
## 459 ctrl_list1
## 460 ctrl_list1
## 461 ctrl_list1
## 462 ctrl_list1
## 463 ctrl_list1
## 464 ctrl_list1
## 465 ctrl_list1
## 466 ctrl_list1
## 467 ctrl_list1
## 468 ctrl_list1
## 469 ctrl_list1
## 470 ctrl_list1
## 471 ctrl_list1
## 472 ctrl_list1
## 473 ctrl_list1
## 474 ctrl_list1
## 475 ctrl_list1
## 476 ctrl_list1
## 477 ctrl_list1
## 478 ctrl_list1
## 479 ctrl_list1
## 480 ctrl_list1
## 481 ctrl_list1
## 482 ctrl_list1
## 483 ctrl_list1
## 484 ctrl_list1
## 485 ctrl_list1
## 486 ctrl_list1
## 487 ctrl_list1
## 488 ctrl_list1
## 489 ctrl_list1
## 490 ctrl_list1
## 491 ctrl_list1
## 492 ctrl_list1
## 493 ctrl_list1
## 494 ctrl_list1
## 495 ctrl_list1
## 496 ctrl_list1
## 497 ctrl_list1
## 498 ctrl_list1
## 499 ctrl_list1
## 500 ctrl_list1
## 501 align1_rep1_list2
## 502 align1_rep1_list2
## 503 align1_rep1_list2
## 504 align1_rep1_list2
## 505 align1_rep1_list2
## 506 align1_rep1_list2
## 507 align1_rep1_list2
## 508 align1_rep1_list2
## 509 align1_rep1_list2
## 510 align1_rep1_list2
## 511 align1_rep1_list2
## 512 align1_rep1_list2
## 513 align1_rep1_list2
## 514 align1_rep1_list2
## 515 align1_rep1_list2
## 516 align1_rep1_list2
## 517 align1_rep1_list2
## 518 align1_rep1_list2
## 519 align1_rep1_list2
## 520 align1_rep1_list2
## 521 align1_rep1_list2
## 522 align1_rep1_list2
## 523 align1_rep1_list2
## 524 align1_rep1_list2
## 525 align1_rep1_list2
## 526 align1_rep1_list2
## 527 align1_rep1_list2
## 528 align1_rep1_list2
## 529 align1_rep1_list2
## 530 align1_rep1_list2
## 531 align1_rep1_list2
## 532 align1_rep1_list2
## 533 align1_rep1_list2
## 534 align1_rep1_list2
## 535 align1_rep1_list2
## 536 align1_rep1_list2
## 537 align1_rep1_list2
## 538 align1_rep1_list2
## 539 align1_rep1_list2
## 540 align1_rep1_list2
## 541 align1_rep1_list2
## 542 align1_rep1_list2
## 543 align1_rep1_list2
## 544 align1_rep1_list2
## 545 align1_rep1_list2
## 546 align1_rep1_list2
## 547 align1_rep1_list2
## 548 align1_rep1_list2
## 549 align1_rep1_list2
## 550 align1_rep1_list2
## 551 align1_rep1_list2
## 552 align1_rep1_list2
## 553 align1_rep1_list2
## 554 align1_rep1_list2
## 555 align1_rep1_list2
## 556 align1_rep1_list2
## 557 align1_rep1_list2
## 558 align1_rep1_list2
## 559 align1_rep1_list2
## 560 align1_rep1_list2
## 561 align1_rep1_list2
## 562 align1_rep1_list2
## 563 align1_rep1_list2
## 564 align1_rep1_list2
## 565 align1_rep1_list2
## 566 align1_rep1_list2
## 567 align1_rep1_list2
## 568 align1_rep1_list2
## 569 align1_rep1_list2
## 570 align1_rep1_list2
## 571 align1_rep1_list2
## 572 align1_rep1_list2
## 573 align1_rep1_list2
## 574 align1_rep1_list2
## 575 align1_rep1_list2
## 576 align1_rep1_list2
## 577 align1_rep1_list2
## 578 align1_rep1_list2
## 579 align1_rep1_list2
## 580 align1_rep1_list2
## 581 align1_rep1_list2
## 582 align1_rep1_list2
## 583 align1_rep1_list2
## 584 align1_rep1_list2
## 585 align1_rep1_list2
## 586 align1_rep1_list2
## 587 align1_rep1_list2
## 588 align1_rep1_list2
## 589 align1_rep1_list2
## 590 align1_rep1_list2
## 591 align1_rep1_list2
## 592 align1_rep1_list2
## 593 align1_rep1_list2
## 594 align1_rep1_list2
## 595 align1_rep1_list2
## 596 align1_rep1_list2
## 597 align1_rep1_list2
## 598 align1_rep1_list2
## 599 align1_rep1_list2
## 600 align1_rep1_list2
## 601 align1_rep2_list2
## 602 align1_rep2_list2
## 603 align1_rep2_list2
## 604 align1_rep2_list2
## 605 align1_rep2_list2
## 606 align1_rep2_list2
## 607 align1_rep2_list2
## 608 align1_rep2_list2
## 609 align1_rep2_list2
## 610 align1_rep2_list2
## 611 align1_rep2_list2
## 612 align1_rep2_list2
## 613 align1_rep2_list2
## 614 align1_rep2_list2
## 615 align1_rep2_list2
## 616 align1_rep2_list2
## 617 align1_rep2_list2
## 618 align1_rep2_list2
## 619 align1_rep2_list2
## 620 align1_rep2_list2
## 621 align1_rep2_list2
## 622 align1_rep2_list2
## 623 align1_rep2_list2
## 624 align1_rep2_list2
## 625 align1_rep2_list2
## 626 align1_rep2_list2
## 627 align1_rep2_list2
## 628 align1_rep2_list2
## 629 align1_rep2_list2
## 630 align1_rep2_list2
## 631 align1_rep2_list2
## 632 align1_rep2_list2
## 633 align1_rep2_list2
## 634 align1_rep2_list2
## 635 align1_rep2_list2
## 636 align1_rep2_list2
## 637 align1_rep2_list2
## 638 align1_rep2_list2
## 639 align1_rep2_list2
## 640 align1_rep2_list2
## 641 align1_rep2_list2
## 642 align1_rep2_list2
## 643 align1_rep2_list2
## 644 align1_rep2_list2
## 645 align1_rep2_list2
## 646 align1_rep2_list2
## 647 align1_rep2_list2
## 648 align1_rep2_list2
## 649 align1_rep2_list2
## 650 align1_rep2_list2
## 651 align1_rep2_list2
## 652 align1_rep2_list2
## 653 align1_rep2_list2
## 654 align1_rep2_list2
## 655 align1_rep2_list2
## 656 align1_rep2_list2
## 657 align1_rep2_list2
## 658 align1_rep2_list2
## 659 align1_rep2_list2
## 660 align1_rep2_list2
## 661 align1_rep2_list2
## 662 align1_rep2_list2
## 663 align1_rep2_list2
## 664 align1_rep2_list2
## 665 align1_rep2_list2
## 666 align1_rep2_list2
## 667 align1_rep2_list2
## 668 align1_rep2_list2
## 669 align1_rep2_list2
## 670 align1_rep2_list2
## 671 align1_rep2_list2
## 672 align1_rep2_list2
## 673 align1_rep2_list2
## 674 align1_rep2_list2
## 675 align1_rep2_list2
## 676 align1_rep2_list2
## 677 align1_rep2_list2
## 678 align1_rep2_list2
## 679 align1_rep2_list2
## 680 align1_rep2_list2
## 681 align1_rep2_list2
## 682 align1_rep2_list2
## 683 align1_rep2_list2
## 684 align1_rep2_list2
## 685 align1_rep2_list2
## 686 align1_rep2_list2
## 687 align1_rep2_list2
## 688 align1_rep2_list2
## 689 align1_rep2_list2
## 690 align1_rep2_list2
## 691 align1_rep2_list2
## 692 align1_rep2_list2
## 693 align1_rep2_list2
## 694 align1_rep2_list2
## 695 align1_rep2_list2
## 696 align1_rep2_list2
## 697 align1_rep2_list2
## 698 align1_rep2_list2
## 699 align1_rep2_list2
## 700 align1_rep2_list2
## 701 align2_rep1_list2
## 702 align2_rep1_list2
## 703 align2_rep1_list2
## 704 align2_rep1_list2
## 705 align2_rep1_list2
## 706 align2_rep1_list2
## 707 align2_rep1_list2
## 708 align2_rep1_list2
## 709 align2_rep1_list2
## 710 align2_rep1_list2
## 711 align2_rep1_list2
## 712 align2_rep1_list2
## 713 align2_rep1_list2
## 714 align2_rep1_list2
## 715 align2_rep1_list2
## 716 align2_rep1_list2
## 717 align2_rep1_list2
## 718 align2_rep1_list2
## 719 align2_rep1_list2
## 720 align2_rep1_list2
## 721 align2_rep1_list2
## 722 align2_rep1_list2
## 723 align2_rep1_list2
## 724 align2_rep1_list2
## 725 align2_rep1_list2
## 726 align2_rep1_list2
## 727 align2_rep1_list2
## 728 align2_rep1_list2
## 729 align2_rep1_list2
## 730 align2_rep1_list2
## 731 align2_rep1_list2
## 732 align2_rep1_list2
## 733 align2_rep1_list2
## 734 align2_rep1_list2
## 735 align2_rep1_list2
## 736 align2_rep1_list2
## 737 align2_rep1_list2
## 738 align2_rep1_list2
## 739 align2_rep1_list2
## 740 align2_rep1_list2
## 741 align2_rep1_list2
## 742 align2_rep1_list2
## 743 align2_rep1_list2
## 744 align2_rep1_list2
## 745 align2_rep1_list2
## 746 align2_rep1_list2
## 747 align2_rep1_list2
## 748 align2_rep1_list2
## 749 align2_rep1_list2
## 750 align2_rep1_list2
## 751 align2_rep1_list2
## 752 align2_rep1_list2
## 753 align2_rep1_list2
## 754 align2_rep1_list2
## 755 align2_rep1_list2
## 756 align2_rep1_list2
## 757 align2_rep1_list2
## 758 align2_rep1_list2
## 759 align2_rep1_list2
## 760 align2_rep1_list2
## 761 align2_rep1_list2
## 762 align2_rep1_list2
## 763 align2_rep1_list2
## 764 align2_rep1_list2
## 765 align2_rep1_list2
## 766 align2_rep1_list2
## 767 align2_rep1_list2
## 768 align2_rep1_list2
## 769 align2_rep1_list2
## 770 align2_rep1_list2
## 771 align2_rep1_list2
## 772 align2_rep1_list2
## 773 align2_rep1_list2
## 774 align2_rep1_list2
## 775 align2_rep1_list2
## 776 align2_rep1_list2
## 777 align2_rep1_list2
## 778 align2_rep1_list2
## 779 align2_rep1_list2
## 780 align2_rep1_list2
## 781 align2_rep1_list2
## 782 align2_rep1_list2
## 783 align2_rep1_list2
## 784 align2_rep1_list2
## 785 align2_rep1_list2
## 786 align2_rep1_list2
## 787 align2_rep1_list2
## 788 align2_rep1_list2
## 789 align2_rep1_list2
## 790 align2_rep1_list2
## 791 align2_rep1_list2
## 792 align2_rep1_list2
## 793 align2_rep1_list2
## 794 align2_rep1_list2
## 795 align2_rep1_list2
## 796 align2_rep1_list2
## 797 align2_rep1_list2
## 798 align2_rep1_list2
## 799 align2_rep1_list2
## 800 align2_rep1_list2
## 801 align2_rep2_list2
## 802 align2_rep2_list2
## 803 align2_rep2_list2
## 804 align2_rep2_list2
## 805 align2_rep2_list2
## 806 align2_rep2_list2
## 807 align2_rep2_list2
## 808 align2_rep2_list2
## 809 align2_rep2_list2
## 810 align2_rep2_list2
## 811 align2_rep2_list2
## 812 align2_rep2_list2
## 813 align2_rep2_list2
## 814 align2_rep2_list2
## 815 align2_rep2_list2
## 816 align2_rep2_list2
## 817 align2_rep2_list2
## 818 align2_rep2_list2
## 819 align2_rep2_list2
## 820 align2_rep2_list2
## 821 align2_rep2_list2
## 822 align2_rep2_list2
## 823 align2_rep2_list2
## 824 align2_rep2_list2
## 825 align2_rep2_list2
## 826 align2_rep2_list2
## 827 align2_rep2_list2
## 828 align2_rep2_list2
## 829 align2_rep2_list2
## 830 align2_rep2_list2
## 831 align2_rep2_list2
## 832 align2_rep2_list2
## 833 align2_rep2_list2
## 834 align2_rep2_list2
## 835 align2_rep2_list2
## 836 align2_rep2_list2
## 837 align2_rep2_list2
## 838 align2_rep2_list2
## 839 align2_rep2_list2
## 840 align2_rep2_list2
## 841 align2_rep2_list2
## 842 align2_rep2_list2
## 843 align2_rep2_list2
## 844 align2_rep2_list2
## 845 align2_rep2_list2
## 846 align2_rep2_list2
## 847 align2_rep2_list2
## 848 align2_rep2_list2
## 849 align2_rep2_list2
## 850 align2_rep2_list2
## 851 align2_rep2_list2
## 852 align2_rep2_list2
## 853 align2_rep2_list2
## 854 align2_rep2_list2
## 855 align2_rep2_list2
## 856 align2_rep2_list2
## 857 align2_rep2_list2
## 858 align2_rep2_list2
## 859 align2_rep2_list2
## 860 align2_rep2_list2
## 861 align2_rep2_list2
## 862 align2_rep2_list2
## 863 align2_rep2_list2
## 864 align2_rep2_list2
## 865 align2_rep2_list2
## 866 align2_rep2_list2
## 867 align2_rep2_list2
## 868 align2_rep2_list2
## 869 align2_rep2_list2
## 870 align2_rep2_list2
## 871 align2_rep2_list2
## 872 align2_rep2_list2
## 873 align2_rep2_list2
## 874 align2_rep2_list2
## 875 align2_rep2_list2
## 876 align2_rep2_list2
## 877 align2_rep2_list2
## 878 align2_rep2_list2
## 879 align2_rep2_list2
## 880 align2_rep2_list2
## 881 align2_rep2_list2
## 882 align2_rep2_list2
## 883 align2_rep2_list2
## 884 align2_rep2_list2
## 885 align2_rep2_list2
## 886 align2_rep2_list2
## 887 align2_rep2_list2
## 888 align2_rep2_list2
## 889 align2_rep2_list2
## 890 align2_rep2_list2
## 891 align2_rep2_list2
## 892 align2_rep2_list2
## 893 align2_rep2_list2
## 894 align2_rep2_list2
## 895 align2_rep2_list2
## 896 align2_rep2_list2
## 897 align2_rep2_list2
## 898 align2_rep2_list2
## 899 align2_rep2_list2
## 900 align2_rep2_list2
## 901 ctrl_list2
## 902 ctrl_list2
## 903 ctrl_list2
## 904 ctrl_list2
## 905 ctrl_list2
## 906 ctrl_list2
## 907 ctrl_list2
## 908 ctrl_list2
## 909 ctrl_list2
## 910 ctrl_list2
## 911 ctrl_list2
## 912 ctrl_list2
## 913 ctrl_list2
## 914 ctrl_list2
## 915 ctrl_list2
## 916 ctrl_list2
## 917 ctrl_list2
## 918 ctrl_list2
## 919 ctrl_list2
## 920 ctrl_list2
## 921 ctrl_list2
## 922 ctrl_list2
## 923 ctrl_list2
## 924 ctrl_list2
## 925 ctrl_list2
## 926 ctrl_list2
## 927 ctrl_list2
## 928 ctrl_list2
## 929 ctrl_list2
## 930 ctrl_list2
## 931 ctrl_list2
## 932 ctrl_list2
## 933 ctrl_list2
## 934 ctrl_list2
## 935 ctrl_list2
## 936 ctrl_list2
## 937 ctrl_list2
## 938 ctrl_list2
## 939 ctrl_list2
## 940 ctrl_list2
## 941 ctrl_list2
## 942 ctrl_list2
## 943 ctrl_list2
## 944 ctrl_list2
## 945 ctrl_list2
## 946 ctrl_list2
## 947 ctrl_list2
## 948 ctrl_list2
## 949 ctrl_list2
## 950 ctrl_list2
## 951 ctrl_list2
## 952 ctrl_list2
## 953 ctrl_list2
## 954 ctrl_list2
## 955 ctrl_list2
## 956 ctrl_list2
## 957 ctrl_list2
## 958 ctrl_list2
## 959 ctrl_list2
## 960 ctrl_list2
## 961 ctrl_list2
## 962 ctrl_list2
## 963 ctrl_list2
## 964 ctrl_list2
## 965 ctrl_list2
## 966 ctrl_list2
## 967 ctrl_list2
## 968 ctrl_list2
## 969 ctrl_list2
## 970 ctrl_list2
## 971 ctrl_list2
## 972 ctrl_list2
## 973 ctrl_list2
## 974 ctrl_list2
## 975 ctrl_list2
## 976 ctrl_list2
## 977 ctrl_list2
## 978 ctrl_list2
## 979 ctrl_list2
## 980 ctrl_list2
## 981 ctrl_list2
## 982 ctrl_list2
## 983 ctrl_list2
## 984 ctrl_list2
## 985 ctrl_list2
## 986 ctrl_list2
## 987 ctrl_list2
## 988 ctrl_list2
## 989 ctrl_list2
## 990 ctrl_list2
## 991 ctrl_list2
## 992 ctrl_list2
## 993 ctrl_list2
## 994 ctrl_list2
## 995 ctrl_list2
## 996 ctrl_list2
## 997 ctrl_list2
## 998 ctrl_list2
## 999 ctrl_list2
## 1000 ctrl_list25.1.4 get_params
The various parameters used during the initialization of the metagene object,
the production of the table and the production of the plot are saved and can
be accessed with the get_params function:
mg <- get_demo_metagene()
mg$get_params()## $padding_size
## [1] 0
##
## $verbose
## [1] FALSE
##
## $bam_files
## align1_rep1
## "/tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/align1_rep1.bam"
## align1_rep2
## "/tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/align1_rep2.bam"
## align2_rep1
## "/tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/align2_rep1.bam"
## align2_rep2
## "/tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/align2_rep2.bam"
## ctrl
## "/tmp/RtmpgWfRmY/Rinst78c362b7aa685/metagene/extdata/ctrl.bam"
##
## $force_seqlevels
## [1] FALSE
##
## $flip_regions
## [1] FALSE
##
## $assay
## [1] "chipseq"
##
## $df_needs_update
## [1] TRUE
##
## $df_arguments
## [1] ""
##
## $table_needs_update
## [1] TRUE5.1.5 get_design
To get the design that was used to produce the last version of the table,
you can use the get_design function:
mg$produce_table(design = get_demo_design())## Warning in normalizePath(bam_files): path[1]="align1_rep1.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align1_rep2.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep1.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep2.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="ctrl.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align1_rep1.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align1_rep2.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align1_rep1.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align1_rep2.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep1.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep2.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep1.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep2.bam": No such file or
## directory## produce data table : ChIP-Seq## Alternatively, it is also possible to add a design without producing the
## table:
mg$add_design(get_demo_design())
mg$get_design()## Samples align1 align2
## 1 align1_rep1.bam 1 0
## 2 align1_rep2.bam 1 0
## 3 align2_rep1.bam 0 1
## 4 align2_rep2.bam 0 1
## 5 ctrl.bam 2 25.1.6 get_bam_count
To get the number of aligned read in a BAM file, you can use the
get_bam_count
function:
mg$get_bam_count(bam_files[1])## [1] 46355.1.7 get_regions
To get all the regions, you can use the get_regions function:
mg$get_regions()## GRangesList object of length 2:
## $list1
## GRanges object with 50 ranges and 0 metadata columns:
## seqnames ranges strand
## <Rle> <IRanges> <Rle>
## [1] chr1 16103663-16105662 *
## [2] chr1 23921318-23923317 *
## [3] chr1 34848977-34850976 *
## [4] chr1 36368182-36370181 *
## [5] chr1 36690488-36692487 *
## ... ... ... ...
## [46] chr1 172081530-172083529 *
## [47] chr1 172081796-172083795 *
## [48] chr1 172147016-172149015 *
## [49] chr1 172205805-172207804 *
## [50] chr1 172260642-172262641 *
## -------
## seqinfo: 1 sequence from an unspecified genome; no seqlengths
##
## $list2
## GRanges object with 50 ranges and 0 metadata columns:
## seqnames ranges strand
## <Rle> <IRanges> <Rle>
## [1] chr1 3670499-3672498 *
## [2] chr1 5916399-5918398 *
## [3] chr1 9699210-9701209 *
## [4] chr1 9907639-9909638 *
## [5] chr1 10718946-10720945 *
## ... ... ... ...
## [46] chr1 81075951-81077950 *
## [47] chr1 85108854-85110853 *
## [48] chr1 85960056-85962055 *
## [49] chr1 86110971-86112970 *
## [50] chr1 87155522-87157521 *
## -------
## seqinfo: 1 sequence from an unspecified genome; no seqlengthsIt is also possible to extract a subset of the regions with the get_regions
function:
mg$get_regions(region_names = c(regions[1]))## GRangesList object of length 1:
## $list1
## GRanges object with 50 ranges and 0 metadata columns:
## seqnames ranges strand
## <Rle> <IRanges> <Rle>
## [1] chr1 16103663-16105662 *
## [2] chr1 23921318-23923317 *
## [3] chr1 34848977-34850976 *
## [4] chr1 36368182-36370181 *
## [5] chr1 36690488-36692487 *
## ... ... ... ...
## [46] chr1 172081530-172083529 *
## [47] chr1 172081796-172083795 *
## [48] chr1 172147016-172149015 *
## [49] chr1 172205805-172207804 *
## [50] chr1 172260642-172262641 *
## -------
## seqinfo: 1 sequence from an unspecified genome; no seqlengths5.1.8 get_raw_coverages
To get the coverages produced during the initialization of the metagene
object, you can use the get_raw_coverages function. Please note that to save
space, metagene will only extract the coverages in the regions provided.
coverages <- mg$get_raw_coverages()
coverages[[1]]## RleList of length 22
## $chr10
## integer-Rle of length 130694993 with 1 run
## Lengths: 130694993
## Values : 0
##
## $chr11
## integer-Rle of length 122082543 with 1 run
## Lengths: 122082543
## Values : 0
##
## $chr12
## integer-Rle of length 120129022 with 1 run
## Lengths: 120129022
## Values : 0
##
## $chr13
## integer-Rle of length 120421639 with 1 run
## Lengths: 120421639
## Values : 0
##
## $chr14
## integer-Rle of length 124902244 with 1 run
## Lengths: 124902244
## Values : 0
##
## ...
## <17 more elements>length(coverages)## [1] 5It is also possible to extract a subset of all the coverages by providing the filenames:
coverages <- mg$get_raw_coverages(filenames = bam_files[1:2])
length(coverages)## [1] 25.1.9 get_normalized_coverages
The get_normalized_coverages function works exactly like the
get_raw_coverages function except that it returns the coverages in read per
million aligned (RPM).
5.2 Chaining functions
Every function of metagene (except for the getters) invisibly return a pointer to itself. This means that the functions can be chained:
rg <- get_demo_regions()
bam <- get_demo_bam_files()
d <- get_demo_design()
title <- "Show chain"
mg <- metagene$new(rg, bam)$produce_table(design = d)$plot(title = title)## Warning in normalizePath(bam_files): path[1]="align1_rep1.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align1_rep2.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep1.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep2.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="ctrl.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align1_rep1.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align1_rep2.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align1_rep1.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align1_rep2.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep1.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep2.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep1.bam": No such file or
## directory## Warning in normalizePath(bam_files): path[1]="align2_rep2.bam": No such file or
## directory## produce data table : ChIP-Seq## produce data frame : ChIP-Seq## Plot : ChIP-Seq
5.3 Copying a metagene object
To copy a metagene object, you have to use the clone function:
mg_copy <- mg$clone()6 Managing large datasets
While metagene try to reduce it’s memory usage, it’s possible to run into
memory limits when working with multiple large datasets (especially when there
is a lot of regions with a large width).
One way to avoid this is to analyse each dataset seperately and then merge just before producing the metagene plot:
mg1 <- metagene$new(bam_files = bam_files, regions = regions[1])
mg1$produce_data_frame()
## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory
## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory
## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory
## produce data table : ChIP-Seq
## produce data frame : ChIP-Seq
mg2 <- metagene$new(bam_files = bam_files, regions = regions[2])
mg2$produce_data_frame()
## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory
## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align1_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align1_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep1": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="align2_rep2": No such file or
## directory
## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory
## Warning in normalizePath(bam_files): path[1]="ctrl": No such file or directory
## produce data table : ChIP-Seq
## produce data frame : ChIP-SeqThen you can extract the data.frames and combine them with rbind:
df1 <- mg1$get_data_frame()
df2 <- mg2$get_data_frame()
df <- rbind(df1, df2)Finally, you can use the plot_metagene function to produce the metagene plot:
p <- plot_metagene(df)## Plot : ChIP-Seqp + ggplot2::ggtitle("Managing large datasets")
7 Comparing profiles with permutations
It is possible to compare two metagene profiles using the permutation_test
function provided with the metagene package. Please note that the permutation
tests functionality is still in development and is expected to change in future
releases.
The first step is to decide which profiles we want to compare and extract the corresponding tables :
tab <- mg$get_table()
tab0 <- tab[which(tab$region == "list1"),]
tab1 <- tab0[which(tab0$design == "align1"),]
tab2 <- tab0[which(tab0$design == "align2"),]Then we defined to function to use to compare the two profiles. For this, a
companion package of metagene named similaRpeak provides
multiple metrics.
For this example, we will prepare a function to calculate the RATIO_NORMALIZED_INTERSECT between two profiles:
library(similaRpeak)
perm_fun <- function(profile1, profile2) {
sim <- similarity(profile1, profile2)
sim[["metrics"]][["RATIO_NORMALIZED_INTERSECT"]]
}We then compare our two profiles using this metric:
ratio_normalized_intersect <-
perm_fun(tab1[, .(moy=mean(value)), by=bin]$moy,
tab2[, .(moy=mean(value)), by=bin]$moy)
ratio_normalized_intersect## [1] 0.7387951To check if this value is significant, we can permute the two tables that were used to produce the profile and calculate their RATIO_NORMALIZED_INTERSECT:
permutation_results <- permutation_test(tab1, tab2, sample_size = 50,
sample_count = 1000, FUN = perm_fun)Finally, we check how often the calculated value is greater than the results of the permutations:
sum(ratio_normalized_intersect >= permutation_results) /
length(permutation_results)## [1] 0.026