start 2026-09-28T08:57:51Z host n3469 pair=NS18_Sequim_Wild__SS18_North_Bay_Wild [persaf::persaf_init] Version of output/06_angsd_structure/saf/NS18_Sequim_Wild.saf.idx is 3 [persaf::persaf_init] Assuming .saf.gz file is output/06_angsd_structure/saf/NS18_Sequim_Wild.saf.gz [persaf::persaf_init] Assuming .saf.pos.gz file is output/06_angsd_structure/saf/NS18_Sequim_Wild.saf.pos.gz [persaf::persaf_init] Version of output/06_angsd_structure/saf/SS18_North_Bay_Wild.saf.idx is 3 [persaf::persaf_init] Assuming .saf.gz file is output/06_angsd_structure/saf/SS18_North_Bay_Wild.saf.gz [persaf::persaf_init] Assuming .saf.pos.gz file is output/06_angsd_structure/saf/SS18_North_Bay_Wild.saf.pos.gz -> args: tole:0.000000 nthreads:16 maxiter:60 nsites(block):0 start:(null) chr:CM196826.1 start:-1 stop:-1 fstout:(null) oldout:0 seed:-1 bootstrap:0 resample_chr:0 whichFst:0 fold:1 ref:(null) anc:(null) [main] Multi SFS is 'still' under development. Please report strange behaviour. -> The choice of -nSites will require atleast: 7693.371582 megabyte memory, that is at least: 0.50% of total memory -> dim(output/06_angsd_structure/saf/NS18_Sequim_Wild.saf.idx):17 -> dim(output/06_angsd_structure/saf/SS18_North_Bay_Wild.saf.idx):17 -> Dimension of parameter space: 289 -> Done reading data from chromosome will prepare next chromosome -> hello Im the master merge part of realSFS. and I'll now do a tripple bypass to find intersect -> 1) Will set iter according to chooseChr and start and stop, and possibly using -sites -> Sites to keep[CM196826.1] from pop0: 79332133 -> Sites to keep[CM196826.1] from pop1: 79332133 -> Will run optimization on nSites: 79332133 ------------ startlik=-375131247.023087 -> Instability detected, accelerated guess is too close to bound or outside will fallback to regular EM for this step likelihood: -12843880.206516 ------------ -> Only read nSites: 0 will therefore prepare next chromosome (or exit) -> NB output is no longer log probs of the frequency spectrum! -> Output is now simply the expected values! -> You can convert to the old format simply with log(norm(x)) -> Please check that it has converged! -> You can add start new optimization by supplying -sfs FILE, where is >FILE from this run -> -maxIter INT -tole FLOAT [persaf::persaf_init] Version of output/06_angsd_structure/saf/NS18_Sequim_Wild.saf.idx is 3 [persaf::persaf_init] Assuming .saf.gz file is output/06_angsd_structure/saf/NS18_Sequim_Wild.saf.gz [persaf::persaf_init] Assuming .saf.pos.gz file is output/06_angsd_structure/saf/NS18_Sequim_Wild.saf.pos.gz [persaf::persaf_init] Version of output/06_angsd_structure/saf/SS18_North_Bay_Wild.saf.idx is 3 [persaf::persaf_init] Assuming .saf.gz file is output/06_angsd_structure/saf/SS18_North_Bay_Wild.saf.gz [persaf::persaf_init] Assuming .saf.pos.gz file is output/06_angsd_structure/saf/SS18_North_Bay_Wild.saf.pos.gz -> args: tole:0.000000 nthreads:16 maxiter:100 nsites(block):0 start:output/06_angsd_structure/fst/NS18_Sequim_Wild__SS18_North_Bay_Wild.ml chr:(null) start:-1 stop:-1 fstout:output/06_angsd_structure/fst/NS18_Sequim_Wild__SS18_North_Bay_Wild oldout:0 seed:-1 bootstrap:0 resample_chr:0 whichFst:0 fold:1 ref:(null) anc:(null) -> nSites: 100000 -> IMPORTANT: please make sure that your saf files haven't been folded with -fold 1 in -doSaf in angsd -> [reynoldFst] sfs1:16 sfs2:16 dimspace:289 -> generating offset remapper lookup -> Reading: output/06_angsd_structure/fst/NS18_Sequim_Wild__SS18_North_Bay_Wild.ml assuming counts (will normalize to probs internally) -> Done reading data from chromosome will prepare next chromosome -> hello Im the master merge part of realSFS. and I'll now do a tripple bypass to find intersect -> 1) Will set iter according to chooseChr and start and stop, and possibly using -sites -> Sites to keep[CM196826.1] from pop0: 79332133 -> Sites to keep[CM196826.1] from pop1: 79332133 -> Only read nSites: 0 will therefore prepare next chromosome (or exit) -> Done reading data from chromosome will prepare next chromosome -> hello Im the master merge part of realSFS. and I'll now do a tripple bypass to find intersect -> 1) Will set iter according to chooseChr and start and stop, and possibly using -sites -> Sites to keep[CM196827.1] from pop0: 78783622 -> Sites to keep[CM196827.1] from pop1: 78783622 -> Only read nSites: 0 will therefore prepare next chromosome (or exit) -> Done reading data from chromosome will prepare next chromosome -> hello Im the master merge part of realSFS. and I'll now do a tripple bypass to find intersect -> 1) Will set iter according to chooseChr and start and stop, and possibly using -sites -> Sites to keep[CM196828.1] from pop0: 78043019 -> Sites to keep[CM196828.1] from pop1: 78043019 -> Only read nSites: 0 will therefore prepare next chromosome (or exit) -> Done reading data from chromosome will prepare next chromosome -> hello Im the master merge part of realSFS. and I'll now do a tripple bypass to find intersect -> 1) Will set iter according to chooseChr and start and stop, and possibly using -sites -> Sites to keep[CM196829.1] from pop0: 68302922 -> Sites to keep[CM196829.1] from pop1: 68302922 -> Only read nSites: 0 will therefore prepare next chromosome (or exit) -> Done reading data from chromosome will prepare next chromosome -> hello Im the master merge part of realSFS. and I'll now do a tripple bypass to find intersect -> 1) Will set iter according to chooseChr and start and stop, and possibly using -sites -> Sites to keep[CM196830.1] from pop0: 68278780 -> Sites to keep[CM196830.1] from pop1: 68278780 -> Only read nSites: 0 will therefore prepare next chromosome (or exit) -> Done reading data from chromosome will prepare next chromosome -> hello Im the master merge part of realSFS. and I'll now do a tripple bypass to find intersect -> 1) Will set iter according to chooseChr and start and stop, and possibly using -sites -> Sites to keep[CM196831.1] from pop0: 41371036 -> Sites to keep[CM196831.1] from pop1: 41371036 -> Only read nSites: 0 will therefore prepare next chromosome (or exit) -> Done reading data from chromosome will prepare next chromosome -> hello Im the master merge part of realSFS. and I'll now do a tripple bypass to find intersect -> 1) Will set iter according to chooseChr and start and stop, and possibly using -sites -> Sites to keep[CM196832.1] from pop0: 54082496 -> Sites to keep[CM196832.1] from pop1: 54082496 -> Only read nSites: 0 will therefore prepare next chromosome (or exit) -> Done reading data from chromosome will prepare next chromosome -> hello Im the master merge part of realSFS. and I'll now do a tripple bypass to find intersect -> 1) Will set iter according to chooseChr and start and stop, and possibly using -sites -> Sites to keep[CM196833.1] from pop0: 44669003 -> Sites to keep[CM196833.1] from pop1: 44669003 -> Only read nSites: 0 will therefore prepare next chromosome (or exit) -> Done reading data from chromosome will prepare next chromosome -> hello Im the master merge part of realSFS. and I'll now do a tripple bypass to find intersect -> 1) Will set iter according to chooseChr and start and stop, and possibly using -sites -> Sites to keep[CM196834.1] from pop0: 30962892 -> Sites to keep[CM196834.1] from pop1: 30962892 -> Only read nSites: 0 will therefore prepare next chromosome (or exit) -> Done reading data from chromosome will prepare next chromosome -> hello Im the master merge part of realSFS. and I'll now do a tripple bypass to find intersect -> 1) Will set iter according to chooseChr and start and stop, and possibly using -sites -> Sites to keep[CM196835.1] from pop0: 28978448 -> Sites to keep[CM196835.1] from pop1: 28978448 -> Only read nSites: 0 will therefore prepare next chromosome (or exit) -> fst index finished with no errors! Example runs: realSFS fst stats output/06_angsd_structure/fst/NS18_Sequim_Wild__SS18_North_Bay_Wild.fst.idx realSFS fst stats2 output/06_angsd_structure/fst/NS18_Sequim_Wild__SS18_North_Bay_Wild.fst.idx end 2026-09-28T10:12:40Z pair=NS18_Sequim_Wild__SS18_North_Bay_Wild fst=0.033844 0.051931