Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 5,378,332 +A intergenic (‑61/‑132) grpE ← / → recN heat shock protein GrpE/DNA repair/recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,378,3311.A97.6% 137.3 / ‑4.7 42intergenic (‑60/‑133)grpE/recNheat shock protein GrpE/DNA repair/recombination protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (19/22);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCGGGAGCCGACCCCATGCTGG  >  NC_002947/5378182‑5378470
                                                                                                                                                       |                                                                                                                                            
cTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCtat                                                                                                                                               >  2:1037978/1‑150 (MQ=255)
        ccAGAACACGGGTATCAGCTGCATTGCACGCACTGGCTTCCTATCTGTCCAGAATGTTCTCGTTCAGCTGCTACTTAGCTATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGc                                                                                                                                       <  1:762766/150‑1 (MQ=255)
                   gTACCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCa                                                                                                                            >  2:561509/1‑150 (MQ=255)
                    tATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAg                                                                                                                           <  1:929766/150‑1 (MQ=255)
                         gCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGTGAGCGTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTCTAATGGCGCCGGAAAAACCAGGTTCa                                                                                                                      >  1:333589/1‑150 (MQ=255)
                         gCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCtataatg                                                                                                                                          >  1:613284/1‑126 (MQ=255)
                         gCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCtataatg                                                                                                                                          <  2:613284/129‑4 (MQ=255)
                           tGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAg                                                                                                                    <  1:1037978/150‑1 (MQ=255)
                                     gCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTa                                                                                                         <  1:326990/150‑1 (MQ=255)
                                           gTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTc                                                                                                   <  1:624752/150‑1 (MQ=255)
                                                ttcCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGcc                                                                                                                                      <  2:580523/111‑1 (MQ=255)
                                                ttcCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGcc                                                                                                                                      >  1:580523/1‑111 (MQ=255)
                                                ttcCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGAc                                                                                              >  1:550891/1‑150 (MQ=255)
                                                   cACGTTGAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGa                                                                                               <  2:567044/146‑1 (MQ=255)
                                                   cACGTTGAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGa                                                                                               >  1:567044/1‑146 (MQ=255)
                                                   cACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGAGCGGTGTAGTTTTCAAGACTGc                                                                                           >  1:322831/1‑150 (MQ=255)
                                                    aCGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAg                                                                                                                           >  1:492245/1‑118 (MQ=255)
                                                    aCGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAg                                                                                                                           <  2:492245/118‑1 (MQ=255)
                                                     cGTTAAGGTCTTTCTCTTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCACAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCaa                                                                                         <  1:305633/150‑1 (MQ=255)
                                                            gTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAg                                                                                                                    >  2:872877/1‑117 (MQ=255)
                                                            gTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAg                                                                                                                    <  1:872877/117‑1 (MQ=255)
                                                                      ttCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAgggg                                                                        >  2:664735/1‑150 (MQ=255)
                                                                        cAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCa                                                                      <  2:888225/150‑1 (MQ=255)
                                                                           cgatctatcaGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGGTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTCGTCAAGACTGCAACATAATCTGCGAAGGGGCATTg                                                                   <  1:612836/144‑1 (MQ=255)
                                                                              ctCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCa                                                                <  1:169250/150‑1 (MQ=255)
                                                                                  tcaGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGGGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGACTATGAGACCTGTACATAATCTGCGAAGGGGCATTggagatcg                                                            >  2:612836/1‑144 (MQ=255)
                                                                                    aGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTg                                                          <  1:664735/150‑1 (MQ=255)
                                                                                      ccATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGaa                                                        >  2:506678/1‑150 (MQ=255)
                                                                                        aTGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGaaaa                                                      >  2:406333/1‑150 (MQ=255)
                                                                                         tGGGGTCTCTCCTGCGCAATTTTGTGGACTACAGGAGAGGCTCGCTATTGATTTCGTCTACAATGGGGCTGGAAAACCAGCTGTAGTGGTAGAAGTGAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTAGCTGGCCTGaaaaa                                                     >  1:504757/1‑150 (MQ=255)
                                                                                           gggTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAg                                                   <  1:561509/150‑1 (MQ=255)
                                                                                              tctctcCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGGTATAATGGGGCCGGGAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTAc                                                >  1:889812/1‑150 (MQ=255)
                                                                                               cgctcCTGCGCAATTTGGTGGACTGCGAGCCAGGCACGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACt                                               <  2:328299/148‑1 (MQ=255)
                                                                                                         cAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTataaataa                                     <  1:506678/150‑1 (MQ=255)
                                                                                                            ttttGTGGACTGTAAGCGAGGCTCGCCTTGGATGTCCGCTATAATATGGACTGCAACAACAAGATCGAGAGGCAAGGGAGTTGTGAAGACTGGTACATAATCTGGTAAGGGGCATTGGCAGGCCTGAAAAAAATAAAATAAAAATAACCa                                  >  2:419669/1‑150 (MQ=255)
                                                                                                               tgtgGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGAc                               <  2:550891/150‑1 (MQ=255)
                                                                                                                      tGCGAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGAttt                        >  1:282730/1‑150 (MQ=255)
                                                                                                                         gAGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCg                     <  2:504757/150‑1 (MQ=255)
                                                                                                                          aGCCAGGCTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCgg                    >  1:116243/1‑150 (MQ=255)
                                                                                                                                 cTCGCCTTGGATGTCGGCTATAATGGGGCCTGAAAAACCAGGTTCAAGGGGCAAGGTAGCTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCGGGAGCCGa             >  2:322094/1‑150 (MQ=255)
                                                                                                                                 cTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAggg                                                                         <  1:639283/90‑1 (MQ=255)
                                                                                                                                 cTCGCCTTGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAggg                                                                         >  2:639283/1‑90 (MQ=255)
                                                                                                                                       ttGGATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCGGGAGCCGACCCCAt       <  2:322831/150‑1 (MQ=255)
                                                                                                                                         ggATGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCGGGAGCCGACCCCATGc     <  1:406333/150‑1 (MQ=255)
                                                                                                                                            tGTCGGCTATAATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCGGGAGCCGACCCCATGCTgg  >  1:645934/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                            
CTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCGGGAGCCGACCCCATGCTGG  >  NC_002947/5378182‑5378470

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: