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

Added an argument ('extended' or 'ext') that outputs the four-base probs in a long format. 

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@485 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
kiran 2009-04-21 22:27:26 +00:00
parent dac76f041b
commit 5abfc7d079
1 changed files with 8 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

24
java/src/org/broadinstitute/sting/gatk/walkers/PileupWalker.java
View File

@ -24,6 +24,9 @@ import java.util.ArrayList;
public class PileupWalker extends LocusWalker<Integer, Integer> {
@Argument(fullName="verbose",required=false,defaultValue="false")
public boolean VERBOSE;
@Argument(fullName="extended",shortName="ext",required=false,defaultValue="false")
public boolean EXTENDED;
public boolean FLAG_UNCOVERED_BASES = true; // todo: how do I make this a command line argument?
@ -47,22 +50,6 @@ public class PileupWalker extends LocusWalker<Integer, Integer> {
if ( VERBOSE ) {
extras += " BQ=" + pileup.getQualsAsInts();
extras += " MQ=" + pileup.getMappingQualsAsInts();
String sqbases = pileup.getSecondaryBasePileup();
String sqquals = pileup.getSecondaryQualPileup();
if (sqbases != null && sqquals != null) {
assert(sqbases.length() == sqquals.length());
extras += " SQ=";
for (int i = 0; i < sqbases.length(); i++) {
extras += sqbases.charAt(i);
extras += sqquals.charAt(i);
if (i < sqbases.length() - 1) {
extras += ',';
}
}
}
}
String rodString = "";
@ -85,6 +72,11 @@ public class PileupWalker extends LocusWalker<Integer, Integer> {
out.printf("%s%s %s%n", pileup.getPileupString(), extras, rodString);
//}
if ( EXTENDED ) {
String probDists = pileup.getProbDistPileup();
System.out.println(probDists);
}
return 1;
}