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.
Files changed (55) hide show
  1. bin/partis-test.py +3 -3
  2. {partis_bcr-1.0.11.post1.dev3.data → partis_bcr-1.0.12.data}/scripts/partis-test.py +3 -3
  3. {partis_bcr-1.0.11.post1.dev3.dist-info → partis_bcr-1.0.12.dist-info}/METADATA +1 -1
  4. {partis_bcr-1.0.11.post1.dev3.dist-info → partis_bcr-1.0.12.dist-info}/RECORD +55 -54
  5. test/new-results/run-times.csv +5 -5
  6. test/new-results/test/simu.yaml +1 -1
  7. test/new-results/test.log +26 -26
  8. test/paired/new-results/partition-new-simu/fasttree/iclust-2/log +1 -1
  9. test/paired/new-results/partition-new-simu/igh+igk/partition-igh/fasttree/iclust-1/log +2 -2
  10. test/paired/new-results/partition-new-simu-annotation-performance/plots/flcount-matrix.svg +107 -107
  11. test/paired/new-results/partition-new-simu-annotation-performance/plots/func-non-func-per-drop-log.svg +83 -83
  12. test/paired/new-results/partition-new-simu-annotation-performance/plots/func-non-func-per-drop.svg +61 -61
  13. test/paired/new-results/partition-new-simu-annotation-performance/plots/paired-seqs-per-seq-after.svg +18 -18
  14. test/paired/new-results/partition-new-simu-annotation-performance/plots/paired-seqs-per-seq-before.svg +24 -24
  15. test/paired/new-results/partition-new-simu-annotation-performance/plots/pseq-matrix.svg +35 -35
  16. test/paired/new-results/partition-new-simu-annotation-performance/plots/seqs-per-droplet.svg +22 -22
  17. test/paired/new-results/partition-new-simu-annotation-performance/plots/true-pair-clean-performance-correct-family.svg +45 -45
  18. test/paired/new-results/partition-new-simu-annotation-performance/plots/true-pair-clean-performance-correct.svg +45 -45
  19. test/paired/new-results/partition-new-simu-annotation-performance/plots/true-pair-clean-performance-mispaired.svg +45 -45
  20. test/paired/new-results/partition-new-simu-annotation-performance/plots/true-pair-clean-performance-multiple.svg +15 -15
  21. test/paired/new-results/partition-new-simu-annotation-performance/plots/true-pair-clean-performance-near-family.svg +15 -15
  22. test/paired/new-results/partition-new-simu-annotation-performance/plots/true-pair-clean-performance-total.svg +47 -47
  23. test/paired/new-results/partition-new-simu-annotation-performance/plots/true-pair-clean-performance-unpaired.svg +48 -48
  24. test/paired/new-results/partition-new-simu-annotation-performance/plots/true-pair-clean-performance.svg +27 -27
  25. test/paired/new-results/run-times.csv +6 -6
  26. test/paired/new-results/subset-partition-new-simu/isub-0/partition.log +9 -5
  27. test/paired/new-results/subset-partition-new-simu/isub-1/partition.log +6 -2
  28. test/paired/new-results/subset-partition-new-simu/merged-partition.log +7 -3
  29. test/paired/new-results/subset-partition-new-simu/merged-subsets/parameters/igh/sw-cache.yaml +1 -1
  30. test/paired/new-results/subset-partition-new-simu/merged-subsets/parameters/igk/sw-cache.yaml +1 -1
  31. test/paired/new-results/subset-partition-new-simu/merged-subsets/parameters/igl/sw-cache.yaml +1 -1
  32. test/paired/new-results/test/simu/all-seqs.fa +168 -212
  33. test/paired/new-results/test/simu/igh+igk/igh.yaml +1 -1
  34. test/paired/new-results/test/simu/igh+igk/igk.yaml +1 -1
  35. test/paired/new-results/test/simu/igh+igl/igh.yaml +1 -1
  36. test/paired/new-results/test/simu/igh+igl/igl.yaml +1 -1
  37. test/paired/new-results/test/simu/igh.yaml +1 -1
  38. test/paired/new-results/test/simu/igk.yaml +1 -1
  39. test/paired/new-results/test/simu/igl.yaml +1 -1
  40. test/paired/new-results/test/simu/meta.yaml +1 -1
  41. test/paired/new-results/test.log +115 -118
  42. test/test.py +3 -3
  43. test/trees.nwk +100 -0
  44. {partis_bcr-1.0.11.post1.dev3.data → partis_bcr-1.0.12.data}/scripts/cf-alleles.py +0 -0
  45. {partis_bcr-1.0.11.post1.dev3.data → partis_bcr-1.0.12.data}/scripts/cf-germlines.py +0 -0
  46. {partis_bcr-1.0.11.post1.dev3.data → partis_bcr-1.0.12.data}/scripts/compare-plotdirs.py +0 -0
  47. {partis_bcr-1.0.11.post1.dev3.data → partis_bcr-1.0.12.data}/scripts/extract-pairing-info.py +0 -0
  48. {partis_bcr-1.0.11.post1.dev3.data → partis_bcr-1.0.12.data}/scripts/gctree-run.py +0 -0
  49. {partis_bcr-1.0.11.post1.dev3.data → partis_bcr-1.0.12.data}/scripts/get-naive-probabilities.py +0 -0
  50. {partis_bcr-1.0.11.post1.dev3.data → partis_bcr-1.0.12.data}/scripts/parse-output.py +0 -0
  51. {partis_bcr-1.0.11.post1.dev3.data → partis_bcr-1.0.12.data}/scripts/split-loci.py +0 -0
  52. {partis_bcr-1.0.11.post1.dev3.dist-info → partis_bcr-1.0.12.dist-info}/WHEEL +0 -0
  53. {partis_bcr-1.0.11.post1.dev3.dist-info → partis_bcr-1.0.12.dist-info}/entry_points.txt +0 -0
  54. {partis_bcr-1.0.11.post1.dev3.dist-info → partis_bcr-1.0.12.dist-info}/licenses/COPYING +0 -0
  55. {partis_bcr-1.0.11.post1.dev3.dist-info → partis_bcr-1.0.12.dist-info}/top_level.txt +0 -0
@@ -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.5 sec)
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.5 sec)
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.1
33
+ hmm step time: 2.2
34
34
  writing hmms removed 3 genes from glfo (leaving v 6 d 15 j 4)
35
- (0.3 sec)
36
- total time: 3.9
35
+ (0.4 sec)
36
+ total time: 4.1
37
37
  cache-parameters igk:
38
38
  note: adding input metafname from default location (test/paired/ref-results/test/simu/meta.yaml)
39
39
  run 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.4 processing 0.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.7
107
- total time: 12.3
106
+ total time: 1.6
107
+ total time: 12.6
108
108
  partition-new-simu 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
  warning ignoring --plot-annotation-performance for paired clustering since it's going to be a bit fiddly to implement
110
110
  partition igh:
@@ -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.0s)
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.5s)
127
+ 50 clusters with 10 procs (0.6s)
128
128
  44 clusters with 7 procs (0.5s)
129
- 36 clusters with 5 procs (0.4s)
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
  note not merging entire cpath history
134
- partition loop time: 2.8
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
  warning 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: 4.9
141
+ total time: 5.0
142
142
  partition igk:
143
143
  note: adding input metafname from default location (test/paired/ref-results/test/simu/meta.yaml)
144
144
  run 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
  annotate igh+igk: igh
232
233
  note: adding input metafname from default location (test/paired/ref-results/test/simu/meta.yaml)
233
234
  run 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
  annotate igh+igk: igk
249
251
  note: adding input metafname from default location (test/paired/ref-results/test/simu/meta.yaml)
250
252
  run 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
  annotate igh+igl: igh
290
293
  note: adding input metafname from default location (test/paired/ref-results/test/simu/meta.yaml)
291
294
  run 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
  annotate igh+igl: igl
307
311
  note: adding input metafname from default location (test/paired/ref-results/test/simu/meta.yaml)
308
312
  run 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
  h 20 --> 23 +3
332
336
  k 6 --> 15 +9
333
337
  l 3 --> 8 +5
334
- total time: 19.9
338
+ total time: 20.5
335
339
  subset-partition-new-simu 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
  run 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: 68.6
372
+ total time: 69.6
366
373
  seed-partition-new-simu 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.6
433
+ total time: 0.7
427
434
  partition igl:
428
435
  note: adding input metafname from default location (test/paired/ref-results/test/simu/meta.yaml)
429
436
  run 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
  annotate igh+igk: igh
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
  annotate igh+igk: igk
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
  h 12 --> 1 -11
554
563
  k 2 --> 1 -1
555
- total time: 9.3
564
+ total time: 9.4
556
565
  get-selection-metrics-new-simu 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
  get-selection-metrics igh+igk: igh
558
567
  run 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.6
685
+ total time: 5.9
677
686
  cache-parameters-data 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
- run extract-pairing-info.py test/paired-data/all-seqs.fa /tmp/runner/hmms/888813/meta.yaml --n-max-queries 100
688
+ run 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
- run split-loci.py test/paired-data/all-seqs.fa --outdir /tmp/runner/hmms/888813 --input-metafname /tmp/runner/hmms/888813/meta.yaml --n-max-queries 100
699
+ run 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/888813/meta.yaml: paired-uids
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 h l
713
722
  1 0.07 h k
714
723
  1 0.07 h h k
715
- writing to /tmp/runner/hmms/888813/
716
- igh: 48 to 888813/igh.fa
717
- igk: 37 to 888813/igk.fa
718
- igl: 15 to 888813/igl.fa
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
  cache-parameters igh:
727
- run partis cache-parameters --locus igh --dont-write-git-info --infname /tmp/runner/hmms/888813/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/888813/meta.yaml
728
- --input-metafnames: added meta info for 48 sequences from /tmp/runner/hmms/888813/meta.yaml: paired-uids
736
+ run 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.5
765
+ total time: 7.6
757
766
  cache-parameters igk:
758
- run partis cache-parameters --locus igk --dont-write-git-info --infname /tmp/runner/hmms/888813/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/888813/meta.yaml
759
- --input-metafnames: added meta info for 37 sequences from /tmp/runner/hmms/888813/meta.yaml: paired-uids
767
+ run 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.1 sec)
795
+ (1.0 sec)
787
796
  total time: 5.7
788
797
  cache-parameters igl:
789
- run partis cache-parameters --locus igl --dont-write-git-info --infname /tmp/runner/hmms/888813/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/888813/meta.yaml
790
- --input-metafnames: added meta info for 15 sequences from /tmp/runner/hmms/888813/meta.yaml: paired-uids
798
+ run 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
- warning expected to remove 1/11 files+dirs that weren't in /tmp/runner/hmms/888813/: failed.fa
826
+ warning expected to remove 1/11 files+dirs that weren't in /tmp/runner/hmms/870150/: failed.fa
818
827
  total time: 22.1
819
- simulate 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
+ simulate 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
  igh+igk: simulating 7 events
821
- generate-trees igh+igk: igh
822
- run 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
+ warning drawing input simulation trees from --input-simulation-treefname *in order* for both heavy and light (rather than randomly, as for single chain)
828
831
  simulate igh+igk: igh
829
- run 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 --input-simulation-treefname test/paired/new-results/test/simu/igh+igk/trees.nwk
832
+ run 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 1 to match trees read from /tmp/runner/hmms/390989/trees-sub-0.nwk
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 12 seqs in 0.1s (0.1s of which was running bppseqgen)
837
- total time: 0.1
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 1 to match trees read from /tmp/runner/hmms/390989/trees-sub-1.nwk
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 3 seqs in 0.2s (0.0s of which was running bppseqgen)
844
- total time: 0.2
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 1 to match trees read from /tmp/runner/hmms/390989/trees-sub-2.nwk
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 3 seqs in 0.0s (0.0s of which was running bppseqgen)
851
- total time: 0.1
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 1 to match trees read from /tmp/runner/hmms/390989/trees-sub-3.nwk
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 2 seqs in 0.1s (0.0s of which was running bppseqgen)
858
- total time: 0.1
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 1 to match trees read from /tmp/runner/hmms/390989/trees-sub-4.nwk
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 2 seqs in 0.0s (0.0s of which was running bppseqgen)
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 1 to match trees read from /tmp/runner/hmms/390989/trees-sub-5.nwk
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 9 seqs in 0.1s (0.0s of which was running bppseqgen)
872
- total time: 0.1
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 1 to match trees read from /tmp/runner/hmms/390989/trees-sub-6.nwk
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 8 seqs in 0.1s (0.0s of which was running bppseqgen)
879
- total time: 0.1
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 39 seqs from 7 .yaml files
882
- total time: 4.1
884
+ read 7 events with 40 seqs from 7 .yaml files
885
+ total time: 4.2
883
886
  simulate igh+igk: igk
884
- run 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 --input-simulation-treefname test/paired/new-results/test/simu/igh+igk/trees.nwk
887
+ run 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 1 to match trees read from /tmp/runner/hmms/780025/trees-sub-0.nwk
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 12 seqs in 0.1s (0.1s of which was running bppseqgen)
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 1 to match trees read from /tmp/runner/hmms/780025/trees-sub-1.nwk
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 3 seqs in 0.0s (0.0s of which was running bppseqgen)
899
- total time: 0.1
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 1 to match trees read from /tmp/runner/hmms/780025/trees-sub-2.nwk
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 3 seqs in 0.0s (0.0s of which was running bppseqgen)
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 1 to match trees read from /tmp/runner/hmms/780025/trees-sub-3.nwk
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 2 seqs in 0.0s (0.0s of which was running bppseqgen)
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 1 to match trees read from /tmp/runner/hmms/780025/trees-sub-4.nwk
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 2 seqs in 0.0s (0.0s of which was running bppseqgen)
920
- total time: 0.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 1 to match trees read from /tmp/runner/hmms/780025/trees-sub-5.nwk
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 9 seqs in 0.1s (0.0s of which was running bppseqgen)
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 1 to match trees read from /tmp/runner/hmms/780025/trees-sub-6.nwk
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 8 seqs in 0.1s (0.0s of which was running bppseqgen)
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 39 seqs from 7 .yaml files
939
+ read 7 events with 40 seqs from 7 .yaml files
937
940
  total time: 4.1
938
941
  igh+igl: simulating 3 events
939
- generate-trees igh+igl: igh
940
- run 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
+ warning drawing input simulation trees from --input-simulation-treefname *in order* for both heavy and light (rather than randomly, as for single chain)
946
943
  simulate igh+igl: igh
947
- run 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 --input-simulation-treefname test/paired/new-results/test/simu/igh+igl/trees.nwk
944
+ run 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 1 to match trees read from /tmp/runner/hmms/99260/trees-sub-0.nwk
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 14 seqs in 0.1s (0.0s of which was running bppseqgen)
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 1 to match trees read from /tmp/runner/hmms/99260/trees-sub-1.nwk
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 3 seqs in 0.0s (0.0s of which was running bppseqgen)
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 1 to match trees read from /tmp/runner/hmms/99260/trees-sub-2.nwk
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 10 seqs in 0.1s (0.0s of which was running bppseqgen)
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 27 seqs from 3 .yaml files
972
- total time: 2.0
968
+ read 3 events with 13 seqs from 3 .yaml files
969
+ total time: 2.1
973
970
  simulate igh+igl: igl
974
- run 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 --input-simulation-treefname test/paired/new-results/test/simu/igh+igl/trees.nwk
971
+ run 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 1 to match trees read from /tmp/runner/hmms/63603/trees-sub-0.nwk
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 14 seqs in 0.1s (0.0s of which was running bppseqgen)
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 1 to match trees read from /tmp/runner/hmms/63603/trees-sub-1.nwk
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 3 seqs in 0.0s (0.0s of which was running bppseqgen)
989
- total time: 0.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 1 to match trees read from /tmp/runner/hmms/63603/trees-sub-2.nwk
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 10 seqs in 0.0s (0.0s of which was running bppseqgen)
996
- total time: 0.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 27 seqs from 3 .yaml files
999
- total time: 2.2
995
+ read 3 events with 13 seqs from 3 .yaml files
996
+ total time: 2.0
1000
997
  igh+igk: synchronizing heavy and light chain simulation trees and rewriting output files in test/paired/new-results/test/simu/igh+igk/
1001
998
  igh+igl: 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
- h starting with 7 annotations (39 seqs)
1004
- k starting with 7 annotations (39 seqs)
1005
- h adding 3 annotations (27 seqs) for total 7 (39)
1006
- l starting with 3 annotations (27 seqs)
1000
+ h starting with 7 annotations (40 seqs)
1001
+ k starting with 7 annotations (40 seqs)
1002
+ h adding 3 annotations (13 seqs) for total 7 (40)
1003
+ l starting with 3 annotations (13 seqs)
1007
1004
  finished with:
1008
- h: 10 annotations (66 seqs)
1009
- apportioned 132 seqs among 52 droplets (mean/2 1.3): 4 2 2 2 2 4 2 2 2 2 2 2 4 2 2 2 4 4 2 2 4 2 2 2 2 2 2 4 2 2 4 2 4 2 2 2 2 2 2 2 2 2 4 4 2 4 2 2 2 4 4 2
1010
- removed 19 / 132 = 0.14 seqs from outfos (leaving 7 / 113 unpaired)
1005
+ h: 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: 22.8
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
- # if args.no_tree_gen:
303
- # argfo['extras'] += ['--input-simulation-treefname', '%s/test/trees.nwk'%utils.get_partis_dir()]
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
- # 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)')
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')