Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • hellohenry
    Junior Member
    • Mar 2012
    • 4

    BWA parameters affects SNP calling with samtools?

    hello everyone, I am new in using BWA and samtools, but I found something really strange: I used two different parameter settings of BWA for mapping, and I got two mapping results. Then I used samtools to call SNP from the two mapping results respectively, got SNP collection A and B. a SNP X presents in A while absent in B. But the real problem is: when I browsered the mapping results with "samtools tview", I found the locus X is a real SNP in both mapping results! So my problem is : the BWA parameter did not affected mapping results, but affected snp calling with samtools!
  • hellohenry
    Junior Member
    • Mar 2012
    • 4

    #2
    appreciaton

    all replies are greatly appreciated!

    Comment

    • dpryan
      Devon Ryan
      • Jul 2011
      • 3478

      #3
      It's unlikely that the parameters you feed to bwa would change the results of SNP calling without also resulting in different mapping. In the case that you're describing, was there a difference in depth at the position in question between the SAM files? That would seem to be the most likely cause.

      Comment

      • hellohenry
        Junior Member
        • Mar 2012
        • 4

        #4
        Originally posted by dpryan View Post
        It's unlikely that the parameters you feed to bwa would change the results of SNP calling without also resulting in different mapping. In the case that you're describing, was there a difference in depth at the position in question between the SAM files? That would seem to be the most likely cause.
        Thanks!~ I have checked the depth at the locus on both mapping results, one is 241, the other 243. So it is unlikely to be the reason....

        Comment

        • hellohenry
          Junior Member
          • Mar 2012
          • 4

          #5
          OK, I got it, finally~~ when I applied -BE when calling variants with the bcftools, the missing locus could be detected again. Inspired by this post: http://seqanswers.com/forums/showthread.php?t=16786.

          Comment

          • wzhjlau2009
            Junior Member
            • Jun 2011
            • 5

            #6
            i also meet the same problem! may i have a chat with you? hellohenry.my email:[email protected] QQ:872664567

            Comment

            Latest Articles

            Collapse

            • 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
            • SEQadmin2
              Cancer Drug Resistance: The Lingering Barrier to Rising Survival
              by SEQadmin2



              Cancer survival rates have significantly increased in the last few decades in the United States, reaching a combined 70% 5-year survival rate by 2021. Behind this number, there are years of research to find new therapies, drug targets, and early detection methods. But there is one core challenge that keeps slowing down these advances, and it’s about drug resistance.

              There is no single reason why many patients don’t respond to treatment as expected. Cancer is...
              07-08-2026, 05:17 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by SEQadmin2, 07-20-2026, 11:10 AM
            0 responses
            20 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 07-13-2026, 10:26 AM
            0 responses
            33 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 07-09-2026, 10:04 AM
            0 responses
            44 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 07-08-2026, 10:08 AM
            0 responses
            30 views
            0 reactions
            Last Post SEQadmin2  
            Working...