zzh
/
gatk-3.8
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Total ploidy does not have to equal to the number of priors

This commit is contained in:
Ron Levine 2017-02-20 20:28:41 -05:00
parent c55c77ab27
commit 3c347c1385
2 changed files with 10 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
protected/gatk-tools-protected/src/main/java/org/broadinstitute/gatk/tools/walkers/genotyper/CustomAFPriorProvider.java
View File

@ -86,8 +86,6 @@ public class CustomAFPriorProvider extends AFPriorProvider {
@Override
protected double[] buildPriors(final int totalPloidy) {
if (totalPloidy != priors.length - 1)
throw new IllegalStateException("requesting an invalid prior total ploidy " + totalPloidy + " != " + (priors.length - 1));
return priors;
}
}

10
protected/gatk-tools-protected/src/test/java/org/broadinstitute/gatk/tools/walkers/variantutils/GenotypeGVCFsIntegrationTest.java
View File

@ -700,4 +700,14 @@ public class GenotypeGVCFsIntegrationTest extends WalkerTest {
spec.disableShadowBCF();
executeTest("testHomRefHighMQ", spec);
}
@Test
public void testInputPrior() {
final WalkerTestSpec spec = new WalkerTestSpec(
baseTestString(" -V " + privateTestDir + "gvcfExample1.vcf -V " + privateTestDir + "gvcfExample2.vcf " +
"-L 1:69485-69791 -inputPrior 0.2 -inputPrior 0.2 -inputPrior 0.2 -inputPrior 0.2", b37KGReference),
Collections.singletonList("860d133262160bbc75ce5849e5fa491f"));
spec.disableShadowBCF();
executeTest("testInputPrior", spec);
}
}