Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,594,115 G→A 10.6% intergenic (‑271/+30) PP_1399 ← / ← kgtP hypothetical protein/alpha‑ketoglutarate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,594,1150GA10.6% 121.4 / 4.8 47intergenic (‑271/+30)PP_1399/kgtPhypothetical protein/alpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base G (20/22);  new base A (2/3);  total (22/25)
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.87e-01

GCGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGG  >  NC_002947/1593966‑1594257
                                                                                                                                                      |                                                                                                                                              
gcgTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATggg                                                                                                                                                >  2:66082/1‑150 (MQ=255)
          tGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCttt                                                                                                                                      <  2:789418/150‑1 (MQ=255)
                aGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGCAc                                                                                                                                <  1:1425132/150‑1 (MQ=255)
                 ggCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACt                                                                                                                               >  1:95754/1‑150 (MQ=255)
                 ggCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCCCCATTGCATGGAGCGCCCCttttctattc                                                                                                                               <  2:593237/150‑7 (MQ=38)
                    cGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATgggg                                                                                                                                               <  2:1210930/131‑1 (MQ=255)
                    cGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATgggg                                                                                                                                               >  1:1210930/1‑131 (MQ=255)
                     gCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTgcca                                                                                                                           <  2:95754/150‑1 (MQ=255)
                       aaGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATCAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTgccagg                                                                                                                         >  1:1229015/1‑150 (MQ=255)
                        aGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCTTTTTAGTACTgccaggc                                                                                                                        <  1:608742/150‑1 (MQ=255)
                                    gAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCa                                                                                                            >  2:226878/1‑150 (MQ=255)
                                       cGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCGATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTg                                                                                                         <  1:66082/150‑1 (MQ=255)
                                         ttCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGt                                                                                                       >  2:541370/1‑150 (MQ=255)
                                                gggCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCcgc                                                                                                <  1:1208211/150‑1 (MQ=255)
                                                  gcgcCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGa                                                                                              <  1:541370/150‑1 (MQ=255)
                                                      cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtgc                                                                                          >  2:335401/1‑150 (MQ=255)
                                                           ttCTGCCAGCCATTGCAGCTTGAAAAGGTGGAATGAATAGGGCGTTAAGAGGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCtt                                                                                     >  2:1378832/1‑150 (MQ=255)
                                                              tGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGgt                                                                                  <  1:226878/150‑1 (MQ=255)
                                                                    ccaTTGCAGCTTGATAAGGGGGAATTAATAGGGCGGTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCttt                                                                            >  1:1177071/1‑150 (MQ=255)
                                                                     caTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTc                                                                           <  2:269604/150‑1 (MQ=255)
                                                                       ttGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGTCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCAt                                                                         <  2:91472/150‑1 (MQ=255)
                                                                             cTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTcgc                                                                   >  1:182952/1‑150 (MQ=255)
                                                                                gATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTa                                                                >  1:822017/1‑150 (MQ=255)
                                                                                    aGGTGGAATTAATAGGGCGGTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGTTCCAGGCAATGGGGAGCCCCTATTAATACTGCCAGGCCAGCGTTACTCGGTGGTAATCCGCGAGTGCTGCCTGGTGTATTTGATGGTCGCGTAGAcg                                                            >  2:540649/1‑149 (MQ=255)
                                                                                       tGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACcagc                                                         >  1:1234786/1‑150 (MQ=255)
                                                                                        ggAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACcagca                                                        <  1:409311/150‑1 (MQ=255)
                                                                                        ggAATTAATACGACGTCAAGAAGGAAAAGTGTGAATATAAAATGGGCCCGCCATACAATGGGGGCCCCCTTTTAGTACTGCCAGGCCAGCGTTACGCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCACGGTCGCGTAAACcagca                                                        <  2:392512/150‑1 (MQ=255)
                                                                                            ttAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTc                                                                               <  2:1060082/123‑1 (MQ=255)
                                                                                            ttAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTc                                                                               >  1:1060082/1‑123 (MQ=255)
                                                                                               aTAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGccc                                                                                                                                          <  2:813121/61‑1 (MQ=255)
                                                                                               aTAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGccc                                                                                                                                          >  1:813121/1‑61 (MQ=255)
                                                                                               aTAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCACGGAGGGTCGCGTAAACCAGCAGCGAGCa                                                 <  1:37321/150‑1 (MQ=255)
                                                                                                       ttAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATAc                                         >  1:1131120/1‑150 (MQ=255)
                                                                                                          agaagGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTcgc                                                                   <  2:1210407/121‑1 (MQ=255)
                                                                                                          agaagGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTcgc                                                                   >  1:1210407/1‑121 (MQ=255)
                                                                                                          agaagGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGCCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGCGTCCTTCATGGCCGCGTGCACCGGCAGCGTGGAGGCGCTACAgg                                      >  2:7058/1‑150 (MQ=255)
                                                                                                           gcagGAAAAGTCTGAAGAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACGCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCATCGAGCAGGCGATACAggg                                     >  2:1096358/3‑149 (MQ=255)
                                                                                                            aagGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCg                                    <  2:1229015/150‑1 (MQ=255)
                                                                                                             agGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATggg                                                                                                                                                >  1:811381/1‑41 (MQ=255)
                                                                                                             agGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATggg                                                                                                                                                <  2:811381/41‑1 (MQ=255)
                                                                                                              ggAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGt                                  <  1:1378832/150‑1 (MQ=255)
                                                                                                                     tCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTAc                           <  1:488271/150‑1 (MQ=255)
                                                                                                                     tCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTAc                           <  2:822017/150‑1 (MQ=255)
                                                                                                                     tCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTAc                           >  2:1205499/1‑150 (MQ=255)
                                                                                                                       tGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCa                         <  1:1205499/150‑1 (MQ=255)
                                                                                                                       tGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGGGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCa                         >  2:754092/1‑150 (MQ=255)
                                                                                                                               aaaGGGGCGCCCCATTGCATGGAGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAgg                                               >  1:950916/1‑120 (MQ=38)
                                                                                                                               aaaGGGGCGCCCCATTGCATGGAGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAgg                                               <  2:950916/120‑1 (MQ=38)
                                                                                                                               aaaGGGGCGCCCCATTGCATGGAGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGt                                  >  1:450539/1‑133 (MQ=39)
                                                                                                                               aaaGGGGCGCCCCATTGCATGGAGCGCACCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGTAGCGAGCAGGCGATACAGGCGGt                                  <  2:450539/133‑1 (MQ=255)
                                                                                                                                              tGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCg   <  2:1131120/150‑1 (MQ=255)
                                                                                                                                               gCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCgg  >  1:1164105/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
GCGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGG  >  NC_002947/1593966‑1594257

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: