My data is a VCF file from exome sequencing variant call. I'm not very familiar with the sequencing process and variant calling process. I noticed that there are some missing genotypes, which is recorded as "./." at the GT field. From googling I learned that they're not homozygous reference genotype ("0/0"), but missing calls due to some sequencing ignorance.
In order to make my data "cleaner", I thought it's better to filter out the loci with missing calls if there're not too many of them. I also checked the corresponding "DP" of the loci with and without missing calls. For example:
chr1:123 GT:DP 0/0:2 1/1:4 ./. chr1:234 GT:DP 0/0:10 1/1:11 1/1:20 chr1:345 GT:DP 0/1:40 1/1:37 0/0:78 chr1:456 GT:DP 0/1:7 0/0:23 ./. chr1:567 GT:DP 0/1:34 1/1:39 0/0:58
In the above toy example, there're 3 people and 5 loci. I checked the mean DP of all 5 loci and I found that the mean DP of loci with missing calls (1st and 4th loci) to be significantly lower than the mean DP of loci without missing calls (2nd, 3rd and 5th loci). Is that a coincidence? Or are there any specific reason that the loci with missing calls have lower coverage than the "normal" loci? Thanks!