-
Notifications
You must be signed in to change notification settings - Fork 10
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
How to control minimap 2 parameters during uLTRA alignment #6
Comments
Hi, There are unfortunately no parameters to do this. I can add it to the next version. You can however do it "manually".
Then I have file
You can manually edit the For instance, if you want to set
|
Thanks for your quick response. It would be great to be able to control the minimap2 flags in a future release of uLTRA. I am a big fan of uLTRA by the way. We have done extensive benchmarking for aligning direct RNA reads and uLTRA outperforms minimap2 in all contexts we have tested. Cheers. |
Ok, I will do that. I'm really glad to hear that! Especially since we did not include dRNA sequencing in our evaluations. |
Thanks, I had one last question. I am wondering what parameters are used in this comparison? Does the algorithm filter for the simplest CIGAR? Or the most aligned bases? |
Good question. The computation is currently simple: Let The, the computation is as follows: I should mention that softclips (at either 5' or 3' ends) are disregarded in the above computations. I observed that minimap2 is very good at softclipping at the right positions and may therefore get a better score according to the computation above, even though they are aligned to identical splice sites (uLTRA extends the alignment out fully in ends). The reason I mention this is that uLTRA prints a table at the end of the alignment stage that states how many alignments were preferred for either aligner. This table is not accurate for low score differences (less than 5-10 score difference) because of the above reason. There is certainly room to explore alternative decision-making on the final alignment. |
Interesting, especially the last part about softclipping at read-ends. When using direct RNA reads, one would expect that there is no need to softclip at the 3' end of the read (except for sometimes a few nt, representing sequencing/basecalling errors) since this represents the bona fide end of the RNA molecule. Cheers! |
Hi,
I understand uLTRA relies on mm2 when aligning to unannotated regions.
How can I control minimap2 mapping parameters?
Cheers,
The text was updated successfully, but these errors were encountered: