Add fraction with 20X coverage VA quality measure

Author: phildarnowsky-broad

graphql-api/src/graphql/resolvers/va.spec.ts

@@ -65,7 +65,7 @@ describe('resolveVACohortAlleleFrequency', () => {
     variant_id: '1-123-G-A',
     exome: exomeEsDocument,
     joint: { fafmax: { faf95_max: 0.234, faf95_max_gen_anc: 'amr' } },
-    coverage: { exome: { mean: 0.345 } },
+    coverage: { exome: { mean: 0.345, over_20: 0.456 } },
   }
 
   test('parses a single CohortAlleleFrequency correctly', async () => {
@@ -95,7 +95,7 @@ describe('resolveVACohortAlleleFrequency', () => {
         subcohortFrequency: [],
         qualityMeasures: {
           meanDepth: 0.345,
-          fractionCoverage20x: null,
+          fractionCoverage20x: 0.456,
           monoallelic: null,
           qcFilters: null,
           lowComplexityRegion: null,

graphql-api/src/graphql/resolvers/va.ts

@@ -197,7 +197,7 @@ type Subset = {
   homozygote_count: number
   grpMax?: GrpMaxFAF95
   jointGrpMax?: GrpMaxFAF95
-  meanDepth: number
+  qualityMeasures: QualityMeasures
 }
 
 const GNOMAD_V4_DERIVATION = {
@@ -271,18 +271,6 @@ const resolveVACohortAlleleFrequency = (
     hemizygotes: subset.hemizygote_count !== undefined ? subset.hemizygote_count : null,
   }
 
-  // const coverage = obj.coverage.exome && obj.coverage.exome
-  const qualityMeasures = {
-    meanDepth: subset.meanDepth,
-    fractionCoverage20x: null, // TK
-    qcFilters: null, // TK
-    monoallelic: null, // TK
-    lowComplexityRegion: null, // TK
-    lowConfidenceLossOfFunctionError: null, // TK
-    lossOfFunctionWarning: null, // TK
-    heterozygousSkewedAlleleCount: null, // TK
-  }
-
   return {
     id,
     label,
@@ -294,7 +282,7 @@ const resolveVACohortAlleleFrequency = (
     alleleFrequency: subset.ac / subset.an,
     cohort,
     ancillaryResults,
-    qualityMeasures,
+    qualityMeasures: subset.qualityMeasures,
   }
 }
 
@@ -397,6 +385,17 @@ export const resolveVACohortAlleleFrequencies = async (
     ? [obj.exome, obj.coverage.exome]
     : [obj.genome, obj.coverage.genome]
 
+  const qualityMeasures = {
+    meanDepth: coverage && coverage.mean ? coverage.mean : null,
+    fractionCoverage20x: coverage && coverage.over_20 ? coverage.over_20 : null,
+    monoallelic: null,
+    qcFilters: null,
+    lowComplexityRegion: null,
+    lowConfidenceLossOfFunctionError: null,
+    lossOfFunctionWarning: null,
+    heterozygousSkewedAlleleCount: null,
+  }
+
   const fullSet: Subset = {
     ac: frequencies.ac,
     an: frequencies.an,
@@ -415,7 +414,7 @@ export const resolveVACohortAlleleFrequencies = async (
             confidenceInterval: 0.95,
           }
         : undefined,
-    meanDepth: coverage && coverage.mean ? coverage.mean : null,
+    qualityMeasures,
   }
   const subsets = [fullSet, ...(frequencies.ancestry_groups as Subset[])]
   const cohortsWithoutSubcohorts = subsets.map((subset) =>