Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Compare Bowtie2 and BWA summary

    How to compare bowtie and bwa summary?

    Here is the output of bowtie2:
    28115453 reads; of these:
    28115453 (100.00%) were paired; of these:
    9177458 (32.64%) aligned concordantly 0 times
    15185809 (54.01%) aligned concordantly exactly 1 time
    3752186 (13.35%) aligned concordantly >1 times
    ----
    9177458 pairs aligned concordantly 0 times; of these:
    2358270 (25.70%) aligned discordantly 1 time
    ----
    6819188 pairs aligned 0 times concordantly or discordantly; of these:
    13638376 mates make up the pairs; of these:
    9978644 (73.17%) aligned 0 times
    1757628 (12.89%) aligned exactly 1 time
    1902104 (13.95%) aligned >1 times
    82.25% overall alignment rate

    How to compare the information between BWA and bowtie2 for the same file.
    I know that

    For multiples hits with BWA, you have XT:A:R flag (R for multiple, you can have XT:A:U for unique reads. For paired-end reads, you might also consider XT:A:M (one-mate recovered) which means that one of the pairs is uniquely mapped and the other isn't)
    For Bowtie2, XS:i flag determined if alignment is unique or not. Only present if the SAM record is for an aligned read and more than one alignment was found for the read.

    I am using paired end sequencing. How to compare and verify if the summary from both bowtie2 and bwa is same

  • #2
    Please forgot the entire concept of a "unique alignment"; it's misleading.

    A proper comparison would be to look at concordance/discordance as a function of MAPQ. See this thread, starting around comment 8 for an example. You could also look at how often paired reads were aligned as singletons ("mate-recovered" in bwa parlance) in one aligner and not the other.

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Current Approaches to Protein Sequencing
      by seqadmin


      Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
      04-04-2024, 04:25 PM
    • 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

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 04-11-2024, 12:08 PM
    0 responses
    30 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 10:19 PM
    0 responses
    32 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 09:21 AM
    0 responses
    28 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-04-2024, 09:00 AM
    0 responses
    53 views
    0 likes
    Last Post seqadmin  
    Working...
    X