minimap2 icon indicating copy to clipboard operation
minimap2 copied to clipboard

Adjusting penalty scores in the miniprot splice model

Open RNA-ROB opened this issue 10 months ago • 1 comments

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!

RNA-ROB avatar Aug 09 '23 20:08 RNA-ROB

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.

lh3 avatar Aug 09 '23 21:08 lh3