-x sr is silently identical to -x adap - #63
Closed
nh13 wants to merge 1 commit into
Closed
Conversation
MB_F_ADAP was set unconditionally, overriding the conditional added in r234. Also drop two statements made dead by the reset in mb_opt_preset().
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
-x srgives identical output to-x adap.mb_opt_presetsetsMB_F_ADAPunconditionally two lines below the gated set that r234 added when it separated the two presets, so both end up at flag0x4048andsris unreachable. The man page hasadapas-xsr --adap=yes, which is not currently true.Deleting the stray line restores
sr. Two no-ops go with it:mb_opt_presetresetsopton entry, so themb_opt_resetinmb_opt_initand the~MB_F_PEin thelrbranch do nothing. Drop those two if you prefer the one-line fix.On 1M HG002 WGS pairs,
-x adapis byte-identical to before;-x srmoves 0.11% of pairs, 69% of them at MAPQ <= 1.If you meant the two presets to stay identical, the gated line should go instead and the man page needs correcting.