error when doing the final epa-ng placement

[AnaMariaCabello]

Hi! I got this error when running: epa-ng -t reference_tree.raxml.bestTree -s reference_aligment.fasta -q query.fasta -w ./ --model reference_tree.raxml.bestModel

Please find attached the .err and the info.log files jobLog_447616.err.txt epa_info.log.txt

It looks like the query.fasta produces after ena-ng --split has a problem, but it looks normal to me.

Thank you so much for your help.

[pierrebarbera]

Hi @AnaMariaCabello !

looks like you're on an older version of epa-ng (0.3.5), could you try again with the current version?

Pierre

[AnaMariaCabello]

oK! I will try that. which one is the current version?

Thank you so much Pierre

[pierrebarbera]

That would be v0.3.8. How did you install the program?

[AnaMariaCabello]

I didn't install it myself, it's installed in a cluster.

[pierrebarbera]

Hi @AnaMariaCabello, have you been able to update the version yet? Alternatively you can install the program locally for your user, following the instructions in the readme in this repository. let me know if you need any help.

[AnaMariaCabello]

Hi Pierre! not yet, I’m waiting for the cluster tech to uodate the version. I’ll let you know how it goes as soon as I try again.

Thanks, Ana

On Nov 24, 2021, at 6:35 PM, Pierre Barbera @.***> wrote:

Hi @AnaMariaCabello https://github.com/AnaMariaCabello, have you been able to update the version yet? Alternatively you can install the program locally for your user, following the instructions in the readme in this repository. let me know if you need any help.

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub https://github.com/Pbdas/epa-ng/issues/40#issuecomment-978089148, or unsubscribe https://github.com/notifications/unsubscribe-auth/APMXFJF6A7INXIV5HF36R4TUNUO55ANCNFSM5ITKR2CA. Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

[AnaMariaCabello]

Hi Pierre, I finally tried to do a placement of short sequences (ca. 400 bps, a total of 39 sequences ) to a reference tree (24 reference sequences, ca. 1800 bps) using version v0.3.8. The papara alignment worked and the epa-ng seemed to work too, but the epa_result.jplace output after running the epa-ng command did not include the short sequences, so I obtained only the sequences in the reference tree (I used itol to visualize the epa_result.jplace tree).

Is there any further step needed?

please find attached the commands and the .out files for both. The error file for epa-ng was empty and the error file for papara only contained this: sweep: 138 -> 0 commands_placeme [epang_jobLog_458414.out.txt](https://github.com/Pbdas/epa-ng/files/7888559/epang_jobLog_458414.out.txt) nt.txt papara_jobLog_458409.out.txt

Thanks a lot! Ana

[pierrebarbera]

Hi Ana,

are you sure the .jplace doesn't contain the query sequences? When visualising in iTOL, placements are displayed as red circles (when clicking "on" in the box that says phylogenetic placements), as opposed to new branches in the tree, could this be the mixup?

Pierre

[AnaMariaCabello]

Dear Pierre, You were right, that was the mixup. Is there a way to visualise the placements as a new branches in the reference tree? this is my goal, but maybe I should have used a different approach. Should I use instead guppy or pplacer? could I use as input file for this the output of epa-ng, so the epa_result.jplace file? Any guidance on this will be appreciated. Thank you so much, Ana

[pierrebarbera]

Dear Ana,

you can generate such a tree using our other software called gappa, which includes a whole host of post-analysis methods that operate on phylogenetic placement results (in the form of .jplace files). The command you are looking for is called graft. This will generate a newick file with the queries mixed in.

Let us know if you run into any other issues! Pierre

[AnaMariaCabello]

OK, great! Thank you so much!