Filter out reads that are over-soft-clipped
Category Read Filters
Overview
Filter out reads where the number of bases without soft-clips (M, I, X, and = CIGAR operators) is lower than a threshold.This filter is intended to filter out reads that are potentially from foreign organisms. From experience with sequencing of human DNA we have found cases of contamination by bacterial organisms; the symptoms of such contamination are a class of reads with only a small number of aligned bases and additionally many soft-clipped bases. This filter is intended to remove such reads.
Note: Consecutive soft-clipped blocks are treated as a single block. For example, 1S2S10M1S2S is treated as 3S10M3S
OverclippedReadFilter 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 Tool Arguments | |||
--dont-require-soft-clips-both-ends |
false | Allow a read to be filtered out based on having only 1 soft-clipped block. By default, both ends must have a soft-clipped block, setting this flag requires only 1 soft-clipped block | |
--filter-too-short |
30 | Minimum number of aligned bases |
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.
--dont-require-soft-clips-both-ends / NA
Allow a read to be filtered out based on having only 1 soft-clipped block. By default, both ends must have a soft-clipped block, setting this flag requires only 1 soft-clipped block
boolean false
--filter-too-short / NA
Minimum number of aligned bases
int 30 [ [ -∞ ∞ ] ]
GATK version 4.1.6.0-SNAPSHOT built at Thu, 2 Apr 2020 14:54:17 -0400.
0 comments
Please sign in to leave a comment.