Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • mygump
    Junior Member
    • Nov 2013
    • 2

    #1

    samtools BAQ difference between 0.1.18 and 1.2 version

    Hi, there
    I compare the mpileup outputs between samtools 1.2 and 0.1.18. And I find
    some discordance of the base quality between the old and new version.
    I extract the pileup file with -B option to disable BAQ of samtools firstly. And I get a position like this:
    gi|57116681|ref|NC_000962.2| 104830 C 23 AAAAAAAAAAAAAAaAAAAAAA^IA BBBB?@@ABBB@=@?C>ABBCB;

    Then I perform the same command except the -B option, but for 0.1.18, I get this:
    gi|57116681|ref|NC_000962.2| 104830 C 23 AAAAAAAAAAAAAAaAAAAAAA^IA <<<<>>:=;<<5=@;:>==;::'
    for 1.2:
    gi|57116681|ref|NC_000962.2| 104830 C 23 AAAAAAAAAAAAAAaAAAAAAA^IA BBBB?@@ABBB@=@?C>ABBCB'

    There is a huge difference of the base quality between these two version. The old version seems change base quality a lot.
    I failed to find such information through google. I am confused about this. Did I miss something?
    Thanks
  • Brian Bushnell
    Super Moderator
    • Jan 2014
    • 2709

    #2
    Different algorithms will give different results; I don't really see what's surprising. Specifically, probability calculations are highly variable because they are subjective; so, it's not unlikely that different versions of a program that does probabilistic calculations will give vastly different results.

    For the record, my company still uses samtools 0.1.19 due to problems with more recent versions.

    Comment

    • mygump
      Junior Member
      • Nov 2013
      • 2

      #3
      Hi, Brian
      Do you mean different algorithms have been used by these two version?
      The different base quality will influence the downstream SNPs analysis. If I perform base quality filter >20 and strand filter with Varscan2. I will get the SNP at this position from old version but not 1.2 version. Is it normal?

      Thanks!

      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
      • SEQadmin2
        Advanced Sequencing Platforms Tackle Neuroscience’s Toughest Genomics Problems
        by SEQadmin2



        Genomics studies in neuroscience face a special challenge due to the brain’s complexity and scarcity of samples. Mapping changes in cell type and state using conventional next-generation sequencing methods remains challenging. Advances in technologies like single-cell sequencing, spatial transcriptomics, and long-read sequencing have opened the door to deeper studies of the brain and diseases like Alzheimer’s, amyotrophic lateral sclerosis (ALS), and schizophrenia.
        ...
        07-09-2026, 11:10 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by SEQadmin2, 08-03-2026, 10:13 AM
      0 responses
      21 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 07-31-2026, 02:55 AM
      0 responses
      35 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 07-24-2026, 12:17 PM
      0 responses
      25 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 07-23-2026, 11:41 AM
      0 responses
      21 views
      0 reactions
      Last Post SEQadmin2  
      Working...