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

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
This commit is contained in:
fromer 2011-01-18 16:09:34 +00:00
parent f8ba76d87c
commit 4bec93e3e4
1 changed files with 8 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
java/src/org/broadinstitute/sting/playground/gatk/walkers/phasing/ReadBasedPhasingValidationWalker.java
View File

@ -175,9 +175,12 @@ public class ReadBasedPhasingValidationWalker extends RodWalker<Integer, Integer
buildReads.updateBases(1, readBases);
List<PhasingRead> reads = new LinkedList<PhasingRead>();
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<String, PhasingRead> 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<PhasingRead> {
class PhasingReadList {
private Map<String, PhasingRead> readsAtSites = null;
private int numSites;
@ -366,8 +369,8 @@ class PhasingReadList implements Iterable<PhasingRead> {
}
}
public Iterator<PhasingRead> iterator() {
return readsAtSites.values().iterator();
public Set<Map.Entry<String, PhasingRead>> entrySet() {
return readsAtSites.entrySet();
}
public int size() {