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

TGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTA  >  REL606/3733647‑3733793
                                                                        |                                                                          
tgtgAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGaa                                                                         <  2:907585/76‑1 (MQ=255)
 gtgAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAag                                                                        <  1:468394/76‑1 (MQ=255)
 gtgAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAag                                                                        <  1:1348093/76‑1 (MQ=255)
    aaCCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGagagg                                                                     >  2:567616/1‑75 (MQ=255)
      ccTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTc                                                                   >  1:1907260/1‑76 (MQ=255)
      ccTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTc                                                                   <  2:1901394/76‑1 (MQ=255)
       cTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGTTCg                                                                  >  1:1520562/1‑76 (MQ=255)
         ggtggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGc                                                                <  1:206700/76‑1 (MQ=255)
          gtggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCg                                                               <  1:1204727/76‑1 (MQ=255)
            ggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGAt                                                             <  2:1213878/76‑1 (MQ=255)
            ggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGAt                                                             <  2:1104811/76‑1 (MQ=255)
              aTCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTTGGCGGTgg                                                           >  2:1402760/1‑75 (MQ=255)
               tCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTc                                                          <  2:1266014/76‑1 (MQ=255)
                 tGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTGGGCGGTGTCgg                                                        >  2:934826/1‑75 (MQ=255)
                  gAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGaa                                                       <  2:447118/76‑1 (MQ=255)
                   aaTATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGAGAAGAGTTCGGCGATGTCGAAg                                                      >  2:106262/1‑76 (MQ=255)
                    atatCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGc                                                     <  1:1896698/76‑1 (MQ=255)
                      atCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGAAGAGGTCGGCGGTGTCGAAGCgc                                                   >  1:1360682/1‑75 (MQ=255)
                         gCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGAttt                                                <  2:1270860/76‑1 (MQ=255)
                          cAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGGGGTTGGCGCTGTTGCAGGGActtt                                               >  2:8017/1‑76 (MQ=255)
                          cAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGGTGGCGAAGGGAtttt                                               >  1:1828000/1‑76 (MQ=255)
                            gACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTcc                                             <  1:25482/76‑1 (MQ=255)
                             aCTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGGTGGCGAAGCGGttttttc                                            >  2:434759/1‑73 (MQ=255)
                             aCTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCa                                            <  2:1803492/76‑1 (MQ=255)
                             aCTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCa                                            <  2:1513395/76‑1 (MQ=255)
                             aCTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCa                                            <  1:333968/76‑1 (MQ=255)
                                 tatGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGAAGAGGTTGGCGGTGTCGGAGGGGTTTTTCACGaa                                        >  1:1388149/1‑76 (MQ=255)
                                   tGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAgg                                      <  2:958865/76‑1 (MQ=255)
                                   tGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAgg                                      <  2:958885/76‑1 (MQ=255)
                                    ggATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGCAGAGTTCGGCGGGGTCGGAGGGCTTTTTCACCAACgc                                     >  1:1684374/1‑76 (MQ=255)
                                      aTAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGGTGGCGGAGCGGTTTTCCACGAAggggg                                   >  1:1826414/1‑73 (MQ=255)
                                      aTAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGgcgc                                   <  1:253592/76‑1 (MQ=255)
                                          ggAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGc                               <  1:870006/76‑1 (MQ=255)
                                           gAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGGTGTCGAAGCGATTTTTCACGAAGGCGCATGct                              >  1:670062/1‑76 (MQ=255)
                                            aaGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGGTGTCGAAGCGATTTTTCACGAAGGCGCATGctc                             >  1:663885/1‑76 (MQ=255)
                                             aGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGctct                            <  1:145791/76‑1 (MQ=255)
                                              gACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCtt                           <  2:1048156/76‑1 (MQ=255)
                                              gACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCtt                           <  1:1267270/76‑1 (MQ=255)
                                               aCTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTc                          <  1:1767946/76‑1 (MQ=255)
                                               aCTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTc                          <  2:1879820/76‑1 (MQ=255)
                                                cTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTTCACGAAGGCGCACGCTCTTcc                         >  2:1743013/1‑76 (MQ=255)
                                                 ttCCTGATCCAGATTATGGCTGGGGGCGAGGGTGGCGATGTGGAGACGAGTTTTTCCAAAAGGGCATGGTCTTccc                        >  1:743918/1‑75 (MQ=255)
                                                 ttCCTGATCCAGATTATGGCTGGCGGAGAGGTCGGCGGAGGCCGAGCGAGTTTTCCCGAAAGCGCCAGCTCTTccc                        >  1:1254177/1‑75 (MQ=255)
                                                     tGATCCAGATTATGGCTGGGGGAGAGGTCGGGGGTGTCGAAGCGGTTTTTCACCAAGGCGCAGGCTCTTCccccac                    >  1:1454193/1‑76 (MQ=255)
                                                     tGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCcaccac                    <  2:653153/76‑1 (MQ=255)
                                                     tGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCcaccac                    <  2:658144/76‑1 (MQ=255)
                                                      gATCCAGATTATGGCTGGGGGAGAGGTCGGCGGTGTCGAAGCGGTTTTTCACGAAGGGGCATGCTCTTCCaccccc                   >  2:1910128/1‑76 (MQ=255)
                                                      gATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCaccacc                   <  1:1449275/76‑1 (MQ=255)
                                                          cAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGt               <  1:1435174/76‑1 (MQ=255)
                                                          cAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGt               <  1:14290/76‑1 (MQ=255)
                                                           aGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTg              <  2:343965/76‑1 (MQ=255)
                                                            gATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTgg             <  2:1445650/76‑1 (MQ=255)
                                                             aTTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTggg            <  2:658423/76‑1 (MQ=255)
                                                                      ggcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGt   <  1:947849/76‑1 (MQ=255)
                                                                      ggcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGt   <  1:442773/76‑1 (MQ=255)
                                                                       gcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATACTCTTCCACCACCGAGTGGGATGGCAAGta  >  2:815560/1‑76 (MQ=255)
                                                                        |                                                                          
TGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTA  >  REL606/3733647‑3733793

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