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

Adjusting penalty scores in the miniprot splice model #1102

Open
RNA-ROB opened this issue Aug 9, 2023 · 1 comment
Open

Adjusting penalty scores in the miniprot splice model #1102

RNA-ROB opened this issue Aug 9, 2023 · 1 comment

Comments

@RNA-ROB
Copy link

RNA-ROB commented Aug 9, 2023

For minimap2 (versions 2.25 and above), it seems like the miniprot splice model is used for aligning RNA-seq reads by default. According to the miniprot paper, there are fixed penalty scores associated with read alignments against different donor/acceptor splice site dinucleotides (see screenshots below):

Screenshot 2023-08-09 163039
Screenshot 2023-08-09 163030

Would it be possible for users to adjust these specific penalty scores when running minimap2? Thanks!

@lh3
Copy link
Owner

lh3 commented Aug 9, 2023

The miniprot scoring is unfortunately hard coded for now. I thought to add an option to rescale the penalty on the command line but I didn't do that for some reason I forgot. That said, you can still use the old model with option -J0.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

2 participants