diff --git a/protected/java/src/org/broadinstitute/sting/gatk/walkers/diagnostics/missing/QualifyMissingIntervals.java b/protected/java/src/org/broadinstitute/sting/gatk/walkers/diagnostics/missing/QualifyMissingIntervals.java index 62716d6d2..d0db3ef98 100644 --- a/protected/java/src/org/broadinstitute/sting/gatk/walkers/diagnostics/missing/QualifyMissingIntervals.java +++ b/protected/java/src/org/broadinstitute/sting/gatk/walkers/diagnostics/missing/QualifyMissingIntervals.java @@ -87,7 +87,7 @@ import java.util.List; * *

Input

*

- * A reference file + * A reference file (for GC content), the input bam file (for base and mapping quality calculation), the missing intervals (in the -L), the baits/targets used to sequence (in the -targets) and a bed file with the coding sequence intervals of the genome (in the -cds) *

* *

Output

@@ -100,6 +100,7 @@ import java.util.List; * java -Xmx2g -jar GenomeAnalysisTK.jar \ * -T QualifyMissingIntervals \ * -R ref.fasta \ + * -I input.bam \ * -o output.grp \ * -L input.intervals \ * -cds cds.intervals \ @@ -117,7 +118,7 @@ public final class QualifyMissingIntervals extends LocusWalker public File targetsFile; @Argument(shortName = "cds", required = false) - public File cdsFile; + public File cdsFile = null; GATKReport simpleReport; GenomeLocSortedSet target; @@ -133,7 +134,8 @@ public final class QualifyMissingIntervals extends LocusWalker target = new GenomeLocSortedSet(parser); cds = new GenomeLocSortedSet(parser); parseFile(targetsFile, target, parser); - parseFile(cdsFile, cds, parser); + if (cdsFile != null) + parseFile(cdsFile, cds, parser); } public Metrics reduceInit() {