[go: nahoru, domu]

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

Can minimap2 really do glocal alignment ? #143

Closed
tseemann opened this issue Apr 2, 2018 · 4 comments
Closed

Can minimap2 really do glocal alignment ? #143

tseemann opened this issue Apr 2, 2018 · 4 comments
Labels

Comments

@tseemann
Copy link
Contributor
tseemann commented Apr 2, 2018

@lh3 In your recent blog post you stated minimap2 does global alignment.

How can I force this?
Like bowtie2 --end-to-end?
I'm currently using minimap2 -x sr -a and still getting soft-clipping eg. 16S53M31S
I can't seem to see anything in the -h or the README.md.

Any help appreciated!

See also lh3/bwa#64

@lh3 lh3 added the question label Apr 3, 2018
@lh3
Copy link
Owner
lh3 commented Apr 3, 2018

Bwa-mem uses global alignment in aligned regions to reconstruct CIGAR at the end. It doesn't support glocal alignment.

Minimap2 implements a similar but not identical algorithm, but this difference might emulate glocal alignment to some extent. To try that, you may use a large "--end-bonus" like

minimap2 -axsr --end-bonus=100

It worked on a toy example I tried. Please beware that you often need to filter out poor glocal alignments as they may be not real.

@tseemann
Copy link
Contributor Author
tseemann commented May 1, 2018

Thanks @lh3 I'll see how it goes.

@tseemann tseemann closed this as completed May 1, 2018
@sdwfrost
Copy link

Dear @tseemann, how did it go?

@tseemann
Copy link
Contributor Author
tseemann commented Feb 18, 2020

@sdwfrost i never really got back to it. in fact i completely lost track of this. thanks for the reminder!
Another lesson in map minimap2 having more info than minimap2 --help

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

No branches or pull requests

3 participants