ReassignMappingQualityFilter

Set the mapping quality of all reads to a given value

Category Read Filters


Overview

If a BAM file contains erroneous or missing mapping qualities (MAPQ), this read transformer will set all your mapping qualities to a given value (see arguments list for default value).

See also

ReassignOneMappingQualityFilter: reassigns a single MAPQ value, as opposed to all those found in the BAM file.

Caveats

Note that due to the order of operations involved in applying filters, it is possible that other read filters (determined either at command-line or internally by the tool you are using) will be applied to your data before this read transformation can be applied. If one of those other filters acts on the read mapping quality (MAPQ), then you may not obtain the expected results. Unfortunately it is currently not possible to change the order of operations from command line. To avoid the problem, we recommend applying this filter separately from any other analysis, using PrintReads.

Input

BAM file(s)

Output

BAM file(s) with the mapping qualities of all reads reassigned to the specified value

Usage example

  java -jar GenomeAnalysisTK.jar \
      -T PrintReads \
      -R reference.fasta \
      -I input.bam \
      -o output.file \
      -rf ReassignMappingQuality \
      -DMQ 35
  

Command-line Arguments

ReassignMappingQualityFilter specific arguments

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
--default_mapping_quality
 -DMQ
60 Default read mapping quality to assign to all reads

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.


--default_mapping_quality / -DMQ

Default read mapping quality to assign to all reads

int  60  [ [ -∞  ∞ ] ]


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.