Make chimeric read filter optional

getzlab · Sep 13, 2021 · f22f927 · f22f927
1 parent df19552
commit f22f927
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 1 deletion.
diff --git a/src/org/broadinstitute/cga/tools/gatk/walkers/cancer/mutect/MuTect.java b/src/org/broadinstitute/cga/tools/gatk/walkers/cancer/mutect/MuTect.java
@@ -791,7 +791,7 @@ private LocusReadPile filterReads(final ReferenceContext ref, final ReadBackedPi
 
             // is this a discordant read pair?
             // TODO: add option to disable this filter
-            if(read.getMateReferenceIndex() != read.getReferenceIndex()) {
+            if(MTAC.EXCLUDE_CHIMERIC_READS && read.getMateReferenceIndex() != read.getReferenceIndex()) {
                 continue;
             }
 

diff --git a/src/org/broadinstitute/cga/tools/gatk/walkers/cancer/mutect/MuTectArgumentCollection.java b/src/org/broadinstitute/cga/tools/gatk/walkers/cancer/mutect/MuTectArgumentCollection.java
@@ -130,6 +130,9 @@ public class MuTectArgumentCollection {
     @Argument(fullName = "heavily_clipped_read_fraction", required = false, doc = "if this fraction or more of the bases in a read are soft/hard clipped, do not use this read for mutation calling")
     public float HEAVILY_CLIPPED_READ_FRACTION = 0.30f;
 
+    @Argument(fullName = "exclude_chimeric_reads", required = false, doc = "exclude chimeric reads from processing")
+    public boolean EXCLUDE_CHIMERIC_READS = false;
+
     @Argument(fullName = "fraction_mapq0_threshold", required = false, doc = "threshold for determining if there is relatedness between the alt and ref allele read piles")
     public float FRACTION_MAPQ0_THRESHOLD = 0.5f;