Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • GUZUMO
    Junior Member
    • Nov 2016
    • 7

    how to count the mapped reads with bbmap?

    hello, I have assembled paired reads with Megahit, and now I want to know how many reads are uesd to make up those contigs(reads%).
    I have tried bbmap and find this in STR:
    Pairing data: pct reads num reads pct bases num bases

    mated pairs: 36.1655% 52614378 37.2210% 15687793904
    bad pairs: 2.9658% 4314738 3.0270% 1275797796
    insert size avg: 338.64


    Read 1 data: pct reads num reads pct bases num bases

    mapped: 46.8503% 68158783 46.8339% 10145704082
    unambiguous: 46.3132% 67377420 46.2987% 10029749277
    ambiguous: 0.5371% 781363 0.5353% 115954805
    low-Q discards: 0.0000% 0 0.0000% 0

    perfect best site: 18.9679% 27594901 19.0008% 4116175782
    semiperfect site: 21.5567% 31361110 21.5960% 4678380777
    rescued: 0.7859% 1143367

    Match Rate: NA NA 95.8928% 9735647890
    Error Rate: 41.6441% 36902108 2.0807% 211247627
    Sub Rate: 41.6137% 36875119 1.9638% 199375874
    Del Rate: 0.7019% 622010 0.0682% 6927772
    Ins Rate: 0.8580% 760297 0.0487% 4943981
    N Rate: 6.7397% 5972234 2.0264% 205736337
    (read2 data omitted)
    Average coverage: 305.47
    Standard deviation: 902.02
    Percent scaffolds with any coverage: 100.00
    Percent of reference bases covered: 99.95

    And I wonder what does this means?
    How to calculate the reads% mapped to contigs?
    Thanks.
    (new to this ,question may be stupid but I have no idea about this.)
  • GenoMax
    Senior Member
    • Feb 2008
    • 7142

    #2
    Just to clarify. You took the original reads and then mapped them back to the assembly done with MegaHit using BBMap?

    Comment

    • GUZUMO
      Junior Member
      • Nov 2016
      • 7

      #3
      yes,is there something wrong?
      thx

      Comment

      • Brian Bushnell
        Super Moderator
        • Jan 2014
        • 2709

        #4
        For paired reads, you can find the total percent mapped by adding the read 1 percent (where it says "mapped: 46.8503%") and read 2 percent, then dividing by 2. The different columns tell you the count/percent of each event. For example, for the "num bases" column and "sub rate" row, the entry indicates there were 199375874 base substitutions (mismatches) compared to the reference, for read 1.

        Comment

        • GUZUMO
          Junior Member
          • Nov 2016
          • 7

          #5
          Thanks,
          So how to understand "Match Rate" and "Error Rate"?
          is there a pdf or something that explain the result of bbmap?

          Comment

          • Brian Bushnell
            Super Moderator
            • Jan 2014
            • 2709

            #6
            Considering the cigar strings from alignment, "Match Rate" is the number of symbols indicating a reference match (=) and error rate is the number indicating substitution, insertion, or deletion (X, I, D).

            BBMap does not have a PDF, but its documentation is included in /bbmap/docs/guides/BBMapGuide.txt

            Comment

            Latest Articles

            Collapse

            • mylaser
              Reply to Advanced Sequencing Platforms Tackle Neuroscience’s Toughest Genomics Problems
              by mylaser
              Kheloyaar: The Complete Guide to Kheloyaar Loginand Kheloyaar ID
              The online gaming industry has transformed the way people enjoy digital entertainment. As technology continues to improve, players are looking for platforms that offer convenience, security, and a seamless user experience. Kheloyaarhas gained attention among users who value an easy-to-use platform, quick account access, and a simple registration process.
              Whether you're exploring Kheloyaar for the first time or want to understand...
              Today, 09:27 PM
            • SEQadmin2
              Advanced Sequencing Platforms Tackle Neuroscience’s Toughest Genomics Problems
              by SEQadmin2



              Genomics studies in neuroscience face a special challenge due to the brain’s complexity and scarcity of samples. Mapping changes in cell type and state using conventional next-generation sequencing methods remains challenging. Advances in technologies like single-cell sequencing, spatial transcriptomics, and long-read sequencing have opened the door to deeper studies of the brain and diseases like Alzheimer’s, amyotrophic lateral sclerosis (ALS), and schizophrenia.
              ...
              Today, 11:10 AM
            • SEQadmin2
              Cancer Drug Resistance: The Lingering Barrier to Rising Survival
              by SEQadmin2



              Cancer survival rates have significantly increased in the last few decades in the United States, reaching a combined 70% 5-year survival rate by 2021. Behind this number, there are years of research to find new therapies, drug targets, and early detection methods. But there is one core challenge that keeps slowing down these advances, and it’s about drug resistance.

              There is no single reason why many patients don’t respond to treatment as expected. Cancer is...
              Yesterday, 05:17 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by SEQadmin2, Today, 10:04 AM
            0 responses
            8 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, Yesterday, 10:08 AM
            0 responses
            7 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 07-07-2026, 11:05 AM
            0 responses
            10 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 07-02-2026, 11:08 AM
            0 responses
            31 views
            0 reactions
            Last Post SEQadmin2  
            Working...