ReassignOneMappingQualityFilter

A read filter (transformer) that changes a given read mapping quality to a different value.

Category Read Filters

Overview

This 'filter' 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.

Input

BAM file(s)

Output

BAM file(s) with one read mapping quality selectively reassigned as desired

Examples

    java
      -jar GenomeAnalysisTK.jar
      -rf ReassignOneMappingQuality
      -RMQF 255
      -RMQT 60

Command-line Arguments

ReassignOneMappingQualityFilter specific arguments

This table summarizes the command-line arguments that are specific to this tool. For details, see the list further down below the table.

Name	Type	Default value	Summary
Optional
--reassign_mapping_quality_to	int	60	Desired mapping quality
-RMQF	int	255	Original mapping quality

Argument details

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.

--reassign_mapping_quality_to / -RMQT ( int with default value 60 )

Desired mapping quality.

-RMQF / --reassign_mapping_quality_from ( int with default value 255 )

Original mapping quality.

GATK version 2.5-2-gdb4546e built at 2013/05/01 09:32:36.