Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • fatakias
    Member
    • Apr 2011
    • 11

    #1

    samtools and IGV question

    I am new to samtools, IGV display and NGS in general. I have a two part question.

    Part 1-
    I would like to know how to interpret the output from "samtools mpileup" command. For example, in the attached txt file, column 5 has symbols like:
    1- "^]," what does this denote?
    2- "t" or "T" what is the difference between lower case and upper case?
    3- What is "-C"?
    4- why is "." or "," is used?
    and so on.

    Is there a comprehensive list where I can find all such information. Samtools FAQ has limited information. Any help is appreciated.
    (a text file is attached herein that I am trying to interpret)

    Part 2-
    Is there a samtools script or easily available script that can parse the information and tell us the allele frequency / total number of A,T,G, C at each mapped position. An icon displays frequency when mouse is moved over mapped position. It would be helpful if I can extract the information from the command line .. what IGV tools function / application would be used for this.

    Thanks in advance.
    Attached Files
    Last edited by fatakias; 01-16-2012, 07:51 AM. Reason: 2nd attachment
  • Robby
    Member
    • Mar 2011
    • 68

    #2
    regarding part 1:

    1. a symbol `^' marks the start of a read segment which is a contiguous subsequence on the read separated by `N/S/H' CIGAR operations. The ASCII of the character following `^' minus 33 gives the mapping quality.
    2. T for a mismatch on the forward strand and t for a mismatch on the reverse strand
    3. -C -> deletion of C
    4. a dot stands for a match to the reference base on the forward strand, a comma for a match on the reverse strand

    Comment

    • swbarnes2
      Senior Member
      • May 2008
      • 910

      #3
      This is the link to the page that explains pileup format:

      http://samtools.sourceforge.net/pileup.shtml

      You can get allele counts from the DP4 values in the vcf, which is simpler than parsing the pileup. And its a much smaller file, too.

      Comment

      Latest Articles

      Collapse

      • SEQadmin2
        Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
        by SEQadmin2



        Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

        Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
        ...
        10-14-2026, 09:55 AM
      • SEQadmin2
        New Genomics Technologies Take Aim at Long-Standing Limits
        by SEQadmin2


        Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

        We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
        ...
        09-28-2026, 10:25 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by SEQadmin2, Yesterday, 11:13 AM
      0 responses
      19 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, Yesterday, 10:48 AM
      0 responses
      18 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-15-2026, 10:52 AM
      0 responses
      36 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-14-2026, 05:39 AM
      0 responses
      42 views
      0 reactions
      Last Post SEQadmin2  
      Working...