code: def split_by_category(biom_cols, mapping, category_id):
      
      columns = defaultdict(list)
      for i, col in enumerate(biom_cols):
          columns[mapping[col['id']][category_id]].append((i, col))
  
      return columns
docs: Split up the column data in a biom table by mapping category value.