Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0020942,058,4080AC53.8% 44.6 / 114.8 78G11G (GGT→GGG) glmUbifunctional N‑acetylglucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (32/4);  new base C (0/42);  total (32/46)
Fisher's exact test for biased strand distribution p-value = 7.53e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-11
Rejected as consensus: Frequency below/above cutoff threshold.

TCCACTTCGGATACACTTTTAAGCACATGCGCCAGCATTGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTACCTAATTATACCATAATCACTGACCA  >  CP002094/2058310‑2058506
                                                                                                  |                                                                                                  
tcCACTTCGGATACACTTTTAAGCACATGCGCCAGCATTGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTa                                                                                                    >  2:1094047/1‑99 (MQ=255)
    cTTCGGATACACTTTTAAGCACATGCGCCAGCATTGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAg                                                                                                >  2:1411442/1‑99 (MQ=255)
    cTTCGGATACACTTTTAAGCACATGCGCCAGCATTGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAg                                                                                                >  1:200295/1‑99 (MQ=255)
        ggATACACTTTTAAGCACATGCGCCAGCATTGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGc                                                                                            >  2:787113/1‑99 (MQ=255)
        gaaTACACTTTTAAGCACAGGCGCCAGCATTGTTTTTCCAGCTACTTGATGCAACCCTTTTGGAAGTGCTGATTTCATTCGTGTTCTTTTCCCAGCCGc                                                                                            <  2:1661370/97‑1 (MQ=255)
         gATACACTTTTAAGCACATGCGCCAGCATTGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCt                                                                                           >  2:752517/1‑99 (MQ=255)
          aTACACTTTTACGCACATGCGCCAGCATTGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTa                                                                                          >  2:338353/1‑99 (MQ=255)
          aTACACTTTTAAGCACATGCGCCAGCATTGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTa                                                                                          >  1:1395347/1‑99 (MQ=255)
           tACACTTTTAAGCACATGCGCCAGCATTGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAg                                                                                         >  1:1331107/1‑99 (MQ=255)
             cacTTTTAAGCACATGCGCCAGCATTGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGa                                                                                       >  2:1619630/1‑99 (MQ=255)
               cTTTTAAGCACATGCGCCAGCATTGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACa                                                                                     >  2:909934/1‑99 (MQ=255)
                ttttAAGCACATGCGCCAGCATTGATTTTCCAGCTACTTGATGCAACTCTTTTGGAGGTGCTTTTTTCATTCGTTTTCCTTTCCCAGCCGCTAGGACaa                                                                                    <  2:118976/99‑1 (MQ=255)
                ttttAAGCACATGCGCCAGCATTGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTTTTTTCATTCGTTTTCCTTTCCCAGCCGCTAGGACaa                                                                                    <  2:352859/99‑1 (MQ=255)
                  ttAAGCACATGCGCCAGCATTGATTTTCCAGCTACTTGATGCAACACTTTGGGAAGTGCTGTTTTCATTCGTGTTCCTTTCCCAGCCGCTAGGACAATa                                                                                  <  2:1416438/99‑1 (MQ=255)
                     aGCACATGCGCCAGCATTGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCg                                                                               >  2:93820/1‑99 (MQ=255)
                     aGCACATGCGCCAGCATTGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCg                                                                               >  1:974194/1‑99 (MQ=255)
                       cacaTGCGCCAGCATTGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGaa                                                                             >  1:993935/1‑99 (MQ=255)
                        acaTGCGCCAGCATTTATTTTCCAGCTACTTGATGCAACCCTTTTGGAAGTGCTTATTTCATTCGTGTTCCTTTCCCAGCCGCTAGGACAATAGCGAAt                                                                            <  2:421576/99‑1 (MQ=255)
                        acaTGCGCCAGCATTGTTTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGTTTTCATTCGTGTTCCTTTCCCAGCCGCTAGGACAATAGCGAAt                                                                            <  2:283125/99‑1 (MQ=255)
                         caTGCGCCAGCATTGTTTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTCCCAGCCGCTAGGACAATAGCGAAtt                                                                           <  2:777297/99‑1 (MQ=255)
                         caTGCCCCAGCATTGTTTTTCCAGCTACTTGATACAACATTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTCCCAGCCGCTAGGACAATAGCGAAtt                                                                           <  2:937758/99‑1 (MQ=255)
                           tGCGCCAGCATTGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCAGCTAGGACAATAGCGAATTTa                                                                         >  1:796637/1‑99 (MQ=255)
                             cgcCAGCATTGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTAtt                                                                       >  2:1384870/1‑99 (MQ=255)
                                 accATTTATTTTCCAGCTCCTTGATGCAACTCTTTTGGAAGTGCTGTTTTCATTCGTGTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATa                                                                   <  2:628494/97‑1 (MQ=255)
                                 aGCATTGTTTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGTTTTCATTCGTGTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATa                                                                   <  2:1100697/99‑1 (MQ=255)
                                  gCATTGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAg                                                                  >  2:1409436/1‑99 (MQ=255)
                                     ttGTTTTTCCAGTTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGAtt                                                               <  2:1174615/99‑1 (MQ=255)
                                     ttGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGTTTTCATTCGTGTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGAtt                                                               <  2:382738/99‑1 (MQ=255)
                                     ttGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGAtt                                                               >  2:1419143/1‑99 (MQ=255)
                                     ttGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGAtt                                                               >  2:734253/1‑99 (MQ=255)
                                      tGATTTTCCAGCTACTTGATGCAACCCTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGAttt                                                              <  2:550453/99‑1 (MQ=255)
                                        tttttCCAGCTACTTGATGCAACTCTTTGGGAAGTGCTGATTTCATTCGTGTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTtct                                                            <  2:738530/98‑1 (MQ=255)
                                        aTTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTtct                                                            <  2:493289/99‑1 (MQ=255)
                                        aTTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTtct                                                            >  2:1200766/1‑99 (MQ=255)
                                        aTTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTtct                                                            >  2:579418/1‑99 (MQ=255)
                                          tttCCAGCTGCTTATTCCAACCCTTTTGGACGTGGTGATTTCATTCGTTTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTcc                                                          <  2:1206063/99‑1 (MQ=255)
                                             ccAGCTACTTGATGCAACACTTTTGGAAGTGCTTATTTCATTCGTGTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCttt                                                       <  1:1644482/99‑1 (MQ=255)
                                              cAGCTATTTGATGCAACTTTTTGGGAAGTGCTGTTTTAGTCCGTTTCCTTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCtttt                                                      <  1:1384163/99‑1 (MQ=255)
                                               aGCTACTTGATGCAACCCTTTTGGAAGTGCTGATTTCATTCGTTTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCttttt                                                     <  2:301792/99‑1 (MQ=255)
                                               aGCTACTGGATCCAACTTTTTTGGAGGTGTTTTTTTCATTCGTTTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCttttt                                                     <  1:1428505/99‑1 (MQ=255)
                                                 cTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTg                                                   >  2:991043/1‑99 (MQ=255)
                                                   aCTTGATGCAACATTTTGGGAAGTGCTGTTTTCTTTCTTTTTCTTTCCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGtt                                                 <  1:1530055/99‑1 (MQ=255)
                                                     ttGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTAc                                               <  2:192974/99‑1 (MQ=255)
                                                       gATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACaa                                             >  2:73807/1‑99 (MQ=255)
                                                        aTGCAACACTTTTGGAATTCTTACTTCCATTCTTGTTCCTTTCCCAGCCGCTAGGACAATAGCGCATTTATTCATAGATTTCTCCTTTTTTGTTACaaa                                            <  2:598772/99‑1 (MQ=255)
                                                          gCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACaaata                                          <  2:1144067/99‑1 (MQ=255)
                                                           cAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACaaataa                                         >  1:1188729/1‑99 (MQ=255)
                                                            aaCACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACaaataaa                                        >  2:673974/1‑99 (MQ=255)
                                                             acctttttGGAGTTTCTGATTTCATTCTTGTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACaaataaaa                                       <  1:752517/95‑1 (MQ=255)
                                                             acacTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACaaataaaa                                       >  1:1408937/1‑99 (MQ=255)
                                                               acTTTTGGAGGGTCTGTTTTCATTCGGTTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACaaataaaata                                     <  1:342757/99‑1 (MQ=255)
                                                                cTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAaataaaataa                                    >  1:967040/1‑99 (MQ=255)
                                                                 ttgtGGAATTGTTTATTTCTTTCGTTTTCCTTCCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAataaaataaa                                   <  1:1395215/96‑1 (MQ=255)
                                                                     ggAAGTGTTTTTTTCTTTCGTTTTCCTTCCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGt                               <  1:1200766/99‑1 (MQ=255)
                                                                     ggAAGTGCTGTTTTCTTTGGTTTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGt                               <  1:387435/99‑1 (MQ=255)
                                                                     ggAAGTGCTGATTTCATTCGTGTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGt                               <  2:1407381/99‑1 (MQ=255)
                                                                      gAAGTTCTGATTTCATTCGTGTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGtt                              <  2:701498/99‑1 (MQ=255)
                                                                        ggTGCTGATTTCATGCGATTTCCTATCCCAGACGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTATAAATAAAATAAAATGTTTa                            <  2:835090/98‑1 (MQ=255)
                                                                        aGTGCTGATTTTATTCGTTTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTa                            <  1:1481343/99‑1 (MQ=255)
                                                                        aGTGCTGATTTCATTCGTGTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTa                            <  2:725913/99‑1 (MQ=255)
                                                                         gTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTAc                           >  1:667810/1‑99 (MQ=255)
                                                                          tGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTAcc                          >  1:1065875/1‑99 (MQ=255)
                                                                          gtcTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTAcc                          <  2:865126/97‑1 (MQ=255)
                                                                            cTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTACCTa                        >  1:463638/1‑99 (MQ=255)
                                                                             ttaTTTCATTCGTGTCCCTTTCCCAGCCGCTAGGACAATAGTGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTACCTaa                       <  1:1655070/97‑1 (MQ=255)
                                                                             tGATTTCATTCGTGTTCTTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTACCTaa                       <  2:819666/99‑1 (MQ=255)
                                                                              gATTTCATGCGTGTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTACCTAAt                      <  1:1640092/99‑1 (MQ=255)
                                                                                tttCATTCGTTTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTACCTAATta                    <  2:813856/99‑1 (MQ=255)
                                                                                     ttCGTGTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTACCTAATTATACCa               <  2:253960/99‑1 (MQ=255)
                                                                                      tCGTTTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTACCTAATTATACCAt              <  1:292445/99‑1 (MQ=255)
                                                                                      tCGTGTTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTACCTAATTATACCAt              <  2:549753/99‑1 (MQ=255)
                                                                                       cGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTACCTAATTATACCATa             >  1:328038/1‑99 (MQ=255)
                                                                                       cGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTACCTAATTATACCATa             >  2:1213902/1‑99 (MQ=255)
                                                                                        gtgtTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTACCTAATTATACCATaa            <  2:1219609/99‑1 (MQ=255)
                                                                                         tgtTCCTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTACCTAATTATACCATAAt           <  1:1204878/99‑1 (MQ=255)
                                                                                         tgtTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTACCTAATTATACCATAAt           >  1:1460191/1‑99 (MQ=255)
                                                                                           ttCTTTTCCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTACCTAATTATACCATAATCa         <  1:22477/99‑1 (MQ=255)
                                                                                           ttCCTTTCCCAGCCGCTAGGACAATAGCGAATTTCTTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTACCTAATTATACCATAATCa         <  1:129997/99‑1 (MQ=255)
                                                                                            tCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTACCTAATTATACCATAATCAc        >  1:137871/1‑99 (MQ=255)
                                                                                                  aCCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTACCTAATTATACCATAATCACTGACCa  >  2:1375052/1‑99 (MQ=255)
                                                                                                  |                                                                                                  
TCCACTTCGGATACACTTTTAAGCACATGCGCCAGCATTGATTTTCCAGCTACTTGATGCAACACTTTTGGAAGTGCTGATTTCATTCGTGTTCCTTTACCAGCCGCTAGGACAATAGCGAATTTATTCATAGATTTCTCCTTTTTTGTTACAAATAAAATAAAATGTTTACCTAATTATACCATAATCACTGACCA  >  CP002094/2058310‑2058506

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