diff --git a/articles/custom_engine.html b/articles/custom_engine.html index 6e19cc3..25f0bbd 100644 --- a/articles/custom_engine.html +++ b/articles/custom_engine.html @@ -159,9 +159,9 @@

Customizing the mutation engine#> [█---------------------------------------] 0% [00m:00s] Saving context index [████████████████████████████████████████] 100% [00m:00s] Context index saved #> done #> Building repeated sequence index... -#> [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█---------------------------------------] 0% [00m:03s] Processing chr. 22 [█---------------------------------------] 0% [00m:04s] Processing chr. 22 [█---------------------------------------] 0% [00m:06s] Processing chr. 22 [█---------------------------------------] 0% [00m:07s] Processing chr. 22 [█---------------------------------------] 0% [00m:08s] Processing chr. 22 [█---------------------------------------] 0% [00m:12s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:12s] RS index built -#> [█---------------------------------------] 0% [00m:00s] Saving RS index [███████████████████████-----------------] 56% [00m:01s] Saving RS indexdone -#> [████████████████████████████████████████] 100% [00m:01s] RS index saved +#> [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█---------------------------------------] 0% [00m:03s] Processing chr. 22 [█---------------------------------------] 0% [00m:04s] Processing chr. 22 [█---------------------------------------] 0% [00m:06s] Processing chr. 22 [█---------------------------------------] 0% [00m:07s] Processing chr. 22 [█---------------------------------------] 0% [00m:08s] Processing chr. 22 [█---------------------------------------] 0% [00m:11s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:12s] RS index built +#> [█---------------------------------------] 0% [00m:00s] Saving RS index [█---------------------------------------] 0% [00m:00s] Saving RS index [███████████████████████████-------------] 65% [00m:02s] Saving RS indexdone +#> [████████████████████████████████████████] 100% [00m:02s] RS index saved #> [█---------------------------------------] 0% [00m:00s] Loading germline [████████████████████████████████████████] 100% [00m:00s] Germline loaded #> [█---------------------------------------] 0% [00m:00s] Saving germline [████████████████████████████████████████] 100% [00m:00s] Germline saved @@ -207,13 +207,13 @@

The Context Index#> Downloading germline mutations... #> Germline mutations downloaded #> Building context index... -#> [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█████████████████-----------------------] 40% [00m:01s] Processing chr. 22 [█████████████████████████████████-------] 81% [00m:02s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:02s] Context index built +#> [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█████████████████-----------------------] 40% [00m:01s] Processing chr. 22 [█████████████████████████████████-------] 81% [00m:02s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:03s] Context index built #> [█---------------------------------------] 0% [00m:00s] Saving context index [████████████████████████████████████████] 100% [00m:00s] Context index saved #> done #> Building repeated sequence index... -#> [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█---------------------------------------] 0% [00m:03s] Processing chr. 22 [█---------------------------------------] 0% [00m:04s] Processing chr. 22 [█---------------------------------------] 0% [00m:06s] Processing chr. 22 [█---------------------------------------] 0% [00m:07s] Processing chr. 22 [█---------------------------------------] 0% [00m:08s] Processing chr. 22 [█---------------------------------------] 0% [00m:12s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:12s] RS index built -#> [█---------------------------------------] 0% [00m:00s] Saving RS index [████████████████------------------------] 38% [00m:01s] Saving RS indexdone -#> [████████████████████████████████████████] 100% [00m:01s] RS index saved +#> [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█---------------------------------------] 0% [00m:03s] Processing chr. 22 [█---------------------------------------] 0% [00m:04s] Processing chr. 22 [█---------------------------------------] 0% [00m:05s] Processing chr. 22 [█---------------------------------------] 0% [00m:07s] Processing chr. 22 [█---------------------------------------] 0% [00m:08s] Processing chr. 22 [█---------------------------------------] 0% [00m:12s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:12s] RS index built +#> [█---------------------------------------] 0% [00m:00s] Saving RS index [█---------------------------------------] 0% [00m:00s] Saving RS index [██████████████████████████--------------] 64% [00m:02s] Saving RS indexdone +#> [████████████████████████████████████████] 100% [00m:02s] RS index saved #> [█---------------------------------------] 0% [00m:00s] Loading germline [████████████████████████████████████████] 100% [00m:00s] Germline loaded #> [█---------------------------------------] 0% [00m:00s] Saving germline [████████████████████████████████████████] 100% [00m:00s] Germline saved @@ -226,10 +226,10 @@

The Context Index# and setting context_sampling to 50 m_engine <- MutationEngine(setup_code = "demo", context_sampling = 50) #> Building context index... -#> [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█████████████████-----------------------] 40% [00m:01s] Processing chr. 22 [█████████████████████████████████-------] 81% [00m:02s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:02s] Context index built +#> [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█████████████████-----------------------] 40% [00m:01s] Processing chr. 22 [██████████████████████████████----------] 73% [00m:02s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:02s] Context index built #> [█---------------------------------------] 0% [00m:00s] Saving context index [████████████████████████████████████████] 100% [00m:00s] Context index saved #> done -#> [█---------------------------------------] 0% [00m:00s] Loading RS index [███████████████-------------------------] 36% [00m:01s] Loading RS index [██████████████████████████████----------] 73% [00m:02s] Loading RS index [████████████████████████████████████████] 100% [00m:02s] RS index loaded +#> [█---------------------------------------] 0% [00m:00s] Loading RS index [████████████████------------------------] 39% [00m:01s] Loading RS index [███████████████████████████████---------] 76% [00m:02s] Loading RS index [████████████████████████████████████████] 100% [00m:02s] RS index loaded #> [█---------------------------------------] 0% [00m:00s] Loading germline [████████████████████████████████████████] 100% [00m:00s] Germline loaded

Above call uses the already downloaded data in the directory "demo" to build the context and produces the file @@ -272,9 +272,9 @@

The Repeated Sequence Index= 10000000) #> [█---------------------------------------] 0% [00m:00s] Loading context index [████████████████████████████████████████] 100% [00m:00s] Context index loaded #> Building repeated sequence index... -#> [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█---------------------------------------] 0% [00m:03s] Processing chr. 22 [█---------------------------------------] 0% [00m:04s] Processing chr. 22 [█---------------------------------------] 0% [00m:06s] Processing chr. 22 [█---------------------------------------] 0% [00m:07s] Processing chr. 22 [█---------------------------------------] 0% [00m:08s] Processing chr. 22 [█---------------------------------------] 0% [00m:10s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:10s] RS index built -#> [█---------------------------------------] 0% [00m:00s] Saving RS index [█---------------------------------------] 1% [00m:01s] Saving RS index [███████████-----------------------------] 25% [00m:02s] Saving RS index [█████████████████████-------------------] 52% [00m:03s] Saving RS index [███████████████████████████████---------] 77% [00m:04s] Saving RS indexdone -#> [████████████████████████████████████████] 100% [00m:05s] RS index saved +#> [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█---------------------------------------] 0% [00m:03s] Processing chr. 22 [█---------------------------------------] 0% [00m:04s] Processing chr. 22 [█---------------------------------------] 0% [00m:05s] Processing chr. 22 [█---------------------------------------] 0% [00m:06s] Processing chr. 22 [█---------------------------------------] 0% [00m:08s] Processing chr. 22 [█---------------------------------------] 0% [00m:09s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:10s] RS index built +#> [█---------------------------------------] 0% [00m:00s] Saving RS index [█---------------------------------------] 0% [00m:00s] Saving RS index [█---------------------------------------] 0% [00m:01s] Saving RS index [█---------------------------------------] 0% [00m:02s] Saving RS index [█---------------------------------------] 0% [00m:03s] Saving RS index [█---------------------------------------] 0% [00m:04s] Saving RS index [█████████-------------------------------] 20% [00m:06s] Saving RS index [█████████████████-----------------------] 40% [00m:07s] Saving RS index [██████████████████████████--------------] 63% [00m:08s] Saving RS index [███████████████████████████████████-----] 86% [00m:09s] Saving RS indexdone +#> [████████████████████████████████████████] 100% [00m:09s] RS index saved #> [█---------------------------------------] 0% [00m:00s] Loading germline [████████████████████████████████████████] 100% [00m:00s] Germline loaded

Above call uses the already downloaded data in the directory "demo" to build the repeated sequence index and produces diff --git a/articles/growth_models.html b/articles/growth_models.html index 270cb56..95dc7b9 100644 --- a/articles/growth_models.html +++ b/articles/growth_models.html @@ -250,7 +250,7 @@

“Border” Growth # Let the simulation evolve until "A+" consists of 1300 cells sim$run_up_to_size("A+", 1300) -#> [████████████████████████████████████████] 100% [00m:00s] Saving snapshot +#> [████████████████████████████████████----] 89% [00m:00s] Cells: 30477 [████████████████████████████████████████] 100% [00m:00s] Saving snapshot bbox_width <- 15 diff --git a/articles/infinite_site_model.html b/articles/infinite_site_model.html index fc26fb7..0bc05fe 100644 --- a/articles/infinite_site_model.html +++ b/articles/infinite_site_model.html @@ -168,7 +168,7 @@ # place the mutations on the same samples forest used above phylo_forest2 <- m_engine$place_mutations(samples_forest, 1000, 500) -#> [█---------------------------------------] 0% [00m:00s] Placing mutations [███████████████████---------------------] 45% [00m:01s] Placing mutations [██████████████████████████████████------] 84% [00m:02s] Placing mutations [████████████████████████████████████----] 89% [00m:03s] Placing mutations [████████████████████████████████████████] 100% [00m:03s] Mutations placed +#> [█---------------------------------------] 0% [00m:00s] Placing mutations [███████████████████---------------------] 46% [00m:01s] Placing mutations [██████████████████████████████████------] 84% [00m:02s] Placing mutations [████████████████████████████████████----] 89% [00m:03s] Placing mutations [████████████████████████████████████████] 100% [00m:03s] Mutations placed # test whether the infinite sites conditions hold in the new forest test_infinite_sites_model(phylo_forest2) diff --git a/articles/mutations.html b/articles/mutations.html index 15a2d3b..e6b19e5 100644 --- a/articles/mutations.html +++ b/articles/mutations.html @@ -137,13 +137,13 @@

Setting Up Mutation Engine#> Downloading germline mutations... #> Germline mutations downloaded #> Building context index... -#> [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█████████████████-----------------------] 40% [00m:01s] Processing chr. 22 [█████████████████████████████████-------] 81% [00m:02s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:03s] Context index built +#> [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█████████████████-----------------------] 40% [00m:01s] Processing chr. 22 [█████████████████████████████████-------] 81% [00m:02s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:02s] Context index built #> [█---------------------------------------] 0% [00m:00s] Saving context index [████████████████████████████████████████] 100% [00m:00s] Context index saved #> done #> Building repeated sequence index... #> [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█---------------------------------------] 0% [00m:03s] Processing chr. 22 [█---------------------------------------] 0% [00m:04s] Processing chr. 22 [█---------------------------------------] 0% [00m:06s] Processing chr. 22 [█---------------------------------------] 0% [00m:07s] Processing chr. 22 [█---------------------------------------] 0% [00m:08s] Processing chr. 22 [█---------------------------------------] 0% [00m:12s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:12s] RS index built -#> [█---------------------------------------] 0% [00m:00s] Saving RS index [█████████████████████████---------------] 60% [00m:01s] Saving RS indexdone -#> [████████████████████████████████████████] 100% [00m:01s] RS index saved +#> [█---------------------------------------] 0% [00m:00s] Saving RS index [█---------------------------------------] 0% [00m:00s] Saving RS index [███████████████████████████-------------] 66% [00m:02s] Saving RS indexdone +#> [████████████████████████████████████████] 100% [00m:02s] RS index saved #> [█---------------------------------------] 0% [00m:00s] Loading germline [████████████████████████████████████████] 100% [00m:00s] Germline loaded #> [█---------------------------------------] 0% [00m:00s] Saving germline [████████████████████████████████████████] 100% [00m:00s] Germline saved @@ -575,7 +575,7 @@

Building Phylogenetic Forests# place mutations on the sample forest assuming 1000 pre-neoplastic SNVs and # 500 indels phylo_forest <- m_engine$place_mutations(samples_forest, 1000, 500) -#> [█---------------------------------------] 0% [00m:00s] Placing mutations [█████████-------------------------------] 20% [00m:01s] Placing mutations [█████████████████-----------------------] 40% [00m:02s] Placing mutations [█████████████████████████---------------] 60% [00m:03s] Placing mutations [█████████████████████████████████-------] 80% [00m:04s] Placing mutations [█████████████████████████████████████---] 92% [00m:05s] Placing mutations [████████████████████████████████████████] 100% [00m:05s] Mutations placed +#> [█---------------------------------------] 0% [00m:00s] Placing mutations [████████--------------------------------] 19% [00m:01s] Placing mutations [████████████████------------------------] 39% [00m:02s] Placing mutations [████████████████████████----------------] 58% [00m:03s] Placing mutations [███████████████████████████████---------] 75% [00m:04s] Placing mutations [████████████████████████████████████----] 89% [00m:05s] Placing mutations [████████████████████████████████████████] 100% [00m:05s] Mutations placed phylo_forest #> PhylogeneticForest @@ -662,11 +662,11 @@

Building Phylogenetic Forests
 phylo_forest$get_samples_info()
-#>    name xmin ymin xmax ymax tumour_cells tumour_cells_in_bbox     time
-#> 1 S_1_1  480  480  530  530         2588                 2588 365.9279
-#> 2 S_1_2  500  500  550  550         1632                 1632 365.9279
-#> 3 S_2_1  370  498  394  522          572                  572 421.2143
-#> 4 S_2_2  370  373  394  397          604                  604 421.2143
+#>    name id xmin ymin xmax ymax tumour_cells tumour_cells_in_bbox     time
+#> 1 S_1_1  0  480  480  530  530         2588                 2588 365.9279
+#> 2 S_1_2  1  500  500  550  550         1632                 1632 365.9279
+#> 3 S_2_1  2  370  498  394  522          572                  572 421.2143
+#> 4 S_2_2  3  370  373  394  397          604                  604 421.2143
 #>   DNA_quantity equivalent_normal_cells
 #> 1 530855879155                5173.573
 #> 2 333071264048                3246.020
@@ -743,12 +743,12 @@ 

Building Phylogenetic Forests# print the cell allelic fragmentation phylo_forest$get_cell_allelic_fragmentation() %>% head() #> cell_id chr begin end major minor -#> 1 2790 22 1 5009998 2 2 -#> 2 2790 22 5009999 5009999 2 0 -#> 3 2790 22 5010000 5209998 2 0 -#> 4 2790 22 5209999 5209999 2 0 -#> 5 2790 22 5210000 10303468 2 2 -#> 6 2790 22 10303469 10303469 2 2

+#> 1 2790 22 1 5009999 2 2 +#> 2 2790 22 5010000 5209999 2 0 +#> 3 2790 22 5210000 10303469 2 2 +#> 4 2790 22 10303470 10503469 4 2 +#> 5 2790 22 10503470 51304566 2 2 +#> 6 5000 22 1 5009999 2 2

The details about the SNV and indel signatures adopted during the evolution are available in the mutation engine and they can be retrieved by using the methods MutationEngine$get_SNV_signatures() @@ -790,11 +790,11 @@

Storing Phylogenetic Forests
 # save the phylogenetic forest in the file "phylo_forest.sff"
 phylo_forest$save("phylo_forest.sff")
-#>  [█---------------------------------------] 0% [00m:00s] Saving forest [█---------------------------------------] 1% [00m:08s] Saving forest [███████---------------------------------] 15% [00m:09s] Saving forest [███████████-----------------------------] 26% [00m:10s] Saving forest [███████████████-------------------------] 37% [00m:11s] Saving forest [████████████████████--------------------] 48% [00m:12s] Saving forest [████████████████████████----------------] 58% [00m:13s] Saving forest [████████████████████████████------------] 69% [00m:15s] Saving forest [████████████████████████████████--------] 79% [00m:16s] Saving forest [███████████████████████████████████-----] 86% [00m:17s] Saving forest [███████████████████████████████████████-] 95% [00m:18s] Saving forest [████████████████████████████████████████] 100% [00m:18s] Forest saved
+#>  [█---------------------------------------] 0% [00m:00s] Saving forest [█---------------------------------------] 0% [00m:00s] Saving forest [█---------------------------------------] 0% [00m:01s] Saving forest [█---------------------------------------] 0% [00m:02s] Saving forest [█---------------------------------------] 0% [00m:03s] Saving forest [█---------------------------------------] 0% [00m:04s] Saving forest [█---------------------------------------] 0% [00m:05s] Saving forest [█---------------------------------------] 0% [00m:06s] Saving forest [█---------------------------------------] 0% [00m:07s] Saving forest [█---------------------------------------] 0% [00m:08s] Saving forest [█---------------------------------------] 0% [00m:09s] Saving forest [█████-----------------------------------] 12% [00m:11s] Saving forest [█████████-------------------------------] 22% [00m:12s] Saving forest [█████████████---------------------------] 32% [00m:13s] Saving forest [█████████████████-----------------------] 42% [00m:14s] Saving forest [█████████████████████-------------------] 52% [00m:15s] Saving forest [█████████████████████████---------------] 62% [00m:16s] Saving forest [█████████████████████████████-----------] 72% [00m:17s] Saving forest [█████████████████████████████████-------] 82% [00m:18s] Saving forest [█████████████████████████████████████---] 92% [00m:19s] Saving forest [████████████████████████████████████████] 100% [00m:20s] Forest saved
 
 # loading the saved forest
 loaded_phylo_forest <- load_phylogenetic_forest("phylo_forest.sff")
-#>  [█---------------------------------------] 0% [00m:00s] Loading forest [█████-----------------------------------] 10% [00m:01s] Loading forest [███████---------------------------------] 17% [00m:02s] Loading forest [██████████------------------------------] 24% [00m:03s] Loading forest [█████████████---------------------------] 31% [00m:04s] Loading forest [████████████████------------------------] 38% [00m:05s] Loading forest [███████████████████---------------------] 45% [00m:06s] Loading forest [█████████████████████-------------------] 52% [00m:07s] Loading forest [████████████████████████----------------] 58% [00m:08s] Loading forest [███████████████████████████-------------] 66% [00m:09s] Loading forest [██████████████████████████████----------] 73% [00m:11s] Loading forest [█████████████████████████████████-------] 80% [00m:12s] Loading forest [███████████████████████████████████-----] 87% [00m:13s] Loading forest [██████████████████████████████████████--] 94% [00m:14s] Loading forest [████████████████████████████████████████] 100% [00m:15s] Forest loaded
+#>  [█---------------------------------------] 0% [00m:00s] Loading forest [█████-----------------------------------] 10% [00m:01s] Loading forest [███████---------------------------------] 17% [00m:02s] Loading forest [██████████------------------------------] 24% [00m:03s] Loading forest [█████████████---------------------------] 31% [00m:04s] Loading forest [████████████████------------------------] 38% [00m:05s] Loading forest [███████████████████---------------------] 45% [00m:06s] Loading forest [█████████████████████-------------------] 52% [00m:07s] Loading forest [████████████████████████----------------] 58% [00m:08s] Loading forest [███████████████████████████-------------] 65% [00m:09s] Loading forest [█████████████████████████████-----------] 72% [00m:11s] Loading forest [████████████████████████████████--------] 79% [00m:12s] Loading forest [███████████████████████████████████-----] 86% [00m:13s] Loading forest [██████████████████████████████████████--] 93% [00m:14s] Loading forest [████████████████████████████████████████] 100% [00m:15s] Forest loaded
 
 loaded_phylo_forest
 #> PhylogeneticForest
diff --git a/articles/sample_partition.html b/articles/sample_partition.html
index 608ad87..5635593 100644
--- a/articles/sample_partition.html
+++ b/articles/sample_partition.html
@@ -127,7 +127,7 @@
 
 seq_results <- simulate_seq(phylo_forest, coverage = 0.5,
                             cell_labelling = epi_labelling)
-#>  [█---------------------------------------] 0% [00m:00s] Found 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:01s] Processing chr. 22 [█---------------------------------------] 0% [00m:06s] Processing chr. 22 [█████████████---------------------------] 32% [00m:07s] Processing chr. 22 [███████████████████---------------------] 47% [00m:08s] Processing chr. 22 [████████████████████--------------------] 48% [00m:09s] Processing chr. 22 [████████████████████████████------------] 69% [00m:10s] Processing chr. 22 [███████████████████████████████---------] 77% [00m:11s] Processing chr. 22 [████████████████████████████████--------] 78% [00m:12s] Processing chr. 22 [███████████████████████████████████-----] 85% [00m:13s] Processing chr. 22 [████████████████████████████████████----] 88% [00m:14s] Processing chr. 22 [████████████████████████████████████----] 88% [00m:15s] Processing chr. 22 [████████████████████████████████████----] 89% [00m:16s] Processing chr. 22 [██████████████████████████████████████--] 93% [00m:17s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:18s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:19s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:23s] Reads simulated
+#>  [█---------------------------------------] 0% [00m:00s] Found 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:01s] Reading 22 [█---------------------------------------] 0% [00m:01s] Processing chr. 22 [█---------------------------------------] 0% [00m:05s] Processing chr. 22 [█████████████---------------------------] 32% [00m:06s] Processing chr. 22 [███████████████████---------------------] 47% [00m:07s] Processing chr. 22 [███████████████████---------------------] 47% [00m:08s] Processing chr. 22 [██████████████████████████--------------] 63% [00m:09s] Processing chr. 22 [███████████████████████████████---------] 77% [00m:10s] Processing chr. 22 [████████████████████████████████--------] 78% [00m:11s] Processing chr. 22 [█████████████████████████████████-------] 82% [00m:12s] Processing chr. 22 [████████████████████████████████████----] 88% [00m:13s] Processing chr. 22 [████████████████████████████████████----] 88% [00m:14s] Processing chr. 22 [████████████████████████████████████----] 89% [00m:15s] Processing chr. 22 [████████████████████████████████████----] 89% [00m:16s] Processing chr. 22 [███████████████████████████████████████-] 95% [00m:17s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:18s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:19s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:22s] Reads simulated
 
 library(dplyr)
 #> 
@@ -199,7 +199,7 @@
 
 seq_results <- simulate_seq(phylo_forest, coverage = 0.5,
                             cell_labelling = mutant_labelling)
-#>  [█---------------------------------------] 0% [00m:00s] Found 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:01s] Processing chr. 22 [█---------------------------------------] 0% [00m:05s] Processing chr. 22 [█████████-------------------------------] 21% [00m:06s] Processing chr. 22 [███████████████████---------------------] 47% [00m:07s] Processing chr. 22 [█████████████████████████---------------] 61% [00m:08s] Processing chr. 22 [████████████████████████████████--------] 78% [00m:09s] Processing chr. 22 [███████████████████████████████████-----] 85% [00m:10s] Processing chr. 22 [█████████████████████████████████████---] 91% [00m:11s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:12s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:13s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:16s] Reads simulated
+#>  [█---------------------------------------] 0% [00m:00s] Found 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:01s] Reading 22 [█---------------------------------------] 0% [00m:01s] Processing chr. 22 [█---------------------------------------] 0% [00m:05s] Processing chr. 22 [█████████████---------------------------] 30% [00m:06s] Processing chr. 22 [██████████████████████------------------] 54% [00m:07s] Processing chr. 22 [██████████████████████████████----------] 74% [00m:08s] Processing chr. 22 [██████████████████████████████████------] 83% [00m:09s] Processing chr. 22 [████████████████████████████████████----] 88% [00m:10s] Processing chr. 22 [██████████████████████████████████████--] 94% [00m:11s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:12s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:13s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:16s] Reads simulated
 
 seq_results %>% head()
 #>   chr  chr_pos ref   alt causes  classes S_1_1_A.occurrences S_1_1_A.coverage
@@ -254,7 +254,7 @@
 
 seq_results <- simulate_seq(phylo_forest, coverage = 0.5,
                             cell_labelling = birth_time_labelling)
-#>  [█---------------------------------------] 0% [00m:00s] Found 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:01s] Reading 22 [█---------------------------------------] 0% [00m:01s] Processing chr. 22 [█---------------------------------------] 0% [00m:05s] Processing chr. 22 [███-------------------------------------] 6% [00m:06s] Processing chr. 22 [█████████-------------------------------] 21% [00m:07s] Processing chr. 22 [████████████████████--------------------] 48% [00m:08s] Processing chr. 22 [█████████████████████-------------------] 51% [00m:09s] Processing chr. 22 [█████████████████████████---------------] 61% [00m:10s] Processing chr. 22 [████████████████████████████████--------] 78% [00m:11s] Processing chr. 22 [████████████████████████████████--------] 78% [00m:12s] Processing chr. 22 [█████████████████████████████████-------] 81% [00m:13s] Processing chr. 22 [████████████████████████████████████----] 88% [00m:14s] Processing chr. 22 [████████████████████████████████████----] 88% [00m:15s] Processing chr. 22 [█████████████████████████████████████---] 92% [00m:16s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:17s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:18s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:19s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:24s] Reads simulated
+#>  [█---------------------------------------] 0% [00m:00s] Found 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:01s] Processing chr. 22 [█---------------------------------------] 0% [00m:05s] Processing chr. 22 [███-------------------------------------] 6% [00m:06s] Processing chr. 22 [█████████-------------------------------] 20% [00m:07s] Processing chr. 22 [███████████████████---------------------] 47% [00m:08s] Processing chr. 22 [█████████████████████-------------------] 51% [00m:09s] Processing chr. 22 [███████████████████████-----------------] 57% [00m:10s] Processing chr. 22 [██████████████████████████████----------] 74% [00m:11s] Processing chr. 22 [████████████████████████████████--------] 78% [00m:12s] Processing chr. 22 [████████████████████████████████--------] 78% [00m:13s] Processing chr. 22 [███████████████████████████████████-----] 85% [00m:14s] Processing chr. 22 [████████████████████████████████████----] 88% [00m:15s] Processing chr. 22 [████████████████████████████████████----] 88% [00m:16s] Processing chr. 22 [██████████████████████████████████████--] 94% [00m:17s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:18s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:19s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:20s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:21s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:26s] Reads simulated
 
 seq_results %>% head()
 #>   chr  chr_pos ref   alt causes  classes S_1_1_MIDDLE_AGED.occurrences
@@ -362,7 +362,7 @@
 
 seq_results <- simulate_seq(phylo_forest, coverage = 0.5,
                             cell_labelling = mutations_labelling)
-#>  [█---------------------------------------] 0% [00m:00s] Found 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:01s] Reading 22 [█---------------------------------------] 0% [00m:01s] Processing chr. 22 [█---------------------------------------] 0% [00m:05s] Processing chr. 22 [███████████-----------------------------] 27% [00m:06s] Processing chr. 22 [███████████████████---------------------] 45% [00m:07s] Processing chr. 22 [███████████████████---------------------] 47% [00m:08s] Processing chr. 22 [███████████████████████-----------------] 55% [00m:09s] Processing chr. 22 [███████████████████████████████---------] 75% [00m:10s] Processing chr. 22 [██████████████████████████████████------] 83% [00m:11s] Processing chr. 22 [█████████████████████████████████████---] 90% [00m:12s] Processing chr. 22 [████████████████████████████████████████] 98% [00m:13s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:14s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:17s] Reads simulated
+#>  [█---------------------------------------] 0% [00m:00s] Found 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:01s] Reading 22 [█---------------------------------------] 0% [00m:01s] Processing chr. 22 [█---------------------------------------] 0% [00m:05s] Processing chr. 22 [█████████████---------------------------] 31% [00m:06s] Processing chr. 22 [███████████████████---------------------] 45% [00m:07s] Processing chr. 22 [███████████████████---------------------] 47% [00m:08s] Processing chr. 22 [██████████████████████████--------------] 64% [00m:09s] Processing chr. 22 [████████████████████████████████--------] 79% [00m:10s] Processing chr. 22 [███████████████████████████████████-----] 86% [00m:11s] Processing chr. 22 [██████████████████████████████████████--] 93% [00m:12s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:13s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:14s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:17s] Reads simulated
 
 seq_results %>% head()
 #>   chr  chr_pos ref   alt causes  classes S_1_1_HAS_MUTATION.occurrences
diff --git a/articles/sampling.html b/articles/sampling.html
index 064a3cc..623d0a5 100644
--- a/articles/sampling.html
+++ b/articles/sampling.html
@@ -119,7 +119,7 @@ 

Custom multi-region samplingsim$place_cell("A", 500, 500) sim$run_up_to_size("A", 60000) -#> [█████████████████████████████████-------] 82% [00m:00s] Cells: 49350 [████████████████████████████████████████] 100% [00m:00s] Saving snapshot

+#> [██████████████████████████████████------] 84% [00m:00s] Cells: 50483 [████████████████████████████████████████] 100% [00m:00s] Saving snapshot @@ -240,9 +240,9 @@

Cell division tree for sampled cel

We can also query the forest about the samples used to build it.

 forest$get_samples_info()
-#>    name xmin ymin xmax ymax tumour_cells tumour_cells_in_bbox     time
-#> 1 S_1_1  400  400  450  450         2588                 2588 397.0568
-#> 2 S_1_2  500  500  550  550         2593                 2593 397.0568
+#> name id xmin ymin xmax ymax tumour_cells tumour_cells_in_bbox time +#> 1 S_1_1 0 400 400 450 450 2588 2588 397.0568 +#> 2 S_1_2 1 500 500 550 550 2593 2593 397.0568

We can visualise the forest. This plot reports the cells and, on the y-axis, their time of birth.

+#> [██████████████████████████████████------] 84% [00m:00s] Cells: 50483 [████████████████████████████████████████] 100% [00m:00s] Saving snapshot

We include a new mutant and let it grow. This new mutant has much higher growth rates than its ancestor.

diff --git a/articles/sequencing.html b/articles/sequencing.html index 4ab667d..ce237b2 100644 --- a/articles/sequencing.html +++ b/articles/sequencing.html @@ -101,7 +101,7 @@

Sequencing Simulationlibrary(rRACES) phylo_forest <- load_phylogenetic_forest("phylo_forest.sff") -#> [█---------------------------------------] 0% [00m:00s] Loading forest [█████-----------------------------------] 10% [00m:01s] Loading forest [███████---------------------------------] 17% [00m:02s] Loading forest [██████████------------------------------] 24% [00m:03s] Loading forest [█████████████---------------------------] 31% [00m:04s] Loading forest [████████████████------------------------] 38% [00m:05s] Loading forest [███████████████████---------------------] 45% [00m:06s] Loading forest [█████████████████████-------------------] 52% [00m:07s] Loading forest [████████████████████████----------------] 58% [00m:08s] Loading forest [███████████████████████████-------------] 66% [00m:09s] Loading forest [█████████████████████████████-----------] 72% [00m:10s] Loading forest [████████████████████████████████--------] 79% [00m:11s] Loading forest [███████████████████████████████████-----] 86% [00m:13s] Loading forest [██████████████████████████████████████--] 93% [00m:14s] Loading forest [████████████████████████████████████████] 100% [00m:15s] Forest loaded +#> [█---------------------------------------] 0% [00m:00s] Loading forest [█████-----------------------------------] 10% [00m:01s] Loading forest [███████---------------------------------] 17% [00m:02s] Loading forest [██████████------------------------------] 23% [00m:03s] Loading forest [█████████████---------------------------] 30% [00m:04s] Loading forest [███████████████-------------------------] 37% [00m:05s] Loading forest [██████████████████----------------------] 44% [00m:06s] Loading forest [█████████████████████-------------------] 51% [00m:07s] Loading forest [████████████████████████----------------] 58% [00m:08s] Loading forest [███████████████████████████-------------] 65% [00m:09s] Loading forest [█████████████████████████████-----------] 72% [00m:10s] Loading forest [████████████████████████████████--------] 79% [00m:11s] Loading forest [███████████████████████████████████-----] 86% [00m:13s] Loading forest [██████████████████████████████████████--] 93% [00m:14s] Loading forest [████████████████████████████████████████] 100% [00m:15s] Forest loaded phylo_forest #> PhylogeneticForest @@ -120,7 +120,7 @@

Sequencing Simulation
 # let us simulate a 50x sequencing of the four sample
 seq_results <- simulate_seq(phylo_forest, coverage = 30)
-#>  [█---------------------------------------] 0% [00m:00s] Found 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:01s] Reading 22 [█---------------------------------------] 0% [00m:01s] Processing chr. 22 [█---------------------------------------] 0% [00m:05s] Processing chr. 22 [█---------------------------------------] 0% [00m:06s] Processing chr. 22 [█---------------------------------------] 1% [00m:07s] Processing chr. 22 [█---------------------------------------] 1% [00m:08s] Processing chr. 22 [█---------------------------------------] 2% [00m:09s] Processing chr. 22 [█---------------------------------------] 2% [00m:10s] Processing chr. 22 [██--------------------------------------] 3% [00m:11s] Processing chr. 22 [██--------------------------------------] 4% [00m:12s] Processing chr. 22 [██--------------------------------------] 4% [00m:13s] Processing chr. 22 [███-------------------------------------] 5% [00m:14s] Processing chr. 22 [███-------------------------------------] 5% [00m:15s] Processing chr. 22 [███-------------------------------------] 6% [00m:16s] Processing chr. 22 [███-------------------------------------] 7% [00m:17s] Processing chr. 22 [███-------------------------------------] 7% [00m:18s] Processing chr. 22 [████------------------------------------] 8% [00m:19s] Processing chr. 22 [████------------------------------------] 8% [00m:20s] Processing chr. 22 [████------------------------------------] 9% [00m:21s] Processing chr. 22 [█████-----------------------------------] 10% [00m:22s] Processing chr. 22 [█████-----------------------------------] 10% [00m:23s] Processing chr. 22 [█████-----------------------------------] 11% [00m:24s] Processing chr. 22 [█████-----------------------------------] 11% [00m:25s] Processing chr. 22 [█████-----------------------------------] 12% [00m:26s] Processing chr. 22 [█████-----------------------------------] 12% [00m:27s] Processing chr. 22 [██████----------------------------------] 13% [00m:28s] Processing chr. 22 [██████----------------------------------] 14% [00m:29s] Processing chr. 22 [██████----------------------------------] 14% [00m:30s] Processing chr. 22 [███████---------------------------------] 15% [00m:31s] Processing chr. 22 [███████---------------------------------] 15% [00m:32s] Processing chr. 22 [███████---------------------------------] 16% [00m:33s] Processing chr. 22 [███████---------------------------------] 17% [00m:34s] Processing chr. 22 [███████---------------------------------] 17% [00m:35s] Processing chr. 22 [████████--------------------------------] 18% [00m:36s] Processing chr. 22 [████████--------------------------------] 18% [00m:37s] Processing chr. 22 [████████--------------------------------] 19% [00m:38s] Processing chr. 22 [█████████-------------------------------] 20% [00m:39s] Processing chr. 22 [█████████-------------------------------] 20% [00m:40s] Processing chr. 22 [█████████-------------------------------] 21% [00m:41s] Processing chr. 22 [█████████-------------------------------] 21% [00m:42s] Processing chr. 22 [█████████-------------------------------] 22% [00m:43s] Processing chr. 22 [██████████------------------------------] 23% [00m:44s] Processing chr. 22 [██████████------------------------------] 23% [00m:45s] Processing chr. 22 [██████████------------------------------] 24% [00m:46s] Processing chr. 22 [██████████------------------------------] 24% [00m:47s] Processing chr. 22 [███████████-----------------------------] 25% [00m:48s] Processing chr. 22 [███████████-----------------------------] 26% [00m:49s] Processing chr. 22 [███████████-----------------------------] 26% [00m:50s] Processing chr. 22 [███████████-----------------------------] 27% [00m:51s] Processing chr. 22 [███████████-----------------------------] 27% [00m:52s] Processing chr. 22 [████████████----------------------------] 28% [00m:53s] Processing chr. 22 [████████████----------------------------] 28% [00m:54s] Processing chr. 22 [████████████----------------------------] 29% [00m:55s] Processing chr. 22 [█████████████---------------------------] 30% [00m:56s] Processing chr. 22 [█████████████---------------------------] 30% [00m:57s] Processing chr. 22 [█████████████---------------------------] 31% [00m:58s] Processing chr. 22 [█████████████---------------------------] 31% [00m:59s] Processing chr. 22 [█████████████---------------------------] 32% [01m:00s] Processing chr. 22 [██████████████--------------------------] 33% [01m:01s] Processing chr. 22 [██████████████--------------------------] 33% [01m:02s] Processing chr. 22 [██████████████--------------------------] 34% [01m:03s] Processing chr. 22 [██████████████--------------------------] 34% [01m:04s] Processing chr. 22 [███████████████-------------------------] 35% [01m:05s] Processing chr. 22 [███████████████-------------------------] 36% [01m:06s] Processing chr. 22 [███████████████-------------------------] 36% [01m:07s] Processing chr. 22 [███████████████-------------------------] 37% [01m:08s] Processing chr. 22 [███████████████-------------------------] 37% [01m:09s] Processing chr. 22 [████████████████------------------------] 38% [01m:10s] Processing chr. 22 [████████████████------------------------] 39% [01m:11s] Processing chr. 22 [████████████████------------------------] 39% [01m:12s] Processing chr. 22 [█████████████████-----------------------] 40% [01m:13s] Processing chr. 22 [█████████████████-----------------------] 40% [01m:14s] Processing chr. 22 [█████████████████-----------------------] 41% [01m:15s] Processing chr. 22 [█████████████████-----------------------] 41% [01m:16s] Processing chr. 22 [█████████████████-----------------------] 42% [01m:17s] Processing chr. 22 [██████████████████----------------------] 43% [01m:18s] Processing chr. 22 [██████████████████----------------------] 43% [01m:19s] Processing chr. 22 [██████████████████----------------------] 44% [01m:20s] Processing chr. 22 [██████████████████----------------------] 44% [01m:21s] Processing chr. 22 [███████████████████---------------------] 45% [01m:22s] Processing chr. 22 [███████████████████---------------------] 46% [01m:23s] Processing chr. 22 [███████████████████---------------------] 46% [01m:24s] Processing chr. 22 [███████████████████---------------------] 47% [01m:25s] Processing chr. 22 [███████████████████---------------------] 47% [01m:26s] Processing chr. 22 [████████████████████--------------------] 48% [01m:27s] Processing chr. 22 [████████████████████--------------------] 48% [01m:28s] Processing chr. 22 [████████████████████--------------------] 48% [01m:29s] Processing chr. 22 [████████████████████--------------------] 49% [01m:30s] Processing chr. 22 [████████████████████--------------------] 49% [01m:31s] Processing chr. 22 [█████████████████████-------------------] 50% [01m:32s] Processing chr. 22 [█████████████████████-------------------] 50% [01m:33s] Processing chr. 22 [█████████████████████-------------------] 50% [01m:34s] Processing chr. 22 [█████████████████████-------------------] 51% [01m:35s] Processing chr. 22 [█████████████████████-------------------] 51% [01m:36s] Processing chr. 22 [█████████████████████-------------------] 51% [01m:37s] Processing chr. 22 [█████████████████████-------------------] 52% [01m:38s] Processing chr. 22 [█████████████████████-------------------] 52% [01m:39s] Processing chr. 22 [██████████████████████------------------] 52% [01m:40s] Processing chr. 22 [██████████████████████------------------] 52% [01m:41s] Processing chr. 22 [██████████████████████------------------] 52% [01m:42s] Processing chr. 22 [██████████████████████------------------] 54% [01m:43s] Processing chr. 22 [██████████████████████------------------] 54% [01m:44s] Processing chr. 22 [██████████████████████------------------] 54% [01m:45s] Processing chr. 22 [███████████████████████-----------------] 55% [01m:46s] Processing chr. 22 [███████████████████████-----------------] 55% [01m:47s] Processing chr. 22 [███████████████████████-----------------] 56% [01m:48s] Processing chr. 22 [███████████████████████-----------------] 56% [01m:49s] Processing chr. 22 [███████████████████████-----------------] 56% [01m:50s] Processing chr. 22 [███████████████████████-----------------] 57% [01m:51s] Processing chr. 22 [███████████████████████-----------------] 57% [01m:52s] Processing chr. 22 [███████████████████████-----------------] 57% [01m:53s] Processing chr. 22 [████████████████████████----------------] 58% [01m:54s] Processing chr. 22 [████████████████████████----------------] 58% [01m:55s] Processing chr. 22 [████████████████████████----------------] 58% [01m:56s] Processing chr. 22 [████████████████████████----------------] 58% [01m:57s] Processing chr. 22 [████████████████████████----------------] 58% [01m:58s] Processing chr. 22 [█████████████████████████---------------] 60% [01m:59s] Processing chr. 22 [█████████████████████████---------------] 60% [02m:00s] Processing chr. 22 [█████████████████████████---------------] 60% [02m:01s] Processing chr. 22 [█████████████████████████---------------] 61% [02m:02s] Processing chr. 22 [█████████████████████████---------------] 61% [02m:03s] Processing chr. 22 [█████████████████████████---------------] 62% [02m:04s] Processing chr. 22 [█████████████████████████---------------] 62% [02m:05s] Processing chr. 22 [█████████████████████████---------------] 62% [02m:06s] Processing chr. 22 [██████████████████████████--------------] 63% [02m:07s] Processing chr. 22 [██████████████████████████--------------] 63% [02m:08s] Processing chr. 22 [██████████████████████████--------------] 63% [02m:09s] Processing chr. 22 [██████████████████████████--------------] 64% [02m:10s] Processing chr. 22 [██████████████████████████--------------] 64% [02m:11s] Processing chr. 22 [██████████████████████████--------------] 64% [02m:12s] Processing chr. 22 [███████████████████████████-------------] 65% [02m:13s] Processing chr. 22 [███████████████████████████-------------] 65% [02m:14s] Processing chr. 22 [███████████████████████████-------------] 66% [02m:15s] Processing chr. 22 [███████████████████████████-------------] 66% [02m:16s] Processing chr. 22 [███████████████████████████-------------] 66% [02m:17s] Processing chr. 22 [███████████████████████████-------------] 67% [02m:18s] Processing chr. 22 [███████████████████████████-------------] 67% [02m:19s] Processing chr. 22 [███████████████████████████-------------] 67% [02m:20s] Processing chr. 22 [████████████████████████████------------] 68% [02m:21s] Processing chr. 22 [████████████████████████████------------] 68% [02m:22s] Processing chr. 22 [████████████████████████████------------] 69% [02m:23s] Processing chr. 22 [████████████████████████████------------] 69% [02m:24s] Processing chr. 22 [████████████████████████████------------] 69% [02m:25s] Processing chr. 22 [█████████████████████████████-----------] 70% [02m:26s] Processing chr. 22 [█████████████████████████████-----------] 70% [02m:27s] Processing chr. 22 [█████████████████████████████-----------] 70% [02m:28s] Processing chr. 22 [█████████████████████████████-----------] 71% [02m:29s] Processing chr. 22 [█████████████████████████████-----------] 71% [02m:30s] Processing chr. 22 [█████████████████████████████-----------] 72% [02m:31s] Processing chr. 22 [█████████████████████████████-----------] 72% [02m:32s] Processing chr. 22 [█████████████████████████████-----------] 72% [02m:33s] Processing chr. 22 [██████████████████████████████----------] 73% [02m:34s] Processing chr. 22 [██████████████████████████████----------] 73% [02m:35s] Processing chr. 22 [██████████████████████████████----------] 73% [02m:36s] Processing chr. 22 [██████████████████████████████----------] 74% [02m:37s] Processing chr. 22 [██████████████████████████████----------] 74% [02m:38s] Processing chr. 22 [███████████████████████████████---------] 75% [02m:39s] Processing chr. 22 [███████████████████████████████---------] 75% [02m:40s] Processing chr. 22 [███████████████████████████████---------] 75% [02m:41s] Processing chr. 22 [███████████████████████████████---------] 76% [02m:42s] Processing chr. 22 [███████████████████████████████---------] 76% [02m:43s] Processing chr. 22 [███████████████████████████████---------] 76% [02m:44s] Processing chr. 22 [███████████████████████████████---------] 77% [02m:45s] Processing chr. 22 [███████████████████████████████---------] 77% [02m:46s] Processing chr. 22 [███████████████████████████████---------] 77% [02m:47s] Processing chr. 22 [████████████████████████████████--------] 78% [02m:48s] Processing chr. 22 [████████████████████████████████--------] 78% [02m:49s] Processing chr. 22 [████████████████████████████████--------] 78% [02m:50s] Processing chr. 22 [████████████████████████████████--------] 78% [02m:51s] Processing chr. 22 [████████████████████████████████--------] 78% [02m:52s] Processing chr. 22 [████████████████████████████████--------] 78% [02m:53s] Processing chr. 22 [████████████████████████████████--------] 78% [02m:54s] Processing chr. 22 [████████████████████████████████--------] 79% [02m:55s] Processing chr. 22 [████████████████████████████████--------] 79% [02m:56s] Processing chr. 22 [████████████████████████████████--------] 79% [02m:57s] Processing chr. 22 [████████████████████████████████--------] 79% [02m:58s] Processing chr. 22 [████████████████████████████████--------] 79% [02m:59s] Processing chr. 22 [████████████████████████████████--------] 79% [03m:00s] Processing chr. 22 [████████████████████████████████--------] 79% [03m:01s] Processing chr. 22 [████████████████████████████████--------] 79% [03m:02s] Processing chr. 22 [█████████████████████████████████-------] 80% [03m:03s] Processing chr. 22 [█████████████████████████████████-------] 80% [03m:04s] Processing chr. 22 [█████████████████████████████████-------] 80% [03m:05s] Processing chr. 22 [█████████████████████████████████-------] 80% [03m:06s] Processing chr. 22 [█████████████████████████████████-------] 80% [03m:07s] Processing chr. 22 [█████████████████████████████████-------] 80% [03m:08s] Processing chr. 22 [█████████████████████████████████-------] 80% [03m:09s] Processing chr. 22 [█████████████████████████████████-------] 80% [03m:10s] Processing chr. 22 [█████████████████████████████████-------] 80% [03m:11s] Processing chr. 22 [█████████████████████████████████-------] 81% [03m:12s] Processing chr. 22 [█████████████████████████████████-------] 81% [03m:13s] Processing chr. 22 [█████████████████████████████████-------] 81% [03m:14s] Processing chr. 22 [█████████████████████████████████-------] 81% [03m:15s] Processing chr. 22 [█████████████████████████████████-------] 81% [03m:16s] Processing chr. 22 [█████████████████████████████████-------] 81% [03m:17s] Processing chr. 22 [█████████████████████████████████-------] 81% [03m:18s] Processing chr. 22 [█████████████████████████████████-------] 81% [03m:19s] Processing chr. 22 [█████████████████████████████████-------] 81% [03m:20s] Processing chr. 22 [█████████████████████████████████-------] 82% [03m:21s] Processing chr. 22 [█████████████████████████████████-------] 82% [03m:22s] Processing chr. 22 [█████████████████████████████████-------] 82% [03m:23s] Processing chr. 22 [█████████████████████████████████-------] 82% [03m:24s] Processing chr. 22 [█████████████████████████████████-------] 82% [03m:25s] Processing chr. 22 [█████████████████████████████████-------] 82% [03m:26s] Processing chr. 22 [█████████████████████████████████-------] 82% [03m:27s] Processing chr. 22 [█████████████████████████████████-------] 82% [03m:28s] Processing chr. 22 [██████████████████████████████████------] 83% [03m:29s] Processing chr. 22 [██████████████████████████████████------] 83% [03m:30s] Processing chr. 22 [██████████████████████████████████------] 83% [03m:31s] Processing chr. 22 [██████████████████████████████████------] 83% [03m:32s] Processing chr. 22 [██████████████████████████████████------] 83% [03m:33s] Processing chr. 22 [██████████████████████████████████------] 83% [03m:34s] Processing chr. 22 [██████████████████████████████████------] 83% [03m:35s] Processing chr. 22 [██████████████████████████████████------] 83% [03m:36s] Processing chr. 22 [██████████████████████████████████------] 84% [03m:37s] Processing chr. 22 [██████████████████████████████████------] 84% [03m:38s] Processing chr. 22 [██████████████████████████████████------] 84% [03m:39s] Processing chr. 22 [██████████████████████████████████------] 84% [03m:40s] Processing chr. 22 [██████████████████████████████████------] 84% [03m:41s] Processing chr. 22 [██████████████████████████████████------] 84% [03m:42s] Processing chr. 22 [██████████████████████████████████------] 84% [03m:43s] Processing chr. 22 [██████████████████████████████████------] 84% [03m:44s] Processing chr. 22 [███████████████████████████████████-----] 85% [03m:45s] Processing chr. 22 [███████████████████████████████████-----] 85% [03m:46s] Processing chr. 22 [███████████████████████████████████-----] 85% [03m:47s] Processing chr. 22 [███████████████████████████████████-----] 85% [03m:48s] Processing chr. 22 [███████████████████████████████████-----] 85% [03m:49s] Processing chr. 22 [███████████████████████████████████-----] 85% [03m:50s] Processing chr. 22 [███████████████████████████████████-----] 85% [03m:51s] Processing chr. 22 [███████████████████████████████████-----] 85% [03m:52s] Processing chr. 22 [███████████████████████████████████-----] 85% [03m:53s] Processing chr. 22 [███████████████████████████████████-----] 86% [03m:54s] Processing chr. 22 [███████████████████████████████████-----] 86% [03m:55s] Processing chr. 22 [███████████████████████████████████-----] 86% [03m:56s] Processing chr. 22 [███████████████████████████████████-----] 86% [03m:57s] Processing chr. 22 [███████████████████████████████████-----] 86% [03m:58s] Processing chr. 22 [███████████████████████████████████-----] 86% [03m:59s] Processing chr. 22 [███████████████████████████████████-----] 86% [04m:00s] Processing chr. 22 [███████████████████████████████████-----] 87% [04m:01s] Processing chr. 22 [███████████████████████████████████-----] 87% [04m:02s] Processing chr. 22 [███████████████████████████████████-----] 87% [04m:03s] Processing chr. 22 [███████████████████████████████████-----] 87% [04m:04s] Processing chr. 22 [███████████████████████████████████-----] 87% [04m:05s] Processing chr. 22 [███████████████████████████████████-----] 87% [04m:06s] Processing chr. 22 [███████████████████████████████████-----] 87% [04m:07s] Processing chr. 22 [███████████████████████████████████-----] 87% [04m:08s] Processing chr. 22 [████████████████████████████████████----] 88% [04m:09s] Processing chr. 22 [████████████████████████████████████----] 88% [04m:10s] Processing chr. 22 [████████████████████████████████████----] 88% [04m:11s] Processing chr. 22 [████████████████████████████████████----] 88% [04m:12s] Processing chr. 22 [████████████████████████████████████----] 88% [04m:13s] Processing chr. 22 [████████████████████████████████████----] 88% [04m:14s] Processing chr. 22 [████████████████████████████████████----] 88% [04m:15s] Processing chr. 22 [████████████████████████████████████----] 88% [04m:16s] Processing chr. 22 [████████████████████████████████████----] 89% [04m:17s] Processing chr. 22 [████████████████████████████████████----] 89% [04m:18s] Processing chr. 22 [████████████████████████████████████----] 89% [04m:19s] Processing chr. 22 [████████████████████████████████████----] 89% [04m:20s] Processing chr. 22 [████████████████████████████████████----] 89% [04m:21s] Processing chr. 22 [████████████████████████████████████----] 89% [04m:22s] Processing chr. 22 [████████████████████████████████████----] 89% [04m:23s] Processing chr. 22 [████████████████████████████████████----] 89% [04m:24s] Processing chr. 22 [█████████████████████████████████████---] 90% [04m:25s] Processing chr. 22 [█████████████████████████████████████---] 90% [04m:26s] Processing chr. 22 [█████████████████████████████████████---] 90% [04m:27s] Processing chr. 22 [█████████████████████████████████████---] 90% [04m:28s] Processing chr. 22 [█████████████████████████████████████---] 90% [04m:29s] Processing chr. 22 [█████████████████████████████████████---] 90% [04m:30s] Processing chr. 22 [█████████████████████████████████████---] 90% [04m:31s] Processing chr. 22 [█████████████████████████████████████---] 91% [04m:32s] Processing chr. 22 [█████████████████████████████████████---] 91% [04m:33s] Processing chr. 22 [█████████████████████████████████████---] 91% [04m:34s] Processing chr. 22 [█████████████████████████████████████---] 91% [04m:35s] Processing chr. 22 [█████████████████████████████████████---] 91% [04m:36s] Processing chr. 22 [█████████████████████████████████████---] 91% [04m:37s] Processing chr. 22 [█████████████████████████████████████---] 91% [04m:38s] Processing chr. 22 [█████████████████████████████████████---] 91% [04m:39s] Processing chr. 22 [█████████████████████████████████████---] 92% [04m:40s] Processing chr. 22 [█████████████████████████████████████---] 92% [04m:41s] Processing chr. 22 [█████████████████████████████████████---] 92% [04m:42s] Processing chr. 22 [█████████████████████████████████████---] 92% [04m:43s] Processing chr. 22 [█████████████████████████████████████---] 92% [04m:44s] Processing chr. 22 [█████████████████████████████████████---] 92% [04m:45s] Processing chr. 22 [█████████████████████████████████████---] 92% [04m:46s] Processing chr. 22 [██████████████████████████████████████--] 93% [04m:47s] Processing chr. 22 [██████████████████████████████████████--] 93% [04m:48s] Processing chr. 22 [██████████████████████████████████████--] 93% [04m:49s] Processing chr. 22 [██████████████████████████████████████--] 93% [04m:50s] Processing chr. 22 [██████████████████████████████████████--] 93% [04m:51s] Processing chr. 22 [██████████████████████████████████████--] 93% [04m:52s] Processing chr. 22 [██████████████████████████████████████--] 93% [04m:53s] Processing chr. 22 [██████████████████████████████████████--] 94% [04m:54s] Processing chr. 22 [██████████████████████████████████████--] 94% [04m:55s] Processing chr. 22 [██████████████████████████████████████--] 94% [04m:56s] Processing chr. 22 [██████████████████████████████████████--] 94% [04m:57s] Processing chr. 22 [██████████████████████████████████████--] 94% [04m:58s] Processing chr. 22 [██████████████████████████████████████--] 94% [04m:59s] Processing chr. 22 [██████████████████████████████████████--] 94% [05m:00s] Processing chr. 22 [██████████████████████████████████████--] 94% [05m:01s] Processing chr. 22 [██████████████████████████████████████--] 94% [05m:02s] Processing chr. 22 [███████████████████████████████████████-] 95% [05m:03s] Processing chr. 22 [███████████████████████████████████████-] 95% [05m:04s] Processing chr. 22 [███████████████████████████████████████-] 95% [05m:05s] Processing chr. 22 [███████████████████████████████████████-] 95% [05m:06s] Processing chr. 22 [███████████████████████████████████████-] 95% [05m:07s] Processing chr. 22 [███████████████████████████████████████-] 95% [05m:08s] Processing chr. 22 [███████████████████████████████████████-] 95% [05m:09s] Processing chr. 22 [███████████████████████████████████████-] 95% [05m:10s] Processing chr. 22 [███████████████████████████████████████-] 96% [05m:11s] Processing chr. 22 [███████████████████████████████████████-] 96% [05m:12s] Processing chr. 22 [███████████████████████████████████████-] 96% [05m:13s] Processing chr. 22 [███████████████████████████████████████-] 96% [05m:14s] Processing chr. 22 [███████████████████████████████████████-] 96% [05m:15s] Processing chr. 22 [███████████████████████████████████████-] 96% [05m:16s] Processing chr. 22 [███████████████████████████████████████-] 96% [05m:17s] Processing chr. 22 [███████████████████████████████████████-] 96% [05m:18s] Processing chr. 22 [███████████████████████████████████████-] 97% [05m:19s] Processing chr. 22 [███████████████████████████████████████-] 97% [05m:20s] Processing chr. 22 [███████████████████████████████████████-] 97% [05m:21s] Processing chr. 22 [███████████████████████████████████████-] 97% [05m:22s] Processing chr. 22 [███████████████████████████████████████-] 97% [05m:23s] Processing chr. 22 [███████████████████████████████████████-] 97% [05m:24s] Processing chr. 22 [███████████████████████████████████████-] 97% [05m:25s] Processing chr. 22 [████████████████████████████████████████] 98% [05m:26s] Processing chr. 22 [████████████████████████████████████████] 98% [05m:27s] Processing chr. 22 [████████████████████████████████████████] 98% [05m:28s] Processing chr. 22 [████████████████████████████████████████] 98% [05m:29s] Processing chr. 22 [████████████████████████████████████████] 98% [05m:30s] Processing chr. 22 [████████████████████████████████████████] 98% [05m:31s] Processing chr. 22 [████████████████████████████████████████] 98% [05m:32s] Processing chr. 22 [████████████████████████████████████████] 98% [05m:33s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:34s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:35s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:36s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:37s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:38s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:39s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:40s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:41s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:42s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:43s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:44s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:45s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:46s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:47s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:48s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:49s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:50s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:51s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:52s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:53s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:54s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:55s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:56s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:57s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:58s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:59s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:00s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:01s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:02s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:03s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:04s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:05s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:06s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:07s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:08s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:09s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:10s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:11s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:12s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:13s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:14s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:15s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:16s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:17s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:18s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:19s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:20s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:21s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:22s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:23s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:24s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:25s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:26s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:27s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:28s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:29s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:30s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:31s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:32s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:33s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:34s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:35s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:36s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:37s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:38s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:39s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:40s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:41s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:42s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:43s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:44s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:45s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:46s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:47s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:48s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:49s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:50s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:51s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:52s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:53s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:54s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:55s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:56s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:57s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:58s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:59s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:00s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:01s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:02s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:03s] Processing chr. 22 [████████████████████████████████████████] 100% [07m:06s] Reads simulated
+#>  [█---------------------------------------] 0% [00m:00s] Found 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:01s] Reading 22 [█---------------------------------------] 0% [00m:01s] Processing chr. 22 [█---------------------------------------] 0% [00m:05s] Processing chr. 22 [█---------------------------------------] 0% [00m:06s] Processing chr. 22 [█---------------------------------------] 1% [00m:07s] Processing chr. 22 [█---------------------------------------] 1% [00m:08s] Processing chr. 22 [█---------------------------------------] 2% [00m:09s] Processing chr. 22 [█---------------------------------------] 2% [00m:10s] Processing chr. 22 [██--------------------------------------] 3% [00m:11s] Processing chr. 22 [██--------------------------------------] 3% [00m:12s] Processing chr. 22 [██--------------------------------------] 4% [00m:13s] Processing chr. 22 [███-------------------------------------] 5% [00m:14s] Processing chr. 22 [███-------------------------------------] 5% [00m:15s] Processing chr. 22 [███-------------------------------------] 6% [00m:16s] Processing chr. 22 [███-------------------------------------] 6% [00m:17s] Processing chr. 22 [███-------------------------------------] 7% [00m:18s] Processing chr. 22 [███-------------------------------------] 7% [00m:19s] Processing chr. 22 [████------------------------------------] 8% [00m:20s] Processing chr. 22 [████------------------------------------] 8% [00m:21s] Processing chr. 22 [████------------------------------------] 9% [00m:22s] Processing chr. 22 [█████-----------------------------------] 10% [00m:23s] Processing chr. 22 [█████-----------------------------------] 10% [00m:24s] Processing chr. 22 [█████-----------------------------------] 11% [00m:25s] Processing chr. 22 [█████-----------------------------------] 11% [00m:26s] Processing chr. 22 [█████-----------------------------------] 12% [00m:27s] Processing chr. 22 [█████-----------------------------------] 12% [00m:28s] Processing chr. 22 [██████----------------------------------] 13% [00m:29s] Processing chr. 22 [██████----------------------------------] 13% [00m:30s] Processing chr. 22 [██████----------------------------------] 14% [00m:31s] Processing chr. 22 [██████----------------------------------] 14% [00m:32s] Processing chr. 22 [███████---------------------------------] 15% [00m:33s] Processing chr. 22 [███████---------------------------------] 16% [00m:34s] Processing chr. 22 [███████---------------------------------] 16% [00m:35s] Processing chr. 22 [███████---------------------------------] 17% [00m:36s] Processing chr. 22 [███████---------------------------------] 17% [00m:37s] Processing chr. 22 [████████--------------------------------] 18% [00m:38s] Processing chr. 22 [████████--------------------------------] 18% [00m:39s] Processing chr. 22 [████████--------------------------------] 19% [00m:40s] Processing chr. 22 [█████████-------------------------------] 20% [00m:41s] Processing chr. 22 [█████████-------------------------------] 20% [00m:42s] Processing chr. 22 [█████████-------------------------------] 21% [00m:43s] Processing chr. 22 [█████████-------------------------------] 21% [00m:44s] Processing chr. 22 [█████████-------------------------------] 22% [00m:45s] Processing chr. 22 [█████████-------------------------------] 22% [00m:46s] Processing chr. 22 [██████████------------------------------] 23% [00m:47s] Processing chr. 22 [██████████------------------------------] 23% [00m:48s] Processing chr. 22 [██████████------------------------------] 24% [00m:49s] Processing chr. 22 [███████████-----------------------------] 25% [00m:50s] Processing chr. 22 [███████████-----------------------------] 25% [00m:51s] Processing chr. 22 [███████████-----------------------------] 26% [00m:52s] Processing chr. 22 [███████████-----------------------------] 26% [00m:53s] Processing chr. 22 [███████████-----------------------------] 27% [00m:54s] Processing chr. 22 [███████████-----------------------------] 27% [00m:55s] Processing chr. 22 [████████████----------------------------] 28% [00m:56s] Processing chr. 22 [████████████----------------------------] 29% [00m:57s] Processing chr. 22 [████████████----------------------------] 29% [00m:58s] Processing chr. 22 [█████████████---------------------------] 30% [00m:59s] Processing chr. 22 [█████████████---------------------------] 30% [01m:00s] Processing chr. 22 [█████████████---------------------------] 31% [01m:01s] Processing chr. 22 [█████████████---------------------------] 31% [01m:02s] Processing chr. 22 [█████████████---------------------------] 32% [01m:03s] Processing chr. 22 [██████████████--------------------------] 33% [01m:04s] Processing chr. 22 [██████████████--------------------------] 33% [01m:05s] Processing chr. 22 [██████████████--------------------------] 34% [01m:06s] Processing chr. 22 [██████████████--------------------------] 34% [01m:07s] Processing chr. 22 [███████████████-------------------------] 35% [01m:08s] Processing chr. 22 [███████████████-------------------------] 35% [01m:09s] Processing chr. 22 [███████████████-------------------------] 36% [01m:10s] Processing chr. 22 [███████████████-------------------------] 36% [01m:11s] Processing chr. 22 [███████████████-------------------------] 37% [01m:12s] Processing chr. 22 [████████████████------------------------] 38% [01m:13s] Processing chr. 22 [████████████████------------------------] 38% [01m:14s] Processing chr. 22 [████████████████------------------------] 39% [01m:15s] Processing chr. 22 [████████████████------------------------] 39% [01m:16s] Processing chr. 22 [█████████████████-----------------------] 40% [01m:17s] Processing chr. 22 [█████████████████-----------------------] 40% [01m:18s] Processing chr. 22 [█████████████████-----------------------] 41% [01m:19s] Processing chr. 22 [█████████████████-----------------------] 42% [01m:20s] Processing chr. 22 [█████████████████-----------------------] 42% [01m:21s] Processing chr. 22 [██████████████████----------------------] 43% [01m:22s] Processing chr. 22 [██████████████████----------------------] 43% [01m:23s] Processing chr. 22 [██████████████████----------------------] 44% [01m:24s] Processing chr. 22 [██████████████████----------------------] 44% [01m:25s] Processing chr. 22 [███████████████████---------------------] 45% [01m:26s] Processing chr. 22 [███████████████████---------------------] 46% [01m:27s] Processing chr. 22 [███████████████████---------------------] 46% [01m:28s] Processing chr. 22 [███████████████████---------------------] 47% [01m:29s] Processing chr. 22 [███████████████████---------------------] 47% [01m:30s] Processing chr. 22 [████████████████████--------------------] 48% [01m:31s] Processing chr. 22 [████████████████████--------------------] 48% [01m:32s] Processing chr. 22 [████████████████████--------------------] 48% [01m:33s] Processing chr. 22 [████████████████████--------------------] 49% [01m:34s] Processing chr. 22 [████████████████████--------------------] 49% [01m:35s] Processing chr. 22 [████████████████████--------------------] 49% [01m:36s] Processing chr. 22 [█████████████████████-------------------] 50% [01m:37s] Processing chr. 22 [█████████████████████-------------------] 50% [01m:38s] Processing chr. 22 [█████████████████████-------------------] 50% [01m:39s] Processing chr. 22 [█████████████████████-------------------] 51% [01m:40s] Processing chr. 22 [█████████████████████-------------------] 51% [01m:41s] Processing chr. 22 [█████████████████████-------------------] 51% [01m:42s] Processing chr. 22 [█████████████████████-------------------] 52% [01m:43s] Processing chr. 22 [█████████████████████-------------------] 52% [01m:44s] Processing chr. 22 [█████████████████████-------------------] 52% [01m:45s] Processing chr. 22 [██████████████████████------------------] 52% [01m:46s] Processing chr. 22 [██████████████████████------------------] 52% [01m:47s] Processing chr. 22 [██████████████████████------------------] 52% [01m:48s] Processing chr. 22 [██████████████████████------------------] 54% [01m:49s] Processing chr. 22 [██████████████████████------------------] 54% [01m:50s] Processing chr. 22 [██████████████████████------------------] 54% [01m:51s] Processing chr. 22 [███████████████████████-----------------] 55% [01m:52s] Processing chr. 22 [███████████████████████-----------------] 55% [01m:53s] Processing chr. 22 [███████████████████████-----------------] 55% [01m:54s] Processing chr. 22 [███████████████████████-----------------] 56% [01m:55s] Processing chr. 22 [███████████████████████-----------------] 56% [01m:56s] Processing chr. 22 [███████████████████████-----------------] 56% [01m:57s] Processing chr. 22 [███████████████████████-----------------] 57% [01m:58s] Processing chr. 22 [███████████████████████-----------------] 57% [01m:59s] Processing chr. 22 [████████████████████████----------------] 58% [02m:00s] Processing chr. 22 [████████████████████████----------------] 58% [02m:01s] Processing chr. 22 [████████████████████████----------------] 58% [02m:02s] Processing chr. 22 [████████████████████████----------------] 58% [02m:03s] Processing chr. 22 [████████████████████████----------------] 58% [02m:04s] Processing chr. 22 [████████████████████████----------------] 58% [02m:05s] Processing chr. 22 [█████████████████████████---------------] 60% [02m:06s] Processing chr. 22 [█████████████████████████---------------] 60% [02m:07s] Processing chr. 22 [█████████████████████████---------------] 60% [02m:08s] Processing chr. 22 [█████████████████████████---------------] 61% [02m:09s] Processing chr. 22 [█████████████████████████---------------] 61% [02m:10s] Processing chr. 22 [█████████████████████████---------------] 61% [02m:11s] Processing chr. 22 [█████████████████████████---------------] 62% [02m:12s] Processing chr. 22 [█████████████████████████---------------] 62% [02m:13s] Processing chr. 22 [█████████████████████████---------------] 62% [02m:14s] Processing chr. 22 [██████████████████████████--------------] 63% [02m:15s] Processing chr. 22 [██████████████████████████--------------] 63% [02m:16s] Processing chr. 22 [██████████████████████████--------------] 63% [02m:17s] Processing chr. 22 [██████████████████████████--------------] 64% [02m:18s] Processing chr. 22 [██████████████████████████--------------] 64% [02m:19s] Processing chr. 22 [██████████████████████████--------------] 64% [02m:20s] Processing chr. 22 [███████████████████████████-------------] 65% [02m:21s] Processing chr. 22 [███████████████████████████-------------] 65% [02m:22s] Processing chr. 22 [███████████████████████████-------------] 65% [02m:23s] Processing chr. 22 [███████████████████████████-------------] 66% [02m:24s] Processing chr. 22 [███████████████████████████-------------] 66% [02m:25s] Processing chr. 22 [███████████████████████████-------------] 66% [02m:26s] Processing chr. 22 [███████████████████████████-------------] 67% [02m:27s] Processing chr. 22 [███████████████████████████-------------] 67% [02m:28s] Processing chr. 22 [███████████████████████████-------------] 67% [02m:29s] Processing chr. 22 [████████████████████████████------------] 68% [02m:30s] Processing chr. 22 [████████████████████████████------------] 68% [02m:31s] Processing chr. 22 [████████████████████████████------------] 68% [02m:32s] Processing chr. 22 [████████████████████████████------------] 69% [02m:33s] Processing chr. 22 [████████████████████████████------------] 69% [02m:34s] Processing chr. 22 [████████████████████████████------------] 69% [02m:35s] Processing chr. 22 [█████████████████████████████-----------] 70% [02m:36s] Processing chr. 22 [█████████████████████████████-----------] 70% [02m:37s] Processing chr. 22 [█████████████████████████████-----------] 70% [02m:38s] Processing chr. 22 [█████████████████████████████-----------] 71% [02m:39s] Processing chr. 22 [█████████████████████████████-----------] 71% [02m:40s] Processing chr. 22 [█████████████████████████████-----------] 71% [02m:41s] Processing chr. 22 [█████████████████████████████-----------] 72% [02m:42s] Processing chr. 22 [█████████████████████████████-----------] 72% [02m:43s] Processing chr. 22 [█████████████████████████████-----------] 72% [02m:44s] Processing chr. 22 [██████████████████████████████----------] 73% [02m:45s] Processing chr. 22 [██████████████████████████████----------] 73% [02m:46s] Processing chr. 22 [██████████████████████████████----------] 73% [02m:47s] Processing chr. 22 [██████████████████████████████----------] 74% [02m:48s] Processing chr. 22 [██████████████████████████████----------] 74% [02m:49s] Processing chr. 22 [███████████████████████████████---------] 75% [02m:50s] Processing chr. 22 [███████████████████████████████---------] 75% [02m:51s] Processing chr. 22 [███████████████████████████████---------] 75% [02m:52s] Processing chr. 22 [███████████████████████████████---------] 76% [02m:53s] Processing chr. 22 [███████████████████████████████---------] 76% [02m:54s] Processing chr. 22 [███████████████████████████████---------] 76% [02m:55s] Processing chr. 22 [███████████████████████████████---------] 77% [02m:56s] Processing chr. 22 [███████████████████████████████---------] 77% [02m:57s] Processing chr. 22 [███████████████████████████████---------] 77% [02m:58s] Processing chr. 22 [████████████████████████████████--------] 78% [02m:59s] Processing chr. 22 [████████████████████████████████--------] 78% [03m:00s] Processing chr. 22 [████████████████████████████████--------] 78% [03m:01s] Processing chr. 22 [████████████████████████████████--------] 78% [03m:02s] Processing chr. 22 [████████████████████████████████--------] 78% [03m:03s] Processing chr. 22 [████████████████████████████████--------] 78% [03m:04s] Processing chr. 22 [████████████████████████████████--------] 78% [03m:05s] Processing chr. 22 [████████████████████████████████--------] 79% [03m:06s] Processing chr. 22 [████████████████████████████████--------] 79% [03m:07s] Processing chr. 22 [████████████████████████████████--------] 79% [03m:08s] Processing chr. 22 [████████████████████████████████--------] 79% [03m:09s] Processing chr. 22 [████████████████████████████████--------] 79% [03m:10s] Processing chr. 22 [████████████████████████████████--------] 79% [03m:11s] Processing chr. 22 [████████████████████████████████--------] 79% [03m:12s] Processing chr. 22 [████████████████████████████████--------] 79% [03m:13s] Processing chr. 22 [█████████████████████████████████-------] 80% [03m:14s] Processing chr. 22 [█████████████████████████████████-------] 80% [03m:15s] Processing chr. 22 [█████████████████████████████████-------] 80% [03m:16s] Processing chr. 22 [█████████████████████████████████-------] 80% [03m:17s] Processing chr. 22 [█████████████████████████████████-------] 80% [03m:18s] Processing chr. 22 [█████████████████████████████████-------] 80% [03m:19s] Processing chr. 22 [█████████████████████████████████-------] 80% [03m:20s] Processing chr. 22 [█████████████████████████████████-------] 80% [03m:21s] Processing chr. 22 [█████████████████████████████████-------] 80% [03m:22s] Processing chr. 22 [█████████████████████████████████-------] 81% [03m:23s] Processing chr. 22 [█████████████████████████████████-------] 81% [03m:24s] Processing chr. 22 [█████████████████████████████████-------] 81% [03m:25s] Processing chr. 22 [█████████████████████████████████-------] 81% [03m:26s] Processing chr. 22 [█████████████████████████████████-------] 81% [03m:27s] Processing chr. 22 [█████████████████████████████████-------] 81% [03m:28s] Processing chr. 22 [█████████████████████████████████-------] 81% [03m:29s] Processing chr. 22 [█████████████████████████████████-------] 81% [03m:30s] Processing chr. 22 [█████████████████████████████████-------] 82% [03m:31s] Processing chr. 22 [█████████████████████████████████-------] 82% [03m:32s] Processing chr. 22 [█████████████████████████████████-------] 82% [03m:33s] Processing chr. 22 [█████████████████████████████████-------] 82% [03m:34s] Processing chr. 22 [█████████████████████████████████-------] 82% [03m:35s] Processing chr. 22 [█████████████████████████████████-------] 82% [03m:36s] Processing chr. 22 [█████████████████████████████████-------] 82% [03m:37s] Processing chr. 22 [█████████████████████████████████-------] 82% [03m:38s] Processing chr. 22 [██████████████████████████████████------] 83% [03m:39s] Processing chr. 22 [██████████████████████████████████------] 83% [03m:40s] Processing chr. 22 [██████████████████████████████████------] 83% [03m:41s] Processing chr. 22 [██████████████████████████████████------] 83% [03m:42s] Processing chr. 22 [██████████████████████████████████------] 83% [03m:43s] Processing chr. 22 [██████████████████████████████████------] 83% [03m:44s] Processing chr. 22 [██████████████████████████████████------] 83% [03m:45s] Processing chr. 22 [██████████████████████████████████------] 83% [03m:46s] Processing chr. 22 [██████████████████████████████████------] 84% [03m:47s] Processing chr. 22 [██████████████████████████████████------] 84% [03m:48s] Processing chr. 22 [██████████████████████████████████------] 84% [03m:49s] Processing chr. 22 [██████████████████████████████████------] 84% [03m:50s] Processing chr. 22 [██████████████████████████████████------] 84% [03m:51s] Processing chr. 22 [██████████████████████████████████------] 84% [03m:52s] Processing chr. 22 [██████████████████████████████████------] 84% [03m:53s] Processing chr. 22 [██████████████████████████████████------] 84% [03m:54s] Processing chr. 22 [██████████████████████████████████------] 84% [03m:55s] Processing chr. 22 [███████████████████████████████████-----] 85% [03m:56s] Processing chr. 22 [███████████████████████████████████-----] 85% [03m:57s] Processing chr. 22 [███████████████████████████████████-----] 85% [03m:58s] Processing chr. 22 [███████████████████████████████████-----] 85% [03m:59s] Processing chr. 22 [███████████████████████████████████-----] 85% [04m:00s] Processing chr. 22 [███████████████████████████████████-----] 85% [04m:01s] Processing chr. 22 [███████████████████████████████████-----] 85% [04m:02s] Processing chr. 22 [███████████████████████████████████-----] 85% [04m:03s] Processing chr. 22 [███████████████████████████████████-----] 86% [04m:04s] Processing chr. 22 [███████████████████████████████████-----] 86% [04m:05s] Processing chr. 22 [███████████████████████████████████-----] 86% [04m:06s] Processing chr. 22 [███████████████████████████████████-----] 86% [04m:07s] Processing chr. 22 [███████████████████████████████████-----] 86% [04m:08s] Processing chr. 22 [███████████████████████████████████-----] 86% [04m:09s] Processing chr. 22 [███████████████████████████████████-----] 86% [04m:10s] Processing chr. 22 [███████████████████████████████████-----] 86% [04m:11s] Processing chr. 22 [███████████████████████████████████-----] 87% [04m:12s] Processing chr. 22 [███████████████████████████████████-----] 87% [04m:13s] Processing chr. 22 [███████████████████████████████████-----] 87% [04m:14s] Processing chr. 22 [███████████████████████████████████-----] 87% [04m:15s] Processing chr. 22 [███████████████████████████████████-----] 87% [04m:16s] Processing chr. 22 [███████████████████████████████████-----] 87% [04m:17s] Processing chr. 22 [███████████████████████████████████-----] 87% [04m:18s] Processing chr. 22 [███████████████████████████████████-----] 87% [04m:19s] Processing chr. 22 [████████████████████████████████████----] 88% [04m:20s] Processing chr. 22 [████████████████████████████████████----] 88% [04m:21s] Processing chr. 22 [████████████████████████████████████----] 88% [04m:22s] Processing chr. 22 [████████████████████████████████████----] 88% [04m:23s] Processing chr. 22 [████████████████████████████████████----] 88% [04m:24s] Processing chr. 22 [████████████████████████████████████----] 88% [04m:25s] Processing chr. 22 [████████████████████████████████████----] 88% [04m:26s] Processing chr. 22 [████████████████████████████████████----] 88% [04m:27s] Processing chr. 22 [████████████████████████████████████----] 89% [04m:28s] Processing chr. 22 [████████████████████████████████████----] 89% [04m:29s] Processing chr. 22 [████████████████████████████████████----] 89% [04m:30s] Processing chr. 22 [████████████████████████████████████----] 89% [04m:31s] Processing chr. 22 [████████████████████████████████████----] 89% [04m:32s] Processing chr. 22 [████████████████████████████████████----] 89% [04m:33s] Processing chr. 22 [████████████████████████████████████----] 89% [04m:34s] Processing chr. 22 [████████████████████████████████████----] 89% [04m:35s] Processing chr. 22 [█████████████████████████████████████---] 90% [04m:36s] Processing chr. 22 [█████████████████████████████████████---] 90% [04m:37s] Processing chr. 22 [█████████████████████████████████████---] 90% [04m:38s] Processing chr. 22 [█████████████████████████████████████---] 90% [04m:39s] Processing chr. 22 [█████████████████████████████████████---] 90% [04m:40s] Processing chr. 22 [█████████████████████████████████████---] 90% [04m:41s] Processing chr. 22 [█████████████████████████████████████---] 90% [04m:42s] Processing chr. 22 [█████████████████████████████████████---] 90% [04m:43s] Processing chr. 22 [█████████████████████████████████████---] 91% [04m:44s] Processing chr. 22 [█████████████████████████████████████---] 91% [04m:45s] Processing chr. 22 [█████████████████████████████████████---] 91% [04m:46s] Processing chr. 22 [█████████████████████████████████████---] 91% [04m:47s] Processing chr. 22 [█████████████████████████████████████---] 91% [04m:48s] Processing chr. 22 [█████████████████████████████████████---] 91% [04m:49s] Processing chr. 22 [█████████████████████████████████████---] 91% [04m:50s] Processing chr. 22 [█████████████████████████████████████---] 91% [04m:51s] Processing chr. 22 [█████████████████████████████████████---] 92% [04m:52s] Processing chr. 22 [█████████████████████████████████████---] 92% [04m:53s] Processing chr. 22 [█████████████████████████████████████---] 92% [04m:54s] Processing chr. 22 [█████████████████████████████████████---] 92% [04m:55s] Processing chr. 22 [█████████████████████████████████████---] 92% [04m:56s] Processing chr. 22 [█████████████████████████████████████---] 92% [04m:57s] Processing chr. 22 [█████████████████████████████████████---] 92% [04m:58s] Processing chr. 22 [██████████████████████████████████████--] 93% [04m:59s] Processing chr. 22 [██████████████████████████████████████--] 93% [05m:00s] Processing chr. 22 [██████████████████████████████████████--] 93% [05m:01s] Processing chr. 22 [██████████████████████████████████████--] 93% [05m:02s] Processing chr. 22 [██████████████████████████████████████--] 93% [05m:03s] Processing chr. 22 [██████████████████████████████████████--] 93% [05m:04s] Processing chr. 22 [██████████████████████████████████████--] 93% [05m:05s] Processing chr. 22 [██████████████████████████████████████--] 93% [05m:06s] Processing chr. 22 [██████████████████████████████████████--] 94% [05m:07s] Processing chr. 22 [██████████████████████████████████████--] 94% [05m:08s] Processing chr. 22 [██████████████████████████████████████--] 94% [05m:09s] Processing chr. 22 [██████████████████████████████████████--] 94% [05m:10s] Processing chr. 22 [██████████████████████████████████████--] 94% [05m:11s] Processing chr. 22 [██████████████████████████████████████--] 94% [05m:12s] Processing chr. 22 [██████████████████████████████████████--] 94% [05m:13s] Processing chr. 22 [██████████████████████████████████████--] 94% [05m:14s] Processing chr. 22 [███████████████████████████████████████-] 95% [05m:15s] Processing chr. 22 [███████████████████████████████████████-] 95% [05m:16s] Processing chr. 22 [███████████████████████████████████████-] 95% [05m:17s] Processing chr. 22 [███████████████████████████████████████-] 95% [05m:18s] Processing chr. 22 [███████████████████████████████████████-] 95% [05m:19s] Processing chr. 22 [███████████████████████████████████████-] 95% [05m:20s] Processing chr. 22 [███████████████████████████████████████-] 95% [05m:21s] Processing chr. 22 [███████████████████████████████████████-] 95% [05m:22s] Processing chr. 22 [███████████████████████████████████████-] 96% [05m:23s] Processing chr. 22 [███████████████████████████████████████-] 96% [05m:24s] Processing chr. 22 [███████████████████████████████████████-] 96% [05m:25s] Processing chr. 22 [███████████████████████████████████████-] 96% [05m:26s] Processing chr. 22 [███████████████████████████████████████-] 96% [05m:27s] Processing chr. 22 [███████████████████████████████████████-] 96% [05m:28s] Processing chr. 22 [███████████████████████████████████████-] 96% [05m:29s] Processing chr. 22 [███████████████████████████████████████-] 96% [05m:30s] Processing chr. 22 [███████████████████████████████████████-] 97% [05m:31s] Processing chr. 22 [███████████████████████████████████████-] 97% [05m:32s] Processing chr. 22 [███████████████████████████████████████-] 97% [05m:33s] Processing chr. 22 [███████████████████████████████████████-] 97% [05m:34s] Processing chr. 22 [███████████████████████████████████████-] 97% [05m:35s] Processing chr. 22 [███████████████████████████████████████-] 97% [05m:36s] Processing chr. 22 [███████████████████████████████████████-] 97% [05m:37s] Processing chr. 22 [████████████████████████████████████████] 98% [05m:38s] Processing chr. 22 [████████████████████████████████████████] 98% [05m:39s] Processing chr. 22 [████████████████████████████████████████] 98% [05m:40s] Processing chr. 22 [████████████████████████████████████████] 98% [05m:41s] Processing chr. 22 [████████████████████████████████████████] 98% [05m:42s] Processing chr. 22 [████████████████████████████████████████] 98% [05m:43s] Processing chr. 22 [████████████████████████████████████████] 98% [05m:44s] Processing chr. 22 [████████████████████████████████████████] 98% [05m:45s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:46s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:47s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:48s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:49s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:50s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:51s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:52s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:53s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:54s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:55s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:56s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:57s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:58s] Processing chr. 22 [████████████████████████████████████████] 99% [05m:59s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:00s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:01s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:02s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:03s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:04s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:05s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:06s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:07s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:08s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:09s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:10s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:11s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:12s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:13s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:14s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:15s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:16s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:17s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:18s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:19s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:20s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:21s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:22s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:23s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:24s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:25s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:26s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:27s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:28s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:29s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:30s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:31s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:32s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:33s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:34s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:35s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:36s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:37s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:38s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:39s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:40s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:41s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:42s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:43s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:44s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:45s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:46s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:47s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:48s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:49s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:50s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:51s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:52s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:53s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:54s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:55s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:56s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:57s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:58s] Processing chr. 22 [████████████████████████████████████████] 99% [06m:59s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:00s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:01s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:02s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:03s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:04s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:05s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:06s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:07s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:08s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:09s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:10s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:11s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:12s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:13s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:14s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:15s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:16s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:17s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:18s] Processing chr. 22 [████████████████████████████████████████] 99% [07m:19s] Processing chr. 22 [████████████████████████████████████████] 100% [07m:22s] Reads simulated
 
 # let us load the dplyr library to filter the `simulate_seq` output
 library(dplyr)
@@ -254,7 +254,7 @@ 

Saving the Simulated Reads
 seq_results <- simulate_seq(phylo_forest, coverage = 2.5, write_SAM = TRUE)
-#>  [█---------------------------------------] 0% [00m:00s] Found 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:01s] Processing chr. 22 [█---------------------------------------] 0% [00m:05s] Processing chr. 22 [██--------------------------------------] 4% [00m:06s] Processing chr. 22 [█████-----------------------------------] 10% [00m:07s] Processing chr. 22 [███████---------------------------------] 15% [00m:08s] Processing chr. 22 [█████████-------------------------------] 20% [00m:09s] Processing chr. 22 [███████████-----------------------------] 25% [00m:10s] Processing chr. 22 [█████████████---------------------------] 30% [00m:11s] Processing chr. 22 [███████████████-------------------------] 35% [00m:12s] Processing chr. 22 [█████████████████-----------------------] 40% [00m:13s] Processing chr. 22 [███████████████████---------------------] 45% [00m:14s] Processing chr. 22 [████████████████████--------------------] 49% [00m:15s] Processing chr. 22 [█████████████████████-------------------] 52% [00m:16s] Processing chr. 22 [███████████████████████-----------------] 55% [00m:17s] Processing chr. 22 [████████████████████████----------------] 58% [00m:18s] Processing chr. 22 [█████████████████████████---------------] 62% [00m:19s] Processing chr. 22 [███████████████████████████-------------] 65% [00m:20s] Processing chr. 22 [████████████████████████████------------] 68% [00m:21s] Processing chr. 22 [█████████████████████████████-----------] 71% [00m:22s] Processing chr. 22 [██████████████████████████████----------] 74% [00m:23s] Processing chr. 22 [███████████████████████████████---------] 77% [00m:24s] Processing chr. 22 [████████████████████████████████--------] 79% [00m:25s] Processing chr. 22 [█████████████████████████████████-------] 80% [00m:26s] Processing chr. 22 [█████████████████████████████████-------] 81% [00m:27s] Processing chr. 22 [█████████████████████████████████-------] 82% [00m:28s] Processing chr. 22 [██████████████████████████████████------] 83% [00m:29s] Processing chr. 22 [██████████████████████████████████------] 84% [00m:30s] Processing chr. 22 [███████████████████████████████████-----] 85% [00m:31s] Processing chr. 22 [███████████████████████████████████-----] 86% [00m:32s] Processing chr. 22 [███████████████████████████████████-----] 87% [00m:33s] Processing chr. 22 [████████████████████████████████████----] 88% [00m:34s] Processing chr. 22 [████████████████████████████████████----] 89% [00m:35s] Processing chr. 22 [█████████████████████████████████████---] 90% [00m:36s] Processing chr. 22 [█████████████████████████████████████---] 91% [00m:37s] Processing chr. 22 [██████████████████████████████████████--] 93% [00m:38s] Processing chr. 22 [██████████████████████████████████████--] 94% [00m:39s] Processing chr. 22 [███████████████████████████████████████-] 95% [00m:40s] Processing chr. 22 [███████████████████████████████████████-] 96% [00m:41s] Processing chr. 22 [███████████████████████████████████████-] 97% [00m:42s] Processing chr. 22 [████████████████████████████████████████] 98% [00m:43s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:44s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:45s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:46s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:47s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:48s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:49s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:50s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:51s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:52s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:53s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:56s] Reads simulated
+#>  [█---------------------------------------] 0% [00m:00s] Found 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:01s] Reading 22 [█---------------------------------------] 0% [00m:01s] Processing chr. 22 [█---------------------------------------] 0% [00m:05s] Processing chr. 22 [███-------------------------------------] 5% [00m:06s] Processing chr. 22 [█████-----------------------------------] 10% [00m:07s] Processing chr. 22 [██████----------------------------------] 14% [00m:08s] Processing chr. 22 [████████--------------------------------] 19% [00m:09s] Processing chr. 22 [██████████------------------------------] 24% [00m:10s] Processing chr. 22 [████████████----------------------------] 28% [00m:11s] Processing chr. 22 [██████████████--------------------------] 33% [00m:12s] Processing chr. 22 [███████████████-------------------------] 37% [00m:13s] Processing chr. 22 [█████████████████-----------------------] 42% [00m:14s] Processing chr. 22 [███████████████████---------------------] 47% [00m:15s] Processing chr. 22 [████████████████████--------------------] 49% [00m:16s] Processing chr. 22 [█████████████████████-------------------] 52% [00m:17s] Processing chr. 22 [███████████████████████-----------------] 55% [00m:18s] Processing chr. 22 [████████████████████████----------------] 58% [00m:19s] Processing chr. 22 [█████████████████████████---------------] 61% [00m:20s] Processing chr. 22 [██████████████████████████--------------] 64% [00m:21s] Processing chr. 22 [███████████████████████████-------------] 67% [00m:22s] Processing chr. 22 [█████████████████████████████-----------] 70% [00m:23s] Processing chr. 22 [██████████████████████████████----------] 73% [00m:24s] Processing chr. 22 [███████████████████████████████---------] 75% [00m:25s] Processing chr. 22 [████████████████████████████████--------] 78% [00m:26s] Processing chr. 22 [████████████████████████████████--------] 79% [00m:27s] Processing chr. 22 [█████████████████████████████████-------] 80% [00m:28s] Processing chr. 22 [█████████████████████████████████-------] 81% [00m:29s] Processing chr. 22 [█████████████████████████████████-------] 82% [00m:30s] Processing chr. 22 [██████████████████████████████████------] 83% [00m:31s] Processing chr. 22 [██████████████████████████████████------] 84% [00m:32s] Processing chr. 22 [███████████████████████████████████-----] 85% [00m:33s] Processing chr. 22 [███████████████████████████████████-----] 86% [00m:34s] Processing chr. 22 [███████████████████████████████████-----] 87% [00m:35s] Processing chr. 22 [████████████████████████████████████----] 88% [00m:36s] Processing chr. 22 [████████████████████████████████████----] 89% [00m:37s] Processing chr. 22 [█████████████████████████████████████---] 90% [00m:38s] Processing chr. 22 [█████████████████████████████████████---] 91% [00m:39s] Processing chr. 22 [█████████████████████████████████████---] 92% [00m:40s] Processing chr. 22 [█████████████████████████████████████---] 92% [00m:41s] Processing chr. 22 [██████████████████████████████████████--] 94% [00m:42s] Processing chr. 22 [███████████████████████████████████████-] 95% [00m:43s] Processing chr. 22 [███████████████████████████████████████-] 96% [00m:44s] Processing chr. 22 [███████████████████████████████████████-] 97% [00m:45s] Processing chr. 22 [████████████████████████████████████████] 98% [00m:46s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:47s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:48s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:49s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:50s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:51s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:52s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:53s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:54s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:55s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:56s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:57s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:59s] Reads simulated
 
 sam_files <- list.files("rRACES_SAM/")
 
@@ -301,7 +301,7 @@ 

Sample Purity# tumour cells seq_results <- simulate_seq(phylo_forest, coverage = 2.5, purity = 0.7, write_SAM = TRUE, output_dir = "SAM_0.7") -#> [█---------------------------------------] 0% [00m:00s] Found 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:01s] Processing chr. 22 [█---------------------------------------] 0% [00m:05s] Processing chr. 22 [███-------------------------------------] 5% [00m:06s] Processing chr. 22 [█████-----------------------------------] 10% [00m:07s] Processing chr. 22 [███████---------------------------------] 15% [00m:08s] Processing chr. 22 [█████████-------------------------------] 20% [00m:09s] Processing chr. 22 [███████████-----------------------------] 25% [00m:10s] Processing chr. 22 [█████████████---------------------------] 30% [00m:11s] Processing chr. 22 [███████████████-------------------------] 35% [00m:12s] Processing chr. 22 [█████████████████-----------------------] 40% [00m:13s] Processing chr. 22 [███████████████████---------------------] 45% [00m:14s] Processing chr. 22 [████████████████████--------------------] 49% [00m:15s] Processing chr. 22 [█████████████████████-------------------] 52% [00m:16s] Processing chr. 22 [███████████████████████-----------------] 55% [00m:17s] Processing chr. 22 [████████████████████████----------------] 58% [00m:18s] Processing chr. 22 [█████████████████████████---------------] 61% [00m:19s] Processing chr. 22 [███████████████████████████-------------] 65% [00m:20s] Processing chr. 22 [████████████████████████████------------] 68% [00m:21s] Processing chr. 22 [█████████████████████████████-----------] 71% [00m:22s] Processing chr. 22 [██████████████████████████████----------] 74% [00m:23s] Processing chr. 22 [███████████████████████████████---------] 77% [00m:24s] Processing chr. 22 [████████████████████████████████--------] 78% [00m:25s] Processing chr. 22 [████████████████████████████████--------] 79% [00m:26s] Processing chr. 22 [█████████████████████████████████-------] 81% [00m:27s] Processing chr. 22 [█████████████████████████████████-------] 82% [00m:28s] Processing chr. 22 [██████████████████████████████████------] 83% [00m:29s] Processing chr. 22 [██████████████████████████████████------] 84% [00m:30s] Processing chr. 22 [███████████████████████████████████-----] 85% [00m:31s] Processing chr. 22 [███████████████████████████████████-----] 86% [00m:32s] Processing chr. 22 [███████████████████████████████████-----] 87% [00m:33s] Processing chr. 22 [████████████████████████████████████----] 88% [00m:34s] Processing chr. 22 [█████████████████████████████████████---] 90% [00m:35s] Processing chr. 22 [█████████████████████████████████████---] 91% [00m:36s] Processing chr. 22 [█████████████████████████████████████---] 92% [00m:37s] Processing chr. 22 [██████████████████████████████████████--] 93% [00m:38s] Processing chr. 22 [██████████████████████████████████████--] 94% [00m:39s] Processing chr. 22 [███████████████████████████████████████-] 95% [00m:40s] Processing chr. 22 [███████████████████████████████████████-] 97% [00m:41s] Processing chr. 22 [████████████████████████████████████████] 98% [00m:42s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:43s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:44s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:45s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:46s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:47s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:48s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:49s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:50s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:51s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:52s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:53s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:56s] Reads simulated +#> [█---------------------------------------] 0% [00m:00s] Found 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:01s] Reading 22 [█---------------------------------------] 0% [00m:01s] Processing chr. 22 [█---------------------------------------] 0% [00m:05s] Processing chr. 22 [██--------------------------------------] 4% [00m:06s] Processing chr. 22 [████------------------------------------] 9% [00m:07s] Processing chr. 22 [██████----------------------------------] 14% [00m:08s] Processing chr. 22 [████████--------------------------------] 19% [00m:09s] Processing chr. 22 [██████████------------------------------] 24% [00m:10s] Processing chr. 22 [████████████----------------------------] 29% [00m:11s] Processing chr. 22 [██████████████--------------------------] 34% [00m:12s] Processing chr. 22 [████████████████------------------------] 39% [00m:13s] Processing chr. 22 [██████████████████----------------------] 44% [00m:14s] Processing chr. 22 [████████████████████--------------------] 48% [00m:15s] Processing chr. 22 [█████████████████████-------------------] 51% [00m:16s] Processing chr. 22 [██████████████████████------------------] 54% [00m:17s] Processing chr. 22 [███████████████████████-----------------] 57% [00m:18s] Processing chr. 22 [█████████████████████████---------------] 60% [00m:19s] Processing chr. 22 [██████████████████████████--------------] 63% [00m:20s] Processing chr. 22 [███████████████████████████-------------] 66% [00m:21s] Processing chr. 22 [████████████████████████████------------] 69% [00m:22s] Processing chr. 22 [█████████████████████████████-----------] 72% [00m:23s] Processing chr. 22 [███████████████████████████████---------] 75% [00m:24s] Processing chr. 22 [████████████████████████████████--------] 78% [00m:25s] Processing chr. 22 [████████████████████████████████--------] 79% [00m:26s] Processing chr. 22 [█████████████████████████████████-------] 80% [00m:27s] Processing chr. 22 [█████████████████████████████████-------] 81% [00m:28s] Processing chr. 22 [█████████████████████████████████-------] 82% [00m:29s] Processing chr. 22 [██████████████████████████████████------] 83% [00m:30s] Processing chr. 22 [██████████████████████████████████------] 84% [00m:31s] Processing chr. 22 [███████████████████████████████████-----] 85% [00m:32s] Processing chr. 22 [███████████████████████████████████-----] 86% [00m:33s] Processing chr. 22 [███████████████████████████████████-----] 87% [00m:34s] Processing chr. 22 [████████████████████████████████████----] 88% [00m:35s] Processing chr. 22 [████████████████████████████████████----] 89% [00m:36s] Processing chr. 22 [█████████████████████████████████████---] 90% [00m:37s] Processing chr. 22 [█████████████████████████████████████---] 91% [00m:38s] Processing chr. 22 [██████████████████████████████████████--] 93% [00m:39s] Processing chr. 22 [██████████████████████████████████████--] 94% [00m:40s] Processing chr. 22 [███████████████████████████████████████-] 95% [00m:41s] Processing chr. 22 [███████████████████████████████████████-] 96% [00m:42s] Processing chr. 22 [███████████████████████████████████████-] 97% [00m:43s] Processing chr. 22 [████████████████████████████████████████] 98% [00m:44s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:45s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:46s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:47s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:48s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:49s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:50s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:51s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:52s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:53s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:54s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:55s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:57s] Reads simulated sam_files <- list.files("SAM_0.7/") @@ -409,7 +409,7 @@

Updating the SAM Output Directory# the output directory seq_results <- simulate_seq(phylo_forest, coverage = 2.5, write_SAM = TRUE, update_SAM = TRUE) -#> [█---------------------------------------] 0% [00m:00s] Found 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:01s] Reading 22 [█---------------------------------------] 0% [00m:01s] Processing chr. 22 [█---------------------------------------] 0% [00m:05s] Processing chr. 22 [███-------------------------------------] 5% [00m:06s] Processing chr. 22 [█████-----------------------------------] 10% [00m:07s] Processing chr. 22 [███████---------------------------------] 15% [00m:08s] Processing chr. 22 [█████████-------------------------------] 20% [00m:09s] Processing chr. 22 [███████████-----------------------------] 26% [00m:10s] Processing chr. 22 [█████████████---------------------------] 31% [00m:11s] Processing chr. 22 [███████████████-------------------------] 36% [00m:12s] Processing chr. 22 [█████████████████-----------------------] 41% [00m:13s] Processing chr. 22 [███████████████████---------------------] 47% [00m:14s] Processing chr. 22 [█████████████████████-------------------] 50% [00m:15s] Processing chr. 22 [██████████████████████------------------] 52% [00m:16s] Processing chr. 22 [███████████████████████-----------------] 56% [00m:17s] Processing chr. 22 [█████████████████████████---------------] 60% [00m:18s] Processing chr. 22 [█████████████████████████---------------] 62% [00m:19s] Processing chr. 22 [███████████████████████████-------------] 66% [00m:20s] Processing chr. 22 [████████████████████████████------------] 69% [00m:21s] Processing chr. 22 [█████████████████████████████-----------] 72% [00m:22s] Processing chr. 22 [███████████████████████████████---------] 75% [00m:23s] Processing chr. 22 [████████████████████████████████--------] 78% [00m:24s] Processing chr. 22 [████████████████████████████████--------] 79% [00m:25s] Processing chr. 22 [█████████████████████████████████-------] 80% [00m:26s] Processing chr. 22 [█████████████████████████████████-------] 81% [00m:27s] Processing chr. 22 [█████████████████████████████████-------] 82% [00m:28s] Processing chr. 22 [██████████████████████████████████------] 84% [00m:29s] Processing chr. 22 [███████████████████████████████████-----] 85% [00m:30s] Processing chr. 22 [███████████████████████████████████-----] 86% [00m:31s] Processing chr. 22 [███████████████████████████████████-----] 87% [00m:32s] Processing chr. 22 [████████████████████████████████████----] 88% [00m:33s] Processing chr. 22 [████████████████████████████████████----] 89% [00m:34s] Processing chr. 22 [█████████████████████████████████████---] 90% [00m:35s] Processing chr. 22 [█████████████████████████████████████---] 92% [00m:36s] Processing chr. 22 [██████████████████████████████████████--] 93% [00m:37s] Processing chr. 22 [██████████████████████████████████████--] 94% [00m:38s] Processing chr. 22 [███████████████████████████████████████-] 95% [00m:39s] Processing chr. 22 [███████████████████████████████████████-] 96% [00m:40s] Processing chr. 22 [████████████████████████████████████████] 98% [00m:41s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:42s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:43s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:44s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:45s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:46s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:47s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:48s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:49s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:50s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:51s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:54s] Reads simulated

+#> [█---------------------------------------] 0% [00m:00s] Found 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:01s] Reading 22 [█---------------------------------------] 0% [00m:01s] Processing chr. 22 [█---------------------------------------] 0% [00m:05s] Processing chr. 22 [██--------------------------------------] 4% [00m:06s] Processing chr. 22 [████------------------------------------] 9% [00m:07s] Processing chr. 22 [██████----------------------------------] 14% [00m:08s] Processing chr. 22 [████████--------------------------------] 18% [00m:09s] Processing chr. 22 [██████████------------------------------] 23% [00m:10s] Processing chr. 22 [████████████----------------------------] 28% [00m:11s] Processing chr. 22 [██████████████--------------------------] 33% [00m:12s] Processing chr. 22 [███████████████-------------------------] 37% [00m:13s] Processing chr. 22 [█████████████████-----------------------] 42% [00m:14s] Processing chr. 22 [███████████████████---------------------] 47% [00m:15s] Processing chr. 22 [████████████████████--------------------] 49% [00m:16s] Processing chr. 22 [█████████████████████-------------------] 52% [00m:17s] Processing chr. 22 [███████████████████████-----------------] 55% [00m:18s] Processing chr. 22 [████████████████████████----------------] 58% [00m:19s] Processing chr. 22 [█████████████████████████---------------] 61% [00m:20s] Processing chr. 22 [██████████████████████████--------------] 64% [00m:21s] Processing chr. 22 [███████████████████████████-------------] 67% [00m:22s] Processing chr. 22 [█████████████████████████████-----------] 70% [00m:23s] Processing chr. 22 [██████████████████████████████----------] 73% [00m:24s] Processing chr. 22 [███████████████████████████████---------] 76% [00m:25s] Processing chr. 22 [████████████████████████████████--------] 78% [00m:26s] Processing chr. 22 [████████████████████████████████--------] 79% [00m:27s] Processing chr. 22 [█████████████████████████████████-------] 80% [00m:28s] Processing chr. 22 [█████████████████████████████████-------] 81% [00m:29s] Processing chr. 22 [█████████████████████████████████-------] 82% [00m:30s] Processing chr. 22 [██████████████████████████████████------] 83% [00m:31s] Processing chr. 22 [██████████████████████████████████------] 84% [00m:32s] Processing chr. 22 [███████████████████████████████████-----] 85% [00m:33s] Processing chr. 22 [███████████████████████████████████-----] 86% [00m:34s] Processing chr. 22 [███████████████████████████████████-----] 87% [00m:35s] Processing chr. 22 [████████████████████████████████████----] 88% [00m:36s] Processing chr. 22 [████████████████████████████████████----] 89% [00m:37s] Processing chr. 22 [█████████████████████████████████████---] 90% [00m:38s] Processing chr. 22 [█████████████████████████████████████---] 92% [00m:39s] Processing chr. 22 [██████████████████████████████████████--] 93% [00m:40s] Processing chr. 22 [██████████████████████████████████████--] 94% [00m:41s] Processing chr. 22 [███████████████████████████████████████-] 95% [00m:42s] Processing chr. 22 [███████████████████████████████████████-] 96% [00m:43s] Processing chr. 22 [███████████████████████████████████████-] 97% [00m:44s] Processing chr. 22 [████████████████████████████████████████] 98% [00m:45s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:46s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:47s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:48s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:49s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:50s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:51s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:52s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:53s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:54s] Processing chr. 22 [████████████████████████████████████████] 99% [00m:55s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:58s] Reads simulated

When the output directory already contains a SAM file associated with a chromosome, the alignments on that chromosome are saved in a file whose name is the first available, having the format diff --git a/articles/time_series.html b/articles/time_series.html index c4c6ac4..23f18a5 100644 --- a/articles/time_series.html +++ b/articles/time_series.html @@ -209,7 +209,7 @@

Custom Time-Seriessim$mutate_progeny(sim$choose_cell_in("A"), "D") sim$run_up_to_size("D+", 1000) -#> [███████████████████---------------------] 47% [00m:00s] Cells: 35337 [████████████████████████████████████████] 100% [00m:00s] Saving snapshot +#> [██████████████████----------------------] 44% [00m:00s] Cells: 33978 [████████████████████████████████████████] 100% [00m:00s] Saving snapshot

The time-series can be plot using plot_timeseries().

 # Time-series plot
@@ -295,7 +295,7 @@ 

Time-Varying Evolutionary Rates # Now D will become larger sim$run_up_to_size("D+", 6000) -#> [███████████████████████████-------------] 66% [00m:00s] Cells: 49181 [████████████████████████████████████████] 100% [00m:00s] Saving snapshot +#> [██████████████████████████--------------] 64% [00m:00s] Cells: 47802 [████████████████████████████████████████] 100% [00m:00s] Saving snapshot # Current state sim diff --git a/pkgdown.yml b/pkgdown.yml index 231a99c..0663205 100644 --- a/pkgdown.yml +++ b/pkgdown.yml @@ -15,7 +15,7 @@ articles: three_sampling: three_sampling.html time_series: time_series.html tissue_simulation: tissue_simulation.html -last_built: 2024-11-01T16:48Z +last_built: 2024-11-14T09:03Z urls: reference: https://caravagnalab.github.io/rRACES/reference article: https://caravagnalab.github.io/rRACES/articles diff --git a/reference/MutationEngine-cash-add_exposure.html b/reference/MutationEngine-cash-add_exposure.html index 1191fe6..3e0c9ff 100644 --- a/reference/MutationEngine-cash-add_exposure.html +++ b/reference/MutationEngine-cash-add_exposure.html @@ -91,15 +91,15 @@

Examples#> Downloading germline mutations... #> Germline mutations downloaded #> Building context index... -#> [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█████████████████-----------------------] 40% [00m:01s] Processing chr. 22 [█████████████████████████████████-------] 81% [00m:02s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:02s] Context index built -#> [█---------------------------------------] 0% [00m:00s] Saving context index [████████████████████████████████████████] 100% [00m:00s] Context index saved +#> [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█████████████████-----------------------] 40% [00m:01s] Processing chr. 22 [█████████████████████████████████-------] 81% [00m:02s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:02s] Context index built +#> [█---------------------------------------] 0% [00m:00s] Saving context index [████████████████████████████████████████] 100% [00m:00s] Context index saved #> done #> Building repeated sequence index... -#> [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█---------------------------------------] 0% [00m:03s] Processing chr. 22 [█---------------------------------------] 0% [00m:04s] Processing chr. 22 [█---------------------------------------] 0% [00m:06s] Processing chr. 22 [█---------------------------------------] 0% [00m:07s] Processing chr. 22 [█---------------------------------------] 0% [00m:08s] Processing chr. 22 [█---------------------------------------] 0% [00m:12s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:12s] RS index built -#> [█---------------------------------------] 0% [00m:00s] Saving RS index [████████████████████████----------------] 58% [00m:01s] Saving RS index done -#> [████████████████████████████████████████] 100% [00m:01s] RS index saved -#> [█---------------------------------------] 0% [00m:00s] Loading germline [████████████████████████████████████████] 100% [00m:00s] Germline loaded -#> [█---------------------------------------] 0% [00m:00s] Saving germline [████████████████████████████████████████] 100% [00m:00s] Germline saved +#> [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█---------------------------------------] 0% [00m:00s] Processing chr. 22 [█---------------------------------------] 0% [00m:03s] Processing chr. 22 [█---------------------------------------] 0% [00m:04s] Processing chr. 22 [█---------------------------------------] 0% [00m:06s] Processing chr. 22 [█---------------------------------------] 0% [00m:07s] Processing chr. 22 [█---------------------------------------] 0% [00m:08s] Processing chr. 22 [█---------------------------------------] 0% [00m:12s] Processing chr. 22 [████████████████████████████████████████] 100% [00m:12s] RS index built +#> [█---------------------------------------] 0% [00m:00s] Saving RS index [█---------------------------------------] 0% [00m:01s] Saving RS index [███████████████████████████-------------] 65% [00m:02s] Saving RS index done +#> [████████████████████████████████████████] 100% [00m:02s] RS index saved +#> [█---------------------------------------] 0% [00m:00s] Loading germline [████████████████████████████████████████] 100% [00m:00s] Germline loaded +#> [█---------------------------------------] 0% [00m:00s] Saving germline [████████████████████████████████████████] 100% [00m:00s] Germline saved # add a default set of coefficients that will be used from simulated # time 0 up to the successive coefficient change. The indel and SNV diff --git a/reference/MutationEngine-cash-add_mutant.html b/reference/MutationEngine-cash-add_mutant.html index c0471a8..0c2ccde 100644 --- a/reference/MutationEngine-cash-add_mutant.html +++ b/reference/MutationEngine-cash-add_mutant.html @@ -79,9 +79,9 @@

DetailsExamples

# create a demostrative mutation engine
 m_engine <- MutationEngine(setup_code = "demo")
-#> 
 [█---------------------------------------] 0% [00m:00s] Loading context index                                

 [████████████████████████████████████████] 100% [00m:00s] Context index loaded                               

-#> 
 [█---------------------------------------] 0% [00m:00s] Loading RS index                                     

 [████████████████------------------------] 38% [00m:01s] Loading RS index                                    

 [███████████████████████████████---------] 75% [00m:02s] Loading RS index                                    

 [████████████████████████████████████████] 100% [00m:02s] RS index loaded                                    

-#> 
 [█---------------------------------------] 0% [00m:00s] Loading germline                                     

 [████████████████████████████████████████] 100% [00m:00s] Germline loaded                                    

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading context index                                                                                                                              

 [████████████████████████████████████████] 100% [00m:00s] Context index loaded                                                                                                                             

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading RS index                                                                                                                                   

 [████████████████------------------------] 39% [00m:01s] Loading RS index                                                                                                                                  

 [███████████████████████████████---------] 77% [00m:02s] Loading RS index                                                                                                                                  

 [████████████████████████████████████████] 100% [00m:02s] RS index loaded                                                                                                                                  

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading germline                                                                                                                                   

 [████████████████████████████████████████] 100% [00m:00s] Germline loaded                                                                                                                                  

 
 # add the mutant "A" characterized by one driver SNV on chromosome 22, an
 # indel on the same chromosome, a whole genome doubling event, and finally
@@ -98,7 +98,7 @@ 

Examples chr_pos = 10303470, len = 200000), CNA("D", "22", 5010000, 200000))) -#> [█---------------------------------------] 0% [00m:00s] Retrieving "A" SNVs [█---------------------------------------] 0% [00m:00s] Found 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:01s] Reading 22 [████████████████████████████████████████] 100% [00m:01s] "A" SNVs retrieved +#> [█---------------------------------------] 0% [00m:00s] Retrieving "A" SNVs [█---------------------------------------] 0% [00m:00s] Found 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:01s] Reading 22 [████████████████████████████████████████] 100% [00m:01s] "A" SNVs retrieved m_engine #> MutationEngine diff --git a/reference/MutationEngine-cash-get_SNV_signatures.html b/reference/MutationEngine-cash-get_SNV_signatures.html index bd780d1..f2c94d1 100644 --- a/reference/MutationEngine-cash-get_SNV_signatures.html +++ b/reference/MutationEngine-cash-get_SNV_signatures.html @@ -69,9 +69,9 @@

DetailsExamples

# build a mutation engine
 m_engine <- MutationEngine(setup_code = "demo")
-#> 
 [█---------------------------------------] 0% [00m:00s] Loading context index                                

 [████████████████████████████████████████] 100% [00m:00s] Context index loaded                               

-#> 
 [█---------------------------------------] 0% [00m:00s] Loading RS index                                     

 [████████████████------------------------] 39% [00m:01s] Loading RS index                                    

 [███████████████████████████████---------] 77% [00m:02s] Loading RS index                                    

 [████████████████████████████████████████] 100% [00m:02s] RS index loaded                                    

-#> 
 [█---------------------------------------] 0% [00m:00s] Loading germline                                     

 [████████████████████████████████████████] 100% [00m:00s] Germline loaded                                    

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading context index                                                                                                                              

 [████████████████████████████████████████] 100% [00m:00s] Context index loaded                                                                                                                             

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading RS index                                                                                                                                   

 [█████████████████-----------------------] 40% [00m:01s] Loading RS index                                                                                                                                  

 [████████████████████████████████--------] 78% [00m:02s] Loading RS index                                                                                                                                  

 [████████████████████████████████████████] 100% [00m:02s] RS index loaded                                                                                                                                  

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading germline                                                                                                                                   

 [████████████████████████████████████████] 100% [00m:00s] Germline loaded                                                                                                                                  

 
 # get the indel dataframe
 head(m_engine$get_SNV_signatures(), 5)
diff --git a/reference/MutationEngine-cash-get_active_germline.html b/reference/MutationEngine-cash-get_active_germline.html
index 39c94b5..0ad10b6 100644
--- a/reference/MutationEngine-cash-get_active_germline.html
+++ b/reference/MutationEngine-cash-get_active_germline.html
@@ -69,9 +69,9 @@ 

See alsoExamples

# build a mutation engine
 m_engine <- MutationEngine(setup_code = "demo")
-#> 
 [█---------------------------------------] 0% [00m:00s] Loading context index                                

 [████████████████████████████████████████] 100% [00m:00s] Context index loaded                               

-#> 
 [█---------------------------------------] 0% [00m:00s] Loading RS index                                     

 [████████████████------------------------] 39% [00m:01s] Loading RS index                                    

 [███████████████████████████████---------] 77% [00m:02s] Loading RS index                                    

 [████████████████████████████████████████] 100% [00m:02s] RS index loaded                                    

-#> 
 [█---------------------------------------] 0% [00m:00s] Loading germline                                     

 [████████████████████████████████████████] 100% [00m:00s] Germline loaded                                    

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading context index                                                                                                                              

 [████████████████████████████████████████] 100% [00m:00s] Context index loaded                                                                                                                             

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading RS index                                                                                                                                   

 [█████████████████-----------------------] 40% [00m:01s] Loading RS index                                                                                                                                  

 [████████████████████████████████--------] 78% [00m:02s] Loading RS index                                                                                                                                  

 [████████████████████████████████████████] 100% [00m:02s] RS index loaded                                                                                                                                  

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading germline                                                                                                                                   

 [████████████████████████████████████████] 100% [00m:00s] Germline loaded                                                                                                                                  

 
 # get the active germline subject dataframe
 head(m_engine$get_active_germline(), 5)
diff --git a/reference/MutationEngine-cash-get_germline_subjects.html b/reference/MutationEngine-cash-get_germline_subjects.html
index 0b0ae57..ec81974 100644
--- a/reference/MutationEngine-cash-get_germline_subjects.html
+++ b/reference/MutationEngine-cash-get_germline_subjects.html
@@ -69,9 +69,9 @@ 

See alsoExamples

# build a mutation engine
 m_engine <- MutationEngine(setup_code = "demo")
-#> 
 [█---------------------------------------] 0% [00m:00s] Loading context index                                

 [████████████████████████████████████████] 100% [00m:00s] Context index loaded                               

-#> 
 [█---------------------------------------] 0% [00m:00s] Loading RS index                                     

 [████████████████------------------------] 39% [00m:01s] Loading RS index                                    

 [███████████████████████████████---------] 76% [00m:02s] Loading RS index                                    

 [████████████████████████████████████████] 100% [00m:02s] RS index loaded                                    

-#> 
 [█---------------------------------------] 0% [00m:00s] Loading germline                                     

 [████████████████████████████████████████] 100% [00m:00s] Germline loaded                                    

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading context index                                                                                                                              

 [████████████████████████████████████████] 100% [00m:00s] Context index loaded                                                                                                                             

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading RS index                                                                                                                                   

 [█████████████████-----------------------] 40% [00m:01s] Loading RS index                                                                                                                                  

 [████████████████████████████████--------] 78% [00m:02s] Loading RS index                                                                                                                                  

 [████████████████████████████████████████] 100% [00m:02s] RS index loaded                                                                                                                                  

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading germline                                                                                                                                   

 [████████████████████████████████████████] 100% [00m:00s] Germline loaded                                                                                                                                  

 
 # get the active germline subject dataframe
 head(m_engine$get_germline_subjects(), 5)
diff --git a/reference/MutationEngine-cash-get_indel_signatures.html b/reference/MutationEngine-cash-get_indel_signatures.html
index a8ab0a2..c3aa532 100644
--- a/reference/MutationEngine-cash-get_indel_signatures.html
+++ b/reference/MutationEngine-cash-get_indel_signatures.html
@@ -68,9 +68,9 @@ 

DetailsExamples

# build a mutation engine
 m_engine <- MutationEngine(setup_code = "demo")
-#> 
 [█---------------------------------------] 0% [00m:00s] Loading context index                                

 [████████████████████████████████████████] 100% [00m:00s] Context index loaded                               

-#> 
 [█---------------------------------------] 0% [00m:00s] Loading RS index                                     

 [████████████████------------------------] 39% [00m:01s] Loading RS index                                    

 [███████████████████████████████---------] 76% [00m:02s] Loading RS index                                    

 [████████████████████████████████████████] 100% [00m:02s] RS index loaded                                    

-#> 
 [█---------------------------------------] 0% [00m:00s] Loading germline                                     

 [████████████████████████████████████████] 100% [00m:00s] Germline loaded                                    

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading context index                                                                                                                              

 [████████████████████████████████████████] 100% [00m:00s] Context index loaded                                                                                                                             

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading RS index                                                                                                                                   

 [████████████████------------------------] 39% [00m:01s] Loading RS index                                                                                                                                  

 [█████████████████████████████-----------] 72% [00m:02s] Loading RS index                                                                                                                                  

 [████████████████████████████████████████] 100% [00m:02s] RS index loaded                                                                                                                                  

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading germline                                                                                                                                   

 [████████████████████████████████████████] 100% [00m:00s] Germline loaded                                                                                                                                  

 
 # get the indel dataframe
 head(m_engine$get_indel_signatures(), 5)
diff --git a/reference/MutationEngine-cash-get_known_drivers.html b/reference/MutationEngine-cash-get_known_drivers.html
index 2eddaef..5abe46e 100644
--- a/reference/MutationEngine-cash-get_known_drivers.html
+++ b/reference/MutationEngine-cash-get_known_drivers.html
@@ -79,9 +79,9 @@ 

See alsoExamples

# build a mutation engine
 m_engine <- MutationEngine(setup_code = "demo")
-#> 
 [█---------------------------------------] 0% [00m:00s] Loading context index                                

 [████████████████████████████████████████] 100% [00m:00s] Context index loaded                               

-#> 
 [█---------------------------------------] 0% [00m:00s] Loading RS index                                     

 [████████████████------------------------] 39% [00m:01s] Loading RS index                                    

 [███████████████████████████████---------] 77% [00m:02s] Loading RS index                                    

 [████████████████████████████████████████] 100% [00m:02s] RS index loaded                                    

-#> 
 [█---------------------------------------] 0% [00m:00s] Loading germline                                     

 [████████████████████████████████████████] 100% [00m:00s] Germline loaded                                    

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading context index                                                                                                                              

 [████████████████████████████████████████] 100% [00m:00s] Context index loaded                                                                                                                             

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading RS index                                                                                                                                   

 [████████████████------------------------] 38% [00m:01s] Loading RS index                                                                                                                                  

 [██████████████████████████████----------] 73% [00m:02s] Loading RS index                                                                                                                                  

 [████████████████████████████████████████] 100% [00m:02s] RS index loaded                                                                                                                                  

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading germline                                                                                                                                   

 [████████████████████████████████████████] 100% [00m:00s] Germline loaded                                                                                                                                  

 
 # get the known driver dataframe
 head(m_engine$get_known_drivers(), 5)
diff --git a/reference/MutationEngine-cash-get_population_descritions.html b/reference/MutationEngine-cash-get_population_descritions.html
index ac862e1..ce9e6f4 100644
--- a/reference/MutationEngine-cash-get_population_descritions.html
+++ b/reference/MutationEngine-cash-get_population_descritions.html
@@ -63,9 +63,9 @@ 

DetailsExamples

# build a mutation engine
 m_engine <- MutationEngine(setup_code = "demo")
-#> 
 [█---------------------------------------] 0% [00m:00s] Loading context index                                

 [████████████████████████████████████████] 100% [00m:00s] Context index loaded                               

-#> 
 [█---------------------------------------] 0% [00m:00s] Loading RS index                                     

 [████████████████------------------------] 39% [00m:01s] Loading RS index                                    

 [███████████████████████████████---------] 75% [00m:02s] Loading RS index                                    

 [████████████████████████████████████████] 100% [00m:02s] RS index loaded                                    

-#> 
 [█---------------------------------------] 0% [00m:00s] Loading germline                                     

 [████████████████████████████████████████] 100% [00m:00s] Germline loaded                                    

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading context index                                                                                                                              

 [████████████████████████████████████████] 100% [00m:00s] Context index loaded                                                                                                                             

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading RS index                                                                                                                                   

 [█████████████████-----------------------] 40% [00m:01s] Loading RS index                                                                                                                                  

 [████████████████████████████████--------] 78% [00m:02s] Loading RS index                                                                                                                                  

 [████████████████████████████████████████] 100% [00m:02s] RS index loaded                                                                                                                                  

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading germline                                                                                                                                   

 [████████████████████████████████████████] 100% [00m:00s] Germline loaded                                                                                                                                  

 
 # get the active germline subject dataframe
 head(m_engine$get_population_descritions(), 5)
diff --git a/reference/MutationEngine-cash-infinite_sites_model.html b/reference/MutationEngine-cash-infinite_sites_model.html
index 2b66c06..12d0c7b 100644
--- a/reference/MutationEngine-cash-infinite_sites_model.html
+++ b/reference/MutationEngine-cash-infinite_sites_model.html
@@ -56,9 +56,9 @@ 

DetailsExamples

# create a demostrative mutation engine
 m_engine <- MutationEngine(setup_code = "demo")
-#> 
 [█---------------------------------------] 0% [00m:00s] Loading context index                                

 [████████████████████████████████████████] 100% [00m:00s] Context index loaded                               

-#> 
 [█---------------------------------------] 0% [00m:00s] Loading RS index                                     

 [████████████████------------------------] 39% [00m:01s] Loading RS index                                    

 [███████████████████████████████---------] 77% [00m:02s] Loading RS index                                    

 [████████████████████████████████████████] 100% [00m:02s] RS index loaded                                    

-#> 
 [█---------------------------------------] 0% [00m:00s] Loading germline                                     

 [████████████████████████████████████████] 100% [00m:00s] Germline loaded                                    

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading context index                                                                                                                              

 [████████████████████████████████████████] 100% [00m:00s] Context index loaded                                                                                                                             

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading RS index                                                                                                                                   

 [█████████████████-----------------------] 40% [00m:01s] Loading RS index                                                                                                                                  

 [████████████████████████████████--------] 78% [00m:02s] Loading RS index                                                                                                                                  

 [████████████████████████████████████████] 100% [00m:02s] RS index loaded                                                                                                                                  

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading germline                                                                                                                                   

 [████████████████████████████████████████] 100% [00m:00s] Germline loaded                                                                                                                                  

 
 # the infinite sites model is enabled by default
 m_engine$infinite_sites_model
diff --git a/reference/MutationEngine-cash-place_mutations.html b/reference/MutationEngine-cash-place_mutations.html
index 2807b6c..5bf795b 100644
--- a/reference/MutationEngine-cash-place_mutations.html
+++ b/reference/MutationEngine-cash-place_mutations.html
@@ -100,7 +100,7 @@ 

Examples sim$death_activation_level <- 100 sim$run_up_to_size(species = "A", num_of_cells = 50000) -#> [████████████████████████████████████████] 100% [00m:00s] Saving snapshot +#> [████████████████████████████████████████] 100% [00m:00s] Saving snapshot # sample the region [450,500]x[475,550] sim$sample_cells("S1", lower_corner = c(450, 475), @@ -111,13 +111,13 @@

Examples # build a mutation engine m_engine <- MutationEngine(setup_code = "demo") -#> [█---------------------------------------] 0% [00m:00s] Loading context index [████████████████████████████████████████] 100% [00m:00s] Context index loaded -#> [█---------------------------------------] 0% [00m:00s] Loading RS index [████████████████------------------------] 39% [00m:01s] Loading RS index [███████████████████████████████---------] 77% [00m:02s] Loading RS index [████████████████████████████████████████] 100% [00m:02s] RS index loaded -#> [█---------------------------------------] 0% [00m:00s] Loading germline [████████████████████████████████████████] 100% [00m:00s] Germline loaded +#> [█---------------------------------------] 0% [00m:00s] Loading context index [████████████████████████████████████████] 100% [00m:00s] Context index loaded +#> [█---------------------------------------] 0% [00m:00s] Loading RS index [█████████████████-----------------------] 40% [00m:01s] Loading RS index [████████████████████████████████--------] 78% [00m:02s] Loading RS index [████████████████████████████████████████] 100% [00m:02s] RS index loaded +#> [█---------------------------------------] 0% [00m:00s] Loading germline [████████████████████████████████████████] 100% [00m:00s] Germline loaded # add the mutant "A" to the engine m_engine$add_mutant("A", c(SNV = 3e-9), list(SNV("22", 12028576, "G"))) -#> [█---------------------------------------] 0% [00m:00s] Retrieving "A" SNVs [█---------------------------------------] 0% [00m:00s] Found 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:01s] Reading 22 [████████████████████████████████████████] 100% [00m:01s] "A" SNVs retrieved +#> [█---------------------------------------] 0% [00m:00s] Retrieving "A" SNVs [█---------------------------------------] 0% [00m:00s] Found 22 [█---------------------------------------] 0% [00m:00s] Reading 22 [█---------------------------------------] 0% [00m:01s] Reading 22 [████████████████████████████████████████] 100% [00m:01s] "A" SNVs retrieved # add the default set of SNV signature coefficients m_engine$add_exposure(c(SBS13 = 0.3, SBS1 = 0.7, ID2 = 0.3, ID21 = 0.5, @@ -126,7 +126,7 @@

Examples# place the mutations on the samples forest assuming 1000 pre-neoplastic # SNVs and 500 indels phylogenetic_forest <- m_engine$place_mutations(samples_forest, 1000, 500) -#> [█---------------------------------------] 0% [00m:00s] Placing mutations [████████████████████████████████████████] 100% [00m:00s] Mutations placed +#> [█---------------------------------------] 0% [00m:00s] Placing mutations [████████████████████████████████████████] 100% [00m:00s] Mutations placed phylogenetic_forest #> PhylogeneticForest diff --git a/reference/MutationEngine-cash-set_germline_subject.html b/reference/MutationEngine-cash-set_germline_subject.html index c8835e5..e807d76 100644 --- a/reference/MutationEngine-cash-set_germline_subject.html +++ b/reference/MutationEngine-cash-set_germline_subject.html @@ -66,13 +66,13 @@

See alsoExamples

# build a mutation engine
 m_engine <- MutationEngine(setup_code = "demo")
-#> 
 [█---------------------------------------] 0% [00m:00s] Loading context index                                

 [████████████████████████████████████████] 100% [00m:00s] Context index loaded                               

-#> 
 [█---------------------------------------] 0% [00m:00s] Loading RS index                                     

 [████████████████------------------------] 39% [00m:01s] Loading RS index                                    

 [█████████████████████████████-----------] 70% [00m:02s] Loading RS index                                    

 [████████████████████████████████████████] 100% [00m:02s] RS index loaded                                    

-#> 
 [█---------------------------------------] 0% [00m:00s] Loading germline                                     

 [████████████████████████████████████████] 100% [00m:00s] Germline loaded                                    

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading context index                                                                                                                              

 [████████████████████████████████████████] 100% [00m:00s] Context index loaded                                                                                                                             

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading RS index                                                                                                                                   

 [█████████████████-----------------------] 40% [00m:01s] Loading RS index                                                                                                                                  

 [████████████████████████████████--------] 78% [00m:02s] Loading RS index                                                                                                                                  

 [████████████████████████████████████████] 100% [00m:02s] RS index loaded                                                                                                                                  

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading germline                                                                                                                                   

 [████████████████████████████████████████] 100% [00m:00s] Germline loaded                                                                                                                                  

 
 # set the active germline subject dataframe
 m_engine$set_germline_subject("NA18941")
-#> 
 [█---------------------------------------] 0% [00m:00s] Loading germline                                     

 [████████████████████████████████████████] 100% [00m:00s] Germline loaded                                    

+#> 
 [█---------------------------------------] 0% [00m:00s] Loading germline                                                                                                                                   

 [████████████████████████████████████████] 100% [00m:00s] Germline loaded