jzook / testHG008curation

This repository is primarily for the Genome in a Bottle Consortium to curate structural variants in HG008
https://www.nist.gov/programs-projects/cancer-genome-bottle
1 stars 0 forks source link

[chr16:89166442-89166520][DEL][SVLEN=78bp] #263

Open jzook opened 1 month ago

jzook commented 1 month ago

chr16 89166442 Minda_51 N <DEL> . PASS SVLEN=78;SVTYPE=DEL;SUPP_VEC=PB_ID_29443_2,ONT_ID_63605_2,PB_severus_DEL2591,PB_Sniffles2.DEL.2CEMF,ONT_Sniffles2.DEL.2D0MF,ONT_severus_DEL2900,PB_d_154

https://v2.genomeribbon.com/?session=https://42basepairs.com/download/s3/giab-data/ribbon-json/ribbon-hg008-cov-bedpe.json&locus=chr16:89166442#ribbon

aganezov commented 1 month ago

Too small of a variant to see in SplitThreader on coverage levels.

Screenshot 2024-10-11 at 09 32 58

Visible in ribbon: Screenshot 2024-10-11 at 09 40 14

Variant location is in the repeat region, making mapping challenging even for long reads. Across consistently mapped deletion signatures the chr16:89,166,442 appears indeed the most consistent.

Most likely del length is indeed 78bp, and mapping inconsistencies within technologies and across them can be explained by the repetitive nature of the underlying genomic context.

Screenshot 2024-10-11 at 09 44 45

From the fraction of affected reads with deletion signature in both PB and ONT data, as well as partial read phasing, it appears that the deletion is clonal (i.e., affecting all reads coming from a single haplotype).