Question: How Does Left-Alignment Work For Indels? Does It Address All Complicated Issues With Indel Representation In Vcf File?
0
gravatar for Duck
4.8 years ago by
Duck20
United States
Duck20 wrote:

I can vaguely guess its meaning but want to know, at least, how it works in different situations. Besides, after left alignment, can we write an indel in VCF file in a unique way?

Thank you!

indel vcf ngs • 2.6k views
ADD COMMENTlink modified 2.1 years ago by Biostar ♦♦ 20 • written 4.8 years ago by Duck20

Did you see indel left/right alignment and http://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_sting_gatk_walkers_variantutils_LeftAlignVariants.html ? "Note that this tool cannot handle anything other than bi-allelic, simple indels. Complex events are written out unchanged." That might answer the second part of your question.

ADD REPLYlink modified 4.8 years ago • written 4.8 years ago by Michael Dondrup44k

Thanks, Mike!

After reading your links, I learned that how left-alignment works in a deletion occurring in a repeat region. However, I guess there are more situations. For example, we could have an insertion in a repeat region: ATGATG(GCG)GCGGCGTAGTAG, where (GCG) is an insertion. How would this be handled by left-alignment?

ADD REPLYlink written 4.8 years ago by Duck20
3
gravatar for Len Trigg
2.4 years ago by
Len Trigg1.2k
New Zealand
Len Trigg1.2k wrote:

In terms of providing an actual answer here for people arriving from google, the short answer is no, left alignment does not address complicated issues with indel representation in VCF files. It does help somewhat, but there are plenty of cases where simple normalization tools do not produce a unique variant representation. (As a side note, left-aligning simple indels is an unstated convention in VCF, and in contrast, HGVS guidelines state that indels should be right-aligned in the context of the genetic transcript).

Since there is no easy answer to representing variants uniquely, the state of the art is to perform variant comparisons using haplotype-aware tools such as rtg vcfeval, hap.py, or vgraph.

ADD COMMENTlink written 2.4 years ago by Len Trigg1.2k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1881 users visited in the last hour