-6,--illumina1.3+
Assume the quality is in the Illumina 1.3+ encoding.
-A,--count-orphans
Do not skip anomalous read pairs in variant calling. Anomalous read pairs are those marked in
the FLAG field as paired in sequencing but without the properly-paired flag set.
-b,--bam-listFILE
List of input BAM files, one file per line [null]
-B,--no-BAQ
Disable base alignment quality (BAQ) computation. See BAQ below.
-C,--adjust-MQINT
Coefficient for downgrading mapping quality for reads containing excessive mismatches.
Mismatches are counted as a proportion of the number of aligned bases ("M", "X" or "=" CIGAR
operations), along with their quality, to derive an upper-bound of the mapping quality.
Original mapping qualities lower than this are left intact, while higher ones are capped at the
new adjusted score.
The exact formula is complex and likely tuned to specific instruments and specific alignment
tools, so this option is disabled by default (indicated as having a zero value). Variables in
the formulae and their meaning are defined below.
VariableMeaning/formula
───────────────────────────────────────────────────────────
M The number of matching CIGAR bases (operation
"M", "X" or "=").
X The number of substitutions with quality >= 13.
SubQ The summed quality of substitution bases
included in X, capped at a maximum of quality
33 per mismatching base.
ClipQ The summed quality of soft-clipped or hard-
clipped bases. This has no minimum or maximum
quality threshold per base. For hard-clipped
bases the per-base quality is taken as 13.
T SubQ - 10 * log10(M^X / X!) + ClipQ/5
Cap MAX(0, INT * sqrt((INT - T) / INT))
Some notes on the impact of this.
○ As the number of mismatches increases, the mapping quality cap reduces, eventually resulting
in discarded alignments.
○ High quality mismatches reduces the cap faster than low quality mismatches.
○ The starting INT value also acts as a hard cap on mapping quality, even when zero mismatches
are observed.
○ Indels have no impact on the mapping quality.
The intent of this option is to work around aligners that compute a mapping quality using a
local alignment without having any regard to the degree of clipping required or consideration
of potential contamination or large scale insertions with respect to the reference. A record
may align uniquely and have no close second match, but having a high number of mismatches may
still imply that the reference is not the correct site.
However we do not recommend use of this parameter unless you fully understand the impact of it
and have determined that it is appropriate for your sequencing technology.
-d,--max-depthINT
At a position, read maximally INT reads per input file. Setting this limit reduces the amount
of memory and time needed to process regions with very high coverage. Passing zero for this
option sets it to the highest possible value, effectively removing the depth limit. [8000]
Note that up to release 1.8, samtools would enforce a minimum value for this option. This no
longer happens and the limit is set exactly as specified.
-E,--redo-BAQ
Recalculate BAQ on the fly, ignore existing BQ tags. See BAQ below.
-f,--fasta-refFILE
The faidx-indexed reference file in the FASTA format. The file can be optionally compressed by
bgzip. [null]
Supplying a reference file will enable base alignment quality calculation for all reads aligned
to a reference in the file. See BAQ below.
-G,--exclude-RGFILE
Exclude reads from read groups listed in FILE (one @RG-ID per line)
-l,--positionsFILE
BED or position list file containing a list of regions or sites where pileup or BCF should be
generated. Position list files contain two columns (chromosome and position) and start counting
from 1. BED files contain at least 3 columns (chromosome, start and end position) and are
0-based half-open.
While it is possible to mix both position-list and BED coordinates in the same file, this is
strongly ill advised due to the differing coordinate systems. [null]
-q,--min-MQINT
Minimum mapping quality for an alignment to be used [0]
-Q,--min-BQINT
Minimum base quality for a base to be considered. [13]
Note base-quality 0 is used as a filtering mechanism for overlap removal which marks bases as
having quality zero and lets the base quality filter remove them. Hence using --min-BQ0 will
make the overlapping bases reappear, albeit with quality zero.
-r,--regionSTR
Only generate pileup in region. Requires the BAM files to be indexed. If used in conjunction
with -l then considers the intersection of the two requests. STR [all sites]
-R,--ignore-RG
Ignore RG tags. Treat all reads in one BAM as one sample.
--rf,--incl-flagsSTR|INT
Required flags: only include reads with any of the mask bits set [null]. Note this is
implemented as a filter-out rule, rejecting reads that have none of the mask bits set. Hence
this does not override the --excl-flags option.
--ff,--excl-flagsSTR|INT
Filter flags: skip reads with any of the mask bits set. This defaults to SECONDARY,QCFAIL,DUP.
The option is not accumulative, so specifying e.g. --ffQCFAIL will reenable output of
secondary and duplicate alignments. Note this does not override the --incl-flags option.
-x,--ignore-overlaps-removal,--disable-overlap-removal
Overlap detection and removal is enabled by default. This option turns it off.
When enabled, where the ends of a read-pair overlap the overlapping region will have one base
selected and the duplicate base nullified by setting its phred score to zero. It will then be
discarded by the --min-BQ option unless this is zero.
The quality values of the retained base within an overlap will be the summation of the two
bases if they agree, or 0.8 times the higher of the two bases if they disagree, with the base
nucleotide also being the higher confident call.
-X Include customized index file as a part of arguments. See EXAMPLES section for sample of usage.
OutputOptions:-o,--outputFILE
Write pileup output to FILE, rather than the default of standard output.
-O,--output-BP
Output base positions on reads in orientation listed in the SAM file (left to right).
--output-BP-5
Output base positions on reads in their original 5' to 3' orientation.
-s,--output-MQ
Output mapping qualities encoded as ASCII characters.
--output-QNAME
Output an extra column containing comma-separated read names. Equivalent to --output-extraQNAME.
--output-extraSTR
Output extra columns containing comma-separated values of read fields or read tags. The names
of the selected fields have to be provided as they are described in the SAM Specification (pag.
6) and will be output by the mpileup command in the same order as in the document (i.e. QNAME,
FLAG, RNAME,...) The names are case sensitive. Currently, only the following fields are
supported:
QNAME,FLAG,RNAME,POS,MAPQ,RNEXT,PNEXT,RLEN
Anything that is not on this list is treated as a potential tag, although only two character
tags are accepted. In the mpileup output, tag columns are displayed in the order they were
provided by the user in the command line. Field and tag names have to be provided in a comma-
separated string to the mpileup command. Tags with type B (byte array) type are not supported.
An absent or unsupported tag will be listed as "*". E.g.
samtoolsmpileup--output-extraFLAG,QNAME,RG,NMin.bam
will display four extra columns in the mpileup output, the first being a list of comma-
separated read names, followed by a list of flag values, a list of RG tag values and a list of
NM tag values. Field values are always displayed before tag values.
--output-sepCHAR
Specify a different separator character for tag value lists, when those values might contain
one or more commas (,), which is the default list separator. This option only affects columns
for two-letter tags like NM; standard fields like FLAG or QNAME will always be separated by
commas.
--output-emptyCHAR
Specify a different 'no value' character for tag list entries corresponding to reads that don't
have a tag requested with the --output-extra option. The default is *.
This option only applies to rows that have at least one read in the pileup, and only to columns
for two-letter tags. Columns for empty rows will always be printed as *.
-M,--output-mods
Adds base modification markup into the sequence column. This uses the Mm and Ml auxiliary tags
(or their uppercase equivalents). Any base in the sequence output may be followed by a series
of strandcodequality strings enclosed within square brackets where strand is "+" or "-", code
is a single character (such as "m" or "h") or a ChEBI numeric in parentheses, and quality is an
optional numeric quality value. For example a "C" base with possible 5mC and 5hmC base
modification may be reported as "C[+m179+h40]".
Quality values are from 0 to 255 inclusive, representing a linear scale of probability 0.0 to
1.0 in 1/256ths increments. If quality values are absent (no Ml tag) these are omitted, giving
an example string of "C[+m+h]".
Note the base modifications may be identified on the reverse strand, either due to the native
ability for this detection by the sequencing instrument or by the sequence subsequently being
reverse complemented. This can lead to modification codes, such as "m" meaning 5mC, being
shown for their complementary bases, such as "G[-m50]".
When --output-mods is selected base modifications can appear on any base in the sequence
output, including during insertions. This may make parsing the string more complex, so also
see the --no-output-ins-mods and --no-output-ins options to simplify this process.
--no-output-ins
Do not output the inserted bases in the sequence column. Usually this is reported as "+lengthsequence", but with this option it becomes simply "+length". For example an insertion of AGT
in a pileup column changes from "CCC+3AGTGCC" to "CCC+3GCC".
Specifying this option twice also removes the "+length" portion, changing the example above to
"CCCGCC".
The purpose of this change is to simplify parsing using basic regular expressions, which
traditionally cannot perform counting operations. It is particularly beneficial when used in
conjunction with --output-mods as the syntax of the inserted sequence is adjusted to also
report possible base modifications, but see also --no-output-ins-mods as an alternative.
--no-output-ins-mods
Outputs the inserted bases in the sequence, but excluding any base modifications. This only
affects output when --output-mods is also used.
--no-output-del
Do not output deleted reference bases in the sequence column. Normally this is reported as
"+lengthsequence", but with this option it becomes simply "+length". For example an deletion
of 3 unknown bases (due to no reference being specified) would normally be seen in a column as
e.g. "CCC-3NNNGCC", but will be reported as "CCC-3GCC" with this option.
Specifying this option twice also removes the "-length" portion, changing the example above to
"CCCGCC".
The purpose of this change is to simplify parsing using basic regular expressions, which
traditionally cannot perform counting operations. See also --no-output-ins.
--no-output-ends
Removes the “^” (with mapping quality) and “$” markup from the sequence column.
--reverse-del
Mark the deletions on the reverse strand with the character #, instead of the usual *.
-a Output all positions, including those with zero depth.
-a-a,-aa
Output absolutely all positions, including unused reference sequences. Note that when used in
conjunction with a BED file the -a option may sometimes operate as if -aa was specified if the
reference sequence has coverage outside of the region specified in the BED file.
BAQ(BaseAlignmentQuality)
BAQ is the Phred-scaled probability of a read base being misaligned. It greatly helps to reduce false
SNPs caused by misalignments. BAQ is calculated using the probabilistic realignment method described in
the paper “Improving SNP discovery by base alignment quality”, Heng Li, Bioinformatics, Volume 27, Issue
8 <https://doi.org/10.1093/bioinformatics/btr076>
BAQ is applied to modify quality values before the -Q filtering happens and before the choice of which
sequence to retain when removing overlaps.
BAQ is turned on when a reference file is supplied using the -f option. To disable it, use the -B
option.
It is possible to store precalculated BAQ values in a SAM BQ:Z tag. Samtools mpileup will use the
precalculated values if it finds them. The -E option can be used to make it ignore the contents of the
BQ:Z tag and force it to recalculate the BAQ scores by making a new alignment.
Install Now
PNG Icons