bwa-meth Overriding hardcoded -B values is possible from command line

Hi @brentp ,

While testing this, I noticed that it's possible to pass -B <N> at the command line and override the -B 2 value set in bwa-meth.py. Is this intended behavior to allow users to override -B, or should it always be fixed to -B 2 for things to properly function?

Thanks, Eric

Feb 12 '24 19:02 edawson

Hi Eric, I don't really remember why that's there, but is it a problem to allow a user to set it? I guess it could mess up other thresholds if other penalties/scores are not modified?

Feb 12 '24 19:02 brentp

I don't know enough about the flag to know for sure; just wanted to clarify the intended behavior. I think it's even possible to pass -B multiple times and override the user override.

It might be clearer to add -B to the argparse arguments with a default of 2 and a single int value as input?

Feb 13 '24 01:02 edawson

is there a "not" missing from your original issue? I don't see how it's possible to set -B: https://github.com/brentp/bwa-meth/blob/master/bwameth.py#L358-L361

are you asking to make that available? I'd rather not unless you have compelling reason especially since bwa-meth has been unchanged for so long.

Feb 13 '24 08:02 brentp