diff --git a/README.md b/README.md index ce4a6b4..8b24e91 100644 --- a/README.md +++ b/README.md @@ -93,11 +93,11 @@ Full list of options: ```bash $ filterBAM reassign --help -usage: filterBAM reassign [-h] --bam BAM [-p STR] [-r FILE] [-t INT] [-i INT] [-s FLOAT] - [-A FLOAT] [-l INT] [-n INT] [--match-reward INT] +usage: filterBAM reassign [-h] --bam BAM [-p STR] [-r FILE] [-t INT] [-i INT] + [-s FLOAT] [-A FLOAT] [-l INT] [-n INT] [--match-reward INT] [--mismatch-penalty INT] [--gap-open-penalty INT] - [--gap-extension-penalty INT] [--lambda FLOAT] [-k FLOAT] [-o [FILE]] - [-m STR] [-N] [--tmp-dir DIR] + [--gap-extension-penalty INT] [--lambda FLOAT] [-k FLOAT] + [-o [FILE]] [-m STR] [-N] [--tmp-dir DIR] [--disable-sort] optional arguments: -h, --help show this help message and exit @@ -132,10 +132,13 @@ Re-assign optional arguments: -o [FILE], --out-bam [FILE] Save a BAM file without multimapping reads (default: None) -m STR, --sort-memory STR - Set maximum memory per thread for sorting; suffix K/M/G recognized - (default: 1G) + Set maximum memory per thread for sorting; suffix K/M/G + recognized (default: 1G) -N, --sort-by-name Sort by read names (default: False) --tmp-dir DIR Temporary directory (default: None) + +miscellaneous arguments: + --disable-sort Disable sorting of the filtered BAM file (default: False) ``` diff --git a/bam_filter/reassign.py b/bam_filter/reassign.py index cd03771..44ec999 100644 --- a/bam_filter/reassign.py +++ b/bam_filter/reassign.py @@ -672,6 +672,7 @@ def reassign_reads( reassign_scale=0.9, sort_memory="4G", sort_by_name=False, + disable_sort=False, ): dt.options.progress.enabled = True dt.options.progress.clear_on_success = True @@ -1112,7 +1113,7 @@ def reassign_reads( entries=entries, sort_memory=sort_memory, min_read_ani=min_read_ani, - disable_sort=args.disable_sort, + disable_sort=disable_sort, ) @@ -1172,5 +1173,6 @@ def reassign(args): gap_extension_penalty=args.gap_extension_penalty, lambda_value=args.lambda_value, K_value=args.K_value, + disable_sort=args.disable_sort, ) log.info("Done!")