@@ -973,17 +973,17 @@ void align_or_map_paired(
973
973
Timer nam_timer;
974
974
auto [nonrepetitive_fraction, nams] = find_nams (query_randstrobes, index );
975
975
976
- if (nams.empty ()) {
977
- query_randstrobes = randstrobes_query_rescue (record.seq , index_parameters);
978
- std::tie (nonrepetitive_fraction, nams) = find_nams (query_randstrobes, index );
979
- }
980
976
statistics.tot_find_nams += nam_timer.duration ();
981
977
982
978
if (map_param.rescue_level > 1 ) {
983
979
Timer rescue_timer;
984
980
if (nams.empty () || nonrepetitive_fraction < 0.7 ) {
985
981
nams = find_nams_rescue (query_randstrobes, index , map_param.rescue_cutoff );
986
982
details[is_revcomp].nam_rescue = true ;
983
+ if (nams.empty ()) {
984
+ query_randstrobes = randstrobes_query_rescue (record.seq , index_parameters);
985
+ std::tie (nonrepetitive_fraction, nams) = find_nams (query_randstrobes, index );
986
+ }
987
987
}
988
988
statistics.tot_time_rescue += rescue_timer.duration ();
989
989
}
@@ -1098,18 +1098,17 @@ void align_or_map_single(
1098
1098
Timer nam_timer;
1099
1099
auto [nonrepetitive_fraction, nams] = find_nams (query_randstrobes, index );
1100
1100
1101
- if (nams.empty ()) {
1102
- query_randstrobes = randstrobes_query_rescue (record.seq , index_parameters);
1103
- std::tie (nonrepetitive_fraction, nams) = find_nams (query_randstrobes, index );
1104
- }
1105
-
1106
1101
statistics.tot_find_nams += nam_timer.duration ();
1107
1102
1108
1103
if (map_param.rescue_level > 1 ) {
1109
1104
Timer rescue_timer;
1110
1105
if (nams.empty () || nonrepetitive_fraction < 0.7 ) {
1111
1106
details.nam_rescue = true ;
1112
1107
nams = find_nams_rescue (query_randstrobes, index , map_param.rescue_cutoff );
1108
+ if (nams.empty ()) {
1109
+ query_randstrobes = randstrobes_query_rescue (record.seq , index_parameters);
1110
+ std::tie (nonrepetitive_fraction, nams) = find_nams (query_randstrobes, index );
1111
+ }
1113
1112
}
1114
1113
statistics.tot_time_rescue += rescue_timer.duration ();
1115
1114
}
0 commit comments