Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,594,292 G→A 8.7% G391G (GGC→GGT)  kgtP ← alpha‑ketoglutarate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,594,2920GA 8.7% 151.2 / 2.6 56G391G (GGC→GGT) kgtPalpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base G (30/21);  new base A (3/2);  total (33/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 = 9.86e-01

TGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGG  >  NC_002947/1594153‑1594430
                                                                                                                                           |                                                                                                                                          
tGGTAATCCGCGAGTGCTGGTTGGTGTCTTTCAGGGTCGCGTAAACCAGCAGCGAGCAGGCGATACCGGGGGGGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGCACCGCCGaata                                                                                                                                  <  1:620298/150‑3 (MQ=255)
    aaTCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGa                                                                                                                              <  2:909460/150‑1 (MQ=255)
                   ctTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATa                                                                                                               >  2:676917/1‑150 (MQ=255)
                           cTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGc                                                                                                       <  2:286492/150‑1 (MQ=255)
                            tttCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGcc                                                                                                      >  2:980949/1‑150 (MQ=255)
                                  ggTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGcc                                                                                                >  2:1048952/1‑150 (MQ=255)
                                              cagcagCGAGCAGGCGATACAGGCGGAGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACTGCGCAACCTTGTTGGTGCCGCTCCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCACCGCCCGGGGCGCGAAt                                                                                    >  1:78996/1‑150 (MQ=255)
                                                          ggCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGaa                                                                        <  2:218743/150‑1 (MQ=255)
                                                           gCGATACAGGCGGTGACGTACCAGTAGATGCAGCCCCCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCCACACCGTCAGAGCATACAGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAAc                                                                       >  1:820805/1‑150 (MQ=255)
                                                           gCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAAc                                                                       <  1:365316/150‑1 (MQ=255)
                                                           gCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAAc                                                                       >  1:555675/1‑150 (MQ=255)
                                                                  aGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCgg                                                                >  1:169646/1‑150 (MQ=255)
                                                                         gACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGAc                                                         >  2:1033919/1‑150 (MQ=255)
                                                                            gTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGGATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCAc                                                      <  2:268784/150‑1 (MQ=255)
                                                                            gTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCAc                                                      <  2:767398/150‑1 (MQ=255)
                                                                                cAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCg                                                  >  2:82168/1‑150 (MQ=255)
                                                                                 aGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGt                                                 >  2:528621/1‑150 (MQ=255)
                                                                                 aGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGt                                                 >  1:1002778/1‑150 (MQ=255)
                                                                                  gTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGtt                                                >  2:736136/1‑150 (MQ=255)
                                                                                  gTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGtt                                                <  1:1033919/150‑1 (MQ=255)
                                                                                  gTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGtt                                                >  2:178665/1‑150 (MQ=255)
                                                                                  gTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCAACGATCGACACTGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGtt                                                >  2:623177/1‑150 (MQ=255)
                                                                                     gAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACg                                                     <  2:108670/142‑1 (MQ=255)
                                                                                     gAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACg                                                     >  1:108670/1‑142 (MQ=255)
                                                                                     gAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGAt                                             >  1:750946/1‑150 (MQ=255)
                                                                                      aaGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGGGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATc                                            >  1:631319/1‑150 (MQ=255)
                                                                                        gCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATc                                            <  1:891875/148‑1 (MQ=255)
                                                                                        gCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATc                                            >  2:891875/1‑148 (MQ=255)
                                                                                        gCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGa                                          >  1:234026/1‑150 (MQ=255)
                                                                                        gCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGa                                          >  1:19379/1‑150 (MQ=255)
                                                                                         ccGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAg                                         <  1:82168/150‑1 (MQ=255)
                                                                                         ccGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAg                                         <  1:736136/150‑1 (MQ=255)
                                                                                         ccGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGCCTTTGACCACGGCGTTGATCGAg                                         <  2:555675/150‑1 (MQ=255)
                                                                                          cGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGGTCGAgg                                        >  2:335069/1‑150 (MQ=255)
                                                                                           gCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGt                                       <  1:1048952/150‑1 (MQ=255)
                                                                                                ccATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGc                                  >  1:396174/1‑150 (MQ=255)
                                                                                                ccATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGc                                  <  1:528621/150‑1 (MQ=255)
                                                                                                 cATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGcc                                 >  1:928630/1‑150 (MQ=255)
                                                                                                 cATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGcc                                 <  1:551905/150‑1 (MQ=255)
                                                                                                  aTCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGc                                  >  2:708644/1‑148 (MQ=255)
                                                                                                  aTCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGc                                  <  1:708644/148‑1 (MQ=255)
                                                                                                       ggCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTAc                           >  2:416567/1‑150 (MQ=255)
                                                                                                         cgcTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTATGACCACGGCGTTGATCGAGGTATAGCCGCTTACga                         <  1:261613/150‑1 (MQ=255)
                                                                                                              ttGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATc                    >  1:888036/1‑150 (MQ=255)
                                                                                                               tGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCa                   <  2:19379/150‑1 (MQ=255)
                                                                                                                   ccACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCACCGCAAAGGGTTCATGTTTCGGTAGGGGATAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGc               <  1:980949/150‑1 (MQ=255)
                                                                                                                    cacaGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCg              >  1:63232/1‑150 (MQ=255)
                                                                                                                    cacaGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCg              >  1:406640/1‑150 (MQ=255)
                                                                                                                    cacaGGGCCACGTATTCGGCGGTGCAGCCGAAGATCGACCCGGTCAGGAGCTACGGCAGGCCGACGCCCAGGGCGCGAATTGCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCg              >  2:971831/1‑150 (MQ=255)
                                                                                                                         ggCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCAt         <  1:676917/150‑1 (MQ=255)
                                                                                                                         ggCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCAt         <  2:206836/150‑1 (MQ=255)
                                                                                                                           ccACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGa       >  2:994489/1‑150 (MQ=255)
                                                                                                                              gatcttCGGCGGCACCGCCGAATATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGc            <  2:555163/138‑1 (MQ=255)
                                                                                                                              gatcttCGGCGGCACCGCCGAATATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGc            >  1:555163/5‑142 (MQ=255)
                                                                                                                               atcttCGGCGGCACCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAg   >  1:726935/4‑150 (MQ=255)
                                                                                                                                tcttCGGCGGCACCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAgg  >  2:403965/3‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                          
TGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGG  >  NC_002947/1594153‑1594430

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