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

To address weird case with all hom-refs, but alt allele is present, skip AS_QD for all-ref sites and remove raw annotations whether or not they can be successfully finalized

This commit is contained in:
Laura Gauthier 2016-03-14 13:17:29 -04:00
parent 4efde50225
commit 31fc64f82c
2 changed files with 4 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/annotator/AS_QualByDepth.java
View File

@ -172,6 +172,8 @@ public class AS_QualByDepth extends InfoFieldAnnotation implements ReducibleAnno
return null;
final List<Integer> standardDepth = getAlleleDepths(genotypes);
if (standardDepth == null) //all no-calls and homRefs
return null;
//Parse the VC's allele-specific qual values
List<Object> alleleQualObjList = vc.getAttributeAsList(GATKVCFConstants.AS_QUAL_KEY);

4
public/gatk-tools-public/src/main/java/org/broadinstitute/gatk/tools/walkers/annotator/VariantAnnotatorEngine.java
View File

@ -357,9 +357,9 @@ public class VariantAnnotatorEngine {
final Map<String, Object> annotationsFromCurrentType = currentASannotation.finalizeRawData(vc, originalVC);
if ( annotationsFromCurrentType != null ) {
infoAnnotations.putAll(annotationsFromCurrentType);
//clean up raw annotation data after annotations are finalized
infoAnnotations.remove(currentASannotation.getRawKeyName());
}
//clean up raw annotation data after annotations are finalized
infoAnnotations.remove(currentASannotation.getRawKeyName());
}
// generate a new annotated VC