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
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,
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!