diff --git a/java/src/org/broadinstitute/sting/playground/gatk/walkers/ReadBackedPhasingWalker.java b/java/src/org/broadinstitute/sting/playground/gatk/walkers/ReadBackedPhasingWalker.java index 132c762c0..029cd7c41 100755 --- a/java/src/org/broadinstitute/sting/playground/gatk/walkers/ReadBackedPhasingWalker.java +++ b/java/src/org/broadinstitute/sting/playground/gatk/walkers/ReadBackedPhasingWalker.java @@ -55,11 +55,12 @@ import static org.broadinstitute.sting.utils.vcf.VCFUtils.getVCFHeadersFromRods; @Requires(value = {DataSource.READS, DataSource.REFERENCE}, referenceMetaData = @RMD(name = "variant", type = ReferenceOrderedDatum.class)) @By(DataSource.READS) -@ReadFilters( {ZeroMappingQualityReadFilter.class} ) // Filter out all reads with zero mapping quality +@ReadFilters({ZeroMappingQualityReadFilter.class}) +// Filter out all reads with zero mapping quality public class ReadBackedPhasingWalker extends RodWalker { - @Output(doc="File to which variants should be written",required=true) + @Output(doc = "File to which variants should be written", required = true) protected VCFWriter writer = null; @Argument(fullName = "cacheWindowSize", shortName = "cacheWindow", doc = "The window size (in bases) to cache variant sites and their reads; [default:20000]", required = false) @@ -68,8 +69,8 @@ public class ReadBackedPhasingWalker extends RodWalker P(error) = 10^(-4.77/10) = 0.33, P(correct) = 0.66, so that we have odds ratio of >= 2 + @Argument(fullName = "phaseQualityThresh", shortName = "phaseThresh", doc = "The minimum phasing quality score required to output phasing; [default:10.0]", required = false) + protected Double phaseQualityThresh = 10.0; // PQ = 10.0 <=> P(error) = 10^(-10/10) = 0.1, P(correct) = 0.9 @Argument(fullName = "variantStatsFilePrefix", shortName = "variantStats", doc = "The prefix of the VCF/phasing statistics files", required = false) protected String variantStatsFilePrefix = null; @@ -78,7 +79,6 @@ public class ReadBackedPhasingWalker extends RodWalker phasedSites = null; // the phased VCs to be emitted, and the alignment bases at these positions private static PreciseNonNegativeDouble ZERO = new PreciseNonNegativeDouble(0.0); - private static boolean DEBUG_DETAILED = true; private LinkedList rodNames = null; private PhasingQualityStatsWriter statsWriter = null; @@ -125,7 +125,7 @@ public class ReadBackedPhasingWalker extends RodWalker discardIrrelevantPhasedSites() { List vcList = new LinkedList(); - + VariantContext nextToPhaseVc = null; if (!unphasedSiteQueue.isEmpty()) nextToPhaseVc = unphasedSiteQueue.peek().variant; @@ -215,6 +215,7 @@ public class ReadBackedPhasingWalker extends RodWalker