Set the mapping quality of reads with a given value to another given 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)
BAM file(s) with one read mapping quality selectively reassigned as desired
java -jar GenomeAnalysisTK.jar \
-T PrintReads \
-R reference.fasta \
-I input.bam \
-o output.file \
-rf ReassignOneMappingQuality \
-RMQF 255 \
-RMQT 60
This table summarizes the command-line arguments that are specific to this tool. For more details on each argument, see the list further down below the table or click on an argument name to jump directly to that entry in the list.
| Argument name(s) | Default value | Summary | |
|---|---|---|---|
| Optional Parameters | |||
| --reassign_mapping_quality_from -RMQF |
255 | Original mapping quality | |
| --reassign_mapping_quality_to -RMQT |
60 | Desired 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.
Original mapping quality
int 255 [ [ -∞ ∞ ] ]
Desired mapping quality
int 60 [ [ -∞ ∞ ] ]