Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,972,6470GA23.5% 117.3 / 35.3 68M42I (ATG→ATA) rffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base G (28/24);  new base A (5/11);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.05e-01

TGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCA  >  NC_000913/3972548‑3972746
                                                                                                   |                                                                                                   
tgCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATg                                                                                                     >  1:765485/1‑100 (MQ=255)
   cGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCg                                                                                                  >  2:249087/1‑100 (MQ=255)
     ggTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCt                                                                                                >  2:535139/1‑100 (MQ=255)
      gTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTg                                                                                               <  2:1134444/100‑1 (MQ=255)
       tttATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTgg                                                                                              <  2:987465/100‑1 (MQ=255)
          aTTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCAc                                                                                           >  2:92730/1‑100 (MQ=255)
          aTTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCAc                                                                                           >  1:524312/1‑100 (MQ=255)
            ttgCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGGCGCGGTGGCAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCAc                                                                                           >  2:519451/3‑98 (MQ=255)
             ggCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCgg                                                                                        <  2:717934/100‑1 (MQ=255)
                  ggCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTcgcg                                                                                   <  1:686650/100‑1 (MQ=255)
                  ggCGCGGGTGCGTTATATCATCAACGAGACGGGCGCCGCGGTGGTAGTGGTCGATAAGCTGTCCTACGCCGGAAACCTGATGTCGCTGGCACCGGTcg                                                                                     <  2:312342/98‑1 (MQ=255)
                        ggTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagc                                                                             <  2:268964/100‑1 (MQ=255)
                        ggTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagc                                                                             >  1:1029871/1‑100 (MQ=255)
                          tGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATATCGCTGGCACCGGTCGCGCAAagcga                                                                           >  2:163550/1‑100 (MQ=255)
                           gCGTTATATCATCAACGAGACGAGCGACGCGGTAGTAGTGTTCGATAAGCTGACCGCCGCCCAAAACCTCATATCGCTCGCACCGCTCGCGCAAag                                                                              <  2:923855/96‑1 (MQ=255)
                           gCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgag                                                                          <  2:909814/100‑1 (MQ=255)
                            cGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgagc                                                                         <  1:214076/100‑1 (MQ=255)
                            cGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgagc                                                                         >  1:438359/1‑100 (MQ=255)
                                 tatCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCttt                                                                    <  1:1054871/100‑1 (MQ=255)
                                 tatCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCttt                                                                    >  2:622694/1‑100 (MQ=255)
                                  atcatcAATGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCCGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTg                                                                   <  2:910435/100‑1 (MQ=255)
                                  atcatcAACTAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTg                                                                   <  2:423735/100‑1 (MQ=255)
                                   tcatcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGc                                                                  <  2:1059431/100‑1 (MQ=255)
                                    catcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGcc                                                                 <  1:455015/100‑1 (MQ=255)
                                    catcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATATCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGcc                                                                 <  1:833209/100‑1 (MQ=255)
                                     atcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCt                                                                >  2:747009/1‑100 (MQ=255)
                                     atcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATATCGCTGGCACCGGTCGCGCAAAGCGCTCGCTTTGcgc                                                                <  1:1047077/100‑3 (MQ=255)
                                     atcaACGAAACGAGCGACGCGGTGGTAGAGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGc                                                                  >  1:1011291/1‑98 (MQ=255)
                                       caACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATATCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCttt                                                              <  1:544683/100‑1 (MQ=255)
                                        aacgaaacgaGCGCTGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGGAACCTGATATCGCTGGCACCGCTCGCGCAAAGCGAGCGATTTGCCtt                                                               >  2:757460/1‑98 (MQ=255)
                                        aacgaaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGTGCGCTTTGCCTTTg                                                             <  2:1194852/100‑1 (MQ=255)
                                        aacgaaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTg                                                             >  1:838856/1‑100 (MQ=255)
                                         acgaaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTga                                                            <  1:304471/100‑1 (MQ=255)
                                             aacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGaaa                                                        >  1:29765/1‑100 (MQ=255)
                                               cgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGt                                                      >  1:1247779/1‑100 (MQ=255)
                                                gaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGtt                                                     >  2:940237/1‑100 (MQ=255)
                                                    gACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGata                                                 >  1:679704/1‑100 (MQ=255)
                                                      cgcgGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATc                                               >  1:513819/1‑100 (MQ=255)
                                                        cgGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTg                                             <  1:92732/100‑1 (MQ=255)
                                                        cgGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATATCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTg                                             <  2:1096120/100‑1 (MQ=255)
                                                                gTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCggg                                     >  1:477028/1‑100 (MQ=255)
                                                                gTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCggg                                     <  2:1120817/100‑1 (MQ=255)
                                                                gTGGTCGATAAGCTGACCTACGCCGGAAACCTGATATCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCggg                                     >  2:1199823/1‑100 (MQ=255)
                                                                gTGGTCGATAAGCTGACCTACGCCGGAAACCTGATATCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCggg                                     <  1:147447/100‑1 (MQ=255)
                                                                   gTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAg                                  <  2:524308/100‑1 (MQ=255)
                                                                   gTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAg                                  <  2:29764/100‑1 (MQ=255)
                                                                    tCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGa                                 >  1:1167565/1‑100 (MQ=255)
                                                                     cGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGATCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGaa                                >  1:709105/1‑100 (MQ=255)
                                                                     cGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGAGATCGGGCAGaa                                >  2:253488/1‑100 (MQ=255)
                                                                          aGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGc                           >  1:254946/1‑100 (MQ=255)
                                                                          aGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGc                           <  2:922405/100‑1 (MQ=255)
                                                                          aGCTGACCTACGCCGGAAACCTGATATCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGc                           <  1:233066/100‑1 (MQ=255)
                                                                             tGACCTACGCCGGAAACCTGATATCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCAcg                        <  2:717731/100‑1 (MQ=255)
                                                                               aCCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCAcgcg                      <  1:315366/100‑1 (MQ=255)
                                                                                 cTACGCCGGAAACCTGATATCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTa                    <  2:1248234/100‑1 (MQ=255)
                                                                                 cTACGCCGGAAACCTGATATCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTa                    >  1:706463/1‑100 (MQ=255)
                                                                                    cGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTc                 <  2:304283/100‑1 (MQ=255)
                                                                                    cGCCGGAAACCTGATATCGCTGGCACCGGTCGCGCAAAGCGAGCGCTATGTATCTGACAAAGTTGATATCTGCGAACTGGCAGAACTGGCACGCCTa                    <  2:133182/97‑1 (MQ=255)
                                                                                      ccGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCAc               >  1:1018996/1‑100 (MQ=255)
                                                                                      ccGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCAc               >  1:762712/1‑100 (MQ=255)
                                                                                      ccGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCAc               >  1:53353/1‑100 (MQ=255)
                                                                                        ggAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTg             <  2:1008677/100‑1 (MQ=255)
                                                                                        ggAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTg             <  2:1008617/100‑1 (MQ=255)
                                                                                           aaCCTGATATCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAg           <  2:493570/99‑1 (MQ=255)
                                                                                             ccTGATATCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCAt        >  1:1279741/1‑100 (MQ=255)
                                                                                              cTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATc       >  1:1290603/1‑100 (MQ=255)
                                                                                               tGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCa      >  1:904681/1‑100 (MQ=255)
                                                                                               tGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCa      >  1:562079/1‑100 (MQ=255)
                                                                                               tGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCa      >  1:195395/1‑100 (MQ=255)
                                                                                                   gTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCa  <  2:1247768/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCA  >  NC_000913/3972548‑3972746

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑