Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • BBRepair quality score changes

    After running BBRepair the log file states: "Changed from ASCII-33 to ASCII-64 on input quality 97 for base C while prescanning." It does this for base scores 66, 95, and 98. I would like to run the files through quality trimming, but am now unsure how the quality scores are encoded (sanger, illumina, etc). Prior to BBRepair, some files had sanger scores, and other had illumina. Any help would be greatly appreciated.

  • #2
    If the encoding was detected correctly, or if you add the flag "changequality=f", the output will retain the same encoding as the input. It's always best to tell it what the encoding is, though, because it is not possible to detect the quality score encoding with 100% accuracy. You can specify encoding with these flags:

    qin=64 qout=33


    That will force the input to be treated as ASCII-64 (old Illumina) and the output will be ASCII-33 (Sanger). It's normally easiest to work with data when all sets have the same encoding...

    Comment


    • #3
      Thank you for your quick response. Even though the log file says the quality encoding changes, when I check the file's quality encoding with testformat.sh in bbmap (from your comment here https://www.biostars.org/p/63225/), the output states that the file is in its original encoding. I was hoping you could clarify whether repair.sh permanently or temporarily changes quality encoding.

      Thank you!!
      Last edited by gstone; 02-27-2017, 08:13 PM.

      Comment


      • #4
        Sorry, I guess the message is unclear. When it says this:

        "Changed from ASCII-33 to ASCII-64 on input quality 97 for base C while prescanning."

        That means it started with the assumption that your data is in ASCII-33, but examined the file, and decided the input was ASCII-64. So, it processed the file assuming the data was ASCII-64. In that case, it does not change the data at all - what changes is the assumption about how to interpret the data. In other words, your data is unchanged.

        Comment


        • #5
          Great, thank you so much for the clarification!

          Comment


          • #6
            If your data is in phred+64 format (only possible if it is pretty old) then remember to use qin=64 option with other bbtools.

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