I have generated a gVCF (genomic VCF) with GATK using this command line:

    gatk HaplotypeCaller -R GRCh38.p13.genome.fa.gz -I bam.bam -D dbsnp_146.hg38.vcf.gz -O genome.g.vcf.gz -ERC GVCF --output-mode EMIT_ALL_ACTIVE_SITES

I'm using GRCh38 patch 13, my BAM file is also based on it.

I want to analyze my genotype for rs4633: https://www.snpedia.com/index.php/Rs4633

As I understand from SNPedia, the reference allele is "C". As I have two variants from my mother and my father, I can have either (C;C), (C;T) or (T;T) 

In my gVCF I have the corresponding line:

    CHROM  POS     ID      REF     ALT     QUAL    FILTER  INFO    FORMAT  bar
    chr22   19962712        rs4633  C       T,<NON_REF>     1117.03 .       DB;DP=29;ExcessHet=3.0103;MLEAC=2,0;MLEAF=1.00,0.00;RAW_MQandDP=95148,29

The reference allele in chromosome 22 at position 19962712 is "C". But what exactly does "T,<NON_REF>" mean in the ALT column? How should I interpret this? How can I extract the two alleles (father/mother) from this line and map them to the ones from SNPedia: (C;C), (C;T) or (T;T)?

To ask the question in a different way: how can I obtain both my alleles from the ALT column?

To add to this confusion. Another example. I want to know my genotype for rs1801131. SNPedia says: https://www.snpedia.com/index.php/rs1801131
Possible alleles: (A;A), (A;C), (C;C)

Why the "letters" are different than rs4633? Is it because of A=T and G=C?

The line in the gVCF is:

    CHROM  POS     ID      REF     ALT     QUAL    FILTER  INFO    FORMAT  bar
    chr1    11794419        rs1801131       T       G,<NON_REF>     1193.60 .       BaseQRankSum=0.306;DB;DP=54;ExcessHet=3.0103;MLEAC=1,0;MLEAF=0.500,0.00;MQRankSum=-4.833;RAW_MQandDP=151742,54;ReadPosRankSum=1.170 GT:AD:DP:GQ:PGT:PID:PL:PS:SB    0|1:19,35,0:54:99:0|1:11794400_G_A:1201,0,583,1258,688,1945:11794400:12,7,20,15

So reference allele "T" and "alternative" is "G,<NON_REF>". So what are my two alleles which I can map to the possible alleles from SNPedia: (A;A), (A;C), (C;C)?

My assumption is that REF column is the allele from GRCh38 and ALT is somehow my own, but there should be two - one from mom and one from dad. What I'm doing wrong?


Two things to help you:

> But what exactly does "T,<\non_ref>" mean in the ALT column? 

You have only made a `g.vcf` file with `HaplotypeCaller`. This is a file that is 'poised' to be genotyped, but that has not yet been performed. You need to use the `g.vcf` file as input in the genotyping step of GATK using the [GenotypeGVCFs][1] tool. This will give you the genotypes that you are expecting to see. 

> Why the "letters" are different than rs4633? Is it because of A=T and G=C?

Yes. the MTHFR gene is transcribed on the reverse strand of DNA and SNPedia is showing the alleles corresponding to the mRNA. In the reference sequence these alleles will be the complimentary alleles. from SNPedia:

> rs1801131 is a SNP in the MTHFR gene, representing an A>C mutation at mRNA position 1298, resulting in a glu429-to-ala (E429A) substitution (hence this SNP is also known as A1298C or E429A).


  [1]: https://gatk.broadinstitute.org/hc/en-us/articles/360036351612-GenotypeGVCFs