Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606232,1090G→T10.8% 1.9 ‑5.5 0.07019dkgB2,5‑diketo‑D‑gluconate reductase B
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (10/0):  ref base (31/29):  total (41/29)
Fisher's exact test strand distribution p-value = 3.91e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TTAAAAGAGGCATATTATGGCTATCCCTGCATTTGGTTTAGGTACTTTCCGTCTGAAAGACGACGTTGTT  >  REL606/232075‑232144
                                  |                                   
ttAAAAGAGGCATATTATGGCTATCCCGTCCTTTgg                                    >  1:6664149/1‑36 (MQ=25)
 tAAAAGAGGCATATTATGGCTATCCCTGCATTTTGt                                   >  1:5487917/1‑36 (MQ=255)
 tAAAAGAGGCATATTATGGCTATCCCTGCATTTGGt                                   <  1:5040447/36‑1 (MQ=255)
 tAAAAGAGGCATATTATGGCTATCCACGCATTTTGt                                   >  1:6984533/1‑36 (MQ=255)
  aaaaGAGGCATATTATGGCTATCCCTGCATTTGGtt                                  <  1:79775/36‑1 (MQ=255)
  aaaaGAGGCATATTATGGCTATCCCTGCATTTGGtt                                  <  1:2998777/36‑1 (MQ=255)
  aaaaGAGGCATATTATGGCTATCCCTGCATTTGGtt                                  <  1:3308636/36‑1 (MQ=255)
  aaaaGAGGCATATTATGGCTATCCCTGCATTTGGtt                                  <  1:1740484/36‑1 (MQ=255)
    aaGAGGCATATTATGGCTATCCCTTCATTTTGTTTa                                >  1:4697326/1‑36 (MQ=37)
    aaGAGGCATATTATGGCTATCCCTGCATTTGGTTTa                                <  1:1545941/36‑1 (MQ=255)
    aaGAGGCATATTATGGCTATCCCTGCATTTGGTTTa                                <  1:5928196/36‑1 (MQ=255)
    aaGAGGCATATTATGGCTATCCCTGCATTTGGTTTa                                <  1:1391443/36‑1 (MQ=255)
    aaGAGGCATATTATGGCTATCCCTGCATTTGGTTTa                                <  1:4490191/36‑1 (MQ=255)
     agagGCATATTATGGCTATTCCTGCATTTTGttttg                               >  1:5199208/1‑34 (MQ=37)
     agagGCATATTATGGCTATCCCTGCATTTGGTTTAg                               <  1:4504513/36‑1 (MQ=255)
      gagGCATATTATGGCTATCCCTGCATTTGGTTTAgg                              >  1:682304/1‑36 (MQ=255)
       agGCATATTATGGCTATCCCTTAATTTGGTTTAGGt                             >  1:7186045/1‑36 (MQ=37)
        ggCATATTATGGCTATCCCTGCTTTTTGTTTAGGTa                            >  1:1232585/1‑36 (MQ=37)
        ggCATATTATGGCTATCCCTGCATTTGGTTTATGTa                            <  1:3222899/36‑1 (MQ=255)
         gCATATTATGGCTAACCCTGCATTTTGTTTTTGTAc                           >  1:6150836/1‑36 (MQ=19)
          cATATTATGGCTATCCCTGCATTTGGTTTAGGTACt                          <  1:41581/36‑1 (MQ=255)
           atatTATGGCTATCCCTGCCTTTTTTTTAGGTACtt                         >  1:5036428/1‑36 (MQ=25)
           atatTATGGCTATCCCTGCATTTTGTTTAGGTACtt                         >  1:4789615/1‑36 (MQ=255)
           atatTATGGCTATCCCTGCATTTGGTTTAGGTACtt                         >  1:442827/1‑36 (MQ=255)
           atatTATGGCTATCCCTGCATTTGGTTTAGGTACtt                         >  1:4865458/1‑36 (MQ=255)
           atatTATGGCTATCCCTGCATTTGGTTTAGGTACtt                         <  1:959086/36‑1 (MQ=255)
            tattatGGCTATCCCTTCATTTTGTTTATGTACttt                        >  1:1244383/1‑36 (MQ=25)
            tattatGGCTATCCCTGCATTTGGTTTAGGTACttt                        <  1:7090789/36‑1 (MQ=255)
            tattatGGCTATCCCTGCATTTGGTTTAGGTACttt                        <  1:5279031/36‑1 (MQ=255)
            tattatGGCTATCCCTGCATTTGGTTTAGGTACttt                        >  1:6484467/1‑36 (MQ=255)
              ttatGGCTATCCCTGCATTTGGTTTAGGTACTTTcc                      >  1:2341258/1‑36 (MQ=255)
              ttatGGCTATCCCTGCATTTGGTTTAGGTACTTTcc                      <  1:6551098/36‑1 (MQ=255)
              ttatGGCTATCCCTGCATTTGGTTTAGGTACTTTcc                      <  1:5548949/36‑1 (MQ=255)
               tatGGCTATCCCTGCCTTTTGTTTAGGTACTTTCCg                     >  1:1402253/1‑36 (MQ=37)
               tatGGCTATCCCTGCATTTGGTTTAGGTATTTTcct                     >  1:4056367/1‑35 (MQ=255)
                atGGCTATCCCTGCATTTGGTTTAGGTACTTTCCGt                    >  1:5209579/1‑36 (MQ=255)
                 tGGCTATCCCTGCATTTGGTTTAGGTACTTTCCGTc                   >  1:7149213/1‑36 (MQ=255)
                 tGGCTATCCCTGCATTTGGTTTAGGTACTTTCCGTc                   >  1:4064417/1‑36 (MQ=255)
                  ggCTATCCCTGCATTTGGTTTAGGTACTTTCCGTCt                  <  1:2690038/36‑1 (MQ=255)
                  ggCTATCCCTGCATTTGGTTTAGGTACTTTCCGTCt                  <  1:1441914/36‑1 (MQ=255)
                  ggCTATCCCTGCATTTGGTTTAGGTACTTTCCGTCt                  >  1:1141553/1‑36 (MQ=255)
                   gCTATCCCTGCATTTGGTTTAGGTACTTTCCGTCTg                 <  1:1959209/36‑1 (MQ=255)
                   gCTATCCCTGCATTTGGTTTAGGTACTTTCCGTCTg                 <  1:938067/36‑1 (MQ=255)
                   gCTATCCCTGCATTTGGTTTAGGTACTTTCCGTCTg                 >  1:2557739/1‑36 (MQ=255)
                    cTATCCCTGCATTTGGTTTAGGTACTTTCCGTCTGa                >  1:7857269/1‑36 (MQ=255)
                    cTATCCCTGCATTTGGTTTAGGTACTTTCCGTCTGa                >  1:2656558/1‑36 (MQ=255)
                     tATCCCTGCATTTGGTTTAGGTACTTTCCGTCTGaa               >  1:7616677/1‑36 (MQ=255)
                      aTCCCTGCATTTGGTTTAGGTACTTTCCGTCTGaaa              >  1:5524978/1‑36 (MQ=255)
                      aTCCCTGCATTTGGTTTAGGTACTTTCCGTCTGaaa              >  1:3947992/1‑36 (MQ=255)
                      aTCCCTGCATTTGGTTTAGGTACTTTCCGTCTGaaa              >  1:3074977/1‑36 (MQ=255)
                       tCCCTGCATTTGGTTTAGGTACTTTCCGTCTGAAAg             <  1:854667/36‑1 (MQ=255)
                       tCCCTGCATTTGGTTTAGGTACTTTCCGTCTGAAAg             >  1:274286/1‑36 (MQ=255)
                        cccTGCATTTGGTTTAGGTACTTTCCGTCTGAAAga            >  1:2206258/1‑36 (MQ=255)
                         ccTGCATTTGGTTTAGGTACTTTCCGTCTGAAAgac           >  1:3758374/1‑36 (MQ=255)
                         ccTGCATTTGGTTTAGGTACTTTCCGTCTGAAAgac           <  1:6305565/36‑1 (MQ=255)
                          cTGCATTTGGTTTAGGTACTTTCCGTCTGAAAgacg          >  1:3939468/1‑36 (MQ=255)
                           tGCATTTGGTTTAGGTACTTTCCGTCTGAAAgacga         >  1:258332/1‑36 (MQ=255)
                           tGCATTTGGTTTAGGTACTTTCCGTCTGAAAgacga         >  1:7551706/1‑36 (MQ=255)
                            gCATTTGGTTTAGGTACTTTCCGTCTGAAAgacgac        <  1:2713444/36‑1 (MQ=255)
                            gCATTTGGTTTAGGTACTTTCCGTCTGAAAgacgac        >  1:531046/1‑36 (MQ=255)
                            gCATTTGGTTTAGGTACTTTCCGTCTGAAAgacgac        <  1:6145709/36‑1 (MQ=255)
                              aTTTGGTTTAGGTACTTTCCGTCTGAAAGACGACgt      >  1:7338176/1‑36 (MQ=255)
                              aTTTGGTTTAGGTACTTTCCGTCTGAAAGACGACgt      <  1:1190038/36‑1 (MQ=255)
                              aTTTGGTTTAGGTACTTTCCGTCTGAAAGACGACgt      <  1:7090711/36‑1 (MQ=255)
                              aTTTGGTTTAGGTACTTTCCGTCTGAAAGACGACgt      >  1:6476443/1‑36 (MQ=255)
                               tttGGTTTAGGTACTTTCCGTCTGAAAGACGACgtt     >  1:1950759/1‑36 (MQ=255)
                               tttGGTTTAGGTACTTTCCGTCTGAAAGACGACgtt     <  1:3468858/36‑1 (MQ=255)
                               tttGGTTTAGGTACTTTCCGTCTGAAAGACGACgtt     >  1:1168722/1‑36 (MQ=255)
                               tttGGTTTAGGTACTTTCCGTCTGAAAGACGACgtt     <  1:2459159/36‑1 (MQ=255)
                                 tGGTTTAGGTACTTTCCGTCTGAAAGACGACgttgt   >  1:3823013/1‑36 (MQ=255)
                                 tGGTTTAGGTACTTTCCGTCTGAAAGACGACgttgt   >  1:780935/1‑36 (MQ=255)
                                  ggTTTAGGTACTTTCCGTCTGAAAGACGACgttgtt  >  1:7665328/1‑36 (MQ=255)
                                  ggTTTAGGTACTTTCCGTCTGAAAGACGACgttgtt  <  1:3813305/36‑1 (MQ=255)
                                  |                                   
TTAAAAGAGGCATATTATGGCTATCCCTGCATTTGGTTTAGGTACTTTCCGTCTGAAAGACGACGTTGTT  >  REL606/232075‑232144

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG