Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 321,743 G→A 9.6% intergenic (+28/+430) aguA → / ← PP_0267 agmatine deiminase/ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947321,7430GA 9.6% 136.6 / 5.9 52intergenic (+28/+430)aguA/PP_0267agmatine deiminase/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (26/21);  new base A (3/2);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01

TGGCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACG  >  NC_002947/321600‑321886
                                                                                                                                               |                                                                                                                                               
tGGCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACgg                                                                                                                                           <  2:767046/150‑1 (MQ=255)
    aaGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCct                                                                                                                                       >  1:561598/1‑150 (MQ=255)
          cTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCttttt                                                                                                                                 >  2:570987/1‑150 (MQ=255)
               ccGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCACATGACgg                                                                                                                                           >  1:129634/1‑135 (MQ=255)
               ccGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCACATGACgg                                                                                                                                           <  2:129634/135‑1 (MQ=255)
                        aGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTATTTTTTGCGCGGGCGACGc                                                                                                                   >  2:525872/1‑150 (MQ=255)
                                     cccGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGCGTATGTTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGc                                                                                                      >  1:647384/1‑150 (MQ=255)
                                            cgcgAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCctc                                                                                                                                      >  1:555678/1‑111 (MQ=255)
                                            cgcgAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCctc                                                                                                                                      <  2:555678/111‑1 (MQ=255)
                                                 gctgctCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGttt                                                                                          <  2:286726/150‑1 (MQ=255)
                                                 gctgctCTTCGGCGGCGGCCACCTTCCCTGCCTGACCCCACACCAGCCGGCGCCGCTCAAGCGCTGAGCCCGCAATGAATGCGAGGCGCGTCATGTAACGGGCCACTTTTTTACGCGGGCGAACCCCGGAAGCCGAGCTGGCCTTCGttt                                                                                          >  2:244786/1‑150 (MQ=255)
                                                  ctgctCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGtttt                                                                                         <  2:649433/150‑1 (MQ=255)
                                                        ttGGGTGGCGGCAACATTCACTGCCTGACCAAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATAATAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTg                                                                                   <  1:450637/150‑1 (MQ=255)
                                                        ttGGGTGCCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGACGTGTTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTg                                                                                   <  1:485068/150‑1 (MQ=255)
                                                          gggTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGtt                                                                                 >  2:534920/1‑150 (MQ=255)
                                                             tGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGttttt                                                                              <  1:525872/150‑1 (MQ=255)
                                                                 ggcAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAg                                                                          >  1:404117/1‑150 (MQ=255)
                                                                 ggcAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAg                                                                          >  2:129177/1‑150 (MQ=255)
                                                                  gcAACATTCACGGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCAGCGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGc                                                                         <  1:175440/150‑1 (MQ=255)
                                                                     aCATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGa                                                                      >  2:498597/1‑150 (MQ=255)
                                                                     aCATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGa                                                                      <  2:404117/150‑1 (MQ=255)
                                                                        ttCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAg                                                                   <  1:534920/150‑1 (MQ=255)
                                                                               ccTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATg                                                            >  2:207754/1‑150 (MQ=255)
                                                                                cTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATgg                                                           >  2:461457/1‑150 (MQ=255)
                                                                                 tGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGc                                                          >  2:47095/1‑150 (MQ=255)
                                                                                   aCCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGa                                                                      <  2:228361/136‑1 (MQ=255)
                                                                                   aCCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGa                                                                      >  1:228361/1‑136 (MQ=255)
                                                                                       aaCAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGAc                                                                                                                     >  2:240271/1‑85 (MQ=255)
                                                                                       aaCAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGAc                                                                                                                     <  1:240271/85‑1 (MQ=255)
                                                                                         cagcagCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAg                                                  >  2:133008/1‑150 (MQ=255)
                                                                                                       ggTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGATGGCATACGTTTTGTTTTGTTTGTCAAGCAGATAAGAAGATGTCACGGTCAGCCGGCGCCTCTGAt                                    >  2:243302/1‑150 (MQ=255)
                                                                                                       ggTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGCGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGTCGCGCCTCTGAt                                    <  1:192909/150‑1 (MQ=255)
                                                                                                         tCAAGCGCTGGGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCt                                  >  1:737859/1‑150 (MQ=255)
                                                                                                              cgcTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGttttt                                                                              >  1:720007/1‑101 (MQ=255)
                                                                                                              cgcTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGttttt                                                                              <  2:720007/101‑1 (MQ=255)
                                                                                                                cTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACaa                           <  1:129177/150‑1 (MQ=255)
                                                                                                                 tGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAAt                          <  1:122939/150‑1 (MQ=255)
                                                                                                                 tGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCCTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAAt                          >  2:214631/1‑150 (MQ=255)
                                                                                                                   aGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGt                                >  2:206912/1‑142 (MQ=255)
                                                                                                                   aGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGt                                <  1:206912/142‑1 (MQ=255)
                                                                                                                   aGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATAc                        >  1:280506/1‑150 (MQ=255)
                                                                                                                         tAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCgg                                                                                                                          <  1:610482/46‑1 (MQ=255)
                                                                                                                         tAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCgg                                                                                                                          >  2:610482/1‑46 (MQ=255)
                                                                                                                                tgaGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAg           <  2:737859/150‑1 (MQ=255)
                                                                                                                                 gagaGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGc          <  1:498597/150‑1 (MQ=255)
                                                                                                                                 aagagGCCCGTCACATGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAg                                                  <  2:248275/109‑1 (MQ=255)
                                                                                                                                 aagagGCCCGTCACATGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAg                                                  >  1:248275/2‑110 (MQ=255)
                                                                                                                                     ggCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCt      <  1:47095/150‑1 (MQ=255)
                                                                                                                                      gCCCGTCACATGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTc     >  2:22273/1‑150 (MQ=255)
                                                                                                                                       cccGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCa    >  2:25377/1‑150 (MQ=255)
                                                                                                                                        ccGTCATGTGACGGGCCTCTTTTTTGCGGGGGCGACGCCCGGAAGCCGAGCTGCCATACGTTTTGTATTGTGTTTCACGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTAATCTGCAACAATCCCTACGTAAAGTAGCCACTcgc   >  1:370857/1‑148 (MQ=255)
                                                                                                                                         cGTCATGTGACGGGCCTCTTTTTTTCTCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACg  >  2:696906/1‑150 (MQ=255)
                                                                                                                                              tgtgACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCctct                                       <  1:781544/108‑1 (MQ=255)
                                                                                                                                              tgtgACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCctct                                       >  2:781544/1‑108 (MQ=255)
                                                                                                                                               |                                                                                                                                               
TGGCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACG  >  NC_002947/321600‑321886

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: