Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • 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!

  • #2
    appreciaton

    all replies are greatly appreciated!

    Comment


    • #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


      • #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


        • #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


          • #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

            • seqadmin
              Advancing Precision Medicine for Rare Diseases in Children
              by seqadmin




              Many organizations study rare diseases, but few have a mission as impactful as Rady Children’s Institute for Genomic Medicine (RCIGM). “We are all about changing outcomes for children,” explained Dr. Stephen Kingsmore, President and CEO of the group. The institute’s initial goal was to provide rapid diagnoses for critically ill children and shorten their diagnostic odyssey, a term used to describe the long and arduous process it takes patients to obtain an accurate...
              12-16-2024, 07:57 AM
            • seqadmin
              Recent Advances in Sequencing Technologies
              by seqadmin



              Innovations in next-generation sequencing technologies and techniques are driving more precise and comprehensive exploration of complex biological systems. Current advancements include improved accessibility for long-read sequencing and significant progress in single-cell and 3D genomics. This article explores some of the most impactful developments in the field over the past year.

              Long-Read Sequencing
              Long-read sequencing has seen remarkable advancements,...
              12-02-2024, 01:49 PM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, 12-17-2024, 10:28 AM
            0 responses
            33 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 12-13-2024, 08:24 AM
            0 responses
            48 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 12-12-2024, 07:41 AM
            0 responses
            34 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 12-11-2024, 07:45 AM
            0 responses
            46 views
            0 likes
            Last Post seqadmin  
            Working...
            X