view format_taxon_list.py @ 2:9916308301da draft default tip

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/checkm commit 3ae2bd72f789518e95ef65b97e0e5ac90165e113
author iuc
date Mon, 02 Sep 2024 13:51:37 +0000
parents 356839cd89d2
children
line wrap: on
line source

from pathlib import Path

from checkm.taxonParser import TaxonParser
from checkm.util.taxonomyUtils import taxonomicRanks

if __name__ == '__main__':
    tool_data_dp = Path('tool-data')

    # get all available marker sets
    taxonParser = TaxonParser()
    taxonMarkerSets = taxonParser.readMarkerSets()

    # create a table per rank
    for rank in taxonomicRanks:
        rank_fp = tool_data_dp / Path("%s.loc.sample" % rank)
        with rank_fp.open('w') as rank_f:
            rank_f.write('# File generated by format_taxon_list.py script\n')
            rank_f.write('# taxon\tdescription\n')
            for taxon in sorted(taxonMarkerSets[rank]):
                markerSet = taxonMarkerSets[rank][taxon]
                numMarkers, numMarkerSets = markerSet.size()
                rank_f.write('{taxon}\t{taxon} ({markerSet.numGenomes} genomes, {numMarkers} marker genes, {numMarkerSets} marker sets)\n')