Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,733,7190CG25.6% 93.4 / 10.9 39G63G (GGC→GGG) gumDADP‑L‑glycero‑D‑mannoheptose‑6‑epimerase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (1/28);  new base G (10/0);  total (11/28)
Fisher's exact test for biased strand distribution p-value = 1.73e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTAT  >  REL606/3733648‑3733794
                                                                       |                                                                           
gtgAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAag                                                                         <  1:1350716/76‑1 (MQ=255)
        ggtggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGc                                                                 <  2:1642681/76‑1 (MQ=255)
           ggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGAt                                                              <  2:1258108/76‑1 (MQ=255)
           ggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGAt                                                              <  2:1138507/76‑1 (MQ=255)
           ggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGAt                                                              <  2:1581286/76‑1 (MQ=255)
           ggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGAt                                                              <  1:330966/76‑1 (MQ=255)
            gATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGAAGAGATTGGCGagg                                                             >  2:104995/1‑74 (MQ=255)
              tCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTc                                                           <  2:1155613/76‑1 (MQ=255)
              tCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTc                                                           <  2:351028/76‑1 (MQ=255)
               cTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGGTGTCg                                                          >  1:933793/1‑76 (MQ=255)
               cTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCg                                                          <  1:32414/76‑1 (MQ=255)
                tGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGa                                                         <  1:954384/76‑1 (MQ=255)
                tGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGa                                                         <  1:53678/76‑1 (MQ=255)
                 gAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGAAGAGGTCGGCGATGGCGaa                                                        >  2:1557304/1‑76 (MQ=255)
                  aaTATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTTGGCGCTGTCGGAg                                                       >  2:1049485/1‑76 (MQ=255)
                   atatCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGATCGGCGATGTCGAAGc                                                      >  2:1691616/1‑76 (MQ=255)
                    tatCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCg                                                     <  1:1156786/76‑1 (MQ=255)
                     atCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGGGGGGTCGAAGCGa                                                    >  1:613162/1‑76 (MQ=255)
                       cGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGAtt                                                  <  2:654641/76‑1 (MQ=255)
                        gCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAAAGGTTGGGGATGTCGGAGCGGttt                                                 >  1:92118/1‑76 (MQ=255)
                        gCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGAttt                                                 <  2:1180986/76‑1 (MQ=255)
                         cAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTTGGCGGTGTCGAAGCGGtttt                                                >  1:524624/1‑76 (MQ=255)
                         cAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGGTGTCGAAGCGAtttt                                                >  2:1193824/1‑76 (MQ=255)
                          aGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGGTGTCGGAGGGGttttt                                               >  1:1054788/1‑75 (MQ=255)
                           gACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTcc                                              <  2:961588/76‑1 (MQ=255)
                            aCTATATGGATAAGCAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCa                                             <  2:810332/76‑1 (MQ=255)
                              tataTGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACg                                           <  1:1058525/76‑1 (MQ=255)
                              tataTGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACg                                           <  1:1035478/76‑1 (MQ=255)
                                    gATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGgcg                                     <  1:1141273/76‑1 (MQ=255)
                                     aTAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGGAGAGGTCGGCGGTGTCGGAGCGGTTTTTCACGAAGgcgc                                    >  1:1504700/1‑76 (MQ=255)
                                     aTAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGgcgc                                    <  1:103696/76‑1 (MQ=255)
                                       aaGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGGTGTCGGAGCGGTTTTTCACGAAAgcgccg                                  >  1:1254063/1‑74 (MQ=255)
                                           aaGACTTCCTGATCCAGATTATGGCTGGGGAAGAGGTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGctc                              >  2:1249204/1‑76 (MQ=255)
                                              aCTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTc                           <  2:122959/76‑1 (MQ=255)
                                              aCTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTc                           <  1:1458018/76‑1 (MQ=255)
                                                ttCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTcca                         <  2:1023903/76‑1 (MQ=255)
                                                 tCCTGATCCAGATTATGGCTGGGGGAGGGGTGGGCGGTGTCGACGCGACTTTCCACCCAAGCGGAAGCTCTTcctc                        >  2:322269/1‑74 (MQ=255)
                                                 tCCTGATCCAGATTATGGCTGGGGGAGAGTTCGGCGGTGTGGAAGCGATTTTTCACCAAGGCGGACGGTCTTtcac                        >  2:1264857/1‑76 (MQ=255)
                                                   cTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTccacca                      <  1:131394/76‑1 (MQ=255)
                                                    tGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCcaccac                     <  2:757519/76‑1 (MQ=255)
                                                     gATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCaccacc                    <  2:595851/76‑1 (MQ=255)
                                                          aGATTATGGCTGGGGGAGAGGTCGGCGATGTCGGAGGGAGTTTTCACGAAAGGGCATGCTCTTCCACCACCGAgcg               >  2:806010/1‑74 (MQ=255)
                                                            aTTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTggg             <  2:608214/76‑1 (MQ=255)
                                                             ttATGGCTGGGGAAGAGGTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGCGTgggg            >  2:1585704/1‑75 (MQ=255)
                                                               aTGGCTGGCGAAGAGTTCGGCGNTGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATg          >  2:230535/1‑76 (MQ=255)
                                                               aTGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATg          <  1:1325244/76‑1 (MQ=255)
                                                                 ggctggcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGc        <  2:555513/76‑1 (MQ=255)
                                                                     ggcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGt    <  1:415732/76‑1 (MQ=255)
                                                                       cGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGtat  <  1:735488/76‑1 (MQ=255)
                                                                       cGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGtat  <  2:256079/76‑1 (MQ=255)
                                                                       cGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGtat  <  2:1107501/76‑1 (MQ=255)
                                                                       |                                                                           
GTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTAT  >  REL606/3733648‑3733794

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑