minigraph icon indicating copy to clipboard operation
minigraph copied to clipboard

Published 20 hours ago • verified publisher icon lh3

Incorrect path alignment

Open tobiasrausch opened this issue 1 year ago • 0 comments

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

tobiasrausch avatar Feb 17 '23 15:02 tobiasrausch