partis-bcr 1.0.11.post1.dev3__cp312-cp312-manylinux1_x86_64.whl → 1.0.12__cp312-cp312-manylinux1_x86_64.whl
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- bin/partis-test.py +3 -3
- {partis_bcr-1.0.11.post1.dev3.data → partis_bcr-1.0.12.data}/scripts/partis-test.py +3 -3
- {partis_bcr-1.0.11.post1.dev3.dist-info → partis_bcr-1.0.12.dist-info}/METADATA +1 -1
- {partis_bcr-1.0.11.post1.dev3.dist-info → partis_bcr-1.0.12.dist-info}/RECORD +55 -54
- test/new-results/run-times.csv +5 -5
- test/new-results/test/simu.yaml +1 -1
- test/new-results/test.log +26 -26
- test/paired/new-results/partition-new-simu/fasttree/iclust-2/log +1 -1
- test/paired/new-results/partition-new-simu/igh+igk/partition-igh/fasttree/iclust-1/log +2 -2
- test/paired/new-results/partition-new-simu-annotation-performance/plots/flcount-matrix.svg +107 -107
- test/paired/new-results/partition-new-simu-annotation-performance/plots/func-non-func-per-drop-log.svg +83 -83
- test/paired/new-results/partition-new-simu-annotation-performance/plots/func-non-func-per-drop.svg +61 -61
- test/paired/new-results/partition-new-simu-annotation-performance/plots/paired-seqs-per-seq-after.svg +18 -18
- test/paired/new-results/partition-new-simu-annotation-performance/plots/paired-seqs-per-seq-before.svg +24 -24
- test/paired/new-results/partition-new-simu-annotation-performance/plots/pseq-matrix.svg +35 -35
- test/paired/new-results/partition-new-simu-annotation-performance/plots/seqs-per-droplet.svg +22 -22
- test/paired/new-results/partition-new-simu-annotation-performance/plots/true-pair-clean-performance-correct-family.svg +45 -45
- test/paired/new-results/partition-new-simu-annotation-performance/plots/true-pair-clean-performance-correct.svg +45 -45
- test/paired/new-results/partition-new-simu-annotation-performance/plots/true-pair-clean-performance-mispaired.svg +45 -45
- test/paired/new-results/partition-new-simu-annotation-performance/plots/true-pair-clean-performance-multiple.svg +15 -15
- test/paired/new-results/partition-new-simu-annotation-performance/plots/true-pair-clean-performance-near-family.svg +15 -15
- test/paired/new-results/partition-new-simu-annotation-performance/plots/true-pair-clean-performance-total.svg +47 -47
- test/paired/new-results/partition-new-simu-annotation-performance/plots/true-pair-clean-performance-unpaired.svg +48 -48
- test/paired/new-results/partition-new-simu-annotation-performance/plots/true-pair-clean-performance.svg +27 -27
- test/paired/new-results/run-times.csv +6 -6
- test/paired/new-results/subset-partition-new-simu/isub-0/partition.log +9 -5
- test/paired/new-results/subset-partition-new-simu/isub-1/partition.log +6 -2
- test/paired/new-results/subset-partition-new-simu/merged-partition.log +7 -3
- test/paired/new-results/subset-partition-new-simu/merged-subsets/parameters/igh/sw-cache.yaml +1 -1
- test/paired/new-results/subset-partition-new-simu/merged-subsets/parameters/igk/sw-cache.yaml +1 -1
- test/paired/new-results/subset-partition-new-simu/merged-subsets/parameters/igl/sw-cache.yaml +1 -1
- test/paired/new-results/test/simu/all-seqs.fa +168 -212
- test/paired/new-results/test/simu/igh+igk/igh.yaml +1 -1
- test/paired/new-results/test/simu/igh+igk/igk.yaml +1 -1
- test/paired/new-results/test/simu/igh+igl/igh.yaml +1 -1
- test/paired/new-results/test/simu/igh+igl/igl.yaml +1 -1
- test/paired/new-results/test/simu/igh.yaml +1 -1
- test/paired/new-results/test/simu/igk.yaml +1 -1
- test/paired/new-results/test/simu/igl.yaml +1 -1
- test/paired/new-results/test/simu/meta.yaml +1 -1
- test/paired/new-results/test.log +115 -118
- test/test.py +3 -3
- test/trees.nwk +100 -0
- {partis_bcr-1.0.11.post1.dev3.data → partis_bcr-1.0.12.data}/scripts/cf-alleles.py +0 -0
- {partis_bcr-1.0.11.post1.dev3.data → partis_bcr-1.0.12.data}/scripts/cf-germlines.py +0 -0
- {partis_bcr-1.0.11.post1.dev3.data → partis_bcr-1.0.12.data}/scripts/compare-plotdirs.py +0 -0
- {partis_bcr-1.0.11.post1.dev3.data → partis_bcr-1.0.12.data}/scripts/extract-pairing-info.py +0 -0
- {partis_bcr-1.0.11.post1.dev3.data → partis_bcr-1.0.12.data}/scripts/gctree-run.py +0 -0
- {partis_bcr-1.0.11.post1.dev3.data → partis_bcr-1.0.12.data}/scripts/get-naive-probabilities.py +0 -0
- {partis_bcr-1.0.11.post1.dev3.data → partis_bcr-1.0.12.data}/scripts/parse-output.py +0 -0
- {partis_bcr-1.0.11.post1.dev3.data → partis_bcr-1.0.12.data}/scripts/split-loci.py +0 -0
- {partis_bcr-1.0.11.post1.dev3.dist-info → partis_bcr-1.0.12.dist-info}/WHEEL +0 -0
- {partis_bcr-1.0.11.post1.dev3.dist-info → partis_bcr-1.0.12.dist-info}/entry_points.txt +0 -0
- {partis_bcr-1.0.11.post1.dev3.dist-info → partis_bcr-1.0.12.dist-info}/licenses/COPYING +0 -0
- {partis_bcr-1.0.11.post1.dev3.dist-info → partis_bcr-1.0.12.dist-info}/top_level.txt +0 -0
test/paired/new-results/test.log
CHANGED
@@ -23,17 +23,17 @@ smith-waterman (writing parameters)
|
|
23
23
|
writing parameters to test/paired/new-results/test/parameters/simu/igh/sw (0.6 sec)
|
24
24
|
water time: 0.7 (ig-sw 0.7 processing 0.0)
|
25
25
|
writing hmms removed 16 genes from glfo (leaving v 6 d 18 j 4)
|
26
|
-
(0.
|
26
|
+
(0.6 sec)
|
27
27
|
hmm
|
28
28
|
reading output
|
29
29
|
read 56 hmm output lines with 56 sequences in 56 events (0 failures)
|
30
30
|
writing parameters to test/paired/new-results/test/parameters/simu/igh/hmm (0.5 sec)
|
31
|
-
writing parameters to test/paired/new-results/test/parameters/simu/igh/true (0.
|
31
|
+
writing parameters to test/paired/new-results/test/parameters/simu/igh/true (0.4 sec)
|
32
32
|
infra time: 1.1
|
33
|
-
hmm step time: 2.
|
33
|
+
hmm step time: 2.2
|
34
34
|
writing hmms removed 3 genes from glfo (leaving v 6 d 15 j 4)
|
35
|
-
(0.
|
36
|
-
total time:
|
35
|
+
(0.4 sec)
|
36
|
+
total time: 4.1
|
37
37
|
[44mcache-parameters[0m [94migk[0m:
|
38
38
|
note: adding input metafname from default location (test/paired/ref-results/test/simu/meta.yaml)
|
39
39
|
[91mrun[0m partis cache-parameters --locus igk --dont-write-git-info --parameter-dir test/paired/new-results/test/parameters/simu/igk --is-simu --random-seed 1 --n-procs 10 --infname test/paired/ref-results/test/simu/igk.yaml --sw-cachefname test/paired/new-results/test/parameters/simu/igk/sw-cache.yaml --input-metafnames test/paired/ref-results/test/simu/meta.yaml
|
@@ -91,7 +91,7 @@ smith-waterman (writing parameters)
|
|
91
91
|
kept 14 (0.667) unproductive
|
92
92
|
writing sw results to test/paired/new-results/test/parameters/simu/igl/sw-cache.yaml
|
93
93
|
writing parameters to test/paired/new-results/test/parameters/simu/igl/sw (0.3 sec)
|
94
|
-
water time: 0.4 (ig-sw 0.
|
94
|
+
water time: 0.4 (ig-sw 0.3 processing 0.0)
|
95
95
|
writing hmms removed 7 genes from glfo (leaving v 3 d 1 j 2)
|
96
96
|
(0.1 sec)
|
97
97
|
hmm
|
@@ -103,8 +103,8 @@ hmm
|
|
103
103
|
hmm step time: 0.9
|
104
104
|
writing hmms removed 0 genes from glfo (leaving v 3 d 1 j 2)
|
105
105
|
(0.1 sec)
|
106
|
-
total time: 1.
|
107
|
-
total time: 12.
|
106
|
+
total time: 1.6
|
107
|
+
total time: 12.6
|
108
108
|
[92mpartition-new-simu[0m partis partition --dont-write-git-info --paired-loci --paired-indir test/paired/ref-results/test/simu --parameter-dir test/paired/new-results/test/parameters/simu --plot-annotation-performance --max-ccf-fail-frac 0.10 --is-simu --plotdir test/paired/new-results/partition-new-simu-annotation-performance --only-csv-plots --no-partition-plots --random-seed 1 --n-procs 10 --paired-outdir test/paired/new-results/partition-new-simu
|
109
109
|
[93mwarning[0m ignoring --plot-annotation-performance for paired clustering since it's going to be a bit fiddly to implement
|
110
110
|
[44mpartition[0m [94migh[0m:
|
@@ -122,23 +122,23 @@ smith-waterman
|
|
122
122
|
(0.0 sec)
|
123
123
|
water time: 0.1
|
124
124
|
hmm
|
125
|
-
caching all 56 naive sequences (1.
|
125
|
+
caching all 56 naive sequences (1.1s)
|
126
126
|
collapsed 56 queries into 50 clusters with identical naive seqs (0.0 sec)
|
127
|
-
50 clusters with 10 procs (0.
|
127
|
+
50 clusters with 10 procs (0.6s)
|
128
128
|
44 clusters with 7 procs (0.5s)
|
129
|
-
36 clusters with 5 procs (0.
|
129
|
+
36 clusters with 5 procs (0.5s)
|
130
130
|
31 clusters with 3 procs (0.3s)
|
131
131
|
29 clusters with 2 procs (0.5s)
|
132
132
|
23 clusters with 1 proc (0.6s)
|
133
133
|
[93mnote[0m not merging entire cpath history
|
134
|
-
partition loop time: 2.
|
134
|
+
partition loop time: 2.9
|
135
135
|
getting annotations for final partition
|
136
136
|
subcluster annotating 20 clusters with steps: (0.7s) (0.2s) plotting performance
|
137
137
|
[93mwarning[0m skipped annotation performance evaluation on 1 queries with different true and inferred net shm indel lengths: db63559ea9-igh
|
138
138
|
(0.0 sec)
|
139
139
|
|
140
140
|
subcluster annotation time 0.9
|
141
|
-
total time:
|
141
|
+
total time: 5.0
|
142
142
|
[44mpartition[0m [94migk[0m:
|
143
143
|
note: adding input metafname from default location (test/paired/ref-results/test/simu/meta.yaml)
|
144
144
|
[91mrun[0m partis partition --locus igk --dont-write-git-info --parameter-dir test/paired/new-results/test/parameters/simu/igk --plot-annotation-performance --max-ccf-fail-frac 0.10 --is-simu --plotdir test/paired/new-results/partition-new-simu-annotation-performance/single-chain/plots/igk --only-csv-plots --no-partition-plots --random-seed 1 --n-procs 10 --refuse-to-cache-parameters --infname test/paired/ref-results/test/simu/igk.yaml --outfname test/paired/new-results/partition-new-simu/single-chain/partition-igk.yaml --sw-cachefname test/paired/new-results/test/parameters/simu/igk/sw-cache.yaml --input-metafnames test/paired/ref-results/test/simu/meta.yaml
|
@@ -228,6 +228,7 @@ WARNING using errors instead of sumw2 with weight != 1.0 in Hist::fill_ibin()
|
|
228
228
|
joint 0.963 0.594
|
229
229
|
merge time 0.0
|
230
230
|
need to get annotations for 7/15 joint igh clusters
|
231
|
+
tag: v1.0.12
|
231
232
|
[44mannotate[0m [94migh+igk[0m:[94m igh[0m
|
232
233
|
note: adding input metafname from default location (test/paired/ref-results/test/simu/meta.yaml)
|
233
234
|
[91mrun[0m partis annotate --locus igh --dont-write-git-info --parameter-dir test/paired/new-results/test/parameters/simu/igh --max-ccf-fail-frac 0.10 --only-csv-plots --no-partition-plots --random-seed 1 --n-procs 10 --refuse-to-cache-parameters --infname test/paired/new-results/partition-new-simu/igh+igk/input-seqs-igh.fa --outfname test/paired/new-results/partition-new-simu/igh+igk/partition-igh.yaml --sw-cachefname test/paired/new-results/test/parameters/simu/igh/sw-cache.yaml --input-metafnames test/paired/new-results/partition-new-simu/igh+igk/meta-igh.yaml --ignore-sw-pair-info --input-partition-fname test/paired/new-results/partition-new-simu/igh+igk/partition-igh-only-partition.yaml
|
@@ -245,6 +246,7 @@ hmm
|
|
245
246
|
subcluster annotation time 0.6
|
246
247
|
total time: 0.6
|
247
248
|
need to get annotations for 13/15 joint igk clusters
|
249
|
+
tag: v1.0.12
|
248
250
|
[44mannotate[0m [94migh+igk[0m:[94m igk[0m
|
249
251
|
note: adding input metafname from default location (test/paired/ref-results/test/simu/meta.yaml)
|
250
252
|
[91mrun[0m partis annotate --locus igk --dont-write-git-info --parameter-dir test/paired/new-results/test/parameters/simu/igk --max-ccf-fail-frac 0.10 --only-csv-plots --no-partition-plots --random-seed 1 --n-procs 10 --refuse-to-cache-parameters --infname test/paired/new-results/partition-new-simu/igh+igk/input-seqs-igk.fa --outfname test/paired/new-results/partition-new-simu/igh+igk/partition-igk.yaml --sw-cachefname test/paired/new-results/test/parameters/simu/igk/sw-cache.yaml --input-metafnames test/paired/new-results/partition-new-simu/igh+igk/meta-igk.yaml --ignore-sw-pair-info --input-partition-fname test/paired/new-results/partition-new-simu/igh+igk/partition-igk-only-partition.yaml
|
@@ -286,6 +288,7 @@ hmm
|
|
286
288
|
joint 1.000 0.611
|
287
289
|
merge time 0.0
|
288
290
|
need to get annotations for 5/8 joint igh clusters
|
291
|
+
tag: v1.0.12
|
289
292
|
[44mannotate[0m [94migh+igl[0m:[94m igh[0m
|
290
293
|
note: adding input metafname from default location (test/paired/ref-results/test/simu/meta.yaml)
|
291
294
|
[91mrun[0m partis annotate --locus igh --dont-write-git-info --parameter-dir test/paired/new-results/test/parameters/simu/igh --max-ccf-fail-frac 0.10 --only-csv-plots --no-partition-plots --random-seed 1 --n-procs 10 --refuse-to-cache-parameters --infname test/paired/new-results/partition-new-simu/igh+igl/input-seqs-igh.fa --outfname test/paired/new-results/partition-new-simu/igh+igl/partition-igh.yaml --sw-cachefname test/paired/new-results/test/parameters/simu/igh/sw-cache.yaml --input-metafnames test/paired/new-results/partition-new-simu/igh+igl/meta-igh.yaml --ignore-sw-pair-info --input-partition-fname test/paired/new-results/partition-new-simu/igh+igl/partition-igh-only-partition.yaml
|
@@ -303,6 +306,7 @@ hmm
|
|
303
306
|
subcluster annotation time 0.4
|
304
307
|
total time: 0.5
|
305
308
|
need to get annotations for 8/8 joint igl clusters
|
309
|
+
tag: v1.0.12
|
306
310
|
[44mannotate[0m [94migh+igl[0m:[94m igl[0m
|
307
311
|
note: adding input metafname from default location (test/paired/ref-results/test/simu/meta.yaml)
|
308
312
|
[91mrun[0m partis annotate --locus igl --dont-write-git-info --parameter-dir test/paired/new-results/test/parameters/simu/igl --max-ccf-fail-frac 0.10 --only-csv-plots --no-partition-plots --random-seed 1 --n-procs 10 --refuse-to-cache-parameters --infname test/paired/new-results/partition-new-simu/igh+igl/input-seqs-igl.fa --outfname test/paired/new-results/partition-new-simu/igh+igl/partition-igl.yaml --sw-cachefname test/paired/new-results/test/parameters/simu/igl/sw-cache.yaml --input-metafnames test/paired/new-results/partition-new-simu/igh+igl/meta-igl.yaml --ignore-sw-pair-info --input-partition-fname test/paired/new-results/partition-new-simu/igh+igl/partition-igl-only-partition.yaml
|
@@ -331,7 +335,7 @@ hmm
|
|
331
335
|
[94mh[0m 20 --> 23 [91m+3[0m
|
332
336
|
[95mk[0m 6 --> 15 [91m+9[0m
|
333
337
|
[92ml[0m 3 --> 8 [91m+5[0m
|
334
|
-
total time:
|
338
|
+
total time: 20.5
|
335
339
|
[92msubset-partition-new-simu[0m partis subset-partition --dont-write-git-info --paired-loci --paired-indir test/paired/ref-results/test/simu --max-ccf-fail-frac 0.15 --is-simu --n-subsets 2 --random-seed 1 --n-procs 10 --paired-outdir test/paired/new-results/subset-partition-new-simu
|
336
340
|
note: --input-metafnames not specified, so setting to default location test/paired/ref-results/test/simu/meta.yaml
|
337
341
|
--input-metafnames: added meta info for 113 sequences from test/paired/ref-results/test/simu/meta.yaml: loci paired-uids
|
@@ -360,9 +364,12 @@ hmm
|
|
360
364
|
writing merged outputs from 2 subsets to test/paired/new-results/subset-partition-new-simu/merged-subsets
|
361
365
|
writing combined h/l info to fasta and meta files in test/paired/new-results/subset-partition-new-simu/merged-subsets
|
362
366
|
merging parameters from 2 subdirs (e.g. test/paired/new-results/subset-partition-new-simu/isub-0) to test/paired/new-results/subset-partition-new-simu/merged-subsets
|
367
|
+
tag: v1.0.12
|
368
|
+
tag: v1.0.12
|
369
|
+
tag: v1.0.12
|
363
370
|
[91mrun[0m partis partition --dont-write-git-info --paired-loci --paired-indir test/paired/ref-results/test/simu --max-ccf-fail-frac 0.15 --is-simu --random-seed 1 --n-procs 10 --paired-outdir test/paired/new-results/subset-partition-new-simu --input-metafnames test/paired/new-results/subset-partition-new-simu/merged-subsets/meta.yaml --input-partition-fname test/paired/new-results/subset-partition-new-simu/merged-subsets --continue-from-input-partition --parameter-dir test/paired/new-results/subset-partition-new-simu/merged-subsets/parameters --ignore-sw-pair-info --refuse-to-cache-parameters --ignore-default-input-metafile
|
364
371
|
log: test/paired/new-results/subset-partition-new-simu/merged-partition.log
|
365
|
-
total time:
|
372
|
+
total time: 69.6
|
366
373
|
[92mseed-partition-new-simu[0m partis partition --dont-write-git-info --paired-loci --paired-indir test/paired/ref-results/test/simu --parameter-dir test/paired/new-results/test/parameters/simu --max-ccf-fail-frac 0.10 --is-simu --seed-unique-id 3542bd38d8-igh:3542bd38d8-igk --seed-loci igh:igk --random-seed 1 --n-procs 10 --paired-outdir test/paired/new-results/seed-partition-new-simu
|
367
374
|
removing seqs very different from seed seq
|
368
375
|
wrong/no removed singly kept missing
|
@@ -423,7 +430,7 @@ hmm
|
|
423
430
|
getting annotations for final partition
|
424
431
|
subcluster annotating 2 clusters with steps: (0.1s) (0.1s) (0.1s)
|
425
432
|
subcluster annotation time 0.2
|
426
|
-
total time: 0.
|
433
|
+
total time: 0.7
|
427
434
|
[44mpartition[0m [94migl[0m:
|
428
435
|
note: adding input metafname from default location (test/paired/ref-results/test/simu/meta.yaml)
|
429
436
|
[91mrun[0m partis partition --locus igl --dont-write-git-info --parameter-dir test/paired/new-results/test/parameters/simu/igl --max-ccf-fail-frac 0.10 --random-seed 1 --n-procs 10 --refuse-to-cache-parameters --infname test/paired/new-results/seed-partition-new-simu/seeds/3542bd38d8-igh+3542bd38d8-igk/single-chain/input-seqs-igl.fa --outfname test/paired/new-results/seed-partition-new-simu/seeds/3542bd38d8-igh+3542bd38d8-igk/single-chain/partition-igl.yaml --sw-cachefname test/paired/new-results/test/parameters/simu/igl/sw-cache.yaml --input-metafnames test/paired/ref-results/test/simu/meta.yaml
|
@@ -485,6 +492,7 @@ WARNING using errors instead of sumw2 with weight != 1.0 in Hist::fill_ibin()
|
|
485
492
|
merge time 0.0
|
486
493
|
removed 7 unseeded clusters (of 8 total, leaving 1) from igh joint partition
|
487
494
|
seed cluster size in best joint partition: 1, excluding seed seq: 0
|
495
|
+
tag: v1.0.12
|
488
496
|
[44mannotate[0m [94migh+igk[0m:[94m igh[0m
|
489
497
|
removed 2 arg strs that were specific to action 'partition': --seed-unique-id 3542bd38d8-igh:3542bd38d8-igk
|
490
498
|
note: adding input metafname from default location (test/paired/ref-results/test/simu/meta.yaml)
|
@@ -506,6 +514,7 @@ hmm
|
|
506
514
|
removed 7 unseeded clusters (of 8 total, leaving 1) from igk joint partition
|
507
515
|
need to get annotations for 1/1 joint igk clusters
|
508
516
|
seed cluster size in best joint partition: 1, excluding seed seq: 0
|
517
|
+
tag: v1.0.12
|
509
518
|
[44mannotate[0m [94migh+igk[0m:[94m igk[0m
|
510
519
|
removed 2 arg strs that were specific to action 'partition': --seed-unique-id 3542bd38d8-igh:3542bd38d8-igk
|
511
520
|
note: adding input metafname from default location (test/paired/ref-results/test/simu/meta.yaml)
|
@@ -552,7 +561,7 @@ hmm
|
|
552
561
|
N clusters with all seqs:
|
553
562
|
[94mh[0m 12 --> 1 [91m-11[0m
|
554
563
|
[95mk[0m 2 --> 1 [91m-1[0m
|
555
|
-
total time: 9.
|
564
|
+
total time: 9.4
|
556
565
|
[92mget-selection-metrics-new-simu[0m partis get-selection-metrics --dont-write-git-info --paired-loci --existing-output-run-cfg paired --min-selection-metric-cluster-size 3 --min-paired-cluster-size-to-read 3 --random-seed 1 --n-procs 10 --paired-outdir test/paired/new-results/partition-new-simu --chosen-ab-fname test/paired/new-results/get-selection-metrics-new-simu-chosen-abs.csv
|
557
566
|
[44mget-selection-metrics[0m [94migh+igk[0m:[94m igh[0m
|
558
567
|
[91mrun[0m partis get-selection-metrics --locus igh --dont-write-git-info --min-selection-metric-cluster-size 3 --min-paired-cluster-size-to-read 3 --random-seed 1 --n-procs 10 --chosen-ab-fname test/paired/new-results/get-selection-metrics-new-simu-chosen-abs.csv --refuse-to-cache-parameters --outfname test/paired/new-results/partition-new-simu/igh+igk/partition-igh.yaml
|
@@ -673,10 +682,10 @@ hmm
|
|
673
682
|
skipped 1/4 clusters that had no uids in common with tree
|
674
683
|
writing selection metrics to test/paired/new-results/partition-new-simu/selection-metrics.yaml
|
675
684
|
writing 10 chosen abs to test/paired/new-results/get-selection-metrics-new-simu-chosen-abs.csv
|
676
|
-
total time: 5.
|
685
|
+
total time: 5.9
|
677
686
|
[92mcache-parameters-data[0m partis cache-parameters --dont-write-git-info --paired-loci --infname test/paired-data/all-seqs.fa --parameter-dir test/paired/new-results/test/parameters/data --n-max-queries 100 --random-seed 1 --n-procs 10
|
678
687
|
note: --paired-outdir is not set, so there will be no persistent record of the results (except the parameter directory).
|
679
|
-
[91mrun[0m extract-pairing-info.py test/paired-data/all-seqs.fa /tmp/runner/hmms/
|
688
|
+
[91mrun[0m extract-pairing-info.py test/paired-data/all-seqs.fa /tmp/runner/hmms/870150/meta.yaml --n-max-queries 100
|
680
689
|
--n-max-queries: chose 48 / 48 droplets (first 100, after random shuffling) which had 100 / 100 seqs
|
681
690
|
extract_pairing_info(): read 100 sequences with 48 droplet ids
|
682
691
|
droplet id separators (set automatically): _ indices: [0]
|
@@ -687,9 +696,9 @@ hmm
|
|
687
696
|
2 42 0.875
|
688
697
|
3 2 0.042
|
689
698
|
4 2 0.042
|
690
|
-
[91mrun[0m split-loci.py test/paired-data/all-seqs.fa --outdir /tmp/runner/hmms/
|
699
|
+
[91mrun[0m split-loci.py test/paired-data/all-seqs.fa --outdir /tmp/runner/hmms/870150 --input-metafname /tmp/runner/hmms/870150/meta.yaml --n-max-queries 100
|
691
700
|
--n-max-queries: chose 48 / 48 droplets (first 100, after random shuffling) which had 100 / 100 seqs
|
692
|
-
--input-metafnames: added meta info for 100 sequences from /tmp/runner/hmms/
|
701
|
+
--input-metafnames: added meta info for 100 sequences from /tmp/runner/hmms/870150/meta.yaml: paired-uids
|
693
702
|
read pairing info for 100 seqs from input meta file
|
694
703
|
running vsearch on 100 sequences:
|
695
704
|
igh: 49 / 100 v annotations (51 failed) with 28 v genes in 0.1 sec
|
@@ -712,10 +721,10 @@ totals: igh 48 igk 37 igl 15
|
|
712
721
|
2 0.13 [94mh[0m [92ml[0m
|
713
722
|
1 0.07 [94mh[0m [95mk[0m
|
714
723
|
1 0.07 [94mh[0m [94mh[0m [95mk[0m
|
715
|
-
writing to /tmp/runner/hmms/
|
716
|
-
igh: 48 to
|
717
|
-
igk: 37 to
|
718
|
-
igl: 15 to
|
724
|
+
writing to /tmp/runner/hmms/870150/
|
725
|
+
igh: 48 to 870150/igh.fa
|
726
|
+
igk: 37 to 870150/igk.fa
|
727
|
+
igl: 15 to 870150/igl.fa
|
719
728
|
writing to paired subdirs
|
720
729
|
igh+igk:
|
721
730
|
igh: 36 / 48 to igh+igk/igh.fa
|
@@ -724,8 +733,8 @@ writing to paired subdirs
|
|
724
733
|
igh: 15 / 48 to igh+igl/igh.fa
|
725
734
|
igl: 15 / 15 to igh+igl/igl.fa
|
726
735
|
[44mcache-parameters[0m [94migh[0m:
|
727
|
-
[91mrun[0m partis cache-parameters --locus igh --dont-write-git-info --infname /tmp/runner/hmms/
|
728
|
-
--input-metafnames: added meta info for 48 sequences from /tmp/runner/hmms/
|
736
|
+
[91mrun[0m partis cache-parameters --locus igh --dont-write-git-info --infname /tmp/runner/hmms/870150/igh.fa --parameter-dir test/paired/new-results/test/parameters/data/igh --random-seed 1 --n-procs 10 --sw-cachefname test/paired/new-results/test/parameters/data/igh/sw-cache.yaml --input-metafnames /tmp/runner/hmms/870150/meta.yaml
|
737
|
+
--input-metafnames: added meta info for 48 sequences from /tmp/runner/hmms/870150/meta.yaml: paired-uids
|
729
738
|
caching parameters
|
730
739
|
vsearch: 48 / 48 v annotations (0 failed) with 27 v genes in 0.0 sec
|
731
740
|
keeping 23 / 198 v genes
|
@@ -753,10 +762,10 @@ hmm
|
|
753
762
|
hmm step time: 3.0
|
754
763
|
writing hmms removed 1 genes from glfo (leaving v 23 d 18 j 6)
|
755
764
|
(1.0 sec)
|
756
|
-
total time: 7.
|
765
|
+
total time: 7.6
|
757
766
|
[44mcache-parameters[0m [94migk[0m:
|
758
|
-
[91mrun[0m partis cache-parameters --locus igk --dont-write-git-info --infname /tmp/runner/hmms/
|
759
|
-
--input-metafnames: added meta info for 37 sequences from /tmp/runner/hmms/
|
767
|
+
[91mrun[0m partis cache-parameters --locus igk --dont-write-git-info --infname /tmp/runner/hmms/870150/igk.fa --parameter-dir test/paired/new-results/test/parameters/data/igk --random-seed 1 --n-procs 10 --sw-cachefname test/paired/new-results/test/parameters/data/igk/sw-cache.yaml --input-metafnames /tmp/runner/hmms/870150/meta.yaml
|
768
|
+
--input-metafnames: added meta info for 37 sequences from /tmp/runner/hmms/870150/meta.yaml: paired-uids
|
760
769
|
caching parameters
|
761
770
|
vsearch: 37 / 37 v annotations (0 failed) with 16 v genes in 0.0 sec
|
762
771
|
keeping 16 / 64 v genes
|
@@ -783,11 +792,11 @@ hmm
|
|
783
792
|
infra time: 1.3
|
784
793
|
hmm step time: 1.9
|
785
794
|
writing hmms removed 0 genes from glfo (leaving v 16 d 1 j 7)
|
786
|
-
(1.
|
795
|
+
(1.0 sec)
|
787
796
|
total time: 5.7
|
788
797
|
[44mcache-parameters[0m [94migl[0m:
|
789
|
-
[91mrun[0m partis cache-parameters --locus igl --dont-write-git-info --infname /tmp/runner/hmms/
|
790
|
-
--input-metafnames: added meta info for 15 sequences from /tmp/runner/hmms/
|
798
|
+
[91mrun[0m partis cache-parameters --locus igl --dont-write-git-info --infname /tmp/runner/hmms/870150/igl.fa --parameter-dir test/paired/new-results/test/parameters/data/igl --random-seed 1 --n-procs 10 --sw-cachefname test/paired/new-results/test/parameters/data/igl/sw-cache.yaml --input-metafnames /tmp/runner/hmms/870150/meta.yaml
|
799
|
+
--input-metafnames: added meta info for 15 sequences from /tmp/runner/hmms/870150/meta.yaml: paired-uids
|
791
800
|
caching parameters
|
792
801
|
vsearch: 15 / 15 v annotations (0 failed) with 12 v genes in 0.0 sec
|
793
802
|
keeping 10 / 80 v genes
|
@@ -814,199 +823,187 @@ hmm
|
|
814
823
|
writing hmms removed 0 genes from glfo (leaving v 10 d 1 j 3)
|
815
824
|
(0.6 sec)
|
816
825
|
total time: 3.5
|
817
|
-
[93mwarning[0m expected to remove 1/11 files+dirs that weren't in /tmp/runner/hmms/
|
826
|
+
[93mwarning[0m expected to remove 1/11 files+dirs that weren't in /tmp/runner/hmms/870150/: failed.fa
|
818
827
|
total time: 22.1
|
819
|
-
[92msimulate[0m partis simulate --dont-write-git-info --paired-loci --parameter-dir test/paired/new-results/test/parameters/data --n-sim-events 10 --n-trees 10 --n-leaf-distribution geometric --n-leaves 5 --min-observations-per-gene 5 --mean-cells-per-droplet 1.25 --constant-cells-per-droplet --fraction-of-reads-to-remove 0.15 --no-per-base-mutation --random-seed 1 --n-procs 10 --paired-outdir test/paired/new-results/test/simu --indel-frequency 0.2
|
828
|
+
[92msimulate[0m partis simulate --dont-write-git-info --paired-loci --parameter-dir test/paired/new-results/test/parameters/data --n-sim-events 10 --n-trees 10 --n-leaf-distribution geometric --n-leaves 5 --min-observations-per-gene 5 --mean-cells-per-droplet 1.25 --constant-cells-per-droplet --fraction-of-reads-to-remove 0.15 --input-simulation-treefname /home/runner/work/partis/partis/test/trees.nwk --no-per-base-mutation --random-seed 1 --n-procs 10 --paired-outdir test/paired/new-results/test/simu --indel-frequency 0.2
|
820
829
|
[94migh+igk[0m: simulating 7 events
|
821
|
-
[
|
822
|
-
[91mrun[0m partis simulate --dont-write-git-info --parameter-dir test/paired/new-results/test/parameters/data/igh --n-trees 7 --n-leaf-distribution geometric --n-leaves 5 --min-observations-per-gene 5 --mean-cells-per-droplet 1.25 --constant-cells-per-droplet --fraction-of-reads-to-remove 0.15 --no-per-base-mutation --random-seed 2003260051 --indel-frequency 0.2 --generate-trees --outfname test/paired/new-results/test/simu/igh+igk/trees.nwk
|
823
|
-
note: reducing --n-procs to 1 (was 4) so it isn't bigger than --n-sim-events
|
824
|
-
generating 7 trees, n-leaves from geometric distribution with parameter 5.0
|
825
|
-
mean over 7 trees: depth 0.02051 leaves 5.57
|
826
|
-
writing trees to test/paired/new-results/test/simu/igh+igk/trees.nwk
|
827
|
-
total time: 1.3
|
830
|
+
[93mwarning[0m drawing input simulation trees from --input-simulation-treefname *in order* for both heavy and light (rather than randomly, as for single chain)
|
828
831
|
[44msimulate[0m [94migh+igk[0m:[94m igh[0m
|
829
|
-
[91mrun[0m partis simulate --locus igh --dont-write-git-info --parameter-dir test/paired/new-results/test/parameters/data/igh --n-sim-events 7 --n-trees 10 --n-leaf-distribution geometric --n-leaves 5 --min-observations-per-gene 5 --mean-cells-per-droplet 1.25 --constant-cells-per-droplet --fraction-of-reads-to-remove 0.15 --no-per-base-mutation --random-seed 2003260051 --n-procs 10 --indel-frequency 0.2 --choose-trees-in-order --outfname test/paired/new-results/test/simu/igh+igk/igh.yaml
|
832
|
+
[91mrun[0m partis simulate --locus igh --dont-write-git-info --parameter-dir test/paired/new-results/test/parameters/data/igh --n-sim-events 7 --n-trees 10 --n-leaf-distribution geometric --n-leaves 5 --min-observations-per-gene 5 --mean-cells-per-droplet 1.25 --constant-cells-per-droplet --fraction-of-reads-to-remove 0.15 --input-simulation-treefname /home/runner/work/partis/partis/test/trees.nwk --no-per-base-mutation --random-seed 2003260051 --n-procs 10 --indel-frequency 0.2 --choose-trees-in-order --outfname test/paired/new-results/test/simu/igh+igk/igh.yaml
|
830
833
|
note: reducing --n-procs to 7 (was 10) so it isn't bigger than --n-sim-events
|
831
834
|
simulating
|
832
835
|
--> proc 0
|
833
836
|
reco params: test/paired/new-results/test/parameters/data/igh/hmm shm params: test/paired/new-results/test/parameters/data/igh/hmm
|
834
|
-
resetting --n-trees from 10 to
|
837
|
+
resetting --n-trees from 10 to 14 to match trees read from /tmp/runner/hmms/938371/trees-sub-0.nwk
|
835
838
|
using bpp uniform mutation with path /home/runner/work/partis/partis/packages/bpp
|
836
|
-
made 1 event with
|
837
|
-
total time: 0.
|
839
|
+
made 1 event with 10 seqs in 0.1s (0.0s of which was running bppseqgen)
|
840
|
+
total time: 0.2
|
838
841
|
|
839
842
|
--> proc 1
|
840
843
|
reco params: test/paired/new-results/test/parameters/data/igh/hmm shm params: test/paired/new-results/test/parameters/data/igh/hmm
|
841
|
-
resetting --n-trees from 10 to
|
844
|
+
resetting --n-trees from 10 to 14 to match trees read from /tmp/runner/hmms/938371/trees-sub-1.nwk
|
842
845
|
using bpp uniform mutation with path /home/runner/work/partis/partis/packages/bpp
|
843
|
-
made 1 event with
|
844
|
-
total time: 0.
|
846
|
+
made 1 event with 8 seqs in 0.2s (0.1s of which was running bppseqgen)
|
847
|
+
total time: 0.3
|
845
848
|
|
846
849
|
--> proc 2
|
847
850
|
reco params: test/paired/new-results/test/parameters/data/igh/hmm shm params: test/paired/new-results/test/parameters/data/igh/hmm
|
848
|
-
resetting --n-trees from 10 to
|
851
|
+
resetting --n-trees from 10 to 14 to match trees read from /tmp/runner/hmms/938371/trees-sub-2.nwk
|
849
852
|
using bpp uniform mutation with path /home/runner/work/partis/partis/packages/bpp
|
850
|
-
made 1 event with
|
851
|
-
total time: 0.
|
853
|
+
made 1 event with 1 seqs in 0.1s (0.0s of which was running bppseqgen)
|
854
|
+
total time: 0.2
|
852
855
|
|
853
856
|
--> proc 3
|
854
857
|
reco params: test/paired/new-results/test/parameters/data/igh/hmm shm params: test/paired/new-results/test/parameters/data/igh/hmm
|
855
|
-
resetting --n-trees from 10 to
|
858
|
+
resetting --n-trees from 10 to 14 to match trees read from /tmp/runner/hmms/938371/trees-sub-3.nwk
|
856
859
|
using bpp uniform mutation with path /home/runner/work/partis/partis/packages/bpp
|
857
|
-
made 1 event with
|
858
|
-
total time: 0.
|
860
|
+
made 1 event with 6 seqs in 0.1s (0.0s of which was running bppseqgen)
|
861
|
+
total time: 0.2
|
859
862
|
|
860
863
|
--> proc 4
|
861
864
|
reco params: test/paired/new-results/test/parameters/data/igh/hmm shm params: test/paired/new-results/test/parameters/data/igh/hmm
|
862
|
-
resetting --n-trees from 10 to
|
865
|
+
resetting --n-trees from 10 to 14 to match trees read from /tmp/runner/hmms/938371/trees-sub-4.nwk
|
863
866
|
using bpp uniform mutation with path /home/runner/work/partis/partis/packages/bpp
|
864
|
-
made 1 event with
|
867
|
+
made 1 event with 1 seqs in 0.0s (0.0s of which was running bppseqgen)
|
865
868
|
total time: 0.1
|
866
869
|
|
867
870
|
--> proc 5
|
868
871
|
reco params: test/paired/new-results/test/parameters/data/igh/hmm shm params: test/paired/new-results/test/parameters/data/igh/hmm
|
869
|
-
resetting --n-trees from 10 to
|
872
|
+
resetting --n-trees from 10 to 14 to match trees read from /tmp/runner/hmms/938371/trees-sub-5.nwk
|
870
873
|
using bpp uniform mutation with path /home/runner/work/partis/partis/packages/bpp
|
871
|
-
made 1 event with
|
872
|
-
total time: 0.
|
874
|
+
made 1 event with 1 seqs in 0.1s (0.0s of which was running bppseqgen)
|
875
|
+
total time: 0.2
|
873
876
|
|
874
877
|
--> proc 6
|
875
878
|
reco params: test/paired/new-results/test/parameters/data/igh/hmm shm params: test/paired/new-results/test/parameters/data/igh/hmm
|
876
|
-
resetting --n-trees from 10 to
|
879
|
+
resetting --n-trees from 10 to 16 to match trees read from /tmp/runner/hmms/938371/trees-sub-6.nwk
|
877
880
|
using bpp uniform mutation with path /home/runner/work/partis/partis/packages/bpp
|
878
|
-
made 1 event with
|
879
|
-
total time: 0.
|
881
|
+
made 1 event with 13 seqs in 0.1s (0.1s of which was running bppseqgen)
|
882
|
+
total time: 0.2
|
880
883
|
|
881
|
-
read 7 events with
|
882
|
-
total time: 4.
|
884
|
+
read 7 events with 40 seqs from 7 .yaml files
|
885
|
+
total time: 4.2
|
883
886
|
[44msimulate[0m [94migh+igk[0m:[94m igk[0m
|
884
|
-
[91mrun[0m partis simulate --locus igk --dont-write-git-info --parameter-dir test/paired/new-results/test/parameters/data/igk --n-sim-events 7 --n-trees 10 --n-leaf-distribution geometric --n-leaves 5 --min-observations-per-gene 5 --mean-cells-per-droplet 1.25 --constant-cells-per-droplet --fraction-of-reads-to-remove 0.15 --no-per-base-mutation --random-seed 2077655082 --n-procs 10 --indel-frequency 0.2 --choose-trees-in-order --outfname test/paired/new-results/test/simu/igh+igk/igk.yaml
|
887
|
+
[91mrun[0m partis simulate --locus igk --dont-write-git-info --parameter-dir test/paired/new-results/test/parameters/data/igk --n-sim-events 7 --n-trees 10 --n-leaf-distribution geometric --n-leaves 5 --min-observations-per-gene 5 --mean-cells-per-droplet 1.25 --constant-cells-per-droplet --fraction-of-reads-to-remove 0.15 --input-simulation-treefname /home/runner/work/partis/partis/test/trees.nwk --no-per-base-mutation --random-seed 2077655082 --n-procs 10 --indel-frequency 0.2 --choose-trees-in-order --outfname test/paired/new-results/test/simu/igh+igk/igk.yaml
|
885
888
|
note: reducing --n-procs to 7 (was 10) so it isn't bigger than --n-sim-events
|
886
889
|
simulating
|
887
890
|
--> proc 0
|
888
891
|
reco params: test/paired/new-results/test/parameters/data/igk/hmm shm params: test/paired/new-results/test/parameters/data/igk/hmm
|
889
|
-
resetting --n-trees from 10 to
|
892
|
+
resetting --n-trees from 10 to 14 to match trees read from /tmp/runner/hmms/129826/trees-sub-0.nwk
|
890
893
|
using bpp uniform mutation with path /home/runner/work/partis/partis/packages/bpp
|
891
|
-
made 1 event with
|
894
|
+
made 1 event with 10 seqs in 0.1s (0.1s of which was running bppseqgen)
|
892
895
|
total time: 0.1
|
893
896
|
|
894
897
|
--> proc 1
|
895
898
|
reco params: test/paired/new-results/test/parameters/data/igk/hmm shm params: test/paired/new-results/test/parameters/data/igk/hmm
|
896
|
-
resetting --n-trees from 10 to
|
899
|
+
resetting --n-trees from 10 to 14 to match trees read from /tmp/runner/hmms/129826/trees-sub-1.nwk
|
897
900
|
using bpp uniform mutation with path /home/runner/work/partis/partis/packages/bpp
|
898
|
-
made 1 event with
|
899
|
-
total time: 0.
|
901
|
+
made 1 event with 8 seqs in 0.1s (0.0s of which was running bppseqgen)
|
902
|
+
total time: 0.2
|
900
903
|
|
901
904
|
--> proc 2
|
902
905
|
reco params: test/paired/new-results/test/parameters/data/igk/hmm shm params: test/paired/new-results/test/parameters/data/igk/hmm
|
903
|
-
resetting --n-trees from 10 to
|
906
|
+
resetting --n-trees from 10 to 14 to match trees read from /tmp/runner/hmms/129826/trees-sub-2.nwk
|
904
907
|
using bpp uniform mutation with path /home/runner/work/partis/partis/packages/bpp
|
905
|
-
made 1 event with
|
908
|
+
made 1 event with 1 seqs in 0.0s (0.0s of which was running bppseqgen)
|
906
909
|
total time: 0.1
|
907
910
|
|
908
911
|
--> proc 3
|
909
912
|
reco params: test/paired/new-results/test/parameters/data/igk/hmm shm params: test/paired/new-results/test/parameters/data/igk/hmm
|
910
|
-
resetting --n-trees from 10 to
|
913
|
+
resetting --n-trees from 10 to 14 to match trees read from /tmp/runner/hmms/129826/trees-sub-3.nwk
|
911
914
|
using bpp uniform mutation with path /home/runner/work/partis/partis/packages/bpp
|
912
|
-
made 1 event with
|
915
|
+
made 1 event with 6 seqs in 0.1s (0.0s of which was running bppseqgen)
|
913
916
|
total time: 0.1
|
914
917
|
|
915
918
|
--> proc 4
|
916
919
|
reco params: test/paired/new-results/test/parameters/data/igk/hmm shm params: test/paired/new-results/test/parameters/data/igk/hmm
|
917
|
-
resetting --n-trees from 10 to
|
920
|
+
resetting --n-trees from 10 to 14 to match trees read from /tmp/runner/hmms/129826/trees-sub-4.nwk
|
918
921
|
using bpp uniform mutation with path /home/runner/work/partis/partis/packages/bpp
|
919
|
-
made 1 event with
|
920
|
-
total time: 0.
|
922
|
+
made 1 event with 1 seqs in 0.1s (0.0s of which was running bppseqgen)
|
923
|
+
total time: 0.2
|
921
924
|
|
922
925
|
--> proc 5
|
923
926
|
reco params: test/paired/new-results/test/parameters/data/igk/hmm shm params: test/paired/new-results/test/parameters/data/igk/hmm
|
924
|
-
resetting --n-trees from 10 to
|
927
|
+
resetting --n-trees from 10 to 14 to match trees read from /tmp/runner/hmms/129826/trees-sub-5.nwk
|
925
928
|
using bpp uniform mutation with path /home/runner/work/partis/partis/packages/bpp
|
926
|
-
made 1 event with
|
929
|
+
made 1 event with 1 seqs in 0.1s (0.0s of which was running bppseqgen)
|
927
930
|
total time: 0.2
|
928
931
|
|
929
932
|
--> proc 6
|
930
933
|
reco params: test/paired/new-results/test/parameters/data/igk/hmm shm params: test/paired/new-results/test/parameters/data/igk/hmm
|
931
|
-
resetting --n-trees from 10 to
|
934
|
+
resetting --n-trees from 10 to 16 to match trees read from /tmp/runner/hmms/129826/trees-sub-6.nwk
|
932
935
|
using bpp uniform mutation with path /home/runner/work/partis/partis/packages/bpp
|
933
|
-
made 1 event with
|
936
|
+
made 1 event with 13 seqs in 0.1s (0.1s of which was running bppseqgen)
|
934
937
|
total time: 0.1
|
935
938
|
|
936
|
-
read 7 events with
|
939
|
+
read 7 events with 40 seqs from 7 .yaml files
|
937
940
|
total time: 4.1
|
938
941
|
[94migh+igl[0m: simulating 3 events
|
939
|
-
[
|
940
|
-
[91mrun[0m partis simulate --dont-write-git-info --parameter-dir test/paired/new-results/test/parameters/data/igh --n-trees 3 --n-leaf-distribution geometric --n-leaves 5 --min-observations-per-gene 5 --mean-cells-per-droplet 1.25 --constant-cells-per-droplet --fraction-of-reads-to-remove 0.15 --no-per-base-mutation --random-seed 982264879 --indel-frequency 0.2 --generate-trees --outfname test/paired/new-results/test/simu/igh+igl/trees.nwk
|
941
|
-
note: reducing --n-procs to 1 (was 4) so it isn't bigger than --n-sim-events
|
942
|
-
generating 3 trees, n-leaves from geometric distribution with parameter 5.0
|
943
|
-
mean over 3 trees: depth 0.05242 leaves 9.00
|
944
|
-
writing trees to test/paired/new-results/test/simu/igh+igl/trees.nwk
|
945
|
-
total time: 1.3
|
942
|
+
[93mwarning[0m drawing input simulation trees from --input-simulation-treefname *in order* for both heavy and light (rather than randomly, as for single chain)
|
946
943
|
[44msimulate[0m [94migh+igl[0m:[94m igh[0m
|
947
|
-
[91mrun[0m partis simulate --locus igh --dont-write-git-info --parameter-dir test/paired/new-results/test/parameters/data/igh --n-sim-events 3 --n-trees 10 --n-leaf-distribution geometric --n-leaves 5 --min-observations-per-gene 5 --mean-cells-per-droplet 1.25 --constant-cells-per-droplet --fraction-of-reads-to-remove 0.15 --no-per-base-mutation --random-seed 982264879 --n-procs 10 --indel-frequency 0.2 --choose-trees-in-order --outfname test/paired/new-results/test/simu/igh+igl/igh.yaml
|
944
|
+
[91mrun[0m partis simulate --locus igh --dont-write-git-info --parameter-dir test/paired/new-results/test/parameters/data/igh --n-sim-events 3 --n-trees 10 --n-leaf-distribution geometric --n-leaves 5 --min-observations-per-gene 5 --mean-cells-per-droplet 1.25 --constant-cells-per-droplet --fraction-of-reads-to-remove 0.15 --input-simulation-treefname /home/runner/work/partis/partis/test/trees.nwk --no-per-base-mutation --random-seed 982264879 --n-procs 10 --indel-frequency 0.2 --choose-trees-in-order --outfname test/paired/new-results/test/simu/igh+igl/igh.yaml
|
948
945
|
note: reducing --n-procs to 3 (was 10) so it isn't bigger than --n-sim-events
|
949
946
|
simulating
|
950
947
|
--> proc 0
|
951
948
|
reco params: test/paired/new-results/test/parameters/data/igh/hmm shm params: test/paired/new-results/test/parameters/data/igh/hmm
|
952
|
-
resetting --n-trees from 10 to
|
949
|
+
resetting --n-trees from 10 to 33 to match trees read from /tmp/runner/hmms/373835/trees-sub-0.nwk
|
953
950
|
using bpp uniform mutation with path /home/runner/work/partis/partis/packages/bpp
|
954
|
-
made 1 event with
|
951
|
+
made 1 event with 10 seqs in 0.1s (0.0s of which was running bppseqgen)
|
955
952
|
total time: 0.1
|
956
953
|
|
957
954
|
--> proc 1
|
958
955
|
reco params: test/paired/new-results/test/parameters/data/igh/hmm shm params: test/paired/new-results/test/parameters/data/igh/hmm
|
959
|
-
resetting --n-trees from 10 to
|
956
|
+
resetting --n-trees from 10 to 33 to match trees read from /tmp/runner/hmms/373835/trees-sub-1.nwk
|
960
957
|
using bpp uniform mutation with path /home/runner/work/partis/partis/packages/bpp
|
961
|
-
made 1 event with
|
958
|
+
made 1 event with 2 seqs in 0.0s (0.0s of which was running bppseqgen)
|
962
959
|
total time: 0.1
|
963
960
|
|
964
961
|
--> proc 2
|
965
962
|
reco params: test/paired/new-results/test/parameters/data/igh/hmm shm params: test/paired/new-results/test/parameters/data/igh/hmm
|
966
|
-
resetting --n-trees from 10 to
|
963
|
+
resetting --n-trees from 10 to 34 to match trees read from /tmp/runner/hmms/373835/trees-sub-2.nwk
|
967
964
|
using bpp uniform mutation with path /home/runner/work/partis/partis/packages/bpp
|
968
|
-
made 1 event with
|
965
|
+
made 1 event with 1 seqs in 0.0s (0.0s of which was running bppseqgen)
|
969
966
|
total time: 0.1
|
970
967
|
|
971
|
-
read 3 events with
|
972
|
-
total time: 2.
|
968
|
+
read 3 events with 13 seqs from 3 .yaml files
|
969
|
+
total time: 2.1
|
973
970
|
[44msimulate[0m [94migh+igl[0m:[94m igl[0m
|
974
|
-
[91mrun[0m partis simulate --locus igl --dont-write-git-info --parameter-dir test/paired/new-results/test/parameters/data/igl --n-sim-events 3 --n-trees 10 --n-leaf-distribution geometric --n-leaves 5 --min-observations-per-gene 5 --mean-cells-per-droplet 1.25 --constant-cells-per-droplet --fraction-of-reads-to-remove 0.15 --no-per-base-mutation --random-seed 1083609742 --n-procs 10 --indel-frequency 0.2 --choose-trees-in-order --outfname test/paired/new-results/test/simu/igh+igl/igl.yaml
|
971
|
+
[91mrun[0m partis simulate --locus igl --dont-write-git-info --parameter-dir test/paired/new-results/test/parameters/data/igl --n-sim-events 3 --n-trees 10 --n-leaf-distribution geometric --n-leaves 5 --min-observations-per-gene 5 --mean-cells-per-droplet 1.25 --constant-cells-per-droplet --fraction-of-reads-to-remove 0.15 --input-simulation-treefname /home/runner/work/partis/partis/test/trees.nwk --no-per-base-mutation --random-seed 1083609742 --n-procs 10 --indel-frequency 0.2 --choose-trees-in-order --outfname test/paired/new-results/test/simu/igh+igl/igl.yaml
|
975
972
|
note: reducing --n-procs to 3 (was 10) so it isn't bigger than --n-sim-events
|
976
973
|
simulating
|
977
974
|
--> proc 0
|
978
975
|
reco params: test/paired/new-results/test/parameters/data/igl/hmm shm params: test/paired/new-results/test/parameters/data/igl/hmm
|
979
|
-
resetting --n-trees from 10 to
|
976
|
+
resetting --n-trees from 10 to 33 to match trees read from /tmp/runner/hmms/635551/trees-sub-0.nwk
|
980
977
|
using bpp uniform mutation with path /home/runner/work/partis/partis/packages/bpp
|
981
|
-
made 1 event with
|
978
|
+
made 1 event with 10 seqs in 0.0s (0.0s of which was running bppseqgen)
|
982
979
|
total time: 0.1
|
983
980
|
|
984
981
|
--> proc 1
|
985
982
|
reco params: test/paired/new-results/test/parameters/data/igl/hmm shm params: test/paired/new-results/test/parameters/data/igl/hmm
|
986
|
-
resetting --n-trees from 10 to
|
983
|
+
resetting --n-trees from 10 to 33 to match trees read from /tmp/runner/hmms/635551/trees-sub-1.nwk
|
987
984
|
using bpp uniform mutation with path /home/runner/work/partis/partis/packages/bpp
|
988
|
-
made 1 event with
|
989
|
-
total time: 0.
|
985
|
+
made 1 event with 2 seqs in 0.0s (0.0s of which was running bppseqgen)
|
986
|
+
total time: 0.1
|
990
987
|
|
991
988
|
--> proc 2
|
992
989
|
reco params: test/paired/new-results/test/parameters/data/igl/hmm shm params: test/paired/new-results/test/parameters/data/igl/hmm
|
993
|
-
resetting --n-trees from 10 to
|
990
|
+
resetting --n-trees from 10 to 34 to match trees read from /tmp/runner/hmms/635551/trees-sub-2.nwk
|
994
991
|
using bpp uniform mutation with path /home/runner/work/partis/partis/packages/bpp
|
995
|
-
made 1 event with
|
996
|
-
total time: 0.
|
992
|
+
made 1 event with 1 seqs in 0.0s (0.0s of which was running bppseqgen)
|
993
|
+
total time: 0.1
|
997
994
|
|
998
|
-
read 3 events with
|
999
|
-
total time: 2.
|
995
|
+
read 3 events with 13 seqs from 3 .yaml files
|
996
|
+
total time: 2.0
|
1000
997
|
[94migh+igk[0m: synchronizing heavy and light chain simulation trees and rewriting output files in test/paired/new-results/test/simu/igh+igk/
|
1001
998
|
[94migh+igl[0m: synchronizing heavy and light chain simulation trees and rewriting output files in test/paired/new-results/test/simu/igh+igl/
|
1002
999
|
concatenating heavy chain loci:
|
1003
|
-
[94mh[0m starting with 7 annotations (
|
1004
|
-
[95mk[0m starting with 7 annotations (
|
1005
|
-
[94mh[0m adding 3 annotations (
|
1006
|
-
[92ml[0m starting with 3 annotations (
|
1000
|
+
[94mh[0m starting with 7 annotations (40 seqs)
|
1001
|
+
[95mk[0m starting with 7 annotations (40 seqs)
|
1002
|
+
[94mh[0m adding 3 annotations (13 seqs) for total 7 (40)
|
1003
|
+
[92ml[0m starting with 3 annotations (13 seqs)
|
1007
1004
|
finished with:
|
1008
|
-
[94mh[0m: 10 annotations (
|
1009
|
-
apportioned
|
1010
|
-
removed
|
1005
|
+
[94mh[0m: 10 annotations (53 seqs)
|
1006
|
+
apportioned 106 seqs among 42 droplets (mean/2 1.3): 2 2 2 4 4 2 4 4 4 4 2 2 2 4 2 2 2 2 2 2 2 2 2 2 4 2 4 2 2 2 2 2 4 4 2 2 2 2 2 2 2 2
|
1007
|
+
removed 15 / 106 = 0.14 seqs from outfos (leaving 9 / 91 unpaired)
|
1011
1008
|
writing combined h/l info to fasta and meta files in test/paired/new-results/test/simu
|
1012
|
-
total time:
|
1009
|
+
total time: 17.6
|
test/test.py
CHANGED
@@ -299,8 +299,8 @@ class Tester(object):
|
|
299
299
|
|
300
300
|
if ptest == 'simulate':
|
301
301
|
argfo['parameter-dir'] = self.paramdir(input_stype, 'data')
|
302
|
-
|
303
|
-
|
302
|
+
if args.no_tree_gen:
|
303
|
+
argfo['extras'] += ['--input-simulation-treefname', '%s/test/trees.nwk'%utils.get_partis_dir()]
|
304
304
|
if args.no_per_base_mutation:
|
305
305
|
argfo['extras'] += ['--no-per-base-mutation']
|
306
306
|
else:
|
@@ -936,7 +936,7 @@ parser.add_argument('--only-bust-current', action='store_true', help='only bust
|
|
936
936
|
parser.add_argument('--paired', action='store_true', help='run paired tests, i.e. with --paired-loci. Note that this doesn\'t test all the things (e.g. seed partitioning) that non-paired does.')
|
937
937
|
parser.add_argument('--run-all', action='store_true', help='run all four combinations of tests: paired/non-paired and slow/non-slow (by default only runs one). *Not* for use with --bust-cache, which runs all of them by default.')
|
938
938
|
parser.add_argument('--no-simu', action='store_true', help='don\'t run simulation, e.g. if using a minimal install')
|
939
|
-
|
939
|
+
parser.add_argument('--no-tree-gen', action='store_true', help='if set, use the pre-generated newick file instead of generating trees with R/TreeSim (to avoid R installation)')
|
940
940
|
parser.add_argument('--no-per-base-mutation', action='store_true', help='use simpler, non-per-base mutation model (to avoid bpp-newlik compilation)')
|
941
941
|
parser.add_argument('--ig-or-tr', default='ig')
|
942
942
|
parser.add_argument('--print-width', type=int, default=300, help='set to 0 for infinite')
|