Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,127,318 A→G 18.0% E25G (GAA→GGA)  frwD → predicted enzyme IIB component of PTS

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,127,3180A→G18.0% 2.2 ‑4.1 0.04074frwDpredicted enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  new base (8/0):  ref base (17/15):  total (25/15)
Fisher's exact test strand distribution p-value = 1.62e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCA  >  REL606/4127284‑4127353
                                  |                                   
gtggcgCATACTTATATGGCGGCGGAACGGCTGGaa                                    >  1:2940717/1‑36 (MQ=255)
  ggcgCATACTTATATGGCGGCGGAACGGCTGGaaaa                                  <  1:2313144/36‑1 (MQ=255)
  ggcgCATACTTATATGGCGGCGGAACGGCTGGaaaa                                  <  1:1602041/36‑1 (MQ=255)
   gcgcATACTTATATGGCGGCGGAACGGCTGGgaacg                                 >  1:6954312/1‑34 (MQ=37)
   gcgcATACTTATATGGCGGCGGAACGGCTGGGAAAg                                 >  1:4721435/1‑36 (MQ=255)
    cgcATACTTATATGGCGGCGGAACGGCTGGGAAAGt                                >  1:6881173/1‑36 (MQ=255)
    cgcATACTTATATGGCGGCGGAACGGCTGGAAAAGt                                <  1:353404/36‑1 (MQ=255)
     gcATACTTATATGGCGGCGGAACGGCTGGGAAAGtt                               >  1:2444817/1‑36 (MQ=255)
      cATACTTATATGGCGGCGGAACGGCTGGGAAAGttt                              >  1:7136638/1‑35 (MQ=255)
      cATACTTATATGGCGGCGGAACGGCTGGAAAAGTtg                              <  1:5945787/36‑1 (MQ=255)
        tACTTATATGGCGNCGGAACGGCTGGGAAAGTtgtg                            >  1:4344086/1‑36 (MQ=38)
        tACTTATATGGCGGCGGAACGGCTGGAAAAGTtttg                            >  1:5796074/1‑36 (MQ=255)
          cTTATATGGCGGCGGAACGGCTGGGAAAGTTGTGcc                          >  1:3106994/1‑36 (MQ=255)
             atatGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGt                       <  1:4207195/36‑1 (MQ=255)
             atatGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGt                       >  1:3428222/1‑36 (MQ=255)
             atatGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGt                       <  1:156757/36‑1 (MQ=255)
              tatGGCGGCGGAACGGCTGGGAAAGGTTTGGCTGtt                      >  1:4571645/1‑36 (MQ=16)
                tGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTag                    <  1:4323944/36‑1 (MQ=255)
                 ggcggcGGAACGGCTGGAAAAGTTGTGCCTGTTcga                   >  1:5524625/1‑36 (MQ=255)
                 ggcggcGGAACGGCTGGAAAAGTTGTGCCTGTTaga                   <  1:7362477/36‑1 (MQ=255)
                 ggcggcGGAACGGCTGGAAAAGTTGTGCCTGTTaga                   >  1:4968951/1‑36 (MQ=255)
                 ggcggcGGAACGGATGGAAAAGTTGTGACTGTTaga                   >  1:3211659/1‑36 (MQ=38)
                  gcggcgGAACGGCTGGAAAAGTTGTGCCTGTTagag                  >  1:7469114/1‑36 (MQ=255)
                   cggcggAACGGCTGGAAAAGTTGTGCCTGTTAgaga                 <  1:3605553/36‑1 (MQ=255)
                   cggcggAACGGCTGGAAAAGTTGTGCCTGTTAgaga                 <  1:2551932/36‑1 (MQ=255)
                    ggcggAACGGCTGGAAAAGTTGTGCCTGTTagagta                >  1:6396075/1‑34 (MQ=255)
                    ggcggAACGGCTGGAAAAGTTGTGCCTGTTAGAGaa                >  1:6630866/1‑36 (MQ=255)
                    ggcggAACGGCTGGAAAAGTTGTGCCTGTTAGAGaa                >  1:437837/1‑36 (MQ=255)
                     gcggAACGGCTGGGAAACTTATGCCTGTTATAGAAg               >  1:6644028/1‑36 (MQ=19)
                     gcggAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAg               <  1:1308313/36‑1 (MQ=255)
                       ggAACGGCTGGAAAAGTTGTGCGTGTTAGAGAAggg             >  1:1003409/1‑34 (MQ=38)
                       ggAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTg             >  1:2374992/1‑36 (MQ=255)
                        gAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTgg            >  1:7330664/1‑36 (MQ=255)
                         aaCGGCTGGAAAAGTTTTGCCTTTTAGAGGAGTggg           >  1:2929080/1‑36 (MQ=25)
                          aCGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTgggg          <  1:3380873/36‑1 (MQ=255)
                          aCGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTgggg          >  1:2560572/1‑36 (MQ=255)
                           cGGCTGGAAAAGTTGTGCCTGTTGGAGAAGTGGGGa         <  1:7198146/36‑1 (MQ=255)
                            ggCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAg        >  1:2866425/1‑36 (MQ=255)
                             gCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGt       <  1:4852229/36‑1 (MQ=255)
                              cTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTc      <  1:1284249/36‑1 (MQ=255)
                               tGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCa     >  1:5937951/1‑36 (MQ=255)
                                ggAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAg    >  1:438750/1‑36 (MQ=255)
                                 gAAAAGTTGTGCCTGTTAGAGAATTGGGGAGTCAGc   <  1:166239/36‑1 (MQ=255)
                                 gAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGc   <  1:6787841/36‑1 (MQ=255)
                                  aaaaGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCa  <  1:736582/36‑1 (MQ=255)
                                  aaaaGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCa  >  1:4606935/1‑36 (MQ=255)
                                  |                                   
GTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCA  >  REL606/4127284‑4127353

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 17 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG