diff --git a/public/java/src/org/broadinstitute/sting/gatk/walkers/annotator/InbreedingCoeff.java b/public/java/src/org/broadinstitute/sting/gatk/walkers/annotator/InbreedingCoeff.java index fa3a7459b..715895526 100755 --- a/public/java/src/org/broadinstitute/sting/gatk/walkers/annotator/InbreedingCoeff.java +++ b/public/java/src/org/broadinstitute/sting/gatk/walkers/annotator/InbreedingCoeff.java @@ -68,6 +68,9 @@ public class InbreedingCoeff extends InfoFieldAnnotation implements StandardAnno if ( g.isNoCall() || !g.hasLikelihoods() ) continue; + if (g.getPloidy() != 2) // only work for diploid samples + continue; + N++; final double[] normalizedLikelihoods = MathUtils.normalizeFromLog10( g.getLikelihoods().getAsVector() ); refCount += normalizedLikelihoods[idxAA]; diff --git a/public/java/src/org/broadinstitute/sting/utils/codecs/beagle/BeagleCodec.java b/public/java/src/org/broadinstitute/sting/utils/codecs/beagle/BeagleCodec.java index 3f72359fa..656b95e7a 100755 --- a/public/java/src/org/broadinstitute/sting/utils/codecs/beagle/BeagleCodec.java +++ b/public/java/src/org/broadinstitute/sting/utils/codecs/beagle/BeagleCodec.java @@ -38,25 +38,41 @@ import java.util.*; import java.util.regex.Pattern; /** - * TODO GUILLERMO DEL ANGEL + * Codec for Beagle imputation engine * *
- * Codec Description + * Reads in tabular files with site markers and genotype posteriors, genotypes and phasing that Beagle produced *
* *
- * See also: @see VCF specification
+ * See also: @see BEAGLE home page
*
- * line 1 - * line 2 - * line 3 + * dummy header + * 20:60251 T T T T T T + * 20:60321 G G G G G G + * 20:60467 G G G G G G ** + *
+ * marker alleleA alleleB NA07056 NA07056 NA07056 + * 20:60251 T C 0.9962 0.0038 0 0.99245 0.00755 0 0.99245 0.00755 0 + * 20:60321 G T 0.98747 0.01253 0 0.99922 0.00078 0 0.99368 0.00632 0 + * 20:60467 G C 0.97475 0.02525 0 0.98718 0.01282 0 0.98718 0.01282 0 + *+ * + *
+ * 20:60251 0.747 + * 20:60321 0.763 + * 20:60467 0.524 + *+ *