Question: What Does Genetype ("0/0", "0/1" Or "1/1") In *.Vcf File Represent?
1
gravatar for lv06025158
3.9 years ago by
lv0602515810
lv0602515810 wrote:

Dear everyone,

Recently I was analyzing some NGS data and genome polymorphysm. Though Pindel, I got the insertions and deletions of NGS against Reference genome. After Pindel2vcf the vcf format files were available. Here is my question: What dose the genetype ("0/0", "0/1" or "1/1") in the vcf file represent separately?

For example,

chr10_irgsp5    2279161    .    CA    C    .    PASS    END=2279162;HOMLEN=9;HOMSEQ=AAAAAAAAA;SVLEN=-1;SVTYPE=DEL    GT:AD    0/0:3,3    0/0:0,0

chr10_irgsp5    2313030    .    CA    C    .    PASS    END=2313031;HOMLEN=10;HOMSEQ=AAAAAAAAAA;SVLEN=-1;SVTYPE=DEL    GT:AD    0/0:1,2    0/0:2,2

chr10_irgsp5    2588340    .    GTA    G    .    PASS    END=2588342;HOMLEN=3;HOMSEQ=TAT;SVLEN=-2;SVTYPE=DEL    GT:AD    0/0:0,1    0/1:8,8

Thanks very much!

Yang

vcf pindel • 7.2k views
ADD COMMENTlink modified 4 weeks ago by alerodriguez40 • written 3.9 years ago by lv0602515810
5
gravatar for matted
3.9 years ago by
matted6.5k
Boston, United States
matted6.5k wrote:

From the VCF (Variant Call Format) version 4.1 page:

GT : genotype, encoded as allele values separated by either of ”/” or “|”. The allele values are 0 for the reference allele (what is in the REF field), 1 for the first allele listed in ALT, 2 for the second allele list in ALT and so on. For diploid calls examples could be 0/1, 1|0, or 1/2, etc.

ADD COMMENTlink written 3.9 years ago by matted6.5k
1

Good call - I was about to clarify that is what the genotype means in general, but I am actually a bit confused in this specific case: I don't know what pindel would give 0/0 indel calls.

Nevertheless, you can confirm if the indel is heterozygous (0/1) or homozygous (1/1) by checking the alignment (in addition to checking the reliability of the call)

ADD REPLYlink written 3.9 years ago by Charles Warden4.9k
3
gravatar for Peixe
3.9 years ago by
Peixe510
Spain
Peixe510 wrote:

Just point out that if the bar is | means that the genotypes have been phased, that is have been assigned to which chromosome comes each allele. If the bar is like / means it has not been phased.

ADD COMMENTlink modified 3.9 years ago • written 3.9 years ago by Peixe510
0
gravatar for Charles Warden
3.9 years ago by
Charles Warden4.9k
Duarte, CA
Charles Warden4.9k wrote:

I don't remember off the top of my head, but I think you should be able to figure it out by looking at the pindel output files (prior to creation of the .vcf file)..

There are alignments for every indel that show every read supporting the indel (which you can compare to the total number of reads).

Alternatively, you can always check your alignment using an indexed .bam file in IGV. This is probably a good idea no matter what: unless you are working with a small indel (which would have been identified using an ordinary variant calling tool like VarScan anyways), structural variant calls are hard to make and will often include false positives.

ADD COMMENTlink written 3.9 years ago by Charles Warden4.9k
0
gravatar for alerodriguez
4 weeks ago by
alerodriguez40
alerodriguez40 wrote:

GT, genotype, encodes alleles as numbers: 0 for the reference allele, 1 for the first allele listed in ALT column, 2 for the second allele listed in ALT and so on. The number of alleles suggests ploidy of the sample and the separator indicates whether the alleles are phased (‘|’) or unphased (‘/’) with respect to other data lines (Fig. 1).

see details:

Bioinformatics. 2011 Aug 1;27(15):2156-8. doi: 10.1093/bioinformatics/btr330. Epub 2011 Jun 7. The variant call format and VCFtools. should help.

ADD COMMENTlink modified 4 weeks ago • written 4 weeks ago by alerodriguez40
0
gravatar for alerodriguez
4 weeks ago by
alerodriguez40
alerodriguez40 wrote:

Bioinformatics. 2011 Aug 1;27(15):2156-8. doi: 10.1093/bioinformatics/btr330. Epub 2011 Jun 7. The variant call format and VCFtools. should help.

ADD COMMENTlink modified 4 weeks ago • written 4 weeks ago by alerodriguez40
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: 1543 users visited in the last hour