Skip to content

Commit

Permalink
now we only take one query sequence file
Browse files Browse the repository at this point in the history
  • Loading branch information
ekg committed Jun 1, 2024
1 parent 729d9d7 commit d9a0880
Showing 1 changed file with 4 additions and 7 deletions.
11 changes: 4 additions & 7 deletions src/interface/parse_args.hpp
Original file line number Diff line number Diff line change
Expand Up @@ -65,8 +65,7 @@ void parse_args(int argc,
args::Positional<std::string> target_sequence_file(mandatory_opts, "target", "alignment target/reference sequence file");

args::Group io_opts(parser, "[ Files IO Options ]");
args::PositionalList<std::string> query_sequence_files(io_opts, "queries", "query sequence file(s)");
//args::ValueFlag<std::string> query_sequence_file_list(io_opts, "queries", "alignment queries files list", {'Q', "query-file-list"});
args::Positional<std::string> query_sequence_file(io_opts, "query", "query sequence file (optional)");

args::Group mapping_opts(parser, "[ Mapping Options ]");
args::ValueFlag<float> map_pct_identity(mapping_opts, "%", "percent identity in the mashmap step [default: 90]", {'p', "map-pct-id"});
Expand Down Expand Up @@ -219,11 +218,9 @@ void parse_args(int argc,
}
map_parameters.referenceSize = skch::CommonFunc::getReferenceSize(map_parameters.refSequences);

if (query_sequence_files) {
for (auto& q : args::get(query_sequence_files)) {
map_parameters.querySequences.push_back(q);
align_parameters.querySequences.push_back(q);
}
if (query_sequence_file) {
map_parameters.querySequences.push_back(args::get(query_sequence_file));
align_parameters.querySequences.push_back(args::get(query_sequence_file));
}

if (target_sequence_file && map_parameters.querySequences.empty()
Expand Down

0 comments on commit d9a0880

Please sign in to comment.