# Call Rate after GT revising

**URL:** <https://discuss.hail.is/t/call-rate-after-gt-revising/2287>\
**Category:** Hail Query & hailctl\
**Created:** [October 7, 2021, 2:18pm UTC](https://discuss.hail.is/t/call-rate-after-gt-revising/2287 "2021-10-07T14:18:39Z")\
**Posts on this page:** 1\
**Page:** 1

<div class="post-metadata">

**Author:** ![zhouhufeng](https://yyz2.discourse-cdn.com/flex036/user_avatar/discuss.hail.is/zhouhufeng/32/535_2.png) [@zhouhufeng](https://discuss.hail.is/u/zhouhufeng)\
**Post date:** [October 7, 2021, 2:18pm UTC](https://discuss.hail.is/t/call-rate-after-gt-revising/2287/1 "2021-10-07T14:18:39Z")

</div>

I found a problem of my hail QC pipeline,  
if I revised GT using the following scripts,

```auto
vds = vds.filter_entries((vds.DP > 400)|(vds.DP < 10), keep = False)
vds = hl.variant_qc(vds)

```

and then calculate the call rate, as below,

```auto
callratelist = [hl.agg.filter(vds.Pop == pop, hl.agg.fraction(hl.is_defined(vds.GT))) for pop in custompops]
vds = vds.annotate_rows(
   callrate=hl.Struct(**dict(zip(custompops, callratelist))),
   lowestcallrate = hl.min(hl.array(callratelist))
)

```

Then all the original missed call genotype will be labeled as NA and will not be count in the call rate calculation. I am wondering if there is any way that we can improve the methods by forcing the call rate calculation by counting both the original missing in original genotype and the revised newly reset missing in the genotype?

Thanks very much.