Skip to content

Commit

Permalink
Added parameter for distributing quantification
Browse files Browse the repository at this point in the history
of each enzyme by the KOs identified for it
  • Loading branch information
iquasere committed Dec 20, 2023
1 parent fffcc88 commit 3d59b5d
Showing 1 changed file with 6 additions and 1 deletion.
7 changes: 6 additions & 1 deletion keggcharter.py
Original file line number Diff line number Diff line change
Expand Up @@ -36,6 +36,9 @@ def get_arguments():
"-mm", "--metabolic-maps", help="IDs of metabolic maps to output",
default=','.join(keggcharter_prokaryotic_maps()))
parser.add_argument("-qcol", "--quantification-columns", help="Names of columns with quantification")
parser.add_argument(
"-dq", "--distribute-quantification", default=False, action="store_true",
help="Quantification of each enzyme is divided by all KOs identified for it.")
parser.add_argument(
"-tls", "--taxa-list", help="List of taxa to represent in genomic potential charts (comma separated)") # TODO - must be tested
parser.add_argument(
Expand Down Expand Up @@ -669,7 +672,9 @@ def main():
cog2ko_file=f'{sys.path[0]}/cog2ko_keggcharter.tsv',
threads=args.threads,
step=args.step)
data = prepare_data_for_charting(data, ko_column='KO (KEGGCharter)', mt_cols=args.quantification_columns)
data = prepare_data_for_charting(
data, ko_column='KO (KEGGCharter)', mt_cols=args.quantification_columns,
distribute_quantification=args.distribute_quantification)
data.to_csv(f'{args.output}/data_for_charting.tsv', sep='\t', index=False)
if not args.input_taxonomy:
taxon_to_mmap_to_orthologs = download_resources(
Expand Down

0 comments on commit 3d59b5d

Please sign in to comment.