Uploaded image for project: 'Barna Package'
  1. Barna Package
  2. BARNA-384

Error rate 100x higher than it should be with 76-bp error model?

    XMLWordPrintable

Details

    • Bug
    • Status: Open
    • Critical
    • Resolution: Unresolved
    • Simulator 1.2.1 (API 1.22)
    • None
    • Simulator
    • None

    Description

      I'm getting that Flux Simulator's substitution rate across all read bases is about 2.5%, far higher than the Phred score plot at http://sammeth.net/confluence/display/SIM/4.5.4+-+Error+Models suggests. My strong suspicion is that somewhere critical in the code, a factor of 100 is missing because of some confusion between working with percentages and not. I already see one place where this confusion occurs: http://sammeth.net/fisheye/browse/barna/barna.simulator/src/main/java/barna/flux/simulator/Sequencer.java, line 301: the variable avg is a proportion of bases that are mutated, not a percentage as claimed by the log.

      I have confirmed for a couple of simulations that the mean error rate computed from quality scores agrees with the substitution rate.

      Attachments

        Activity

          People

            thasso Thasso Griebel
            verve Abhinav Nellore
            Votes:
            0 Vote for this issue
            Watchers:
            1 Start watching this issue

            Dates

              Created:
              Updated:

              Time Tracking

                Estimated:
                Original Estimate - 1 hour
                1h
                Remaining:
                Remaining Estimate - 1 hour
                1h
                Logged:
                Time Spent - Not Specified
                Not Specified