ReassignOriginalMQAfterIndelRealignmentFilter

Revert the MQ of reads that were modified by IndelRealigner

Category Read Filters


Overview

IndelRealigner systematically adds +10 to the MQ of the reads it realigns. In some cases, that brings the resulting MQ to a value higher than MQ 60, which is the normal cap for MQ values. Since many downstream tools assume that MQ is <= 60, this may potentially cause problems.

This read filter makes it possible to revert the MQ values of all the reads touched by IndelRealigner. It works by subtracting 10 from the MQ of all reads that have an "OC" tag, which stands for "original CIGAR" and is added by IndelRealigner to any read that it realigns.

Usage example

Enable the filter

     java -jar GenomeAnalysisTk.jar \
         -T ToolName \
         -R reference.fasta \
         -I input.bam \
         -o output.file \
         -rf ReassignOriginalMQAfterIndelRealignmentFilter
 

Caveat

There is currently no way to tell programmatically that a file has already been processed with this filter, so you should check the header manually before running this tool. Running it multiple times on the same BAM file would levy an unjustified penalty on realigned reads.


Return to top


See also GATK Documentation Index | Tool Docs Index | Support Forum

GATK version 3.8-0-ge9d806836 built at 2017/07/29 01:40:22.