Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,645,521 G→A 7.2% G14S (GGC→AGC)  fusB → elongation factor G 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,645,5210GA 7.2% 150.6 / 2.8 56G14S (GGC→AGC) fusBelongation factor G 2
Reads supporting (aligned to +/- strand):  ref base G (28/23);  new base A (2/2);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.21e-01

GGCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTTGGTGGAATCCGTCCACCCACAGGTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATGGACTGGATGGCCCAGGAGCAGGAACGCGGCATCACCAT  >  NC_002947/4645376‑4645669
                                                                                                                                                 |                                                                                                                                                    
ggCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTTGGTGGAATCCGTCCACCCACAGGTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCAt                                                                                                                                                  >  1:511086/1‑150 (MQ=255)
        gATGCTGCAAGCCCAGACCGACTGACGGCCCGGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTTGGTGGAATCCGTCCACCCACAGGTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGccc                                                                                                                                          >  1:683455/1‑150 (MQ=255)
         aTGCTGCAAGCCCAGACCGACTGACGGCCCGGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTTGGTGGAATCCGTCCACCCACAGGGAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCa                                                                                                                                         >  1:605853/1‑150 (MQ=255)
                        aCCGACTGACGGCCCAGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTTGGTGGAATCCGTCCACCCACAGGTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCaa                                                                                                                          >  1:359795/1‑150 (MQ=255)
                                   gCCCGGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTTGGTGGAATCCGTCCACCCACAGGTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAgaccacgacca                                                                                                               >  2:621076/1‑150 (MQ=255)
                                          ggcagggcagGATTCCACCGGTCTGCTAACTTGGTGGAATCCGTCCACCCACAGGTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCg                                                                                                        <  2:511086/150‑1 (MQ=255)
                                                    gATTCCACCGGTCTGCTAACTTGGTGGAATCCGTCCACCCACAGGTAGTCCGCCATGGCCCGCACAACGCCCATCAAGCTGTACCGCAATATCGGCATCGTCCCCCACGTGCATGCCTGCAAGACCACGACCCCCGAGAGGATccggcgc                                                                                              >  2:53589/1‑145 (MQ=255)
                                                             tgTCTGCTAACTTGGTGGAATCCGTCCACCCATAGGTAGTCCGCCATGGCCCGCACAACGCCTATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCggg                                                                                     <  2:78680/149‑1 (MQ=255)
                                                               tCTGCTAACTTGGTGGAATCCGTCCACCCACAGGTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGt                                                                                   >  1:609168/1‑150 (MQ=255)
                                                                 tGCTAACTTGGTGGAATCCGTCCACCCACAGGTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCa                                                                                 >  1:551461/1‑150 (MQ=255)
                                                                           gtggAATCCGTCCACCCACAGGTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGTGGATCCTGTTCTACACCGGGGTCAACCACAAGAt                                                                       <  2:483517/150‑1 (MQ=255)
                                                                            tggAATCCGTCCACCCACAGGTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATg                                                                      >  2:188789/1‑150 (MQ=255)
                                                                                    gTCCACCCACAGGTAGTCCGCCATGGCCCGGACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGt                                                              >  2:9940/1‑150 (MQ=255)
                                                                                       cacccacAGGTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCa                                                           >  1:139460/1‑150 (MQ=255)
                                                                                        acccacAGGTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCAc                                                          >  1:207427/1‑150 (MQ=255)
                                                                                          ccacAGGTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGa                                                        <  2:271348/150‑1 (MQ=255)
                                                                                            acaGGTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGGCGACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATg                                                      >  1:90085/1‑150 (MQ=255)
                                                                                              aGGTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGc                                                            >  1:134873/1‑142 (MQ=255)
                                                                                              aGGTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGc                                                            <  2:134873/142‑1 (MQ=255)
                                                                                                gTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGa                                                        <  1:159407/144‑1 (MQ=255)
                                                                                                gTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGa                                                        >  2:159407/1‑144 (MQ=255)
                                                                                                gTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGgcgc                                                  >  2:695175/1‑150 (MQ=255)
                                                                                                    tCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCgcggcg                                              >  2:337414/1‑150 (MQ=255)
                                                                                                      cGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGc                                                                                                                                            >  1:579611/1‑54 (MQ=255)
                                                                                                      cGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGc                                                                                                                                            <  2:579611/54‑1 (MQ=255)
                                                                                                      cGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGAc                                            <  1:442403/150‑1 (MQ=255)
                                                                                                        ccATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATggg                                                                    <  1:560373/124‑1 (MQ=255)
                                                                                                        ccATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATggg                                                                    >  2:560373/1‑124 (MQ=255)
                                                                                                           tGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATgg                                       <  1:621076/150‑1 (MQ=255)
                                                                                                             gCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGTAACGGCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATGGAc                                     <  1:196672/150‑1 (MQ=255)
                                                                                                             gCCCGCACAACGCCCATCGAGCTGTACCGCAATATCAGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCggg                                                                                     <  2:472253/102‑1 (MQ=255)
                                                                                                             gCCCGCACAACGCCCATCGAGCTGTACCGCAATATCAGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCggg                                                                                     <  2:472242/102‑1 (MQ=255)
                                                                                                             gCCCGCACAACGCCCATCGAGCTGTACCGCAATATCAGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCggg                                                                                     >  1:472253/1‑102 (MQ=255)
                                                                                                             gCCCGCACAACGCCCATCGAGCTGTACCGCAATATCAGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCggg                                                                                     >  1:472242/1‑102 (MQ=255)
                                                                                                                   acaACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATGGACTGGATg                               >  2:390198/1‑150 (MQ=255)
                                                                                                                      aCGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATGGACTGGATGGcc                            <  1:9940/150‑1 (MQ=255)
                                                                                                                       tctccATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATGGACTGGATGGccc                           <  1:724931/147‑1 (MQ=255)
                                                                                                                        gCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATGGACTGGATGGCCCa                          <  2:207427/150‑1 (MQ=255)
                                                                                                                          ccATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATGGACTGGATGGCCCAgg                        >  1:458981/1‑150 (MQ=255)
                                                                                                                           cATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGtt                                                                                               >  2:187204/1‑78 (MQ=255)
                                                                                                                           cATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGtt                                                                                               <  1:187204/78‑1 (MQ=255)
                                                                                                                           cATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATGGACTGGATGGCCCAGGa                       <  1:337414/150‑1 (MQ=255)
                                                                                                                           cATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGAGCGAAGTGAACCCTGGCGCCGCGACCATGGACCGGATGGCCCaata                       >  2:724931/1‑147 (MQ=38)
                                                                                                                            aTCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATGGACTGGATGGCCCAGGAg                      >  1:207328/1‑150 (MQ=255)
                                                                                                                              cGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATGGACTGGATGGCCCAGGAGCa                    <  1:493832/150‑1 (MQ=255)
                                                                                                                               gAGCTGCGCCACAATATCTGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATGGACTGGATGGCCCAGGAGCAg                   <  2:605853/150‑1 (MQ=255)
                                                                                                                                   tGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATGGACTGGATGGCCCAGGAGCAGGAAc               <  2:609168/150‑1 (MQ=255)
                                                                                                                                       ccGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATGGACTGGATGGCCCAGGAGCAGGAACGCgg           <  1:695175/150‑1 (MQ=255)
                                                                                                                                       ccGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGGCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATGGACTGGATGGCCCAGGAGCAGGAACCCgg           >  2:596126/1‑150 (MQ=255)
                                                                                                                                        cGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATGGACTGGATGGCCCAGGAGTAGGAACGCGGc          >  1:173787/1‑150 (MQ=255)
                                                                                                                                        cGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATGGACTGGATGGCCCAGGAGCAGGAACGCGGc          <  2:139460/150‑1 (MQ=255)
                                                                                                                                        cGCAATATCGGCATCATCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCACGGACTGGATGGCCCAg                         >  1:378966/1‑135 (MQ=255)
                                                                                                                                        cGCAATATCGGCATCATCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCACGGACTGGATGGCCCAg                         <  2:378966/135‑1 (MQ=255)
                                                                                                                                              atCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATGGACTGGATGGCCCAGGAGCAGGAAcgc             <  1:15971/141‑1 (MQ=255)
                                                                                                                                              atCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATGGACTGGATGGCCCAGGAGCAGGAAcgc             >  2:15971/1‑141 (MQ=255)
                                                                                                                                                cGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATGGACTGGATGGCCCAGGAGCAGGAACGCGGcatcaccat  <  2:683455/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTTGGTGGAATCCGTCCACCCACAGGTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACCGCAATATCGGCATCGTCGCCCACGTGGATGCCGGCAAGACCACGACCACCGAGCGGATCCTGTTCTACACCGGGGTCAACCACAAGATGGGCGAAGTGCACGATGGCGCGGCGACCATGGACTGGATGGCCCAGGAGCAGGAACGCGGCATCACCAT  >  NC_002947/4645376‑4645669

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: