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

Support for multiple records at the same site. The first record gets chr:start, and subsequent records get chr:start_2, chr:start_3, etc.

This commit is contained in:
Mark DePristo 2011-07-18 15:43:52 -04:00
parent 782a05e9b5
commit c05451047c
1 changed files with 7 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
public/java/src/org/broadinstitute/sting/gatk/walkers/diffengine/VCFDiffableReader.java
View File

@ -72,13 +72,19 @@ public class VCFDiffableReader implements DiffableReader {
} }
String line = lineReader.readLine(); String line = lineReader.readLine();
int count = 0; int count = 0, nRecordsAtPos = 1;
String prevName = "";
while ( line != null ) { while ( line != null ) {
if ( count++ > maxElementsToRead && maxElementsToRead != -1) if ( count++ > maxElementsToRead && maxElementsToRead != -1)
break; break;
VariantContext vc = (VariantContext)vcfCodec.decode(line); VariantContext vc = (VariantContext)vcfCodec.decode(line);
String name = vc.getChr() + ":" + vc.getStart(); String name = vc.getChr() + ":" + vc.getStart();
if ( name.equals(prevName) ) {
name += "_" + ++nRecordsAtPos;
} else {
prevName = name;
}
DiffNode vcRoot = DiffNode.empty(name, root); DiffNode vcRoot = DiffNode.empty(name, root);
// add fields // add fields