1 .TH samtools 1 "21 May 2009" "samtools-0.1.4" "Bioinformatics tools"
4 samtools - Utilities for the Sequence Alignment/Map (SAM) format
7 samtools view -bt ref_list.txt -o aln.bam aln.sam.gz
9 samtools sort aln.bam aln.sorted
11 samtools index aln.sorted.bam
13 samtools view aln.sorted.bam chr2:20,100,000-20,200,000
15 samtools merge out.bam in1.bam in2.bam in3.bam
17 samtools faidx ref.fasta
19 samtools pileup -f ref.fasta aln.sorted.bam
21 samtools tview aln.sorted.bam ref.fasta
25 Samtools is a set of utilities that manipulate alignments in the BAM
26 format. It imports from and exports to the SAM (Sequence
27 Alignment/Map) format, does sorting, merging and indexing, and
28 allows to retrieve reads in any regions swiftly.
30 .SH COMMANDS AND OPTIONS
33 samtools import <in.ref_list> <in.sam> <out.bam>
35 Since 0.1.4, this command is an alias of:
37 samtools view -bt <in.ref_list> -o <out.bam> <in.sam>
41 samtools sort [-n] [-m maxMem] <in.bam> <out.prefix>
43 Sort alignments by leftmost coordinates. File
45 will be created. This command may also create temporary files
46 .I <out.prefix>.%d.bam
47 when the whole alignment cannot be fitted into memory (controlled by
54 Sort by read names rather than by chromosomal coordinates
57 Approximately the maximum required memory. [500000000]
62 samtools merge [-n] <out.bam> <in1.bam> <in2.bam> [...]
64 Merge multiple sorted alignments. The header of
68 and the headers of other files will be ignored.
74 The input alignments are sorted by read names rather than by chromosomal
80 samtools index <aln.bam>
82 Index sorted alignment for fast random access. Index file
88 samtools view [-bhHS] [-t in.refList] [-o output] [-f reqFlag] [-F
89 skipFlag] [-q minMapQ] <in.bam> [region1 [...]]
91 Extract/print all or sub alignments in SAM or BAM format. If no region
92 is specified, all the alignments will be printed; otherwise only
93 alignments overlapping with the specified regions will be output. An
94 alignment may be given multiple times if it is overlapping several
95 regions. A region can be presented, for example, in the following
96 format: `chr2', `chr2:1000000' or `chr2:1,000,000-2,000,000'.
102 Output in the BAM format.
105 Include the header in the output.
108 Output the header only.
111 Input is in SAM. If @SQ header lines are absent, the
116 This file is TAB-delimited. Each line must contain the reference name
117 and the length of the reference, one line for each distinct reference;
118 additional fields are ignored. This file also defines the order of the
119 reference sequences in sorting. If you run `samtools faidx <ref.fa>',
120 the resultant index file
130 Only output alignments with all bits in INT present in the FLAG
131 field. INT can be in hex in the format of /^0x[0-9A-F]+/ [0]
134 Skip alignments with bits present in INT [0]
137 Skip alignments with MAPQ smaller than INT [0]
142 samtools faidx <ref.fasta> [region1 [...]]
144 Index reference sequence in the FASTA format or extract subsequence from
145 indexed reference sequence. If no region is specified,
147 will index the file and create
149 on the disk. If regions are speficified, the subsequences will be
150 retrieved and printed to stdout in the FASTA format. The input file can
157 samtools pileup [-f in.ref.fasta] [-t in.ref_list] [-l in.site_list]
158 [-iscg] [-T theta] [-N nHap] [-r pairDiffRate] <in.alignment>
160 Print the alignment in the pileup format. In the pileup format, each
161 line represents a genomic position, consisting of chromosome name,
162 coordinate, reference base, read bases, read qualities and alignment
163 mapping qualities. Information on match, mismatch, indel, strand,
164 mapping quality and start and end of a read are all encoded at the read
165 base column. At this column, a dot stands for a match to the reference
166 base on the forward strand, a comma for a match on the reverse strand,
167 `ACGTN' for a mismatch on the forward strand and `acgtn' for a mismatch
168 on the reverse strand. A pattern `\\+[0-9]+[ACGTNacgtn]+' indicates
169 there is an insertion between this reference position and the next
170 reference position. The length of the insertion is given by the integer
171 in the pattern, followed by the inserted sequence. Similarly, a pattern
172 `-[0-9]+[ACGTNacgtn]+' represents a deletion from the reference. The
173 deleted bases will be presented as `*' in the following lines. Also at
174 the read base column, a symbol `^' marks the start of a read segment
175 which is a contiguous subsequence on the read separated by `N/S/H' CIGAR
176 operations. The ASCII of the character following `^' minus 33 gives the
177 mapping quality. A symbol `$' marks the end of a read segment.
181 is applied, the consensus base, consensus quality, SNP quality and RMS
182 mapping quality of the reads covering the site will be inserted between
183 the `reference base' and the `read bases' columns. An indel occupies an
184 additional line. Each indel line consists of chromosome name,
185 coordinate, a star, the genotype, consensus quality, SNP quality, RMS
186 mapping quality, # covering reads, the first alllele, the second allele,
187 # reads supporting the first allele, # reads supporting the second
188 allele and # reads containing indels different from the top two alleles.
195 Print the mapping quality as the last column. This option makes the
196 output easier to parse, although this format is not space efficient.
200 Only output pileup lines containing indels.
204 The reference sequence in the FASTA format. Index file
211 List of reference names ane sequence lengths, in the format described
214 command. If this option is present, samtools assumes the input
216 is in SAM format; otherwise it assumes in BAM format.
220 List of sites at which pileup is output. This file is space
221 delimited. The first two columns are required to be chromosome and
222 1-based coordinate. Additional columns are ignored. It is
223 recommended to use option
227 as in the default format we may not know the mapping quality.
231 Call the consensus sequence using MAQ consensus model. Options
236 are only effective when
242 Generate genotype likelihood in the binary GLFv3 format. This option
243 suppresses -c, -i and -s.
247 The theta parameter (error dependency coefficient) in the maq consensus
252 Number of haplotypes in the sample (>=2) [2]
256 Expected fraction of differences between a pair of haplotypes [0.001]
262 samtools tview <in.sorted.bam> [ref.fasta]
264 Text alignment viewer (based on the ncurses library). In the viewer,
265 press `?' for help and press `g' to check the alignment start from a
266 region in the format like `chr10:10,000,000'. Note that if the region
267 showed on the screen contains no mapped reads, a blank screen will be
268 seen. This is a known issue and will be improved later.
274 samtools fixmate <in.nameSrt.bam> <out.bam>
276 Fill in mate coordinates, ISIZE and mate related flags from a
277 name-sorted alignment.
281 samtools rmdup <input.srt.bam> <out.bam>
283 Remove potential PCR duplicates: if multiple read pairs have identical
284 external coordinates, only retain the pair with highest mapping quality.
287 works with FR orientation and requires ISIZE is correctly set.
293 samtools rmdupse <input.srt.bam> <out.bam>
295 Remove potential duplicates for single-ended reads. This command will
296 treat all reads as single-ended even if they are paired in fact.
302 samtools fillmd [-e] <aln.bam> <ref.fasta>
304 Generate the MD tag. If the MD tag is already present, this command will
305 give a warning if the MD tag generated is different from the existing
312 Convert a the read base to = if it is identical to the aligned reference
313 base. Indel caller does not support the = bases at the moment.
320 SAM is TAB-delimited. Apart from the header lines, which are started
321 with the `@' symbol, each alignment line consists of:
327 Col Field Description
329 1 QNAME Query (pair) NAME
331 3 RNAME Reference sequence NAME
332 4 POS 1-based leftmost POSition/coordinate of clipped sequence
333 5 MAPQ MAPping Quality (Phred-scaled)
334 6 CIAGR extended CIGAR string
335 7 MRNM Mate Reference sequence NaMe (`=' if same as RNAME)
336 8 MPOS 1-based Mate POSistion
337 9 ISIZE Inferred insert SIZE
338 10 SEQ query SEQuence on the same strand as the reference
339 11 QUAL query QUALity (ASCII-33 gives the Phred base quality)
340 12 OPT variable OPTional fields in the format TAG:VTYPE:VALUE
344 Each bit in the FLAG field is defined as:
352 0x0001 the read is paired in sequencing
353 0x0002 the read is mapped in a proper pair
354 0x0004 the query sequence itself is unmapped
355 0x0008 the mate is unmapped
356 0x0010 strand of the query (1 for reverse)
357 0x0020 strand of the mate
358 0x0040 the read is the first read in a pair
359 0x0080 the read is the second read in a pair
360 0x0100 the alignment is not primary
361 0x0200 the read fails platform/vendor quality checks
362 0x0400 the read is either a PCR or an optical duplicate
368 Unaligned words used in bam_import.c, bam_endian.h, bam.c and bam_aux.c.
370 CIGAR operation P is not properly handled at the moment.
372 The text viewer mysteriously crashes in a very rare case.
376 Heng Li from the Sanger Institute wrote the C version of samtools. Bob
377 Handsaker from the Broad Institute implemented the BGZF library and Jue
378 Ruan from Beijing Genomics Institute wrote the RAZF library. Various
379 people in the 1000Genomes Project contributed to the SAM format
384 Samtools website: http://samtools.sourceforge.net