Pileup format explained
Pileup format is a text-based format for summarizing the base calls of aligned reads to a reference sequence. This format facilitates visual display of SNP/indel calling and alignment. It was first used byTony Cox and Zemin Ning at the Wellcome Trust Sanger Institute, and became widely known through its implementation within the SAMtools software suite.[1]
Format
Example
Sequence | Position | Reference Base | Read Count | Read Results | Quality |
---|
seq1 | 272 | T | 24 | ,.$.....,,.,.,...,,,.,..^+. | |
seq1 | 273 | T | 23 | ,.....,,.,.,...,,,.,..A | |
seq1 | 274 | T | 23 | ,.$....,,.,.,...,,,.,... | |
seq1 | 275 | A | 23 | ,$....,,.,.,...,,,.,...^l. | |
seq1 | 276 | G | 22 | ...T,,.,.,...,,,.,.... | |
seq1 | 277 | T | 22 | ....,,.,.,.C.,,,.,..G. | |
seq1 | 278 | G | 23 | ....,,.,.,...,,,.,....^k. | |
seq1 | 279 | C | 23 | A..T,,.,.,...,,,.,..... | | |
The columns
Each line consists of 5 (or optionally 6) tab-separated columns:
- Sequence identifier
- Position in sequence (starting from 1)
- Reference nucleotide at that position
- Number of aligned reads covering that position (depth of coverage)
- Bases at that position from aligned reads
- Phred Quality of those bases, represented in ASCII with -33 offset (OPTIONAL)
Column 5: The bases string
- . (dot) means a base that matched the reference on the forward strand
- , (comma) means a base that matched the reference on the reverse strand
- > (less-/greater-than sign) denotes a reference skip. This occurs, for example, if a base in the reference genome is intronic and a read maps to two flanking exons. If quality scores are given in a sixth column, they refer to the quality of the read and not the specific base.
- AGTCN (upper case) denotes a base that did not match the reference on the forward strand
- agtcn (lower case) denotes a base that did not match the reference on the reverse strand
- A sequence matching the regular expression denotes an insertion of one or more bases starting from the next position. For example, +2AG means insertion of AG in the forward strand
- A sequence matching the regular expression denotes a deletion of one or more bases starting from the next position. For example, -2ct means deletion of CT in the reverse strand
- ^ (caret) marks the start of a read segment and the ASCII of the character following `^' minus 33 gives the mapping quality
- $ (dollar) marks the end of a read segment
- * (asterisk) is a placeholder for a deleted base in a multiple basepair deletion that was mentioned in a previous line by the notation
Column 6: The base quality string
This is an optional column. If present, the ASCII value of the character minus 33 gives the mapping Phred quality of each of the bases in the previous column 5. This is similar to quality encoding in the FASTQ format.
File extension
There is no standard file extension for a Pileup file, but .msf (multiple sequence file), .pup[2] and .pileup[3] [4] are used.
See also
External links
Notes and References
- 10.1093/bioinformatics/btp352 . 2009 . Bioinformatics . 25 . 2078–2079 . The Sequence alignment/map (SAM) format and SAMtools. Li H. . Handsaker B. . Wysoker A. . Fennell T. . Ruan J.. Homer N. . Marth G. . Abecasis G. . Durbin R . 1000 Genome Project Data Processing Subgroup (2009). 16 . 19505943 . 2723002 .
- Web site: QUANTA: Protein Design. 3. Reading and Writing Sequence Data Files . . 1998-10-02 . . 2020-03-27.
- Glez-Peña . Daniel . Gómez-López . Gonzalo . Reboiro-Jato . Miguel . Fdez-Riverola . Florentino . Pisano . David G . 2011-01-24 . PileLine: a toolbox to handle genome position information in next-generation sequencing studies . BMC Bioinformatics . 12 . 31 . 10.1186/1471-2105-12-31 . 1471-2105 . 3037855 . 21261974 . free .
- Web site: Chisom . Halimat . 2023-03-31 . File Formats Every Bioinformatician — Established or Upcoming — Must Know (and then some) . 2023-11-11 . Medium . en.