Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 705,972:1 (T)9→10 100% intergenic (+82/‑121) nagE → / → glnS N‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components/glutamyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913705,9631.T50.9% 93.1 / 102.0 55intergenic (+73/‑130)nagE/glnSN‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base . (15/12);  new base T (16/12);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑TTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGC  >  NC_000913/705867‑706058
                                                                                                 |                                                                                                   
ccGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑ttt                                                                                                >  1:715033/1‑100 (MQ=255)
ccGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑ttt                                                                                                >  2:1489250/1‑100 (MQ=255)
 cGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑tttt                                                                                               >  1:727129/1‑100 (MQ=255)
   cggTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑tttttn                                                                                             <  2:1073466/98‑2 (MQ=255)
   cTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑tttttt                                                                                            <  1:2398961/100‑1 (MQ=255)
    tGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑ttttttt                                                                                           <  2:2069810/100‑1 (MQ=255)
     gTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑tttttttt                                                                                         >  1:2493693/1‑100 (MQ=255)
      tATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGCGATCTTCTCTGCCGC‑‑ttttttttt                                                                                        >  1:409056/1‑100 (MQ=255)
      tATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑tttttttt                                                                                         <  1:81708/99‑1 (MQ=255)
       aTGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑tttttttttt                                                                                       >  2:318643/1‑99 (MQ=255)
        tGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑TTTTTTTTtca                                                                                      <  1:147748/100‑1 (MQ=255)
          aaaTCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑ttttttttttcat                                                                                    >  2:469316/1‑96 (MQ=255)
          aaaTCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑TTTTTTTTtcatc                                                                                    >  1:237389/1‑100 (MQ=255)
          aaaTCAAAAAGTAATCTGCTTTATGCCTGATGCAACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTAGAGAGCTTCTCTGCCGC‑‑TTTTTTTTtcat                                                                                     >  1:386655/1‑99 (MQ=255)
           aaTCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑ttttttttttcatc                                                                                   <  2:2253890/100‑6 (MQ=255)
           aaTCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑TTTTTTTTtcatc                                                                                    >  2:211323/1‑99 (MQ=255)
            aTCAAAAAGTAATCTGCTTTATGCCTGATGCGCCGCTTGAGCGTCGCATCCAACAATGACAAGCTGTAGAGATCGTCGCTGCCGC‑‑gtttttattgtaaca                                                                                  >  2:2367225/1‑94 (MQ=255)
            aTCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTtcatca                                                                                   >  2:2283339/1‑100 (MQ=255)
             tCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑tttttttttttcatca                                                                                 >  1:611623/1‑93 (MQ=255)
              cAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑ttttttttttttttttt                                                                                >  1:2048165/1‑92 (MQ=255)
               aaaaaGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCatcaatc                                                                                >  2:1219505/1‑100 (MQ=255)
               aaaaaGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑ttttttttttttttttt                                                                                >  1:554140/1‑91 (MQ=255)
               aaaaaGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑TTTTTTTTTCatcaatca                                                                               >  2:1016066/1‑100 (MQ=255)
                aaaaGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCatcaatca                                                                               >  2:1894903/1‑100 (MQ=255)
                        tCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑TTTTTTTTTCATCAATCATCCCCATaa                                                                      >  2:1599396/1‑100 (MQ=255)
                          tGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑TTTTTTTTTCATCAATCATCCCCATAAt                                                                     >  2:1030450/1‑99 (MQ=255)
                           gCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGTAGATCTTGGCGGCCGC‑‑CTTTTCTTTAATCAATAAAGCCCATAAtac                                                                   >  1:2499312/1‑98 (MQ=255)
                           gCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATc                                                                    >  2:889194/1‑100 (MQ=255)
                            cTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAAtgc                                                                   >  2:1334134/1‑98 (MQ=255)
                              ttATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑TTTTTTTTTCATCAATCATCCCCATAATCCTTg                                                                <  1:236792/100‑1 (MQ=255)
                                 tGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGtt                                                              >  2:1562541/1‑100 (MQ=255)
                                 tGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGtt                                                              <  1:2272148/100‑1 (MQ=255)
                                 tGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑TTTTTTTTTCATCAATCATCCCCATAATCCTTGTTa                                                             >  1:2187976/1‑100 (MQ=255)
                                   ccTGATACGAAGCTAAAAAGTCGCATCCAACAATGCCAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTg                                                                <  2:549875/96‑1 (MQ=255)
                                      gATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGAtt                                                         <  1:1025759/100‑1 (MQ=255)
                                      gATGCGACGCTCAAGCGTCGCATCCAACAATGACAAGCGGTGGGGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGAtt                                                         >  2:909614/1‑100 (MQ=255)
                                         gCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑TTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCa                                                     >  2:437099/1‑100 (MQ=255)
                                         gCGACGCTCAAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTa                                                        >  2:991089/1‑98 (MQ=255)
                                          cGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCCTCCCCATAATCCTTGTTAGATTAtaa                                                     >  2:2342540/1‑98 (MQ=255)
                                          cGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATcc                                                     >  1:433312/1‑99 (MQ=255)
                                          cGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑TTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCaa                                                    >  2:140089/1‑100 (MQ=255)
                                          cGACGCTCAAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCa                                                     >  2:307610/1‑100 (MQ=255)
                                                   aGCGTCGCATCCAAAAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑TTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTaaaaa                                           <  2:1082666/100‑1 (MQ=255)
                                                     cGTCGCATCCAACAATGACAAGCGGAGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTaaaaaa                                          <  1:217437/100‑1 (MQ=255)
                                                        cGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTa                                       >  2:639311/1‑100 (MQ=255)
                                                           agaaaacaaTGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAAc                                     <  1:414660/95‑1 (MQ=255)
                                                             ccAACAATGACAGCGGTGGAGATCTTCTCTGCCGC‑‑TTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTg                                >  2:2432308/1‑100 (MQ=255)
                                                              caacaaTGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAg                                   >  1:1117569/1‑98 (MQ=255)
                                                               aacaaTGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTg                                <  1:473734/100‑1 (MQ=255)
                                                                 caaTGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTc                              >  2:2454854/1‑100 (MQ=255)
                                                                 caaTGACAAGCGGTGGAGAACTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTc                              <  2:282892/100‑1 (MQ=255)
                                                                  aaTGACAAGCGGTGGAGATCTTGTCTGCCGC‑‑TGTTTTTTTACAAAAAAAATCCCCATAATCCTTGTTAGATTCTCAATTTTAAAAAACTAACa                                    >  1:1735090/1‑93 (MQ=255)
                                                                  aaTGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTc                              <  1:154407/99‑1 (MQ=255)
                                                                   aTGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑TTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCa                             >  1:35378/1‑98 (MQ=255)
                                                                     gACTAGCGGTGGAGAGATTCTCTGCCTCTTTTTTTTTTCATCAATCATCACCATAATCCTTGTTAGATTATCAATTTTAAGAAACTAACAGTTGTCAGcc                          <  1:1176249/100‑1 (MQ=255)
                                                                             gTGGAGATCTTCTCTGCCGC‑‑TTTTTATTTTCATCAATCAACCCCATCATCCTTGTTAGATTAACAAAAAAAAAAAACTAACAGTTGTCAGCCTGTCCCGc                  <  1:1265035/100‑1 (MQ=255)
                                                                             gTGGAGATCTTCTCTGCCGC‑‑TCTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCt                 <  1:2540739/100‑1 (MQ=255)
                                                                                 agaTCTTCTCTGCCGC‑‑TTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTtata             <  1:2589480/100‑1 (MQ=255)
                                                                                  gaTCTTCTCTGCCGCTTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTtat              >  2:475291/1‑100 (MQ=255)
                                                                                  gaTCTTCTCTGCCGC‑‑TATTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTta               <  2:8077/98‑1 (MQ=255)
                                                                                     cttctCTGCAGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAg           <  2:609969/99‑1 (MQ=255)
                                                                                      ttctCTGCCGC‑‑TTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATc        <  1:1016063/100‑1 (MQ=255)
                                                                                        ctctGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCa       >  2:2367729/1‑100 (MQ=255)
                                                                                        ctctGCCGC‑‑TTTTTTTATTAATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGAt         <  2:2330451/98‑1 (MQ=255)
                                                                                        ctccGCCGC‑‑TATTATTTTTCATCAATCATCCCCATAAGCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCa       <  2:1518495/100‑1 (MQ=255)
                                                                                          ctGCCGC‑‑TTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATAc    >  1:2567566/1‑100 (MQ=255)
                                                                                          ctGCCGC‑‑TTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATAc    <  2:727132/100‑1 (MQ=255)
                                                                                          ctGCCGC‑‑TTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATAc    <  2:1423052/100‑1 (MQ=255)
                                                                                           tGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATAc    <  2:1344449/100‑1 (MQ=255)
                                                                                            gCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATaa            <  1:450613/91‑1 (MQ=255)
                                                                                            gCCGC‑‑TTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGc  >  2:1863599/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
CCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑TTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGC  >  NC_000913/705867‑706058

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