math domain error

[Y4nkey]

Hi @Vini2,

I am running metacoag on the latest version and receive this error:

2024-08-29 15:15:36,323 - INFO - Welcome to MetaCoAG: Binning Metagenomic Contigs via Composition, Coverage and Assembly Graphs. 2024-08-29 15:15:36,323 - INFO - Input arguments: 2024-08-29 15:15:36,323 - INFO - Assembler used: spades 2024-08-29 15:15:36,323 - INFO - Contigs file: contigs.fasta 2024-08-29 15:15:36,323 - INFO - Assembly graph file: assembly_graph_after_simplification.gfa 2024-08-29 15:15:36,323 - INFO - Contig paths file: contigs.paths 2024-08-29 15:15:36,323 - INFO - Abundance file: /srv/scratch/z5363929/single_assembly/kris_data_bp_deconstruction/284/trimmed/metaspades_binning/284_metaspades_abundance.tsv 2024-08-29 15:15:36,323 - INFO - Final binning output file: /srv/scratch/z5363929/single_assembly/kris_data_bp_deconstruction/284/trimmed/metaspades_binning/284_metacoag_metaspades_output 2024-08-29 15:15:36,323 - INFO - Marker gene file hmm: auxiliary/marker.hmm 2024-08-29 15:15:36,323 - INFO - Minimum length of contigs to consider: 1000 2024-08-29 15:15:36,323 - INFO - Depth to consider for label propagation: 10 2024-08-29 15:15:36,323 - INFO - p_intra: 0.1 2024-08-29 15:15:36,323 - INFO - p_inter: 0.01 2024-08-29 15:15:36,323 - INFO - Do not use --cut_tc: False 2024-08-29 15:15:36,323 - INFO - mg_threshold: 0.5 2024-08-29 15:15:36,323 - INFO - bin_mg_threshold: 0.33333 2024-08-29 15:15:36,323 - INFO - min_bin_size: 200000 base pairs 2024-08-29 15:15:36,323 - INFO - d_limit: 20 2024-08-29 15:15:36,323 - INFO - Number of threads: 8 2024-08-29 15:15:36,323 - INFO - MetaCoAG started 2024-08-29 15:16:12,468 - INFO - Total number of contigs available: 2464476 2024-08-29 15:16:35,313 - INFO - Total number of edges in the assembly graph: 475216 2024-08-29 15:16:36,582 - INFO - Total isolated contigs in the assembly graph: 2116227 2024-08-29 15:16:36,582 - INFO - Obtaining lengths and coverage values of contigs 2024-08-29 15:16:59,029 - INFO - Total long contigs: 180672 2024-08-29 15:16:59,030 - INFO - Total isolated long contigs in the assembly graph: 88740 2024-08-29 15:16:59,030 - INFO - Obtaining tetranucleotide frequencies of contigs 2024-08-29 15:25:46,681 - INFO - Scanning for single-copy marker genes 2024-08-29 15:25:47,384 - INFO - Obtaining hmmout file 2024-08-29 15:25:47,384 - INFO - Using marker file: /home/z5363929/miniconda3/envs/metacoag/lib/python3.9/site-packages/metacoag/metacoag_utils/auxiliary/marker.hmm 2024-08-29 15:54:56,204 - INFO - Obtaining contigs with single-copy marker genes 2024-08-29 15:54:56,369 - INFO - Number of contigs containing single-copy marker genes: 6783 2024-08-29 15:54:56,369 - INFO - Determining contig counts for each single-copy marker gene 2024-08-29 15:54:56,372 - INFO - Initialising bins 2024-08-29 15:54:56,372 - INFO - Matching and assigning contigs with single-copy marker genes to bins Traceback (most recent call last): File "/home/z5363929/miniconda3/envs/metacoag/bin/metacoag", line 10, in <module> sys.exit(main()) File "/home/z5363929/miniconda3/envs/metacoag/lib/python3.9/site-packages/click/core.py", line 1157, in __call__ return self.main(*args, **kwargs) File "/home/z5363929/miniconda3/envs/metacoag/lib/python3.9/site-packages/click/core.py", line 1078, in main rv = self.invoke(ctx) File "/home/z5363929/miniconda3/envs/metacoag/lib/python3.9/site-packages/click/core.py", line 1434, in invoke return ctx.invoke(self.callback, **ctx.params) File "/home/z5363929/miniconda3/envs/metacoag/lib/python3.9/site-packages/click/core.py", line 783, in invoke return __callback(*args, **kwargs) File "/home/z5363929/miniconda3/envs/metacoag/lib/python3.9/site-packages/metacoag/cli.py", line 266, in main metacoag_runner.main(args) File "/home/z5363929/miniconda3/envs/metacoag/lib/python3.9/site-packages/metacoag/metacoag_runner.py", line 1116, in main run(args) File "/home/z5363929/miniconda3/envs/metacoag/lib/python3.9/site-packages/metacoag/metacoag_runner.py", line 577, in run ) = matching_utils.match_contigs( File "/home/z5363929/miniconda3/envs/metacoag/lib/python3.9/site-packages/metacoag/metacoag_utils/matching_utils.py", line 133, in match_contigs prob_cov = get_cov_probability( File "/home/z5363929/miniconda3/envs/metacoag/lib/python3.9/site-packages/metacoag/metacoag_utils/matching_utils.py", line 53, in get_cov_probability (cov2[i] * math.log(cov1[i])) - math.lgamma(cov2[i] + 1.0) - cov1[i] ValueError: math domain error