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:



      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
        Beyond CRISPR/Cas9: Understand, Choose, and Use the Right Genome Editing Tool
        by SEQadmin2



        CRISPR/Cas9 sparked the gene editing revolution for both research and therapeutics.1 But this system still showed severe issues that limited its applications. The most prominent were the heavy reliance on PAM sequences, delivery limitations, double-stranded breaks that prompt unintended edits and cell death, and editing inefficiency (both in targeting and in knock-in reliability).

        Despite this, “CRISPR helped turn genome editing from a specialized technique into
        ...
        07-31-2026, 11:01 AM
      • SEQadmin2
        Proteomic Platforms: How to Choose the Right Analytical Strategy to Improve Detection and Clinical Applications
        by SEQadmin2


        Proteomics platforms are evolving rapidly, with advances in mass spectrometry and affinity-based approaches expanding what researchers can detect and at what scale. As the field moves toward deeper proteome coverage and clinical applications, scientists face an increasingly complex landscape of tools. This article will explore how researchers are navigating these choices to find the right platform for their work.

        The systematic characterization of the human proteome has
        ...
        07-20-2026, 11:48 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by SEQadmin2, 08-13-2026, 12:22 PM
      0 responses
      25 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 08-11-2026, 10:35 AM
      0 responses
      20 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 08-06-2026, 07:41 AM
      0 responses
      34 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 08-03-2026, 10:13 AM
      0 responses
      51 views
      0 reactions
      Last Post SEQadmin2  
      Working...