Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • bwa quality trimming and samtools rmdup

    I am trying to use bwa and samtools to map 76 bp reads from multiple bacterial strains back onto a reference sequence, with the ultimate goal of extracting SNP frequencies. To obtain accurate variant frequencies, it is important to me to remove PCR duplicates. It occurred to me that quality trimming reads with BWA using the "-q" flag during alignment could affect how well rmdup works downstream. Say 2 reads are PCR duplicates, but one is rather low-quality and is trimmed to a different length than the other. The alignment start and stop positions would no longer be the same for these duplicates, and they would not be filtered using samtools rmdup, which requires identical external coordinates.

    Is this right? Does BWA do hard trimming of reads with the "-q" flag? Or does it ignore low quality bases when calculating the alignment, but still use them to determine alignment coordinates?

    While I'm at it, could anyone explain how the BWA quality trimmer works? Please don't say "read the man page," I did, and was very confused by the explanation there:
    Code:
    -q INT 	
    Parameter for read trimming. BWA trims a read down to argmax_x{\sum_{i=x+1}^l(INT-q_i)} if q_l<INT
    where l is the original read length. [0]
    Noob question #3: Does rmdup work for single read data or doesn't it? The samtools man page states:
    Code:
    Samtools’ rmdup does not work for single-end data and does not remove duplicates across chromosomes. Picard is better.
    However, there is an option "-s" to use rmdup on single read data, and when I apply it to an alignment, it looks to reduce the size of the resulting bam file. The man page and the command itself don't agree!
    Last edited by greigite; 02-16-2010, 03:33 PM. Reason: new question

  • #2
    did you get the answer? nobody know the answer?

    Comment


    • #3
      I did not get an answer- now I use Mosaik instead of bwa, it has much better documentation.

      Comment


      • #4
        Explanation of BWA read trimming

        The BWA trimming feature seems to be explained a little more clearly here: http://seqanswers.com/forums/showthread.php?t=6251 . The real C source code is in the function bwa_trim_read() in the file bwaseqio.c, but I found the comments and variable names of the Perl example referenced in the other thread more clear.

        Comment


        • #5
          Also see the SolexaQA FAQ for an enlightening discussion of the bwa algorithm vs. SolexaQA algorithm.

          Comment


          • #6
            Originally posted by greigite View Post
            I did not get an answer- now I use Mosaik instead of bwa, it has much better documentation.
            Where is the real bottleneck ;-)
            Homepage: Dan Bolser
            MetaBase the database of biological databases.

            Comment


            • #7
              Duplicate Marking &amp; Trimming

              Originally posted by greigite View Post
              I am trying to use bwa and samtools to map 76 bp reads from multiple bacterial strains back onto a reference sequence, with the ultimate goal of extracting SNP frequencies. To obtain accurate variant frequencies, it is important to me to remove PCR duplicates. It occurred to me that quality trimming reads with BWA using the "-q" flag during alignment could affect how well rmdup works downstream. Say 2 reads are PCR duplicates, but one is rather low-quality and is trimmed to a different length than the other. The alignment start and stop positions would no longer be the same for these duplicates, and they would not be filtered using samtools rmdup, which requires identical external coordinates.

              Is this right? Does BWA do hard trimming of reads with the "-q" flag? Or does it ignore low quality bases when calculating the alignment, but still use them to determine alignment coordinates?

              While I'm at it, could anyone explain how the BWA quality trimmer works? Please don't say "read the man page," I did, and was very confused by the explanation there:
              Code:
              -q INT 	
              Parameter for read trimming. BWA trims a read down to argmax_x{\sum_{i=x+1}^l(INT-q_i)} if q_l<INT
              where l is the original read length. [0]
              Noob question #3: Does rmdup work for single read data or doesn't it? The samtools man page states:
              Code:
              Samtools’ rmdup does not work for single-end data and does not remove duplicates across chromosomes. Picard is better.
              However, there is an option "-s" to use rmdup on single read data, and when I apply it to an alignment, it looks to reduce the size of the resulting bam file. The man page and the command itself don't agree!
              I have the same question. I know it has been a long time since you wrote this post, but I am wondering what you ended up doing. I am trying to develop a good pipeline for doing similar analyses.

              Comment


              • #8
                Originally posted by mrood View Post
                I have the same question. I know it has been a long time since you wrote this post, but I am wondering what you ended up doing. I am trying to develop a good pipeline for doing similar analyses.
                Which question? He posted 3. I found a good description of the protocol here:
                Download SAM tools for free. SAM (Sequence Alignment/Map) is a flexible generic format for storing nucleotide sequence alignment. SAMtools provide efficient utilities on manipulating alignments in the SAM format.
                Homepage: Dan Bolser
                MetaBase the database of biological databases.

                Comment

                Latest Articles

                Collapse

                • seqadmin
                  Strategies for Sequencing Challenging Samples
                  by seqadmin


                  Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                  03-22-2024, 06:39 AM
                • seqadmin
                  Techniques and Challenges in Conservation Genomics
                  by seqadmin



                  The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                  Avian Conservation
                  Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                  03-08-2024, 10:41 AM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, Yesterday, 06:37 PM
                0 responses
                10 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, Yesterday, 06:07 PM
                0 responses
                9 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-22-2024, 10:03 AM
                0 responses
                51 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-21-2024, 07:32 AM
                0 responses
                67 views
                0 likes
                Last Post seqadmin  
                Working...
                X