Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,594,103 T→C 8.7% intergenic (‑259/+42) 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,1030TC 8.7% 145.5 / 3.4 55intergenic (‑259/+42)PP_1399/kgtPhypothetical protein/alpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base T (25/24);  new base C (2/3);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

TTTGGGCGCGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCA  >  NC_002947/1593959‑1594234
                                                                                                                                                |                                                                                                                                    
tttGGGCGCGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATg                                                                                                                                >  1:789418/1‑150 (MQ=255)
   gggCGCGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCCGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCaa                                                                                                                             <  2:1107762/150‑1 (MQ=255)
    ggCGCGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAAt                                                                                                                            >  2:1208211/1‑150 (MQ=255)
       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)
                                                                                                       tAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATgg                                                                                                                          <  2:596337/53‑1 (MQ=255)
                                                                                                       tAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATgg                                                                                                                          >  1:596337/1‑53 (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    >  2:1205499/1‑150 (MQ=255)
                                                                                                                            tCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTAc    <  2:822017/150‑1 (MQ=255)
                                                                                                                            tCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTAc    <  1:488271/150‑1 (MQ=255)
                                                                                                                              tGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCa  <  1:1205499/150‑1 (MQ=255)
                                                                                                                              tGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGGGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCa  >  2:754092/1‑150 (MQ=255)
                                                                                                                                      aaaGGGGCGCCCCATTGCATGGAGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAgg                        <  2:950916/120‑1 (MQ=38)
                                                                                                                                      aaaGGGGCGCCCCATTGCATGGAGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAgg                        >  1:950916/1‑120 (MQ=38)
                                                                                                                                      aaaGGGGCGCCCCATTGCATGGAGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGt           >  1:450539/1‑133 (MQ=39)
                                                                                                                                      aaaGGGGCGCCCCATTGCATGGAGCGCACCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGTAGCGAGCAGGCGATACAGGCGGt           <  2:450539/133‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                    
TTTGGGCGCGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCA  >  NC_002947/1593959‑1594234

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: