Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,378,308 A→T 7.1% intergenic (‑37/‑156) 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,3080AT 7.1% 171.2 / 3.7 56intergenic (‑37/‑156)grpE/recNheat shock protein GrpE/DNA repair/recombination protein
Reads supporting (aligned to +/- strand):  ref base A (22/30);  new base T (2/2);  total (24/32)
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.76e-01

TCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCA  >  NC_002947/5378159‑5378438
                                                                                                                                                     |                                                                                                                                  
tCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCa                                                                                                                                    >  1:1307770/1‑150 (MQ=255)
 ccTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAg                                                                                                                                   <  1:283544/150‑1 (MQ=255)
     ggCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTc                                                                                                                               <  2:834452/150‑1 (MQ=255)
       cTGCGGCCAGCTTTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTTCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGc                                                                                                                             <  1:646186/150‑1 (MQ=255)
          cGGCCCGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCAGCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCCGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCtt                                                                                                                          <  2:122412/150‑1 (MQ=255)
                       cTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCtat                                                                                                             >  2:62439/1‑150 (MQ=255)
                              tCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGc                                                                                                      <  2:858708/150‑1 (MQ=255)
                              tCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGc                                                                                                      >  1:198590/1‑150 (MQ=255)
                                  gAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTGATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGa                                                                                                  >  2:192611/1‑150 (MQ=255)
                                  gAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGa                                                                                                  >  2:438884/1‑150 (MQ=255)
                                      acGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGaaaaa                                                                                              >  1:1053410/1‑150 (MQ=255)
                                          gTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAg                                                                                          <  1:62439/150‑1 (MQ=255)
                                          gTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAg                                                                                          <  1:213738/150‑1 (MQ=255)
                                            aTCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGGGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGTGCCGGAAAAACCAGGt                                                                                        >  2:647222/1‑150 (MQ=255)
                                                 cTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAg                                                                                   <  1:192611/150‑1 (MQ=255)
                                                   gCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAggg                                                                                 <  1:401746/150‑1 (MQ=255)
                                                        gCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCCTCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAg                                                                            <  1:1379024/150‑1 (MQ=255)
                                                           tGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTa                                                                         >  1:180949/1‑150 (MQ=255)
                                                            gCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAg                                                                        >  1:300998/1‑150 (MQ=255)
                                                                cGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCgg                                                                                                                 >  2:1323785/1‑105 (MQ=255)
                                                                cGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCgg                                                                                                                 <  1:1323785/105‑1 (MQ=255)
                                                                 ggTTTCTTCCACGTTAAGGTCTTGCTCGTTCAGCTGCTCATCAGCCATGGTGTCGCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTc                                                                   >  2:765139/1‑150 (MQ=255)
                                                                     tcttcCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGa                                                               <  2:180949/150‑1 (MQ=255)
                                                                       ttcCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACt                                                             >  2:19198/1‑150 (MQ=255)
                                                                       ttcCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACt                                                             >  1:807235/1‑150 (MQ=255)
                                                                           aCGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCgg                                                                                                                 <  1:250776/94‑1 (MQ=255)
                                                                           aCGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCgg                                                                                                                 >  2:250776/1‑94 (MQ=255)
                                                                                aaGGTCTTTCTCGTTGAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCTTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATaa                                                    <  2:417788/150‑1 (MQ=255)
                                                                                 aGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAAt                                                   <  2:1307770/150‑1 (MQ=255)
                                                                                         ctcGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAg                                           <  2:198590/150‑1 (MQ=255)
                                                                                           cGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAggg                                         <  1:19198/150‑1 (MQ=255)
                                                                                             ttCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGc                                       >  2:452161/1‑150 (MQ=255)
                                                                                                aGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCAtt                                    >  1:427723/1‑150 (MQ=255)
                                                                                                  ctgctCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTgg                                  <  2:427723/150‑1 (MQ=255)
                                                                                                     ctCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGaa                                                                                                 >  2:161075/1‑84 (MQ=255)
                                                                                                     ctCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGaa                                                                                                 <  1:161075/84‑1 (MQ=255)
                                                                                                             ccATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGaaa                       >  2:1131902/1‑150 (MQ=255)
                                                                                                                tGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGaaaaaa                    <  2:807235/150‑1 (MQ=255)
                                                                                                                 ggggTCTCTCCTGTGCAATTTTGTGGACTGCGAGCAAGGCTCGCCTTGCATGTCGGCTATATGGCGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAg                   <  1:647222/150‑1 (MQ=255)
                                                                                                                     tctctcCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACt               <  1:438884/150‑1 (MQ=255)
                                                                                                                       tctcCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAAACTGCGAAgg                                          >  1:233863/1‑121 (MQ=255)
                                                                                                                       tctcCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTAGATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAAACTGCGAAgg                                          <  2:233863/121‑1 (MQ=255)
                                                                                                                               gcAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTataaataa     <  1:765139/150‑1 (MQ=255)
                                                                                                                                cAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAAc    <  2:300998/150‑1 (MQ=255)
                                                                                                                                 aaTTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAgg                                                                                  <  1:1099378/71‑1 (MQ=255)
                                                                                                                                 aaTTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAgg                                                                                  >  2:1099378/1‑71 (MQ=255)
                                                                                                                                 aaTTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAAcc   <  1:1111438/150‑1 (MQ=255)
                                                                                                                                   ttttGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGc                             >  2:29414/1‑122 (MQ=255)
                                                                                                                                   ttttGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGc                             <  1:29414/122‑1 (MQ=255)
                                                                                                                                    tttGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATa                                                     <  2:161880/97‑1 (MQ=255)
                                                                                                                                    tttGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATa                                                     >  1:161880/1‑97 (MQ=255)
                                                                                                                                           aaggCGAGCCTGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCa                                                                  >  1:18397/4‑77 (MQ=255)
                                                                                                                                           aaggCGAGCCTGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCa                                                                  <  2:18397/74‑1 (MQ=255)
                                                                                                                                           aaggCGAGCCTGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCa  >  1:1164727/4‑141 (MQ=255)
                                                                                                                                           aaggCGAGCCTGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCa  <  2:1164727/138‑1 (MQ=255)
                                                                                                                                                gAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGt                  <  2:119479/120‑1 (MQ=255)
                                                                                                                                                gAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGt                  >  1:119479/1‑120 (MQ=255)
                                                                                                                                                     |                                                                                                                                  
TCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCA  >  NC_002947/5378159‑5378438

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: ‑