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

Incorrect path alignment #91

Open
tobiasrausch opened this issue Feb 17, 2023 · 0 comments
Open

Incorrect path alignment #91

tobiasrausch opened this issue Feb 17, 2023 · 0 comments

Comments

@tobiasrausch
Copy link

Hi Heng,

Extracting a chm13 subsequence and aligning it back to the pre-built graph results in an incorrect path alignment (deviating from rank 0 nodes) with lots of small InDels:

gfatools gfa2fa -s -0 chm13-90c.r518.gfa.gz > chm13.fa
samtools faidx chm13.fa chr1:1025000-1035000 > read.fa
minigraph --vc -k 15 -w 10 -t 24 -cx lr chm13-90c.r518.gfa.gz read.fa 

This subsequence is uniquely mappable in chm13 using minimap2 at MAPQ 60.

Thanks, Tobias

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

No branches or pull requests

1 participant