Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • the number of aligned reads by maq and bwa

    We have applied maq and bwa to align the same production data, and we see an obvious difference on the number of aligned reads between maq and bwa, especially a decrease in the number with longer read length.

    This issue is critical, is it a drawback of bwa that aligns much less reads than maq? I am using the default options for bwa running, anyone has experience to increase the aligned reads dramatically by using different bwa options?

    I don't know if it is a bottleneck of BWT algorithm even it is much faster than maq. If so, I might have to think again if I am really in favor of bwa. Now bwa is pushed forward, anyone has comments of this issue?

    thanks

  • #2
    have you analysed your data further downstream ?

    you are better off with 40% mapped sequences and good results than with 60% mapped rubbish.

    what is the type of experiment you are doing at the moment? You might consider tweaking the options a bit in function of your experiment. If you are mRNA expression counting, you might allow lower quality reads to be mapped because both the control and sample of interest will have the same amount of 'rubbish match'. When you are SNP hunting you might use more stringent parameters to get rid of false positives etc.

    You can always set filters to have almost everything mapped, but I do not think that is the way you want to go?

    Comment


    • #3
      BWA even allows gaps, so essentially more reads should be mapped than those from gapless MAQ alignment.
      --
      bioinfosm

      Comment


      • #4
        maq is more tolerant with the mismatches toward the 3'-end. Given typical sequencing error rate, the mapping ratio of maq and bwa is quite similar to each other. However, when the sequencing error rate is high, maq maps more reads. At the same time, one may consider to run bwa with "bwa aln -q15" to enable base trimming. This only works when the base quality is not rubbish.

        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
        8 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, Yesterday, 06:07 PM
        0 responses
        8 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-22-2024, 10:03 AM
        0 responses
        49 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