Skip to content

Clarify genome size error and update usage documentation #478

@adthrasher

Description

@adthrasher

Description of the bug

When running with --genome GRCh38, I get the following failure.

WARN: =============================================================================
--macs_gsize parameter has not been provided.
It will be auto-calculated by 'khmer unique-kmers.py' using the '--read_length' parameter.
Explicitly provide '--macs_gsize macs3_genome_size' to change this behaviour.

Both '--read_length' and '--macs_gsize' not specified! Please specify either to infer MACS3 genome size for peak calling.

The issue is that --macs_gsize is rejected when using --genome. So the user has to specify --read_length instead. The warning message should be updated. The usage documentation should also be updated to include documentation on both parameters.

Command used and terminal output

Relevant files

No response

System information

No response

Metadata

Metadata

Assignees

No one assigned

    Labels

    bugSomething isn't working

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions