From 4bec93e3e45f628861b1c38b9b4adf929acc8587 Mon Sep 17 00:00:00 2001 From: fromer Date: Tue, 18 Jan 2011 16:09:34 +0000 Subject: [PATCH] Permit retrieval of read names for debugging purposes git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@5011 348d0f76-0448-11de-a6fe-93d51630548a --- .../phasing/ReadBasedPhasingValidationWalker.java | 13 ++++++++----- 1 file changed, 8 insertions(+), 5 deletions(-) diff --git a/java/src/org/broadinstitute/sting/playground/gatk/walkers/phasing/ReadBasedPhasingValidationWalker.java b/java/src/org/broadinstitute/sting/playground/gatk/walkers/phasing/ReadBasedPhasingValidationWalker.java index 6227a1f8a..e263a7217 100755 --- a/java/src/org/broadinstitute/sting/playground/gatk/walkers/phasing/ReadBasedPhasingValidationWalker.java +++ b/java/src/org/broadinstitute/sting/playground/gatk/walkers/phasing/ReadBasedPhasingValidationWalker.java @@ -175,9 +175,12 @@ public class ReadBasedPhasingValidationWalker extends RodWalker reads = new LinkedList(); - for (PhasingRead rd : buildReads) { - if (rd.getNonNullIndices().length == NUM_IN_PAIR) // only want reads with BOTH bases called [possibly as deleted ("D")] + for (Map.Entry readEntry : buildReads.entrySet()) { + PhasingRead rd = readEntry.getValue(); + if (rd.getNonNullIndices().length == NUM_IN_PAIR) { // only want reads with BOTH bases called [possibly as deleted ("D")] reads.add(rd); + logger.debug("Read name: " + readEntry.getKey() + "\trd: " + rd); + } } // Count the occurence of each "haplotype": @@ -339,7 +342,7 @@ class SiteGenotypeAndReads { } } -class PhasingReadList implements Iterable { +class PhasingReadList { private Map readsAtSites = null; private int numSites; @@ -366,8 +369,8 @@ class PhasingReadList implements Iterable { } } - public Iterator iterator() { - return readsAtSites.values().iterator(); + public Set> entrySet() { + return readsAtSites.entrySet(); } public int size() {