readMultiTopology err

[mos]

Dera all,

I'm getting the following error, when trying to run a bootstrap analysis on the snaq best network.

ERROR: Could not find or open  Dendroscope: (Sarcoptes_scabiei,Varroa_destructor,((Acalitus_phloeocoptes_X6,(Calepitrimerus_cf_scutulatus_Z352,(Tetra_qingyuanensis_Z333,(Calepitrimerus_flexuosus_Z342,(Tegolophus_suffruticosae_Z382,(Litaculus_sp_Z392,((Apodiptacus_sp_Z313,Rhinophytoptus_broussonetiae_Z401)1:0.4872546611169829,Keiferella_guanegouensis_QH1)1:0.2523692626914002)1:0.1492247079837156)0.58:0.03723301317997846)0.91:0.049385270622364544)1:0.1842758831788743)1:3.635895063719693)1:3.846803985797675,(Fragariocoptes_setiger,Trisetacus_juniperinus_HEN46)1:0.6584329897379736):5.378431493554643)1; file

...

in expression starting at /gss1/home/szk20200608/temp/snaq/bootSNaQ.jl:21

It seems like something went wrong with "readMultiTopology", it try to open a file named after the input tree, rather than input the tree itself. But I have no idea why this happens, perhaps I used the wrong argument or wrong input format?

Script and input files are linked below.

Thanks a lot!

[Cécile Ané]

The files "net0.out" and "net2.out" do not have multiple topologies: they start with 1 topology, then give other information. So the command asking to read multiple topologies has a problem with the second line, which starts with "Dendroscope..." and is not a topology.

The error should go away if you replace these lines:

net0 = readMultiTopology("net0.out")[1]
net2 = readMultiTopology("net2.out")[1]

by this instead, to ask to read only 1 topology (assumed to be on the first line):

net0 = readTopology("net0.out")
net2 = readTopology("net2.out")

Cecile.