Notice:
If you happen to see a question you know the answer to, please do chime in and help your fellow community members. We appreciate your help!

Test-drive the GATK tools and Best Practices pipelines on Terra


Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.

Position of Indel event based on the REF

rfuentesrfuentes PhilippinesMember

Hi,
How do I know based on the REF and ALT column of a VCF file the actual position where an indel event happened?
I usually see an event that occur in the 2nd base. For example,
REF ALT
GGCGTGGCGT G,GCGCGTGGCGT --deletion; insertion of "C"
ATTT A,ATT --deletions

But I saw a VCF format documentation(http://samtools.github.io/hts-specs/VCFv4.2.pdf) allowing a different case:
GTC G,GTCT --deletion; insertion of "T" at the end

How is UnifiedGenotyper formatting the indels? Does it differ
from the one used in the 1000Genome Project? Thank you!

Roven

Best Answer

Answers

Sign In or Register to comment.