Scripts for parsing CAFE5 output?

[Hanna Dort]

Hello,

I'm currently using CAFE5 for the first time, and I'm having a bit of trouble parsing the output files. With CAFE4, I had a python script written by Gregg Thomas (report_analysis.py) that would take a .cafe file as input and produce a number of files, including a "-fams.txt" file that would show (if using the --rapid option) gene families that underwent a significant expansion or contraction on the branch leading to a given node.

Example line: Pluxyl<0>:      25[+49*],42[+33*],46[+60*],51[+36*],66[+33*],84[+38*],88[+12*]

This python script unfortunately doesn't work with the .cafe file in my CAFE5 output. Does a similar script exist somewhere for CAFE5? Ultimately, I want to be able to compare specific nodes in my tree and see if the same gene families are expanding or contracting.

Thanks in advance!

Hanna