A read filter (transformer) that changes a given read mapping quality to a different value.
This read transformer will change a certain read mapping quality to a different value without affecting reads that have other mapping qualities. This is intended primarily for users of RNA-Seq data handling programs such as TopHat, which use MAPQ = 255 to designate uniquely aligned reads. According to convention, 255 normally designates "unknown" quality, and most GATK tools automatically ignore such reads. By reassigning a different mapping quality to those specific reads, users of TopHat and other tools can circumvent this problem without affecting the rest of their dataset.
This differs from the ReassignMappingQuality filter by its selectivity -- only one mapping quality is targeted. ReassignMappingQuality will change ALL mapping qualities to a single one, and is typically used for datasets that have no assigned mapping qualities.
BAM file(s) with one read mapping quality selectively reassigned as desired
java -jar GenomeAnalysisTK.jar -T PrintReads -rf ReassignOneMappingQuality -RMQF 255 -RMQT 60
This table summarizes the command-line arguments that are specific to this tool. For details, see the list further down below the table.
|--reassign_mapping_quality_to||int||60||Desired mapping quality|
|-RMQF||int||255||Original mapping quality|
Arguments in this list are specific to this tool. Keep in mind that other arguments are available that are shared with other tools (e.g. command-line GATK arguments); see Inherited arguments above.
Desired mapping quality.
Original mapping quality.
GATK version 2.8-1-g2a26ec9 built at 2013/12/06 16:54:02.