Skip to content

Commit

Permalink
compatible with GATK 2.7.1 (includes bugfix for MAPQ > 127 causing ov…
Browse files Browse the repository at this point in the history
…erflows in MuTect)
  • Loading branch information
kcibul committed Aug 22, 2013
1 parent f07a1d1 commit 8f5ae96
Showing 1 changed file with 2 additions and 2 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -575,10 +575,10 @@ public Integer map(final RefMetaDataTracker tracker, final ReferenceContext ref,
final LocusReadPile refPile = new LocusReadPile(referencePileup, altAllele, 0, 0, MTAC.ENABLE_QSCORE_OUTPUT);

// Set the maximum observed mapping quality score for the reference and alternate alleles
byte[] rmq = referencePileup.getMappingQuals();
int[] rmq = referencePileup.getMappingQuals();
candidate.setTumorRefMaxMapQ((rmq.length==0)?0:NumberUtils.max(rmq));

byte[] amq = mutantPileup.getMappingQuals();
int[] amq = mutantPileup.getMappingQuals();
candidate.setTumorAltMaxMapQ((amq.length==0)?0:NumberUtils.max(amq));

candidate.setStrandContingencyTable(SequenceUtils.getStrandContingencyTable(refPile, mutantPile));
Expand Down

0 comments on commit 8f5ae96

Please sign in to comment.