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

update this walker so any variants can be passed in 

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@1426 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
ebanks 2009-08-17 16:30:39 +00:00
parent 8403618846
commit df5744bcd3
1 changed files with 3 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
java/src/org/broadinstitute/sting/gatk/walkers/indels/SNPClusterWalker.java
View File

@ -9,9 +9,7 @@ import org.broadinstitute.sting.gatk.walkers.*;
import org.broadinstitute.sting.utils.cmdLine.Argument;
@WalkerName("SNPClusters")
@By(DataSource.REFERENCE)
@Requires(DataSource.REFERENCE)
@Allows(DataSource.REFERENCE)
@Requires(value={DataSource.REFERENCE},referenceMetaData={@RMD(name="snps",type=AllelicVariant.class)})
public class SNPClusterWalker extends RefWalker<GenomeLoc, GenomeLoc> {
@Argument(fullName="windowSize", shortName="window", doc="window size for calculating clusters", required=false)
int windowSize = 10;
@ -22,10 +20,8 @@ public class SNPClusterWalker extends RefWalker<GenomeLoc, GenomeLoc> {
}
public GenomeLoc map(RefMetaDataTracker tracker, ReferenceContext ref, AlignmentContext context) {
AllelicVariant eval = (AllelicVariant)tracker.lookup("eval", null);
if ( eval instanceof SNPCallFromGenotypes )
return context.getLocation();
return null;
AllelicVariant snp = (AllelicVariant)tracker.lookup("snps", null);
return (snp != null && snp.isSNP()) ? context.getLocation() : null;
}
public void onTraversalDone(GenomeLoc sum) {