Question: 2/50 inidividuals show missing data at almost all positions in final vcf generated using HaplotypeCaller
0
gravatar for shreyasibiswas88
11 months ago by
United States
shreyasibiswas8830 wrote:

Hello all,

I used HaplotypeCaller from GATK to generate VCF from 50 individuals. My reference is Rattus norvegicus and 48/50 individuals are Rattus norvegicus. The remaining two individuals are Rattus rattus species. The WGS for these two individuals have been mapped on the same reference, Rattus norvegicus, as the other 48 individuals. In my final VCF output almost all positions in these two Rattus rattus individuals is reported as ./. I checked using IGV and I find read coverage at all of those sites (avg coverage is around 15x) I have not applied any filter while calling the variants. in fact, variants are called from the same positions when only these two individuals are being used as input files with the same reference.

I have posted this on GATK support. But wondering if anybody in this community has experienced the same and can point out why I am I seeing ./. for only these two samples?

Thank you.

Edit: I calculated the percentage of missing data and found 96% for one individual and 99% for the other individual indicating that the samples are being read.

ADD COMMENTlink modified 11 months ago • written 11 months ago by shreyasibiswas8830

Did you get a response yet on the GATK Support Forum? Can you provide the link to there?

Can you paste some commands that you've used prior to the generation of the VCFs?

Thanks!

ADD REPLYlink modified 11 months ago • written 11 months ago by Kevin Blighe42k
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: 908 users visited in the last hour