Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,083,482 T→C Q611R (CAA→CGA)  pgaB ← poly‑beta‑1,6‑N‑acetyl‑D‑glucosamine N‑deacetylase PgaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,083,4820TC96.4% 186.8 / ‑1.3 55Q611R (CAA→CGA) pgaBpoly‑beta‑1,6‑N‑acetyl‑D‑glucosamine N‑deacetylase PgaB
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (26/27);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.24e-01

TTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTTGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTGGTCAGCCGATTTTTCTGCGACACCTTCCAGATAAGGCATAGCCATAATAGCGGT  >  NZ_CP009273/1083348‑1083615
                                                                                                                                      |                                                                                                                                     
tttCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTcgt                                                                                                                                     <  2:159690/137‑1 (MQ=255)
tttCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTcgt                                                                                                                                     <  1:209815/137‑1 (MQ=255)
          gttATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTc                                                                                                                          >  2:347658/1‑138 (MQ=255)
          gttATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTc                                                                                                                          >  1:438647/1‑138 (MQ=255)
                   aaaaTTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTTGTGCCTGTAATTCTAAAATAGAt                                                                                                                <  1:344851/139‑1 (MQ=255)
                    aaaTTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGAtt                                                                                                               >  2:1242394/1‑139 (MQ=255)
                    aaaTTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGAtt                                                                                                               >  2:1169552/1‑139 (MQ=255)
                    aaaTTGTCGGGATAATAACCATAGGTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGAtt                                                                                                               <  2:1040531/139‑1 (MQ=255)
                        tGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTc                                                                                                           >  1:365030/1‑139 (MQ=255)
                                        aTAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATg                                                                                           >  2:893576/1‑139 (MQ=255)
                                          aGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGAtt                                                                                                               <  2:897345/117‑1 (MQ=255)
                                          aGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGAtt                                                                                                               >  1:897345/1‑117 (MQ=255)
                                                 cACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAAtt                                                                                  >  2:1254146/1‑139 (MQ=255)
                                                      cATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGa                                                                               <  1:893576/137‑1 (MQ=255)
                                                          cAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGt                                                                          <  2:446856/138‑1 (MQ=255)
                                                            gTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCaa                                                                       <  2:438647/139‑1 (MQ=255)
                                                            gTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCaa                                                                       <  2:279260/139‑1 (MQ=255)
                                                            gTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCaa                                                                       <  2:270123/139‑1 (MQ=255)
                                                            gTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCaa                                                                       <  1:656032/139‑1 (MQ=255)
                                                            gTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCaa                                                                       <  1:422238/139‑1 (MQ=255)
                                                                tAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTg                                                                    >  1:920342/1‑138 (MQ=255)
                                                                   tAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTAtt                                                                <  2:920342/139‑1 (MQ=255)
                                                                       cTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGAtt                                                                                                               >  2:948388/1‑88 (MQ=255)
                                                                       cTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGAtt                                                                                                               <  1:948388/88‑1 (MQ=255)
                                                                          aTCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTTGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCAtt                                                         >  2:1044606/1‑139 (MQ=255)
                                                                          aTCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGt                                                                          >  1:910532/1‑122 (MQ=255)
                                                                          aTCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGt                                                                          <  2:910532/122‑1 (MQ=255)
                                                                          aTCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCaa                                                                       >  1:1223715/1‑125 (MQ=255)
                                                                          aTCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCaa                                                                       <  2:1223715/125‑1 (MQ=255)
                                                                          aTCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCAtt                                                         >  1:1135779/1‑139 (MQ=255)
                                                                          aTCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCAtt                                                         >  1:1272320/1‑139 (MQ=255)
                                                                          aTCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCAtt                                                         >  2:428491/1‑139 (MQ=255)
                                                                                  aGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGa                                                                                             >  2:814570/1‑95 (MQ=255)
                                                                                  aGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGa                                                                                             <  1:814570/95‑1 (MQ=255)
                                                                                  aGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTa                                                                  >  2:605952/1‑122 (MQ=255)
                                                                                  aGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTa                                                                  <  1:605952/122‑1 (MQ=255)
                                                                                    cGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTGGTCAGCCGa                                               >  2:243343/1‑139 (MQ=255)
                                                                                     gAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTGGTCAGCCGAt                                              >  2:988513/1‑139 (MQ=255)
                                                                                            tgCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAggg                                                                                              >  1:621279/1‑84 (MQ=255)
                                                                                            tgCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAggg                                                                                              <  2:621279/84‑1 (MQ=255)
                                                                                                aagaaATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTGGTCAGCCGATTTTTCTGCGac                                   <  1:1242394/139‑1 (MQ=255)
                                                                                                aagaaATAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTGGTCAGCCGATTTTTCTGCGac                                   <  1:794292/139‑1 (MQ=255)
                                                                                                     aTAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTGGTCAGCCGATTTTTCTGCGACACCtt                              <  1:311748/139‑1 (MQ=255)
                                                                                                     aTAGCCTGATGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTGGTCAGCCGATTTTTCTGCGACACCtt                              <  2:502307/139‑1 (MQ=255)
                                                                                                             aTGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCaa                                                                       >  2:372085/1‑90 (MQ=255)
                                                                                                             aTGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCaa                                                                       <  1:372085/90‑1 (MQ=255)
                                                                                                              tGCTGACCATTTTTCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTGGTCAGCCGATTTTTCTGCGACACCTTCCAGCTAAg                     >  2:193230/1‑139 (MQ=255)
                                                                                                                       tttttCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCaa                                                                       >  1:893859/1‑80 (MQ=255)
                                                                                                                       tttttCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCaa                                                                       <  2:893859/80‑1 (MQ=255)
                                                                                                                       tttttCTGCCAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTGGTCAGCCGATTTTTCTGCGACACCTTCCAGATAAGGCATAGCCa            >  2:230097/1‑139 (MQ=255)
                                                                                                                               ccAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAAt                                                                      <  1:870333/73‑1 (MQ=255)
                                                                                                                               ccAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAAt                                                                      >  2:870333/1‑73 (MQ=255)
                                                                                                                               ccAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTGGTCAGCCGATTTTTCTGCGACACCTTCCAGATAAGGCATAGCCataata       <  1:1254146/136‑1 (MQ=255)
                                                                                                                               ccAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTGGTCAGCCGATTTTTCTGCGACACCTTCCAGATAAGGCATAGCCataata       <  1:1169552/136‑1 (MQ=255)
                                                                                                                               ccAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTGGTCAGCCGATTTTTCTGCGACACCTTCCAGATAAGGCATAGCCATAATAg      >  2:1290310/1‑137 (MQ=255)
                                                                                                                               ccAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTGGTCAGCCGATTTTTCTGCGACACCTTCCAGATAAGGCATAGCCATAATAg      <  1:1290310/137‑1 (MQ=255)
                                                                                                                                cAGTTTCGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTGGTCAGCCGATTTTTCTGCGACACCTTCCAGATAAGGCTAGCCATAATAGCGGt  <  2:365030/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                     
TTTCAGGTTGGTTATGCAGAAAATTGTCGGGATAATAACCATAGTTTTTCACTCCATTCAGTTGTAATAGGCTCATCCAGTGAGCGAGTTGTTGCGAAGAAATAGCCTGATGCTGACCATTTTTCTGCCAGTTTTGTGCCTGTAATTCTAAAATAGATTTGTCTTTAGCCTGAGGGATGTTTTTAATTTGATTGGTCAATTGTATTAACCATTGGTCAGCCGATTTTTCTGCGACACCTTCCAGATAAGGCATAGCCATAATAGCGGT  >  NZ_CP009273/1083348‑1083615

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