Unconfigured Ad

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

    #1

    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
              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, 10-22-2026, 11:13 AM
            0 responses
            23 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 10-22-2026, 10:48 AM
            0 responses
            20 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...