Question: Interpretation of mixed VCF record
0
gravatar for rutujagtap08
2.9 years ago by
rutujagtap0810
India/Pune/Savitribai Phule Pune University
rutujagtap0810 wrote:

Hi I am confused about interpreting mixed VCF record. I am having 2 nucleotides in the ALT field of the VCF which are separated by comma. The GT field shows that it is homozygous. How should I interpret this? I am pasting the row in the file just as example. Thanks in advance.

chrI    509 .   G   A,T 119 .   DP=60;VDB=8.840780e-04;AF1=1;AC1=2;DP4=0,0,54,3;MQ=20;FQ=-178   GT:PL:GQ    1/1:152,151,0,172,87,151:99
snp • 884 views
ADD COMMENTlink modified 2.9 years ago • written 2.9 years ago by rutujagtap0810

rather a question then a tutorial.

ADD REPLYlink written 2.9 years ago by dago2.5k

Thank you Wouter. Would you please share any written material/link about interpreting genotypes.

ADD REPLYlink written 2.9 years ago by rutujagtap0810

Please use ADD COMMENT to answer to earlier posts, as such this thread remains logically structured and easy to follow.

Among similar resources, the samtools vcf specifications will put you on the right track. You could have found this yourself by googling.

ADD REPLYlink written 2.9 years ago by WouterDeCoster41k

https://samtools.github.io/hts-specs/VCFv4.2.pdf

ADD REPLYlink written 2.9 years ago by Pierre Lindenbaum123k

Our links are particularly similar ;-)

ADD REPLYlink written 2.9 years ago by WouterDeCoster41k
0
gravatar for WouterDeCoster
2.9 years ago by
Belgium
WouterDeCoster41k wrote:

The individual shown here is homozygous (1/1) for the first alternative allele and as such has a AA genotype for this position. The alleles for the second alternative allele would be presented with a '2'. If the individual would have a AT genotype it would be 1/2, if the individual would be heterozygous GT the genotype would be 0/2.

ADD COMMENTlink written 2.9 years ago by WouterDeCoster41k

However if the position is homozygous then why do we get a mixed VCF record

ADD REPLYlink written 2.9 years ago by rutujagtap0810

It was probably part of a larger vcf which was filtered? I don't know how you obtained this file or how it was generated.

ADD REPLYlink written 2.9 years ago by WouterDeCoster41k

Yes the line posted above is just one of the records of my vcf file "S2_vcf". Below is work-flow I used to obtain the file.

bwa mem genome.fsa S2_1.fq S2_2.fq > S2.sam -t 4
samtools view -S S2.sam -b -o S2_pe
samtools sort S2_pe Sorted_S2_pe
samtools mpileup -g -u -f genome.fsa Sorted_S2_pe.bam > S2_var
bcftools view -c -v -g S2_var > S2_vcf
ADD REPLYlink modified 2.9 years ago by WouterDeCoster41k • written 2.9 years ago by rutujagtap0810

yes I am using the above commands to generate the vcf file but what is the reason to have mixed vcf record for homozygous position. Or is there an issue with the workflow.

ADD REPLYlink written 2.9 years ago by rutujagtap0810

I'm not familiar with using samtools mpileup for variant calling (since I use GATK), so I'm not sure.

ADD REPLYlink written 2.9 years ago by WouterDeCoster41k
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: 2246 users visited in the last hour