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

      • kinpict.com
        Reply to 16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
        by kinpict.com
        The part about this topic feels surprisingly grounded, and the details on the page is the bit that made me slow down. I’d keep <a href="https://kinpict.com/">ai family photo generator</a> bookmarked for later.
        Today, 06:37 PM
      • SEQadmin2
        16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
        by SEQadmin2





        NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification,...
        10-30-2026, 06:46 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by SEQadmin2, 11-12-2026, 07:59 AM
      0 responses
      27 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 11-05-2026, 09:25 AM
      0 responses
      35 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 11-03-2026, 11:54 AM
      2 responses
      56 views
      0 reactions
      Last Post kinpict.com  
      Started by SEQadmin2, 10-28-2026, 09:57 AM
      0 responses
      175 views
      0 reactions
      Last Post SEQadmin2  
      Working...