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

GTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGT  >  REL606/3733644‑3733792
                                                                           |                                                                         
gTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTAtggctggc                                                                           <  1:1358897/76‑1 (MQ=255)
     tgAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAaga                                                                      <  2:1479962/76‑1 (MQ=255)
      gAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAagag                                                                     <  1:1572963/76‑1 (MQ=255)
       aaCCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGagagg                                                                    >  1:1687957/1‑75 (MQ=255)
         ccTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTc                                                                  >  1:1100450/1‑76 (MQ=255)
          cTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCg                                                                 <  2:967292/76‑1 (MQ=255)
             gtggGTCCGAATATCGCAGTCTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCg                                                              <  2:254944/76‑1 (MQ=255)
             gtggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCg                                                              <  1:1447171/76‑1 (MQ=255)
              tggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGCAGAGGTCGGCgg                                                             >  1:1067705/1‑75 (MQ=255)
                   cTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCg                                                        <  1:817555/76‑1 (MQ=255)
                    tGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGGAGAGGTCGGGGCTGTCGa                                                       >  2:1556238/1‑76 (MQ=255)
                     gAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTTGGCGGTGTCGaa                                                      >  1:1628003/1‑76 (MQ=255)
                     gAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGAGGTCgca                                                      >  1:324449/1‑74 (MQ=255)
                     gAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGaa                                                      <  2:834701/76‑1 (MQ=255)
                      aaTATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGAAGGGGTCGGGGATGGCGAAg                                                     >  2:1144762/1‑76 (MQ=255)
                         atCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGa                                                  <  2:1414670/76‑1 (MQ=255)
                            gCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGAttt                                               <  2:1377068/76‑1 (MQ=255)
                             cAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGGGGTGTCGAAGCGAtttt                                              >  1:144971/1‑76 (MQ=255)
                              aGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGAAGAGGTTGGGGAGGGCGGAGGGGTTTTc                                             >  1:412944/1‑76 (MQ=255)
                                  tataTGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACg                                         <  1:1169263/76‑1 (MQ=255)
                                   atatGGATAAGGAAGATTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGa                                        <  2:472532/76‑1 (MQ=255)
                                   atatGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTTGGCGGGGTCGGAGCGGTTTTTCACGa                                        >  2:1611328/1‑76 (MQ=255)
                                     atGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTTGGCGGTGTCGGAGGGATTTTCCACGAAg                                      >  1:544340/1‑76 (MQ=255)
                                        gATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGGTGTCGAAGGGATTTTTCACCAAGgcg                                   >  2:1610380/1‑76 (MQ=255)
                                        gATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGgcg                                   <  2:1667802/76‑1 (MQ=255)
                                         aTAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGATGTCGAAGCGGTTTTTCACGAAAgcgc                                  >  2:1022334/1‑76 (MQ=255)
                                               aaGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTTGGCGATGGCCAAGGGGTTTTTCCCGAAAGGGCCTGctc                            >  2:163253/1‑76 (MQ=255)
                                               aaGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGGGGTGTCGGAGGGGTTTTTCACGAAGGGGCATGctc                            >  1:1417787/1‑76 (MQ=255)
                                               aaGACTTCCTGATCCAGAATATGGCTGGGGGAGAGGTCGGCGATGTCGAAGGGGTTTTTCACGAAGGCGGAGGctc                            >  2:792003/1‑76 (MQ=255)
                                               aaGACGTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGctc                            <  2:387654/76‑1 (MQ=255)
                                                    ttCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTcca                       <  2:929458/76‑1 (MQ=255)
                                                    ttCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTcca                       <  2:291913/76‑1 (MQ=255)
                                                     tCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTccac                      <  2:1081506/76‑1 (MQ=255)
                                                         gATCCAGATTATGGCTGGGGGAGAGGTTGGCGGTGTCGGAGCGATTTTTCACGAAGGCGCATGCTCTTcccccccc                  >  2:808718/1‑76 (MQ=255)
                                                         gATCCAGATTATGGCTGGGGGAGAGATCGGGGGTGTCGAAGGGGTTTTCCACGAAGGCGCAGGCTCTTcccccccc                  >  2:1741746/1‑76 (MQ=255)
                                                          aTCCAGATTATGGCTGGGGAAGAGTTCGGCGGTGTCGAAGCGATTTTTCACGAAAGCGCATGCTCTTCCaccaccc                 >  2:1589169/1‑75 (MQ=255)
                                                           tCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGa                <  1:1506314/76‑1 (MQ=255)
                                                             cAGATTATGGCTGGGGGAGAGTTGGGCGGAGTCGAAGCGAGTTTTCACGAAAGCGCATGGTCTTCCCCCCCCCAgg              >  2:1301654/1‑75 (MQ=255)
                                                             cAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGt              <  1:980892/76‑1 (MQ=255)
                                                             cAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGt              <  2:439108/76‑1 (MQ=255)
                                                                  tATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGAt         <  2:1650316/76‑1 (MQ=255)
                                                                   aTGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATg        <  2:228653/76‑1 (MQ=255)
                                                                      gctggcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCa     <  1:1131563/76‑1 (MQ=255)
                                                                      gctggcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCa     <  2:1362892/76‑1 (MQ=255)
                                                                      gctggcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCa     <  2:1562630/76‑1 (MQ=255)
                                                                       ctggcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCaa    <  2:1169408/76‑1 (MQ=255)
                                                                        tggcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAg   >  1:102710/1‑76 (MQ=255)
                                                                         ggcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGt  <  1:936977/76‑1 (MQ=255)
                                                                           |                                                                         
GTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGT  >  REL606/3733644‑3733792

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