Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • 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.)

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

    Comment


    • #3
      yes,is there something wrong?
      thx

      Comment


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


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

          Comment


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

            • 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