From 9de06560df2fa46c5aba1642688b598568c4fc44 Mon Sep 17 00:00:00 2001 From: Menachem Fromer Date: Thu, 11 Aug 2011 17:54:16 -0400 Subject: [PATCH] Update to new RodBinding system --- .../gatk/walkers/phasing/MergeMNPsWalker.java | 12 +++++++----- .../MergeSegregatingAlternateAllelesWalker.java | 14 ++++++-------- .../walkers/phasing/MergeMNPsIntegrationTest.java | 2 +- ...SegregatingAlternateAllelesIntegrationTest.java | 2 +- 4 files changed, 15 insertions(+), 15 deletions(-) diff --git a/public/java/src/org/broadinstitute/sting/gatk/walkers/phasing/MergeMNPsWalker.java b/public/java/src/org/broadinstitute/sting/gatk/walkers/phasing/MergeMNPsWalker.java index 14f54ec7d..809772c05 100644 --- a/public/java/src/org/broadinstitute/sting/gatk/walkers/phasing/MergeMNPsWalker.java +++ b/public/java/src/org/broadinstitute/sting/gatk/walkers/phasing/MergeMNPsWalker.java @@ -25,11 +25,12 @@ package org.broadinstitute.sting.gatk.walkers.phasing; import org.broadinstitute.sting.commandline.Argument; +import org.broadinstitute.sting.commandline.Input; import org.broadinstitute.sting.commandline.Output; +import org.broadinstitute.sting.commandline.RodBinding; import org.broadinstitute.sting.gatk.contexts.AlignmentContext; import org.broadinstitute.sting.gatk.contexts.ReferenceContext; import org.broadinstitute.sting.gatk.refdata.RefMetaDataTracker; -import org.broadinstitute.sting.gatk.refdata.ReferenceOrderedDatum; import org.broadinstitute.sting.gatk.walkers.*; import org.broadinstitute.sting.utils.codecs.vcf.VCFHeader; import org.broadinstitute.sting.utils.codecs.vcf.VCFHeaderLine; @@ -58,7 +59,8 @@ public class MergeMNPsWalker extends RodWalker { @Argument(fullName = "maxGenomicDistanceForMNP", shortName = "maxDistMNP", doc = "The maximum reference-genome distance between consecutive heterozygous sites to permit merging phased VCF records into a MNP record; [default:1]", required = false) protected int maxGenomicDistanceForMNP = 1; - private String rodName = "variant"; + @Input(fullName="variant", shortName = "V", doc="Select variants from this VCF file", required=true) + public RodBinding variants; public void initialize() { initializeVcfWriter(); @@ -74,8 +76,8 @@ public class MergeMNPsWalker extends RodWalker { hInfo.addAll(VCFUtils.getHeaderFields(getToolkit())); hInfo.add(new VCFHeaderLine("reference", getToolkit().getArguments().referenceFile.getName())); - Map rodNameToHeader = getVCFHeadersFromRods(getToolkit(), Arrays.asList(rodName)); - vcMergerWriter.writeHeader(new VCFHeader(hInfo, new TreeSet(rodNameToHeader.get(rodName).getGenotypeSamples()))); + Map rodNameToHeader = getVCFHeadersFromRods(getToolkit(), Arrays.asList(variants.getName())); + vcMergerWriter.writeHeader(new VCFHeader(hInfo, new TreeSet(rodNameToHeader.get(variants.getName()).getGenotypeSamples()))); } public boolean generateExtendedEvents() { @@ -98,7 +100,7 @@ public class MergeMNPsWalker extends RodWalker { if (tracker == null) return null; - for (VariantContext vc : tracker.getValues(VariantContext.class, rodName, context.getLocation())) + for (VariantContext vc : tracker.getValues(variants, context.getLocation())) writeVCF(vc); return 0; diff --git a/public/java/src/org/broadinstitute/sting/gatk/walkers/phasing/MergeSegregatingAlternateAllelesWalker.java b/public/java/src/org/broadinstitute/sting/gatk/walkers/phasing/MergeSegregatingAlternateAllelesWalker.java index e1be5e5c5..96d5c471f 100644 --- a/public/java/src/org/broadinstitute/sting/gatk/walkers/phasing/MergeSegregatingAlternateAllelesWalker.java +++ b/public/java/src/org/broadinstitute/sting/gatk/walkers/phasing/MergeSegregatingAlternateAllelesWalker.java @@ -24,13 +24,10 @@ package org.broadinstitute.sting.gatk.walkers.phasing; -import org.broadinstitute.sting.commandline.Argument; -import org.broadinstitute.sting.commandline.Hidden; -import org.broadinstitute.sting.commandline.Output; +import org.broadinstitute.sting.commandline.*; import org.broadinstitute.sting.gatk.contexts.AlignmentContext; import org.broadinstitute.sting.gatk.contexts.ReferenceContext; import org.broadinstitute.sting.gatk.refdata.RefMetaDataTracker; -import org.broadinstitute.sting.gatk.refdata.ReferenceOrderedDatum; import org.broadinstitute.sting.gatk.walkers.*; import org.broadinstitute.sting.utils.GenomeLocParser; import org.broadinstitute.sting.utils.codecs.vcf.VCFHeader; @@ -81,7 +78,8 @@ public class MergeSegregatingAlternateAllelesWalker extends RodWalker variants; public void initialize() { initializeVcfWriter(); @@ -111,8 +109,8 @@ public class MergeSegregatingAlternateAllelesWalker extends RodWalker rodNameToHeader = getVCFHeadersFromRods(getToolkit(), Arrays.asList(rodName)); - vcMergerWriter.writeHeader(new VCFHeader(hInfo, new TreeSet(rodNameToHeader.get(rodName).getGenotypeSamples()))); + Map rodNameToHeader = getVCFHeadersFromRods(getToolkit(), Arrays.asList(variants.getName())); + vcMergerWriter.writeHeader(new VCFHeader(hInfo, new TreeSet(rodNameToHeader.get(variants.getName()).getGenotypeSamples()))); } public boolean generateExtendedEvents() { @@ -135,7 +133,7 @@ public class MergeSegregatingAlternateAllelesWalker extends RodWalker