Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,966,363 A→G 7.3% intergenic (‑172/+30) flhB ← / ← cheZ flagellar biosynthesis protein FlhB/chemotaxis protein CheZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,966,3630AG 7.3% 161.4 / 3.2 55intergenic (‑172/+30)flhB/cheZflagellar biosynthesis protein FlhB/chemotaxis protein CheZ
Reads supporting (aligned to +/- strand):  ref base A (21/30);  new base G (2/2);  total (23/32)
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.31e-01

TAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCG  >  NC_000913/1966215‑1966511
                                                                                                                                                    |                                                                                                                                                    
tAAGCGTGAATGATTCCAGTGCGCTAAGCGTTCAATGGTTTGAGTAAGTGGCTATACAGGCTGGATTTAGGGGTTTTGCTGGCTCTCATCTAGGATAGTGTGCGTTTGTCGGATGCGGCGTCTTCGCCTTATCAGACCGCCTGATATGAc                                                                                                                                                     <  1:664044/150‑1 (MQ=255)
    cGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTgg                                                                                                                                                 >  1:1121073/1‑150 (MQ=255)
          tgatgCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGc                                                                                                                                           >  1:830832/1‑150 (MQ=255)
                 agagCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTTGTCACGCCACATCa                                                                                                                                    <  2:539254/150‑1 (MQ=255)
                  gagCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAg                                                                                                                                   <  1:389002/150‑1 (MQ=255)
                      gcAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATGTGGCGTGACCACGTCATATCAGGcgg                                                                                                                               <  1:141151/150‑3 (MQ=255)
                       cAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATGTGGCGTGACCACGTCATATCAGGcggt                                                                                                                              <  2:573527/150‑4 (MQ=255)
                       cAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAAt                                                                                                                              <  1:712902/150‑1 (MQ=255)
                        aaaGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATa                                                                                                                             <  2:167161/150‑1 (MQ=255)
                            cGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACaaa                                                                                                                         <  2:1324697/150‑1 (MQ=255)
                              ttCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATc                                                                                                                       <  2:51894/150‑1 (MQ=255)
                                      tttGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCa                                                                                                               >  2:1168241/1‑150 (MQ=255)
                                        tGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAg                                                                                                             <  2:1091380/150‑1 (MQ=255)
                                         gAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGa                                                                                                            <  1:651835/150‑1 (MQ=255)
                                            tAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTa                                                                                                         <  2:196515/150‑1 (MQ=255)
                                            tAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTa                                                                                                         >  2:1357750/1‑150 (MQ=255)
                                              aGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACAAATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATc                                                                                                       >  2:891827/1‑150 (MQ=255)
                                                 ggCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCcaa                                                                                                    >  1:1299520/1‑150 (MQ=255)
                                                        cAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCg                                                                                             <  1:1357750/150‑1 (MQ=255)
                                                           gCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTcc                                                                                          <  1:1065326/150‑1 (MQ=255)
                                                             gggATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCAc                                                                                        >  1:10718/1‑150 (MQ=255)
                                                             gggATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCAc                                                                                        <  2:769184/150‑1 (MQ=255)
                                                                   tAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATc                                                                                  <  2:185934/150‑1 (MQ=255)
                                                                    aGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTTTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTAACCAACAAATCGTACACCTGATcc                                                                                 >  2:1093408/1‑150 (MQ=255)
                                                                     gggCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCcac                                                                                                                                        >  1:608837/1‑94 (MQ=255)
                                                                     gggCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCcac                                                                                                                                        <  2:608837/94‑1 (MQ=255)
                                                                       gCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACAATCCAACAAATCGTCCACCTGATCCTGa                                                                              >  2:719493/1‑150 (MQ=255)
                                                                             gctgcCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCt                                                                        >  2:1282681/1‑150 (MQ=255)
                                                                              ctgcCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTa                                                                       <  2:566362/150‑1 (MQ=255)
                                                                              ctgcCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTa                                                                       <  2:143683/150‑1 (MQ=255)
                                                                                     acaTCAAGCATGGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACcacaac                                                                >  2:1180266/1‑148 (MQ=255)
                                                                                         cAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCt                                                            >  1:806307/1‑150 (MQ=255)
                                                                                          aaGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCtt                                                           >  2:524729/1‑150 (MQ=255)
                                                                                              aTAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCt                                                       >  1:1299432/1‑150 (MQ=255)
                                                                                                aGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTgg                                                     <  1:1368164/150‑1 (MQ=255)
                                                                                                  ggtgCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTa                                                   <  2:778096/149‑1 (MQ=255)
                                                                                                   gtgCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTAt                                                  <  1:1168241/150‑1 (MQ=255)
                                                                                                   gtgCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTAt                                                  >  2:1184287/1‑150 (MQ=255)
                                                                                                   gtgCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTAt                                                  >  2:1189253/1‑150 (MQ=255)
                                                                                                        tttGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGAcc                                             <  1:1189253/150‑1 (MQ=255)
                                                                                                        tttGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGAcc                                             <  1:1184287/150‑1 (MQ=255)
                                                                                                          tGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTg                                           <  2:288389/150‑1 (MQ=255)
                                                                                                          tGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTg                                           >  1:532155/1‑150 (MQ=255)
                                                                                                             cGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAgg                                        >  2:773522/1‑150 (MQ=255)
                                                                                                                       gTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGc                              <  2:996492/150‑1 (MQ=255)
                                                                                                                              ccTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAActct                       <  1:524729/150‑1 (MQ=255)
                                                                                                                                     gtattgCCTGATGTGGCGTGACCACGTCATATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTc                >  1:202840/6‑150 (MQ=38)
                                                                                                                                     aGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTc                <  2:10718/150‑1 (MQ=255)
                                                                                                                                      gACCGCCTGATATGACGTTGTAACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCAAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCa               <  1:746232/150‑1 (MQ=255)
                                                                                                                                       attgCCTGATGTGGCGTGACCACGTCATATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCAc              >  1:861165/4‑150 (MQ=255)
                                                                                                                                             tGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTg        <  1:773522/150‑1 (MQ=255)
                                                                                                                                             tGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAACAAACTCTGGTTTTCACGTTTTg        <  2:1100365/150‑1 (MQ=255)
                                                                                                                                              gATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTgg       <  2:532155/150‑1 (MQ=255)
                                                                                                                                              gATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTgg       >  2:190696/1‑150 (MQ=255)
                                                                                                                                                   gACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGAcgcg  >  1:587494/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCG  >  NC_000913/1966215‑1966511

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