Skip to content

1.0.2

Compare
Choose a tag to compare
@hoelzer hoelzer released this 23 Aug 13:11
· 31 commits to main since this release
0131fc2

What's Changed

  • Hot fix for issue #54, skipping tree calc if no core gene MSA by @hoelzer in #55

Full Changelog: 1.0.1...1.0.2

For diverse input genomes, no core gene may exist for each input species. By default, RIBAP only considers core genes found in all input species for the (optional) --tree calculation. With this update, a new parameter --core_perc is introduced that can be used to lower the threshold (default: 1.0 == 100%). See example below and the RIBAP --help:

    --core_perc         Define how many species are required so that a gene is considered a core gene for tree calculation.
                        Per default, RIBAP will only consider genes that were found in all input genomes (100%).
                        However, this can cause tree calculation to stop when there are no such core genes. 
                        You can lower the threshold to include more homologous genes into the tree caclulation. 
                        The total input genome number will be multiplied by this value and rounded down when the 
                        results is <= x.5 (e.g., 28*0.9=25.2 --> 25) and up otherwise (e.g., 28*0.95=26.6 --> 27).
                        All RIBAP groups that are composed of genes from different species equal or greater this number 
                        will be considered in the core gene MSA and tree [default: $params.core_perc]