You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
@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.
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.
@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
@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 theREADME.md
.Any help appreciated!
See also lh3/bwa#64
The text was updated successfully, but these errors were encountered: