Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

MAPQ calculation #63

Draft
wants to merge 2 commits into
base: main
Choose a base branch
from
Draft

MAPQ calculation #63

wants to merge 2 commits into from

Conversation

jdidion
Copy link
Collaborator

@jdidion jdidion commented Oct 3, 2023

Implementation of single-end MAPQ algorithm from BWA.

  • Constrains MAPQ between 0-60
  • Sets MAPQ to 0 when the best score is <= the suboptimal score
  • Sets MAPQ to 60 when there is no suboptimal score
  • Otherwise uses a formula based on the ratio between sub-optimal and optimal scores, with various scaling constants.

It seems to work in that it generates different values in the expected range. I need to do more work to compare with Minimap2 to see that the results at least make sense in a relative sense.

@jdidion jdidion requested a review from nh13 October 3, 2023 21:39
@jdidion jdidion marked this pull request as draft October 30, 2023 17:19
@jdidion jdidion self-assigned this Oct 30, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant