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

Cleanup for new SampleUtils function 

-- getVCFHeadersFromRods(rods) is now available so that you don't have getVCFHeadersFromRods(rods, null) throughout the codebase
This commit is contained in:
Mark DePristo 2011-11-10 13:27:13 -05:00
parent 35fe9c8a06
commit 67b022c34b
5 changed files with 11 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
public/java/src/org/broadinstitute/sting/gatk/walkers/varianteval/evaluators/CountVariants.java
View File

@ -41,6 +41,9 @@ public class CountVariants extends VariantEvaluator implements StandardEval {
public long nDeletions = 0;
@DataPoint(description = "Number of complex indels")
public long nComplex = 0;
@DataPoint(description = "Number of symbolic events")
public long nSymbolic = 0;
@DataPoint(description = "Number of mixed loci (loci that can't be classified as a SNP, Indel or MNP)")
public long nMixed = 0;
@ -131,8 +134,7 @@ public class CountVariants extends VariantEvaluator implements StandardEval {
nMixed++;
break;
case SYMBOLIC:
// ignore symbolic alleles, but don't fail
// todo - consistent way of treating symbolic alleles thgoughout codebase?
nSymbolic++;
break;
default:
throw new ReviewedStingException("Unexpected VariantContext type " + vc1.getType());

2
public/java/src/org/broadinstitute/sting/gatk/walkers/varianteval/evaluators/G1KPhaseITable.java
View File

@ -103,7 +103,7 @@ public class G1KPhaseITable extends VariantEvaluator {
}
public String update2(VariantContext eval, VariantContext comp, RefMetaDataTracker tracker, ReferenceContext ref, AlignmentContext context) {
if ( eval == null ) return null;
if ( eval == null || eval.isMonomorphic() ) return null;
switch (eval.getType()) {
// case NO_VARIATION:

2
public/java/src/org/broadinstitute/sting/gatk/walkers/variantutils/CombineVariants.java
View File

@ -162,7 +162,7 @@ public class CombineVariants extends RodWalker<Integer, Integer> {
private boolean sitesOnlyVCF = false;
public void initialize() {
Map<String, VCFHeader> vcfRods = VCFUtils.getVCFHeadersFromRods(getToolkit(), null);
Map<String, VCFHeader> vcfRods = VCFUtils.getVCFHeadersFromRods(getToolkit());
if ( PRIORITY_STRING == null ) {
PRIORITY_STRING = Utils.join(",", vcfRods.keySet());

2
public/java/src/org/broadinstitute/sting/utils/SampleUtils.java
View File

@ -150,7 +150,7 @@ public class SampleUtils {
// iterate to get all of the sample names
for ( Map.Entry<String, VCFHeader> pair : VCFUtils.getVCFHeadersFromRods(toolkit, null).entrySet() ) {
for ( Map.Entry<String, VCFHeader> pair : VCFUtils.getVCFHeadersFromRods(toolkit).entrySet() ) {
Set<String> vcfSamples = pair.getValue().getGenotypeSamples();
for ( String sample : vcfSamples )
addUniqueSample(samples, sampleOverlapMap, rodNamesToSampleNames, sample, pair.getKey());

4
public/java/src/org/broadinstitute/sting/utils/codecs/vcf/VCFUtils.java
View File

@ -51,6 +51,10 @@ public class VCFUtils {
return getVCFHeadersFromRods(toolkit, names);
}
public static Map<String, VCFHeader> getVCFHeadersFromRods(GenomeAnalysisEngine toolkit) {
return getVCFHeadersFromRods(toolkit, (Collection<String>)null);
}
public static Map<String, VCFHeader> getVCFHeadersFromRods(GenomeAnalysisEngine toolkit, Collection<String> rodNames) {
Map<String, VCFHeader> data = new HashMap<String, VCFHeader>();