Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,756,5790T→G13.6% 1.6 ‑14.7 0.098290ygbJpredicted dehydrogenase, with NAD(P)‑binding Rossmann‑fold domain
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (19/0):  ref base (29/48):  total (50/48)
Fisher's exact test strand distribution p-value = 4.11e-07
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

AGGCGTTGCACAACATCTGAAACCCGGTACGGCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGATGCGC  >  REL606/2756544‑2756614
                                   |                                    
aGGCGTTGCACAACATCTGAAACCCGGTACGGCAGt                                      <  1:2337668/36‑1 (MQ=255)
aGGCGTTGCACAACATCTGAAACCCGGTACGGCAGt                                      <  1:1836274/36‑1 (MQ=255)
aGGCGTTGCACAACATCTGAAACCCGGTACGGCAGt                                      <  1:5430657/36‑1 (MQ=255)
aGGCGTTGCACAACATCTGAAACCCGGTACGGCAGt                                      <  1:1450955/36‑1 (MQ=255)
aGGCGTTGCACAACATCTGAAACCCGGTACGGCAGt                                      <  1:4032666/36‑1 (MQ=255)
aGGCGTTGCACAACATCTGAAAACCGGTACGGCAGt                                      >  1:4303581/1‑36 (MQ=255)
 ggCGTTGCACAACATCTGAAACCCGGTACGGCAGTa                                     <  1:3888751/36‑1 (MQ=255)
 ggCGTTGCACAACATCTGAAACCCGGTACGGCAGTa                                     <  1:3060697/36‑1 (MQ=255)
   cGTTGCACAACATCTGAAACCCGGTACGGCAGTAA.t                                  >  1:3783508/1‑36 (MQ=255)
   cGTTGCACAACATCTGAAACCCGGTACGGCAGTAA.t                                  <  1:1835460/36‑1 (MQ=255)
   cGTTGCACAACATCTGAAACCCGGTACGGCAGTAA.t                                  <  1:3478341/36‑1 (MQ=255)
    gTTGCACAACATCTGAAACCCGGTACGGCAGTAA.Tg                                 <  1:316548/36‑1 (MQ=255)
    gTTGCACAACATCTGAAACCCGGTACGGCAGTAA.Tg                                 >  1:1369178/1‑36 (MQ=255)
    gTTGCACAACATCTGAAACCCGGTACGGCAGGAA.Tg                                 >  1:4695324/1‑36 (MQ=255)
    gTTGCACAACATCTGAAACACGGTACGGCAGGAA.Tg                                 >  1:4633967/1‑36 (MQ=37)
     ttGCACAACATCTGAAACCCGGTACGGCAGTAA.Tgg                                >  1:2144382/1‑36 (MQ=255)
     ttGCACAACATCTGAAACCCGGTACGGCAGGAA.Tgg                                >  1:2872440/1‑36 (MQ=255)
     ttGCACAACATCTGAAACCCGGTACGGCAGGAA.Tgg                                >  1:4739031/1‑36 (MQ=255)
       gCACAACATCTGAAACCCGGTACGGCAGTTA.TGGtt                              >  1:1320367/1‑36 (MQ=255)
       gCACAACATCTGAAACCCGGTACGGCAGGAA.Tgggt                              >  1:3785889/1‑34 (MQ=37)
        cacaACATCTGAAACCCGGTACGGCAGTAA.TTGttt                             >  1:1610466/1‑36 (MQ=255)
        cacaACATCTGAAACCCGGTACGGCAGTAA.TGGttt                             <  1:2185160/36‑1 (MQ=255)
        cacaACATCTGAAACCCGGTACGGCAGTAA.TGGttt                             <  1:3576221/36‑1 (MQ=255)
         acaacaTCTGAAACCCGGTACGGCAGTCA.TGGtttt                            >  1:210411/1‑35 (MQ=255)
         acaacaTCTGAAACCCGGTACGGCAGTAA.TGGTttc                            <  1:2966140/36‑1 (MQ=255)
         acaacaTCTGAAACCCGGTACGGCAGTAA.TGGTttc                            >  1:1619198/1‑36 (MQ=255)
         acaacaTCTGAAACCCGGTACGGCAGGGA.TGGTttc                            >  1:4060988/1‑36 (MQ=37)
         acaacaTCTGAAACCCGGTACGGCAGGCA.TTGTttc                            >  1:4870364/1‑36 (MQ=25)
         acaacaTCTGAAACCCGGTACGGCAGAGA.TGGtttt                            >  1:2307706/1‑35 (MQ=37)
         acaacaTCTGAAACCCGGGACGACAGGAA.TGGtttg                            >  1:5006609/1‑35 (MQ=255)
          caacaTCTGAAACCCGGTACGGCAGTAA.TGGTttct                           <  1:2583155/36‑1 (MQ=255)
          caacaTCTGAAACCCGGTACGGCAGTAA.TGGTttct                           <  1:5480769/36‑1 (MQ=255)
           aacaTCTGAACCCCGGTACGGCAGTAA.TGGTttctt                          <  1:3374522/36‑1 (MQ=255)
           aacaTCTGAAACCCGGTACGGCAGTTA.TTGTttttt                          >  1:2491616/1‑36 (MQ=25)
           aacaTCTGAAACCCGGTACGGCAGTAA.TGGTttctt                          <  1:1567986/36‑1 (MQ=255)
           aacaTCTGAAACCCGGTACGGCAGTAA.TGGTttctt                          <  1:291067/36‑1 (MQ=255)
           aacaTCTGAAACCCGGTACGGCAGCAA.TTGTttttt                          >  1:928114/1‑36 (MQ=25)
            acaTCTGAAGCCCGTTACGGCCGTAA.TGGTttcttc                         <  1:1573884/36‑1 (MQ=25)
            acaTCTGAAACCCGGTTCGGCAGTAA.TGGTttcttc                         <  1:5442449/36‑1 (MQ=255)
            acaTCTGAAACCCGGTACGGCAGTTA.CGGTttcctt                         >  1:4460162/1‑33 (MQ=25)
            acaTCTGAAACCCGGTACGGCAGTAA.TGGTttcttc                         <  1:5883705/36‑1 (MQ=255)
            acaTCTGAAACCCGGTACGGCAGGAA.TGGGttcttc                         >  1:1332120/1‑36 (MQ=37)
             caTCTGAACCCCGGTACGGCAGTAA.TGGTTTCTTcc                        <  1:3575310/36‑1 (MQ=255)
             caTCTGAAACCCGGTACGGCAGTAA.TGGTTTCTTcc                        <  1:1432108/36‑1 (MQ=255)
             caTCTGAAACCCGGTACGGCAGTAA.TGGTTTCTTcc                        <  1:3304065/36‑1 (MQ=255)
             caTCTGAAACCCGGTACGGCAGTAA.TGGTTTCTTcc                        <  1:5859887/36‑1 (MQ=255)
             caTCTGAAACCCGGTACGGCAGTAA.TGGGttctttc                        >  1:4199695/1‑34 (MQ=37)
              aTCTGAAGCCCGGTACGGCAGTAA.TGGTTTCTTCCa                       <  1:508648/36‑1 (MQ=255)
              aTCTGAAACCCGGTACGGCAGGAA.TTGTttttttcc                       >  1:2861463/1‑33 (MQ=16)
              aTCTGAAACCCGGTACGGCAGGAA.TGGTttttttcc                       >  1:4045622/1‑33 (MQ=25)
               tCTGAAACCCGGTACGGCAGTAA.TGGTTTTTTcccc                      >  1:5247699/1‑34 (MQ=38)
               tCTGAAACCCGGTACGGCAGTAA.TGGTTTCTTCCAc                      >  1:1427369/1‑36 (MQ=255)
               tCTGAAACCCGGTACGGCAGGAAGTGGTTTTTTccc                       >  1:1297591/1‑35 (MQ=16)
               tCTGAAACCCGGTACGGCAGGAA.AGGTtttttctcc                      >  1:4791044/1‑33 (MQ=16)
                cTGAAACCCGGTACGGCAGTAA.TGGTTTCTTCCACt                     >  1:4968746/1‑36 (MQ=255)
                cTGAAACCCGGTACGGCAGTAA.TGGTTTCTTCCACt                     <  1:1897044/36‑1 (MQ=255)
                 tGAAACCCGGTACGGCAGTAA.TGGTTTTTTCCCCTa                    >  1:749265/1‑36 (MQ=37)
                 tGAAACCCGGTACGGCAGTAA.TGGTTTCTTCCACTa                    >  1:5869778/1‑36 (MQ=255)
                 tGAAACCCGGTACGGCAGGAA.TGGGTTTTTTCACtt                    >  1:1306062/1‑35 (MQ=16)
                  gAAACCCGGTACGGCAGTAA.TGGTTTCTTCCACTAt                   <  1:4079393/36‑1 (MQ=255)
                  gAAACCCGGTACGGCAGTAA.TGGTTTCTTCCACTAt                   <  1:2513019/36‑1 (MQ=255)
                  gAAACCCGGTACGGCAGTAA.TGGTTTCTTCCACTAt                   <  1:3043542/36‑1 (MQ=255)
                   caaCCAGGTCCGGCAGTAA.TGGTTTCTTCCACTATc                  <  1:5374225/35‑1 (MQ=37)
                   aaaCCCGGTACGGCAGTAA.TGGTTTCTTCCACTATc                  <  1:53801/36‑1 (MQ=255)
                    aaCCCGGTACGGCAGTAA.TGGTTTTTTTCACTAttc                 >  1:670746/1‑34 (MQ=37)
                    aaCCCGGTACGGCAGTAA.TGGTTTCTTCCCCTATCg                 >  1:4541233/1‑36 (MQ=255)
                    aaCCCGGTACGGCAGTAA.TGGTTTCTTCCACTATCg                 <  1:2164838/36‑1 (MQ=255)
                    aaCCCGGTACGGCAGGAA.AGGGTTTTTTCACTAttt                 >  1:3092859/1‑34 (MQ=12)
                     aCCCGGTACGGCAGTAA.TGGTTTCTTCCTCTAtttc                >  1:3108847/1‑33 (MQ=37)
                     aCCCGGTACGGCAGTAA.TGGTTTCTTCCACTATCGc                <  1:3912907/36‑1 (MQ=255)
                     aCCCGGTACGGCAGTAA.TGGTTTCTTCCACTATCGc                <  1:2555216/36‑1 (MQ=255)
                     aCCCGGTACGGCAGTAA.TGGGTTCTTACACTATCgg                >  1:4130291/1‑35 (MQ=37)
                      cccGGTACGGCAGTAA.TGGTTTCTTCCACTATCGCt               <  1:1894954/36‑1 (MQ=255)
                      cccGGTACGGCAGTAA.TGGTTTCTTCCACTATCGCt               <  1:344390/36‑1 (MQ=255)
                      cccGGTACGGCAGTAA.TGGTTTCTTCCACTATCGCt               <  1:3067282/36‑1 (MQ=255)
                      cccGGTACGGCAGGTA.TGGGTTTTTTCACCATCCCt               >  1:4211464/1‑36 (MQ=0)
                       ccGGTACGTCAGTAA.TGGTTTCTTCCACTATCGCTa              <  1:3719914/36‑1 (MQ=255)
                         ggTACGGCAGTTA.TGGTTTTTTTCCCTATCGCTTGt            >  1:5539446/1‑36 (MQ=16)
                          gTACGGCAGTAA.TGGTTTCTTCCACTATCGCTAGTg           <  1:3650622/36‑1 (MQ=255)
                            aCGGCAGTAA.TGGTTTCTTCCACTATCTCTAGTGCt         >  1:1861362/1‑36 (MQ=255)
                            aCGGCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCt         <  1:681700/36‑1 (MQ=255)
                             cGGCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTg        <  1:4451598/36‑1 (MQ=255)
                             cGGCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTg        <  1:4016012/36‑1 (MQ=255)
                             cGGCAGGAA.TGGTTTCTTCCACTATCGCTAGTGCtc        >  1:5365123/1‑35 (MQ=255)
                              gtcAGTAA.TGGTTTCTTCCACCATCGCTAGTGCTGa       <  1:5419059/34‑1 (MQ=25)
                              ggCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCtta       >  1:4246322/1‑34 (MQ=255)
                              ggCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGa       <  1:5084423/36‑1 (MQ=255)
                              ggCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGa       >  1:1151245/1‑36 (MQ=255)
                               gCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGAt      <  1:79208/36‑1 (MQ=255)
                               gCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGAt      >  1:5905124/1‑36 (MQ=255)
                               gCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGAt      >  1:3346845/1‑36 (MQ=255)
                                cAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGATg     <  1:598804/36‑1 (MQ=255)
                                cAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGATg     <  1:2741762/36‑1 (MQ=255)
                                cAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGATg     <  1:3965434/36‑1 (MQ=255)
                                 aGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGAttc    >  1:1087336/1‑34 (MQ=255)
                                 aGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGATgg    >  1:4741808/1‑35 (MQ=255)
                                 aGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGATgc    <  1:4590110/36‑1 (MQ=255)
                                   tAA.TGGTTTCTTCCACTCTCGCTAGTGCTGATgcgc  <  1:246287/36‑1 (MQ=255)
                                   tAA.TGGTTTCTTCCACTATCGCTAGTGCTGATgcgc  <  1:2637939/36‑1 (MQ=255)
                                   tAA.TGGTTTCTTCCACTATCGCTAGTGCTGATgcgc  <  1:5647171/36‑1 (MQ=255)
                                   |                                    
AGGCGTTGCACAACATCTGAAACCCGGTACGGCAGTAA.TGGTTTCTTCCACTATCGCTAGTGCTGATGCGC  >  REL606/2756544‑2756614

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 14 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG