From 2ba3a15d9f49a675f38641e652418667aabc6161 Mon Sep 17 00:00:00 2001 From: meganshand Date: Thu, 12 Sep 2024 15:53:44 -0400 Subject: [PATCH] Adds QD and AS_QD emission from VariantAnnotator on GVS input (#8978) * Updating VariantAnnotator to emit QD and AS_QD on GVS input * clean up * addressing comments --- .../tools/walkers/annotator/QualByDepth.java | 14 +- .../allelespecific/AS_QualByDepth.java | 6 +- .../VariantAnnotatorIntegrationTest.java | 45 +- .../VariantAnnotator/gvsCallsetWithNoQual.vcf | 3394 +++++++++++++++++ .../gvsCallsetWithNoQual.vcf.idx | Bin 0 -> 114136 bytes 5 files changed, 3450 insertions(+), 9 deletions(-) create mode 100644 src/test/resources/org/broadinstitute/hellbender/tools/walkers/annotator/VariantAnnotator/gvsCallsetWithNoQual.vcf create mode 100644 src/test/resources/org/broadinstitute/hellbender/tools/walkers/annotator/VariantAnnotator/gvsCallsetWithNoQual.vcf.idx diff --git a/src/main/java/org/broadinstitute/hellbender/tools/walkers/annotator/QualByDepth.java b/src/main/java/org/broadinstitute/hellbender/tools/walkers/annotator/QualByDepth.java index 8646a63e21d..ab9d8fc6b71 100644 --- a/src/main/java/org/broadinstitute/hellbender/tools/walkers/annotator/QualByDepth.java +++ b/src/main/java/org/broadinstitute/hellbender/tools/walkers/annotator/QualByDepth.java @@ -7,6 +7,7 @@ import htsjdk.variant.variantcontext.VariantContext; import org.broadinstitute.barclay.help.DocumentedFeature; import org.broadinstitute.hellbender.engine.ReferenceContext; +import org.broadinstitute.hellbender.exceptions.GATKException; import org.broadinstitute.hellbender.utils.MathUtils; import org.broadinstitute.hellbender.utils.Utils; import org.broadinstitute.hellbender.utils.genotyper.AlleleLikelihoods; @@ -57,7 +58,7 @@ public Map annotate(final ReferenceContext ref, final VariantContext vc, final AlleleLikelihoods likelihoods) { Utils.nonNull(vc); - if ( !vc.hasLog10PError() ) { + if ( !(vc.hasLog10PError() || vc.hasAttribute(GATKVCFConstants.RAW_QUAL_APPROX_KEY)) ) { return Collections.emptyMap(); } @@ -72,7 +73,16 @@ public Map annotate(final ReferenceContext ref, return Collections.emptyMap(); } - final double qual = -10.0 * vc.getLog10PError(); + final double qual; + if (vc.hasLog10PError()) { + qual = -10.0 * vc.getLog10PError(); + } else { + try { + qual = vc.getAttributeAsInt(GATKVCFConstants.RAW_QUAL_APPROX_KEY, 0); + } catch (NumberFormatException e) { + throw new GATKException("Error at: " + vc.getContig() + ":" + vc.getStart() + " when parsing " + GATKVCFConstants.RAW_QUAL_APPROX_KEY + ": " + e.getMessage()); + } + } double QD = qual / depth; // Hack: see note in the fixTooHighQD method below diff --git a/src/main/java/org/broadinstitute/hellbender/tools/walkers/annotator/allelespecific/AS_QualByDepth.java b/src/main/java/org/broadinstitute/hellbender/tools/walkers/annotator/allelespecific/AS_QualByDepth.java index f8bb3b3cea9..d0e15fa5ac0 100644 --- a/src/main/java/org/broadinstitute/hellbender/tools/walkers/annotator/allelespecific/AS_QualByDepth.java +++ b/src/main/java/org/broadinstitute/hellbender/tools/walkers/annotator/allelespecific/AS_QualByDepth.java @@ -5,8 +5,6 @@ import htsjdk.variant.variantcontext.GenotypesContext; import htsjdk.variant.variantcontext.VariantContext; import htsjdk.variant.vcf.VCFCompoundHeaderLine; -import htsjdk.variant.vcf.VCFHeaderLine; -import htsjdk.variant.vcf.VCFInfoHeaderLine; import org.apache.commons.lang3.StringUtils; import org.broadinstitute.barclay.help.DocumentedFeature; import org.broadinstitute.hellbender.engine.ReferenceContext; @@ -78,7 +76,9 @@ public List getRawDescriptions() { public Map annotate(final ReferenceContext ref, final VariantContext vc, final AlleleLikelihoods likelihoods ) { - return Collections.emptyMap(); + // first vc is used for the annotation and the second vc here is used just to get the alleles, so in this case we can pass the same vc for both + Map annotation = finalizeRawData(vc, vc); + return (annotation == null ? Collections.emptyMap() : Collections.singletonMap(getKeyNames().get(0), annotation.get(getKeyNames().get(0)))); } /** diff --git a/src/test/java/org/broadinstitute/hellbender/tools/walkers/annotator/VariantAnnotatorIntegrationTest.java b/src/test/java/org/broadinstitute/hellbender/tools/walkers/annotator/VariantAnnotatorIntegrationTest.java index 6f3a87c60e7..7a5268a7d90 100644 --- a/src/test/java/org/broadinstitute/hellbender/tools/walkers/annotator/VariantAnnotatorIntegrationTest.java +++ b/src/test/java/org/broadinstitute/hellbender/tools/walkers/annotator/VariantAnnotatorIntegrationTest.java @@ -29,10 +29,7 @@ import java.io.File; import java.io.IOException; -import java.util.ArrayList; -import java.util.List; -import java.util.Map; -import java.util.Set; +import java.util.*; import java.util.stream.Collectors; import java.util.stream.Stream; @@ -364,6 +361,46 @@ public void testDeNovo() { Assert.assertEquals(lociWithHighConfidenceDeNovo, new int[] {10130767, 10197999}); } + //GVS vcfs have QUALapprox, but no QUAL. This tests that we can still calculate QD and AS_QD without QUAL (using QUALapprox and AS_QUALapprox) + @Test + public void addQdToGvsVcf() { + final File inputVCF = getTestFile("gvsCallsetWithNoQual.vcf"); + final File outputVCF = createTempFile("output", ".vcf"); + final ArgumentsBuilder args = new ArgumentsBuilder() + .addVCF(inputVCF) + .addOutput(outputVCF) + .addInterval("chr1:66506-66507") + .add("A", "QualByDepth") + .add("A", "AS_QualByDepth"); + runCommandLine(args.getArgsList()); + + Map> expectedAnnotations = new HashMap<>(); + Map expectedAnnotations66506 = new HashMap<>(); + Map expectedAnnotations66507 = new HashMap<>(); + + // bi-allelic site + expectedAnnotations66506.put("QD", "4.23"); + expectedAnnotations66506.put("AS_QD", "4.23"); + //make sure that the AS_QUALapprox doesn't change from the input + expectedAnnotations66506.put("AS_QUALapprox", "0|110"); + + //multi-allelic site + expectedAnnotations66507.put("QD", "6.29"); + expectedAnnotations66507.put("AS_QD", "[9.00, 9.26]"); + expectedAnnotations66507.put("AS_QUALapprox", "0|396|537"); + + expectedAnnotations.put(66506, expectedAnnotations66506); + expectedAnnotations.put(66507, expectedAnnotations66507); + + + final List outputVCs = VariantContextTestUtils.getVariantContexts(outputVCF); + for (final VariantContext outputVC : outputVCs) { + Assert.assertEquals(outputVC.getAttribute("QD"), expectedAnnotations.get(outputVC.getStart()).get("QD")); + Assert.assertEquals(outputVC.getAttribute("AS_QD").toString(), expectedAnnotations.get(outputVC.getStart()).get("AS_QD")); + Assert.assertEquals(outputVC.getAttribute("AS_QUALapprox"), expectedAnnotations.get(outputVC.getStart()).get("AS_QUALapprox")); + } + } + private static String keyForVariant( final VariantContext variant ) { return String.format("%s:%d-%d %s", variant.getContig(), variant.getStart(), variant.getEnd(), variant.getAlleles()); } diff --git a/src/test/resources/org/broadinstitute/hellbender/tools/walkers/annotator/VariantAnnotator/gvsCallsetWithNoQual.vcf b/src/test/resources/org/broadinstitute/hellbender/tools/walkers/annotator/VariantAnnotator/gvsCallsetWithNoQual.vcf new file mode 100644 index 00000000000..f9f63bcff26 --- /dev/null +++ b/src/test/resources/org/broadinstitute/hellbender/tools/walkers/annotator/VariantAnnotator/gvsCallsetWithNoQual.vcf @@ -0,0 +1,3394 @@ +##fileformat=VCFv4.2 +##FILTER= +##FILTER= +##FILTER= +##FILTER= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##GATKCommandLine= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##high_CALIBRATION_SENSITIVITY_INDEL=Sample Genotype FT filter value indicating that the genotyped allele failed INDEL model calibration sensitivity cutoff (0.99) +##high_CALIBRATION_SENSITIVITY_SNP=Sample Genotype FT filter value indicating that the genotyped allele failed SNP model calibration sensitivity cutoff (0.997) +##source=SelectVariants +#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT ERS4367795 ERS4367796 ERS4367797 +chr1 66506 . T TAA . NO_HQ_GENOTYPES AC=1;AF=0.167;AN=6;AS_QUALapprox=0|110;CALIBRATION_SENSITIVITY=.;QUALapprox=110;SCORE=. GT:AD:GQ:RGQ 0/1:21,5:99:110 0/0:.:40 0/0:.:40 +chr1 66507 . T TATAATATATATTATATATTTATAGAATATAATATATATTTTATA,A . . AC=1,2;AF=0.167,0.333;AN=6;AS_QUALapprox=0|396|537;CALIBRATION_SENSITIVITY=0.9884,0.5557;QUALapprox=396;SCORE=-0.6019,-0.4314 GT:AD:GQ:RGQ 0/1:18,5,0:99:396 0/2:9,0,5:99:132 0/2:12,0,14:99:405 diff --git a/src/test/resources/org/broadinstitute/hellbender/tools/walkers/annotator/VariantAnnotator/gvsCallsetWithNoQual.vcf.idx b/src/test/resources/org/broadinstitute/hellbender/tools/walkers/annotator/VariantAnnotator/gvsCallsetWithNoQual.vcf.idx new file mode 100644 index 0000000000000000000000000000000000000000..861cbce998263f1b88546d923fe5c6b28efc9b24 GIT binary patch literal 114136 zcma)_%g#5)ao%4i>)eQdm~%dGr&BCBvJ5%0VXp{UW-P(77?6~MY;y-$1z83`E-bf^ zRi2)i?)uG4Jzx0Y6Gi^JrVrKC)%{deSO5CwKl$4~`S|$w*ZSXo{N3-r`uwxcKKtu$ zzxw9e&;I$_-+uXrfB5Xre)rvP|M>Tx{o|M4{nKaPe)F5pzWeIi?>_tHt8c&l<2S$g z>f6u0{^lP)`}=Rc{_-Dw_lIx4``vec{O+sIe*4w$fB*Mi{ox?#u7K{^qm4`SP3Jeffv)e%OBa$A9_FfBefoe)pgL_}!oC-#&iz zzyFo`^YQUN{_B7LKmV(77l+i&`hVcPb2nTDbLHt4sZ-^Y2_ z+i&AXzi-1n&vX1X>9=m^C;iyy=8t~s=W!cAw|w+#w{`Qr#<%N7KkwtD-^TXs_My7^ zaa+cv{kDJThjE|Qer>DF!<*B(D=G+n>%eZS%Jpf=K(~L)+tRN~A0r&P zk7d%Ieu)td{l~Pd-M(+Z9Suvv>$ho%4*spUxO7ve8EozT-N(ER!#efx(_eh_%e-I+ zU;O)@537c|{7e7EHxkaT|5;NscbLb%lB_{84&4}1(kQl|=!bRL5|TYgmStPI*4bJ~ z-v!0gi9v6Q9u&hM8E7JWb+sRYLN#bE-y~y@Ox?;Etd&fl=qYtWrI>@FUwTPqYo?Kq zud^gISU-O8XTRvWP6I$K^E~b``O`_#b=$uTb#{hsmpHVB=BpNwv!Z~Jc^S7JODdW|HDl3_jeS9a9{208QFdEP z5E&_<5d)v=pg5J&CcPa^EtLk7*>9^fA1Yd_)Xc%4$d5q$$apni`>Uqb<$a$G9&0n~PJ~JajC;38j8E^`**Tncn@F z`h&qbK)pIL^_p9q9&m(;3MCnZi|=WfiXQzKM(zJq=Gs@oIhlzFOnn*Hzzr|HI>t4n zZ50vNx2&l{0=DfvDbsKaYh04Arqa}?^^z4X$2<>X?8&W*%hB4J2(ZpTSCwyejOD9+ z;0dYKI2HyMjnpM9UP__1gm=+=$mzh7Q>}9D%D9IShr?FXpVPU$$|g-mbhRaA^6n)} z>*01NMje&jjcvy=+i_56ZEc5!(>X5rs;Cbki)6mmGOke@NcLGNlU~dwb1$t42P7iA zl|dI2m2{-6Pwk8ERXuNoA4@+cyx$8y#ty)y4SFvlo7eq4vhkQqZ9H$@kEe79w7hya zj>n9Ri3p5i>Bg`mGd3QYpzm-Tsy!4+(^tJV`gdDpXEjCRu{676g$Z?*ulhtkkINxA zf_(3`mP2*xq?HsMN>jq06qRlWRwhO3F|KsThfIkZNn6{YJESt-Z5E4kapV`(ER7!FnI7j7n_l0KVD`OAohU zyQ_^86`zNub(gK__2qD!jryhE?i&5V(}(JiE=qO@*+NR2GS9$#A#Pq@9*mv^x8F_n z2hTJPz@It|s;#i5?=nG=C8#thtehO9Wp_~FB%{)MF4n1ZXIZA?II&L0Qrx$7d{4)E zOvf&3S61n<^5>0ku#ICrHQ}};QnSt{>trjr&E*NR=OS%`PHste=E9w2*prjJKwufa zENe+u=fazDS%FIu0$o$e0+qv06J8vF3{}fB)UbI2Yqy6RlRRp&aGIybRDMtz46Gl=t1PQ`7r< zvFR+|R2=;8CDl@Ub2*e=<--K&ZRL`Pvs@Sj8*VMvVjGU~aglt7I%)kAOONuo8=xp# zv1b1$mhF`5qn6LL`AOpkF0AE?#O~~#z8-^p)vT~K_KXWraVxLw%kYd%trgea4}-oB z?$z?8c2CF0o2~ZZ?qSO+v#Y&N=B_q7_4}fR&+h9agrid{Z?h+hTOa7LmAB#Ra0uIO zEiV}6tiSJMB}V&y_KR6&dZDdz_V%S8V?hkH`>lK0@}Vzk^=0oltC#9rVw)+Nb zSUmOnwe@e6nLpUm$}9V1LmO~Lt>tAezt}r=kiJ@8&R9qJ(9NTK3>B-_^0KX8{E;W+ z-F`V>oaNgz;XGcAf#o!I)-Ut|;80AL$M|diWy~P|>i6ZycGlm^#1o42GJVQXam|mM z5wzvB?`!=+1Usxpt^HR2GR-)lwC|gE7w1Xvr$1!8xH?Ld)6h9N!@M z4*6{5jsNRTt4J$v^SjS0tg*$cUidrK@3%qb+BUtVc%O>fsZEbCSPZXK-=m zmoIB7R_;T4Ed8)Q8aUdo^Dpm{RscG1?z?;CkF3d3a45g}-^RBV7aPj2@7w&ypn>z- zq5&4?giZu!2W5cKu`tx{WoF*G*1=VqLv)~;q4zq$y^ewaW~G5+fH8rSmT2o>tG^Cl z6L7}v)5_}%=xC3lsfvFjhoszTmBE{|VaR{$UX$NluEAIH>t5wwJM#*E)@AtKx;DW= zzYu^{UWYHYjh0v&UMp|&H*eIh#Xx_sr}bY)B)9f6RbOjQ|6X2);Ye7$CbQ z6kPkh3Wx@<19NSs0+$!EF2G6KcMPzV*ZCQ=g23QbUguZhFHdaiUoCI@t8C(PR)1xH z>}LZYSNp!UufxFv0k7r7ZhNjkox-7@T6vXk2|vw&t-MaKyvs6=z5Cz3*GWE(2krC5 ze(6`tZ!53zF?REY_20_t^z@^4$XK3xm7mr<9l$YtwfyiZ6a}>QblSAyIy>@z4iMi~ z-nc6_UZ}sz7pMhyF<0l#4y&i{uI?51b8>&(2RMX8aOfwwwkpBaPtgBSLojG39E=)T z6WgXENSmPkI$+AIay20Ur%6l*>J;{?+>kNKsG)L6$1VEb8$mI29mAAmUcnD+;#}qh z&YgS^M(E-gq#FMGX9|zd*4An$>yR~E1g!=+FzGCRSr7!)KD!YqQbkHx-s z@=-2WLYaD;ZQGSbAblNdY6Nw%OK=I)Ve5my>a`;aOkux^J~WBcfr>NHnyAvC0F4nR z`lz9<3PmO?v86|O<>Ra)4Y+WfDYv#}wI#@`5QD4)#YlcilWdUCnSCWF_KK6Bcq}2& zji8M|e9Zb#%&A8Ws$UfEKD3jOe3dxYS5LT3yR`xZKmb?-K@}FkhVmzxxXy+OKB%ME z6u^Q$G=j=Vg#$@kKyjVVbY#1#m{vI{urM1zWyG*Z1B2Ic+&bi9uFzV2Y*=ht;H35hIqeI*FfhnSjkiDD!s3)M- z2Zj1kkGAqUmx{3DrdtnKYi(fM=7lod2&!^Wh}6J%$3{@6vX^QYAkvMXO_dPpaw9l* z-Quu67_IFLo`$}Sp-rWnVqDp!g##tm3S_Td;p`deF;G^TOXdd!&;$q&!PSK(urhLB zT{eOy3>pIV)F^KPD;pvXc|)~w%8zz-K^hQBSmkFAil^n2wnbj`STaFSXG(^^9SSg7 zs90aoFFk_%8i72=XBTAohx<6n`^qJ4m{<|fg}R$TJYPKj;J1wmCD zIs*=v$snlnP;8k%J3&yVTKnaKZ9NF8upEc&Sob10m$fW{Cl|CHAg)n9R4y#HsI9-f zCaU<~kBJ&Kc}_xiweqGmCdJrI_BZ0v3xr>LQ_vDg&=TZ5oYr2-Lo#A{8|sg&Yac#f?E&Ge{q?34|&IQ`Eh zp&EGXxMrewNUtTMGalrzme(aE!w>@w>S;f0doqq8+@>B4Tw32q`%7SA3j@49=r^3@O0`ZvcjS4no6uk%(Dn7zMUXQ2y5gC4*NA_hT z+PJu+5Z_uuuNY{km3pk&OYG1&QF6TTv~clhKVrU*HndGI(67W>XAn7t6IgZT?{p}ZLaQOt8Pai zR-@>?A1%cjQb(Z{fSE4Mq(ke~fFd{rkfV#&ij>N(1l=s~nGn>eC~F$>Sha%iVjY9f z6L#VtXj8R!d!jaEgse4)kCeF9$vQ~&VWmH>5)|kCDGf5%z+~I{&`f?Z8Oo^j9%cz{ zJ8Rfh-AzddDko*rmXD|rT$)4sO07U#f>rL)95O{+<4NK?m?687BF)dDiEv-F4Ub=_0~bGi+qEu;*|Mjbt=1gMsP%J z7m**j?glLDR2?DP=bGqfyl+zl+u5x9E$=>*V6Ahi%1(IOOiu@c*2H!6Qx6bx!f*ceIGr6(h;$&-~xENtWZB5kK)~QT}f;NcOM1^c}5f;?f2-;+Y z*-6WeUf|Y5TX%~>bVGNvCh93(0wpaQ!;B{CUL!ZJy!6qD)>5u+Kp%n@R zieZ(Yj9*HFh#zY8AHMp{*Z)Nk;P@s8%D_(oxoEQ*O0yS%0`^+JqK2ZfP6Bz15`trY zx(H}TYBd~t!$rWg7zBm47l96q?>-#LzxDx;J8CEk;UtiwIqPH*iw)i7MM^eQ3E)+;}%>HfRh3;@7X+*_t}JY z;$$uk<$aA}&t~&b-WS1yeYW-CP~I28^H#~RHV@@|5oilZItiKN(J-0_PRfIwjUr@H zvz%(DBzX)iG3IfQaG~0Wahq3$W1EVaWhLFjd{WJ(5 zld{<-O@Na;;YJPQBo0ME5S%O*c#|eVCS~VISbg|sE=H}m#nhTO4AB<>&s9MXoRk$L zSvZG5{;EO2cJEG>V8W^&Hk{+5iQuFnev(m?u)?bbZRd%|9hT~=54wPw^Z}d{@=XZ9 z$#h*DL2qXhgMy^qH5>-klR!t7L>>-X_DLXopA7S`hF?waW@KYjiBhkeROD*Xhv1~p zXEIskx=!V!?kFcUgiI>>DI*9@=8N~RpI)O7AW52lp492?O>o$e&l(iMmkjeT9bQeW zFg`X$9ai8wp{NF^-VR zn+GNcz)M@TU*SuJYxWs&QU1OKq8fgaqZ%R;<4Wn zb6D3ZL8eL%^V!9y=MAi5_dc#Yp9DHkCZjkEY$t&ZjLB>s7PgB(ar{Xq$;oVO9r0CH zA7rFSY6wmava-Dh`Y2&QCZsn$lIDB-LRT;7O6jXs1- zP782y5QI!Bl||M>aIyp%CwvH*R2)=RLvT`gFy5+En6C{}X1MpERR+dJP%dY``*7IM zszd>cw_y_HJTKf7nOqMO*I9#5YBE{nDsWxhU7~8`Uup!y1|%JQ2u}9Ty9{R^f|Ck4 zN?2X)rqn(t>hjG8pq!;k)`g6qt%hFuKB*yOQlm&%eS}tBtX5=bQbWk3?iXhS!AV8l zXMF%C6?~O65u7aH5eY$X@@{@^lX2Ln&SA)_mdG4nj?Av_ZBL;{jg_*`B%L74(H-)4 zF9CBDM4BuNz?>4pk<<{Z)Fq8<6kw$cRY?uOO5M!OMiH!3m}D|rfH}J5n9Mg}j)XE9 zM!6hy$@T3jMw>0*oJr=O_b`sYqDv|jaG4N^TPkJI*$m9vk6_QPG;+f zr8@}*Wtm9&5Ud=Z_jJ!DVtq>wCX)f2qogkhAIe>+sv+jY%pv<(BT#r$G78`vxo#3p z9x*-FFcmrb9wu;(076C(LOIE+p41RRdAt13ChM>WUVV_yIjJEysUt?RFo1J(Jt)yJ zz&XWn^KJq-M+V|#vVe1REKWFi#HyU5kncZP7r;4E=UEfMN$tG}tB-3k7o!-8(l`m6 zqt%^o@`zPAYbf!cZ$8{{PQobQoH2vea?LfF)P2Nc;Q;4May=$}h}EszhY72X7}9eT z3PDIj4k%~Pg_fj=kVzR@;ecwaEH|RBwO1~UBvnu-c_Xh?0i1ISH4B1R-{oS!yA9wR z1^s0Nai%Ouev%r%OM1bg51}cw`z8`mE=pH*uJECxiQuKKStiUb_p`4WW*q z2r~k3Qde)YP6R8JV&px}^QO(ELv$7A)d^i6Y?Co4v>+jf1wJr1r;3!KgB9{x58T}Dhe3816Yc*!siXa2P_(%IjH7VaE5 z(z6T73cws)E={HzFh{o&lXU@@qX4)>>;ZF1Dw{-Vfpa8t zSsP&GGt1VYgO0HOb24V8-fNRVoWnT1Ab?PgnTi3<85v2_Y9P*$ujDNauCZD-dz0l0m?I65FzU#Da7{Opakj|<%+bL)nG3jc24TkcWWk-Ivsl&#>^YeM za0ASFWd%Ft;>dw;jzU*-6R8Hw(bCHZf|D{CCtB(Ta}uMo@N}Dtl5zUY2f!R<(n}Zx zn4`Oq$LZ}$()=8cA5--FCoTKeMV?-=(9gMSK1S{nLc-Nq-6Ri`9 z5PC0+Gidsn&M~o!ok7tjfqV#A4FPkMQzj#b?R(Uf-q&oMG0#^MnKwor+BOPV=il>v z#zdbrteIG!5z8lm!Zi|JD?m5e;2o>FI?u4`tBMSKpCQc`fg)w!;{?vp{qKZRz&T3r zkW43Vj=b-0T`Nf6Mj6!bgs1Ri`ZJpGW zy0?aqm^A`9QxZl2=iKrDwN3)(+%9kjK?tP|TkrLKMlD~nsoVSsr;f1Xi$WHFH!Bgh zrX2WhC36NHSCPt=FM-aGoZpoeuW?Gw^@LLZI&wNEeIU^3vdI9@QQUK~d~ZM}0elJ? zk4`eJPn!$?ok?cMHzxsflsF=p3jm#xS0Je&gi_v^gaANCcl+L(?+$;CZeeCM1o&3| z!;By_rBF&TisKshDG#r#!^g@2(2+!Cos8wJsFY-J66j>!nEkl>amnN(uip_utSg8y zn|)#k(cRA=$ZT)AVig2>0^&Vc1Ui{vYzNR$tVJT#1Ukz3moyPFDFZ;lD8#KvQ6{qq zoTEcXwzLRzp7~~v>15!V!f)Qa1kh3L?PM4LI=YOW^pZel&a^dfj*^;WdJaHGmjn}` z1JKcIC5jM0=aqo*U<5om%7~hs5rA}bc`WNhh^1}?BohjxBTgmE2GUXbsO*>lsH01z z3AX`tbRNqHV%t&HkE|0Rmy1y2TO{Dqku^18_K~FRTEMy%n$!>*xWXkef?%c67$s5% zs3U7i!fQYs9YYd65b8XWe;*P-sIw&2k{wV-3pi1xh+NabHPI1(I=UvE5#T5!mH*}> zp-vWezSE<_OK5H91JqH%x`fq$Ix@#3tOnGPE%e6fKmYZ=`uBPC)7;ba{P>6J^;W&_pvN;((!R?~wIvi2R`1CNv$J_yx2+y0*0huUp&@~rY zMHhX|+ySA;cjt>0-p7F*7)>M;#0i(8{E=qPY|RvB+y;#22V#VruA{&;hwT-}DKN z1HF_)t)@?*Rn(eMVn8v|cz05uxsdT&jl0-Bak>Tz!0|`eC5VFKvYQ&!gGjdWN?JZG zfIdo->xe>+>!yd{6A+MO&PXOv5VOwQm$?Ir5RmtT+sq%T7Zf~7}T1jU*ePj`v)j!OfA$cV3W zDT#7sDn2!cn7o$qB)blYT%Vv-9?#ogbkHi=I#Z6z1{=*uiy0ztoJqg$9S6sChb^xO zm{S3TjtCrQ0*+gAz?^y=Dpx&nC}4C3y~^|%L=Wp=FcX+l#?*XqDqpxBT4{7 z;P}G_4%#+83Bp9yPPyAv#bODGi^1IdvP$ za2(z<5xG8vjEl&z03269T+Sz8PG#OV(*n#XBEte;P7#^81LjnQ62k&uP7%49Ama`r zNKM9SaTq~4bxGu!BgYLY)6Fu}WKOKo?;vuz04_JZg=Z!hm{XCpMu`D)>VRPsD=?=5 z?2M=&*ib}9#^Hq%kvTsBbL#S?kpc!DDiv5RzQCNS$&3z|Q~8CA_yTi^$Z#B(^Y+Zh zZEpeQ)SV17I(q5WySFm*6ELUlf4Q0<z?`~!>4>m%Dlc3NP+XARUYxoAd1tW?(*d_t~g`VWPqsxLr9dfK#@8ma?L@HFGj`( z#5q0_8DVtf_&st3MC7=!&?wn)bX=kO=>*GSN}^oo6sSO=e0x>WViM)+PI0gxGR7BO z;U-b8=@gSsA~ScuoVqaOSOAXe2CLx{IPQtaadRG5Gy^(*j|qqg2If?DLDO+yP7#@D z0p@(>DMWMR_&sg~MCABAcO;0QRc2$@SBO4XV2Q|agIknX5jB~!ihwM zgy|Dx+=YFV=Q->B`a#(YO-+z-y*+P;$Z@BAi8&wF8bY3kIrZjU-W)NfQb8CZa{L}d zjpoR4<8o5eKlpTy8AVOxxI54gbLys^={PxVoP0Vtqff@VrZC+glJzBTj+j$c1w#ZG z7w9xZ;P@lEPVkAC^N}bFMBq4|3b`!=VovSUrcaRZM@9fN2aXR0DQeVUrU~E^Xa|uK z6=F`EZt{*3bMk$P8w-dzb-*(f z7Sn=Uh{6j-j}vq1nqS&+UFAS?=By%{IEaijLd?F{j>x$>)xk(-VQ?k1J4U z4jg}k_JarU0gC2~ zXIZXP5V>^+Kc_UbSqMl+BBFd#91&Y9qWt_Mzcl*fRxn~~WkGU%BF7Ed7>PMMmFq3?bvhoG#pIBIeZfXVV-w z?ymUAEgdY-9qGI#Vopy4t)d5|O-OeTtj!060*Tu~JyG82s0TE*Gs_5wXZ zf|^Y9gsvcy$kjv{H`b7enDdc`COA&axwu^lv9)@v)l4uk=i+t*Vop7}n~MrD=i+t* zVon#nHxY9_^7VrS3ZnA6<} z)P;SF&RiUvaGz72Iiu`}IUlbqVszxV+Y#XBRC8tzCgya(i4**sO2KP{j+pb26*)Ri z%=yR%2qI{e+3g6#oJxq1F9c#v<+?CL&?>VrzD)3QK9fSDOf4F{dYj zjC&&J@!6Ph6lWLW`$%ZtMyK=^B+B=e8GcUjDW8^^m~%FJi{42>b7n`-qZ1^`5ztJ` zIh*^O%J9LO^3^gEb1J`KDo;xMgqn=WMA7q|(Von!*I>XPYNAMh<$Z@l6 z%*32}ZzE>`F{jILH4}4sBFebgTa<7Q%rtw8^1P5J*S^XaLn0&N@N+6Pmtz67iX|e) z4Y-+!IbA5PvM;0KCQfW7=6qyhk7Z8GsYRGiF#Mc)#={Yj<7SuA``GBX@pI1bb1Da< znLA=mZCR${;P@lAH&_6zVu^@3^%Pm&95JW%e?tTr*Ga+hULAm*G5+L(zsg$45F=;wSqZ5UEOKc@=~pXujx zq2Wqv3}(8mjF|J0Qx<)KR&hb*Gcl(pBF9bQDs6G-xFy2Q`FON5I!??v83d$fjnJGW zqKuoMU_A}Wn#{T*=A6tWy_uL(CZt@Rh&dn6w}A!ZxbaBn4O?{Fc)N9EB9R##F{kX7 zdB=%4JrOzXwpU`#$y`Cy`4b&CaXmW8lF01A^m9I*KM0N!bILHCvw)b>6OrS#ChBqH zH&>oKjLs5~sh`S-Dt46&_9XtX<3Ll-PIL6 zCCZv|W-jz|D#=+o!3!~`OX?#72%0khUZ#}qLj;by>of~7r@a5Bkc1A~QXt zz<5j-`Z-5)QC?mF)MSaMRV)!XZq9TI{hW`NG-85@ITeANMxjdm{U&)l?CM9it=IHo3Y+?}6@IhBRQ)I^S(1YZj=r^{8Tj6mSHv1BQ?(DxytjJx9+ zG3Ve;V#J*8rm=!u*e5f0^mDpQoC?5UO@`z2b1Jn)J}tzYgE7Ia#GIap95To6I9LyxsB?8XN0M%<0nPuJm(yBFeZqwCmv>%$+47$BlPV$?Z^+iR)R3 zIR_Khvl4T<`~@rhoVvQ6PYW@pC!&nonkeH2NUZd8Dv_7z6FF{9_A4>xre!eG5|QI3q;MtXRFI$P6J^|_#?+NN_Q}W- zF{fg{OmmcROGJ)45J=3~8zNmOM90lVuZ@0AWx!86u1k)r$q`Yjm?Q@qF{e`Zo94)I z6J@dybM}^q95=TJbYm2wGniA)5|PNTfSB{~(r!B!Y{ZZy=i`la)I^S(?0Xyi zoO&WW=QuH^EN6y@TEzgajhIufhM1bjapNV?TYccTBO>PP&2{9Bm{SkZ=Y1mP^hD&i zt%+L25>dvD=Vqgyvp4u|Bj)r(PXMAYL31nTuxOt1lg8!@LI2g>J; z_c`^vv?022+!0a64V>H-;J14Z#J1eQN|q+F{iDGv9+Fv957s^&(_*%`pK5p#MXaNOmm)BC3Eli@fq zXJ=CSZ^WFPaTIODoSpI9Y{Z-{iQPs&XJ=xym6ijeGpKAM=5*=sHu^a|5jk#~Bgai> zxE@+z$IZ0R&)FFexG}cY6XC4V84#$qX4t3vFsgUuL1ggVM$Fln_@|AS(>;Q&M{?Ml z;S*zPyVVhq;|AYt#GIW8eBX#UI|GxH><1mUMC7<7BFBxBZzJY(k9#YX5&M*zaW-O3 zUlV2A;Jb~O)7Jz&-kDSv8}D;^q8rD}>5#Ftoe3t}h&ek0OqHvO%rt0&nA6up88^4) zHe$}scyTxSIXeTE^%fR7Zi%pSb_Ogf@fT|{@_n9~vk$G6eM)9vtc$~%!CMt5M&Wpi6TFsE*`8pR6CsVs_)C^)WkqK+syuInvE z(*biT6>N_0c3@83H!&;#=2Svpa|{9I)U}P=q`L!ij_zLm4nL>f!f`c)R#ED6qg8-8 zm2BQH6PQ!?16)nuxQl$zyTVP=>3)D=CNQU-L2yK&$Mv8>ZmQgYIrS1m?&{s)=Ty+Q znLA)k5gB+LwrT_&G&nBomlZ56U|h1jqGKjp;Zrr=p8Z$ALNRv#C3B z29v@7|vxvZP6K*9U$atTd1$!O7g5%18mG=plQ(=c@bika7D0DS}<9)8|)hBS= zxWPmOjvH^1h`@2-sKiQ5rN~L2psPmGl4l3PM7MZz60hIk&#SbP6>e{0>`@? z0jW7~+&rZrqTslD07FEb zp;h$Kg3&6(oaW%Cnu6oXXP+uo5e3JUVBTmIVoq84^TkQbDgU@30>`)9_o_aHj4Ow= znHFMB#en3^5p%A&Yp&Nxwy7Rv!!{i!<}|6GRTDVAx)n^!x#p&C)dU$gXI~LP#@GC) zEh2DS*BlI=z;Ru@Fhr2?H4mHZWibtA>gG^h6EUZ5r5d8xIhE7XObaomuK1doAmb)o zr#cRfn{=Hb0>>A(ZP3rDYwV^^;P{f8anu|*Zi$FF&6C`HN6e{w4tbv#TdS-dh6pll z@iU``!+%tD}_bIJp{)hBS=6M^IA(5{*wvI8CRC#d|H4xdy`L0H9^Lw+>Sq)^dqyF_~Ia6SRs+!lG9~$?@FPD{mh+#&+U2e*QWIfp!$Of^BqWjxHKfSA+dFObYZ z#`Vl~UK25=c_v&nf#XA-NlHYJag)D5MC77MGB{m8%qcUZX^tGv@0+VRa@-P8#&hGaY9hyTmy?Ldai_%S=N!!8 zK}T84oxxH%6oV)a4VTI!$IZDw*``^O(JI88gF)1~AB386zm)m}j+;xLA_B+FeNGX9 z<1!oNqQcl(bHP$If#b3W>)n`@D(2^=>s%8Q5`H^={; zn6o$baowhi2{u>j)Eu;msx{0c$8%ST_G!qthMLzz%-K8rM9kS6w9ylDdLnY%SVr~y z82Ds~DC3rh95>c~-56zaj!)#caS_P!h?;Wmq)tYZaZ5xQcSlHKPIGHc&5`47+aTub zjRjeY7Ng5ub83zpcRKJkd(jhr+ z4#*1PU~~CKJ@n+b*%9O=WlcHSkiQi~#;2%+F(h&c~v$I6xc>eO4`b3WBx5Pz6j+a z&?mDg3PF;{@QIjHj{SUsi8)OuuPi3ycpmR7B68fAmInGc6>nksL>V{GZ6M|}0l4ZD zIi3gLiijLH=B|O5(?p!fXhaz|hlhcfQ;|S93y3*gV1!Z}LJIQ0QuT>4Zi&coO9UBT z6?JcV9xPX2q$3KJZ&tGcbSg!PIfnpru5120jqG%xR%Ruo%sY{yVb? z19ZwYXA~4bryhAR2Pl9}5xG7E$0wc3%;Y@=;Hso0>w%NAqHVIvyNXz^~ai7+mW@^JJlI zOL!DN@;(7{j*?|V6f&;V;%40ebSiY%)f8%FG?_ePunx5{3Imxwf#VKLDLZhmKwflH z6F6=zbgEAw<4T!rSOCzee9@-kaC8oPo={BCH#n|;H&Ou5DI&vhfKIW)h%Z2=xbH?6 z92b!rU2t4}1S50+osxIgr{K7xH&@Ye#Iy+>78A%i0(1_GSyIY5j+)GoRURY|8H0g* zQ6S2r%;ic68CM?UoKGVhoeG=HRkWOF(VQN=Ol4eg^B^*tup;n-Xv(kn%drQJPbSDu zs}?mG)uZ(fB7;nop&vx%Xg>mU>O7fC<_OTKVC{T+mD>`{nM2(ON2iVG9sxRq17oyCJ_FtC}F=CJS)!m7^vTlr;i$$}V6=2hb@CgW(fE=O>83aalKWW&(7|>R}WH zK&Ol>MhXBrMdU^Yj?4RE`UH;4<(Crybjo7oh@e&E*ETgl#ueXTh@e#zPGX42aZ3al z_ccMrO;#T{vLWLp!HS5$@y=~!0G$#{Gj{-;A~K=^&?$qf(JBC)A~Mqg(D?}>aNJxF z?iHpJi__fL6%jdZYyl%2ow8n=#R*5}Cx}AEWh^%}0dtDTtYBbH+2+l;0GQKWNEm@R zKS2~6x3?BXU{2-tb$tqs%Nt^(0GLyLmHeDJ($T4eGx_3VaP5{qlc8XZ(BsmEra55F zPY{JxQC42V0$@%N85RI@%KdCq4=|_Pf@X0txYiSeR19?D z<#RUybLw!I%abCQQIoNZPQaY1DW70HPRp7MGl4l(lOu}Lp$@Eu zZgjx}%&9<|T=h)AoFdAPYXYp`ld+5{*d0VhVaP&DA~P+(oI2PWJ^^#~dayKKEp&A1 zeN;n488=u;FGr)}Mq%huKoIG@UDF)SDtZRl5P{=*P}vZH<9g)T5P{>f(<;E63X(F+ z1m+ZxSqQ+KisdsL2j&!!s|g&Jb0Z%eFsI@R4WEEHMPwENFsFzNGkKwNHm1snj!sJi z%&F)}Bja>*&gNX8_vj&X=A5b9F(fk6LPzK9mN^}rvs>m2uAL2Zn}9j{*?5g7U``Pk zQ32*u^hD6( z?rz;g%sClzmoj{y<4*MubBbxXIwazNYG(_aMJ7ExWI`{KL%;~Q3 zO~jm@2s`KG^b;}XWWe$SKWFa@`x7yzC&KA)GB9Z(=JZ75xUC60=j6^$#GIapGM;B- z>L+4OohNfqA?EZ%)GEgGJrQ&2jGH${%sDwdPR!}9VNURK>avU3fQUJTbWES9$4xkc z6mIAO!A4UPWLyB$5J9WB)IZuz*j%m-wK0LnDEpb1b2Q3+CgyavK43xYswMOOw2i&U22A(vsfGRjuUgbJFU9a zk3N|VNT(qZ85R(8j%Jsdi8(zHId1$2N^F2WSt4@WoT6u9PKD#-%p~R<%~4xNOZ3Sg z>KT4c`Be-vi8)6T6)+QXDzL|NoLa>ak>i$#95;axN=gA17=yJA^WTREJEyV>=A$F# z98E^1nV3^42u)3taZ3a}ZW0~}s*vMm=M*vm3rtjiFcyjO?R6&RboY^!Py{uZC{E!% z5*a=bb1HZxXC^VH!e9&$Ic^Y;q6ni;3R^NYf#dG_^GwXC+=-?paNON&o{2dJ<13qq zIXw|L?k+9Q#GDFS$~#WX=>j8UM*uTT^u$cesn{md964^T^3BAY?y9je>!44@-y%~A ziHtlEbGnF*nV8ds6wbt)o(Maqc~V~1By!x?dSr6K=-dW`*xC{Uo=*!grzaxEO<;tq zZs?OEBIZ=AeBLKwP8S{{b0L~@8xZ}RgWG_JIn9gXGCY#wCcteb<}`_GRTK2MCnCp< zC2J<;G>?fF+bj5F%w046oaUWy)kKaPo0tLw*_=`KywEwAI~lT3vnI2G>E|?AfYotw z+{8l4ZjR6LXr(#JZGDjvF|slx66+F)+@=oaW(fHAjw{fTWq2 zb1<<+GyR+{;%ug$(-V>7=J2o(bDFoN`-Pa(r5#v^ITeeN&fP-H>4~7nJrOwW?jtP3 zoaPy6jSe#I!uJ;XIn6WDstFud2F;uW#GJY_ZisMJ@kEeum%v~l=2TdKX^tE>$B>1X zQyyAV6V579*oKH4Hyiar%&EQI)I^ROV7d@i$#GH&j-$Z?B486VJs*jlC2$ZT#iu!$s$Z=Z}c21>z z$VW#%r}D%YB68epiVOXmE&x}d7Z{y^-3$Gk%6n(}M2;Ioy%2M{^b8CAoSq0AcR3su zVopy4tu4@d#GIap95t3N?h#xtfSM z&2y0oS|rDfk9r~I?A#%Xn6ooca0@Y~c`H$UqKuoc*@c+XJgBIe$Z=z#Q95xj)7VuO zV$RMIQO1q4VIk&pxh@uBP9?6-6$UY9XN)fkF{gQyPszy0ad%uJ=JYj@R@J8I$)=ccxLVZJ9X@IHHD1pl8>t?6pIjwnP9+9K`TcT zGA@8^mJD~nWLu6zJHMP!5yuv2$H49DT> zT=Wn?-W*`3Z0fn0Rm;C=66#0U&+ZD| z3RkCYZWy6Mc&)tlhU0*pv(i~Qe8myP=Bc#jr7(CZp99I;x??x9|MNi$h=0dAXdcrxEallT!V&jNH#ueD_h=SwO<`f2C zr-J3p=m0wvz;Bj0V5hDpn7swCQ+dbCk^=11+cT~va9km)hMC~FdxB^K>{NVwzLjl& zoqEN_5rN~n+?e)hgR4^k>-p&P4qa1XBFZIGM~@&H%;nGxuv7O4b3Q30Th!zdm2ZHZ ziq_6GonF~PP3HK!0e0%5sNoY}r!Lgy%X|au)YV+WC%{fcYa12-c8)sf=W_>Fry_n` z$3u@Roq^#KV5h=@9Z|@*h+M~m<03L02kcaC2vZYWou42I86V|D&iMq`IWA^5-wj$`OI%O3<25Ffgb3ZI}tnsi0>^1Q}PlFw-13?gDZ)U{1Ze?rMToQIZPN z9AsSKUO5pkr_!()W&(3+R?I>G<`j`@F3#g3GP(emQ$$9q0CRqV2po4#Q^CwwoGyj_ z23M!@>^VMx<1WeM2F$5IY%_Oobt(|s5rvHF9V#Q4z?>q=4}Kdkry!?mE@WIU-8-V- zxQGm&fI0OxmXRl5P7xWc0?ertpGH)GIrUh&5jtQ_5t-2eb1E69k#V>>ZN^TeJZeHW zY;IbBIhC=~bsQX*D$bWVFsJg$8YuwgR1m$9Ctywm(VOOgId!TvQUJ_3$b^)49Ij3g z87Tnf9NY`Z8!)GyEOX64#$`D(eFDeji#0@{Rb-$sG7eX#h#a3nt0)VjX%3iE@0yvp zLwK$35f~Q0)hQy=Ctyxx9X6u_<{Vt=kPVnqR|H&5;J7n%2x-TXk~t)w7KGQz@L(1K zT%DF^2j*0A*L)%Dz?`2Tf*zMSB<=VP%&DyOt~qEGmv~?Y=G5JWe1dmiPMwTgpTP0S z$oLM-sTBXNCdjx0@bVqR1W)E9wgYo2m#k?Hm~(Ka=pC3-31@T3+<`eoWS9xeDWZG@ z%eaow863O=b1HSM>rdJ|!37Atvxr$bL>lmH8w!8y#swTrH zU`}0L$u->$SEqsoe9ioQz7rboJuHc zTjin;@H4(v>NPLyLd7#;&3#cQ7q(aT0Ut*$gu} zXqDc?-6{M8%~>LH+#RxrIrT7n-X~&C_sHQ+%;|#BcVf=YoqC8lJ9mmE=5&$RvH=C3 zJP|nF85hBh@LF9$G;>GH>4L0xgxA^>rbx``9_-qQIa{9cpZ(&8|MrJJ==#s=-w*j1F17S58-DqrhYzl^ zCkw~z7auaqUwc}V4C$|9+Qg*Fn<+M`*!635`?O#CUc~SAw5iG{dSR^~UVixz7S8sw z!be^;59C74;&FY^$_9%1XWJ-@CQ36zGdf&%_rhP&3Nb zbt<4JnB5nPD_Bi|P990idOR-wdDj5ycHX*E-nEy*`8qP;=Q-=V>`tPtunp9Q0NiD6 z_`D?7(CppibohKM68h#?Nf;&SyLb+tPh*kida`d#7UO_oB%jJ-m@#jon#yCZT=YMnf}r%YVx#-&PH z4f38Rv|3JT?aquMnde`z0!m>>E`-oL`EQO@TW4S9g^`ow#Xo2T? zA3qqmQ=J6OW*1ppPP`0SVl6T)f=zMuMz~|VRCS|X5%%7rQldz zw@UQuxD|%ixVy6-!*ufLR3vv{h_a$HACeM9U1HceHd)U=w38j#U?S@xZG9Rds@zV8mILxvvZ zfVxlH^!ZjKyd_M$YA(vZjQ_cn3`O~$jJYc$6&yWpht#Vg{j9z`w3mqGz1Fh2$>y-1 zbzjxrv_zC0n&7@JO5wGkiW_%b8)f3(yG-||DFsvZYhR8*?KbQ+y#Dht=*Jov>%g8o zDg}i{H5%%wBv$*QD=N4a&by#&ACGzpUryHC#c8}0LFxWwNVr-)P9IO6O=p{rtz~h- zssc6l!>iGR+RA1O}bsK{AqaT!lIulvJ(Ezk;`=TgXis4zHExQGGoeWR|WygB3 zQNQZwrHpsC;rC5_Y4LD|O2O63cfnfLN%a|Cjb8ihIpS*SdX{PH(;9_9YYgqfvQn!( zDn7hAai8Pz(Wwib<59t}-X#C`@u=WD_UM>-8|eIbdCx9F)>_l&0q4ar(O>iWuBi6u zIl}(b`Ry)Cv{CwW@T}e7ewH0;8OMx+W&Nq}-eW0w#Uci#)io3@OD0}y91f-=gFN=r z^2L)|LK6D*E}--03S>)Ywd$qcpA&JQuvfh@)~N=lP7U;0`e)y1M_^g$&%3N>fNCVX z*r3{y#U;HUc@|*5FO2PT$zoFRt`mufv7`Lx2Zf2J)M;5zxJrRIfAlaZ&;dg&+~&4? zakTp)zb2~2*Pde%x!#T~wLIx%iIYmYu~Wo{3V8i0b@{AeECO@>EoHZQqU`uCDbtpz zc&3}=l#i;~dwo&Y&20+G7o`^(N4;rt4!RT^^8XZD6&!ji>v4gR!VWiC6n**iIlRK@ z(yj__3VJ*X4lQ~g!X>Br;;t-56{K2J1+vcR0p^>g7>OVBcnU=+$_jXFu9Ax*?ZYDA zk3O7bDoAV8AuBJCTxL~WGjiWtOXsBGj$QZ8vR&Hs;AQuaAr+L3@U@Yss)I8%S$7+t zHg8G!v-WZORgv=>tETIWR0?h~coiH{dsASg;GCzH>hx^46zGilSV-m9b~1k~y(lYV ztFa)hwu>7rlW`R2=%Y1r(<`b=VFXF*>($xn;ApiJJ79A%z?zeZvsOYCs(Yuk7&eJW4Y#|Hdpc+jR1)JQVVC-4yCvC0r+m+s<1Gik^S}T7>TLYn%`6$JJQ{m+;Xn+cK_S zM;Y0lf{<6-x|E04Jy2uceNHny@iGnW|K{SSUGtdi3E{L)riL8Z&T>8&pruf;w8PT9!#JH=0>lcFf}IyGuanH9kKI_2$YB*1js#ccm8uSFJ*c z!A=XvwCs)+;$#P-^va~Px_uocTDq9If_+?KgY#`Ljewv+=B3 zCI@M?JIANP4E?P6+Hq6fwz7n`Dl7e@fm9vRuVuq1+`;~R6ez%V*fra9Z)M5EE;=^9S)yOdzW-pq@N%OU@o!ugeRF_+sy45Zb-o7jSB}%Cu zmE|;MwP|~cYoe@5F^Aa~!LoO4XIbeRBZ|&0oMnYM?z_UG8dn(|Ilr@gDT{=z!`4`l z`mCiA+maNN73m5N8Nl)2Y$I);B>vXS++|;XyI zHGSBr-oGB=b=RS)Avurf*F@jySY4tDy2oa=9R`_OaF?b;ziL)HI%kew%bxO3&b=5_ z{=P<13JRk##TjMCceiy=%8J&$A20Ol;rYDp7v!e54L&&5ww*Gt`w_oyWo2xxtnm9?reE9rhI8puXAyVm zpZw~FeGPLRK9r;#6V%GeY!tq$7|QTbe0+@+{3ssF*S>gCalIg3dfUFPYr%Yd7MHfu zuM5XN(ety|_WTew)A0JCr)Lso30f(Y8LSO;BtoOg!Gc55S#U_Y9udF>aTdhDX#!|v zrPFCa`7SNH23v}i2=UckD;_>{hxgS_6|2d1_a^43FzNiOYH7nbyV{ENwVq9ky|fie zGTI9aw`#9X^Rq`WnULxX!rg}Ej-w-gUn6zkHTw^rBP*V(dOOY z)A?~LA5=cwG{^D1R^-uR`Au#b{U~f3hbmu62=1fuR<6n;XXGI}$M@%Pt+cFX0-qo6 z+RgRZ&G~Tj`S5C;1JZNJkI#;9iEFg_{@$UhurzpsbU}{%sp*?%kBws^Ow8^Va ztz3N2Rk;_FTKQPMXQc0={>he~@~GtRub!*&#uLSM!FtmwKg*BK5ksBjb##7q<*z@6k3aeN`1s9lzv(|d{_n@fztVq{AIhKd@Bj0Em4C+n L|LT9A|NZ{}WG>qZ literal 0 HcmV?d00001