coordinate overlaps between pairs
hrrsjeong opened this issue · comments
I found that there are SD pairs that overlap genomic coordinates. Also, there are some SDs pairs that have exactly the same coordinate (self-alignment) (examples below)
h2tg000001l 3733143 3748938 h2tg000001l 3743164 3751506 HG00733_hap1.masked:HG00733_hap1.masked 54.2 + + 13607 13796 232M3D287M17I1M20N46M1D191M29I25N456M9I273M6I16M2I99M3I105M20N4D24S174M1I93M1D45M6D1M7S84M9I7S6I22N125M28S1M28N7M24S24N31M1D13M1I43M2I116M3D183M3D19S19N68M101I39S663I19S458I40S268I23S463I43S199I20S155I30S42I183S399I7S51I23S499I23S103I22S115I29S142I24S150I22S199I22S415I197M2I104M19S19N90M22S22N6M26S1M28D31M317N2M10I312S37M1I67M23I2M23N64M1I29M11D21N4D1M34S227M29S29N150M1I18M68D65S66M7D7S13M5I38M26S1M26N139M20S20N11M1I10M2D315N546M4I210M5D7N8D1M20S238M25S2M23N2D125M3D106M716I20S408I22S126I28S7I24S393I19S369I19S90I22S66I317S107I1M23S104M21S21N232M29S29N149M1D50M23I2M23N136M1I26S27N139M20D20S23M327N4D315S123M24I24N399M4D142M X=2.9;ID=51.3
h2tg000001l 1336074 1337493 h2tg000001l 1336074 1337493 HG00733_hap1.masked:HG00733_hap1.masked 45.7 + - 1282 1367 73M1I125M18I31M7D63M13I26M107S25M107N74M17I43M5D58M16I47M8I64M8D49M16D42M23S15M5I22M23N13M7S10M17D15M7N105M13D66M7I30M18D126M1D75M X=33.3;ID=12.4
Yes, this is an issue but unfortunately I won't be able to fix it anytime soon (it's not major issue in our pipelines).
You can maybe add simple post-processing to filter these out.