Skip to content

Commit a8356a2

Browse files
committedMay 14, 2024·
Clarify what -U does
See #424
1 parent ec49586 commit a8356a2

File tree

2 files changed

+3
-2
lines changed

2 files changed

+3
-2
lines changed
 

‎README.md

+2-1
Original file line numberDiff line numberDiff line change
@@ -158,7 +158,8 @@ options. Some important ones are:
158158
specified multiple times. Example: `--rg-id=1 --rg=SM:mysamle --rg=LB:mylibrary`.
159159
* `-N INT`: Output up to INT secondary alignments. By default, no secondary
160160
alignments are output.
161-
* `-U`: Suppress output of unmapped reads.
161+
* `-U`: Suppress output of unmapped single-end reads and pairs in which both
162+
reads are unmapped.
162163
* `--use-index`: Use a pre-generated index instead of generating a new one.
163164
* `--create-index`, `-i`: Generate a strobemer index file (`.sti`) and write it
164165
to disk next to the input reference FASTA. Do not map reads. If read files are

‎src/cmdline.cpp

+1-1
Original file line numberDiff line numberDiff line change
@@ -36,7 +36,7 @@ CommandLineOptions parse_command_line_arguments(int argc, char **argv) {
3636
args::Group sam(parser, "SAM output:");
3737
args::Flag eqx(parser, "eqx", "Emit =/X instead of M CIGAR operations", {"eqx"});
3838
args::Flag no_pg(parser, "no-PG", "Do not output PG header", {"no-PG"});
39-
args::Flag U(parser, "U", "Suppress output of unmapped reads", {'U'});
39+
args::Flag U(parser, "U", "Do not output unmapped single-end reads. Do not output pairs where both reads are unmapped", {'U'});
4040
args::ValueFlag<std::string> rgid(parser, "ID", "Read group ID", {"rg-id"});
4141
args::ValueFlagList<std::string> rg(parser, "TAG:VALUE", "Add read group metadata to SAM header (can be specified multiple times). Example: SM:samplename", {"rg"});
4242
args::Flag details(parser, "details", "Add debugging details to SAM records", {"details"});

0 commit comments

Comments
 (0)
Please sign in to comment.