Skip to content

"No nesting structure in HOGs defined in OrthoXML" #51

@nromashchenko

Description

@nromashchenko

Version: 2.1.0

I am trying to make a database from orthoxml from the fastoma's example data:

omamer mkdb --db testdata.h5 --log_level debug --orthoxml expected_output/FastOMA_HOGs.orthoxml --sequences in_folder/proteome/*.fa 
 --species_tree in_folder/species_tree.nwk 

this breaks with

...

WARNING: No nesting structure in HOGs defined in OrthoXML.

...

    raise tables.NoSuchNodeError(
tables.exceptions.NoSuchNodeError: Path not found in database (/ChildrenHOG)
...

And in the code, the warning in question has a todo note to check things more carefully if there is no nesting structure. However there are subHOGs in the example file, for example:

<orthologGroup id="HOG:0000004_1" taxonId="1">
  <score id="CompletenessScore" value="1.0" />
  <property name="OMAmerRootHOG" value="HOG:E0990677" />
  <property name="TaxRange" value="inter2" />
  <geneRef id="1000000002" />
  <orthologGroup id="HOG:0000004_2" taxonId="2">
    <score id="CompletenessScore" value="1.0" />
    <property name="TaxRange" value="inter1" />
    <geneRef id="1002000004" />
    <geneRef id="1001000004" />
  </orthologGroup>
</orthologGroup>

I also see the same error if I try toy files from the orthoxml standard (I tried ex1 and ex2. They require changing "geneId" to "protId" which is a separate issue btw).

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions