I ran samtools mpileup & then vcfutils varFilter to call snps on a set of 10 illumina samples against a number of reference genomes. I'm now looking at the results, and I am having trouble understanding the information I've been given per sample.
The VCF file lists my 10 samples as the last 10 columns, and has them in the format GT:PL:GQ. I understand that those show me the genotype, phred-scaled likelihood of the genotypes, and the overall quality of the genotype assignment. But I am confused by the per-sample results. Here are some specific questions I have:
1) I do not expect every sample to have coverage over every loci called. What should I see in the field for a sample that does not contain that loci? I suspect it looks like this:
0/1:0,0,0:3
since that seems to appear a lot in my output. But I want to be sure.
2) The VCF documentation says that for diploid calls (my samples are actually metagenomic, but the caller seems to consider them diploid so I am going with it), the 'PL' portion of the sample info fields is giving me the phred-scaled likelihood of the genotypes in the order:
AA AB BB
Where A is the ref allele & B is the alternate allele. But I am seeing data that looks like this:
1/1:145,21,0:26
From the documentation, I had thought a GT value of 1/1 meant that it was seeing, for that given sample, both alleles being the alternate allele. And from the PL documentation, I would have expected the large PL value (145) to be in the THIRD position. I thought the first position in that string was reserved for REF/REF genotype. Could somebody clarify the order of the PL values for me?
Thanks,
John Martin
p.s. The documentation I am referring to is the VCF format v4.1 manual page at the 1000 genomes website (not sure if its OK to give URLs in this forum, but its pretty easy to find via google if you want to see it).
The VCF file lists my 10 samples as the last 10 columns, and has them in the format GT:PL:GQ. I understand that those show me the genotype, phred-scaled likelihood of the genotypes, and the overall quality of the genotype assignment. But I am confused by the per-sample results. Here are some specific questions I have:
1) I do not expect every sample to have coverage over every loci called. What should I see in the field for a sample that does not contain that loci? I suspect it looks like this:
0/1:0,0,0:3
since that seems to appear a lot in my output. But I want to be sure.
2) The VCF documentation says that for diploid calls (my samples are actually metagenomic, but the caller seems to consider them diploid so I am going with it), the 'PL' portion of the sample info fields is giving me the phred-scaled likelihood of the genotypes in the order:
AA AB BB
Where A is the ref allele & B is the alternate allele. But I am seeing data that looks like this:
1/1:145,21,0:26
From the documentation, I had thought a GT value of 1/1 meant that it was seeing, for that given sample, both alleles being the alternate allele. And from the PL documentation, I would have expected the large PL value (145) to be in the THIRD position. I thought the first position in that string was reserved for REF/REF genotype. Could somebody clarify the order of the PL values for me?
Thanks,
John Martin
p.s. The documentation I am referring to is the VCF format v4.1 manual page at the 1000 genomes website (not sure if its OK to give URLs in this forum, but its pretty easy to find via google if you want to see it).
Comment