This is a test version of Biostars. For the public version, visit https://www.biostars.org.
My Pileup File Doesn'T Match The Documented Format? What Are These Columns?

Hi everyone,

I'm new to Bioinformatics and I'm playing around with some existing code and I'm trying to make sense of the pileup file it's generating.

I'm using Samtools version 0.1.7 and I'm running the pileup command like this:

> samtools pileup -f [fasta file] [another file I'm not sure what it is] -c > my.pileup

My file looks like this:

2L    4043    T    T    108    0    37    27    ...........................    CCCCCCCCCCCCBCCCCCCCCCCCCCC
2L    4044    T    T    108    0    37    27    ...........................    CCCCCCCCCCCCCCCCCCCCCCCCCCC
2L    4045    T    T    81    0    37    27    .................C.........    CCCCCCCCCCCCCCCCC<CCCCCCCCC
2L    4046    C    C    108    0    37    27    ...........................    CCCCCCCCCCCCCCCCCCCCCCCCCAC

But the documentation says it will look like this:

seq1 272 T 24  ,.$.....,,.,.,...,,,.,..^+. <<<+;<<<<<<<<<<<=<;<;7<&
seq1 273 T 23  ,.....,,.,.,...,,,.,..A <<<;<<<<<<<<<3<=<<<;<<+
seq1 274 T 23  ,.$....,,.,.,...,,,.,...    7<7;<;<<<<<<<<<=<;<;<<6
seq1 275 A 23  ,$....,,.,.,...,,,.,...^l.  <+;9*<<<<<<<<<=<<:;<<<<
seq1 276 G 22  ...T,,.,.,...,,,.,....  33;+<<7=7<<7<&<<1;<<6<
seq1 277 T 22  ....,,.,.,.C.,,,.,..G.  +7<;<<<<<<<&<=<<:;<<&<
seq1 278 G 23  ....,,.,.,...,,,.,....^k.   %38*<<;<7<<7<=<<<;<<<<<
seq1 279 C 23  A..T,,.,.,...,,,.,..... ;75&<<<<<<<<<=<<<9<<:<<

Would anyone know why there are differences? And more importantly what are the additional columns?

samtools pileup mpileup
• 2,821 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by Applet • Applet

1 answer

You have used the -c option: it changes the output format. See the (old) FAQ.

Also, pileup is deprecated you should now use mpileup

• 1 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this
  • Unusual samtools mpileup output?
    written by quachtina96 • quachtina96
    Hi, I am using samtools mpileup to generate a pileup file, and my output doesn't match the example pileup files I've seen (e.g. the example …
  • Imitate Pileup Format
    written by Lee Katz Lee Katz
    <p>Hi, I am fixing a script that uses an older version of <a href='http://samtools.sourceforge.net/'>samtools</a>. I just want to fix the system calls to <a href='http://samtools.sourceforge.net/'>samtools</a>. …
  • Creating A Simple Matrix With Perl
    written by loicatraile • loicatraile
    <p>Hi! I´m trying to build a 4x4 nucleotide Count Matrix of a sequence alignment and im having some trouble to display it correctly as a …
  • Split Fastq Paired-Read Using Awk
    written by thiagomafra thiagomafra
    <p>Hi guys,</p> <p>I have a fastq file that contains paired reads. I know that F3 and R3 tag beginning with G and T, respectively in …
  • Heterogeneous Snp
    written by polangxin • polangxin
    <p>i'm using <a href='http://bio-bwa.sourceforge.net/'>bwa</a> to mapping NGS data to reference, and want to find position with heterogeneous data: for example:</p> <pre><code>seq1 274 T 23 AAAAAAACCCCCCC …
  • Calculating Coverage From Pileup File To Find Gene Duplication Events
    written by thecuriousbiologist thecuriousbiologist
    <p>Hello, </p> <p>I have a pileup file like below :</p> <pre><code>seq1 272 T 24 ,.$.....,,.,.,...,,,.,..^+. &lt;&lt;&lt;+;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;=&lt;;&lt;;7&lt;&amp; seq1 273 T 23 ,.....,,.,.,...,,,.,..A &lt;&lt;&lt;;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;3&lt;=&lt;&lt;&lt;;&lt;&lt;+ seq1 274 T …
  • Samtools Pileup Got Really Slow After I Switched To Stampy For Mapping
    written by Applet • Applet
    <p>So here is my perplexing question.</p> <p>When I used to run my pipeline on the same input files but using BWA for the mapping, the …
  • Samtools - Calmd Errors - Different Nm For Read
    written by Applet • Applet
    <p>I'm trying to put MD tags into my BAM file by calling samtools calmd like so: </p> <pre><code>samtools calmd -u update_reads-aligned-parent2.sorted.bam parent2.fa &gt; update_reads-aligned-parent2.sorted.calmd.bam </code></pre> …
  • Help Understanding Pileup Format
    written by Pi • Pi
    <p>Hello</p> <p>I have been given some data in pileup format. Whilst it is beyond the scope of my work to go into the details of …
  • How To Process Raw Pileup
    written by Rm Rm
    <blockquote> <p>Raw pileup file some times before indel (*) have multiple insert (+4caaa) or deletion (-13GGCGCGCGTGCGC) strings in the read colum (column9) </p> <p>how to …