Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,050,872 G→A 9.0% P578P (CCG→CCA)  spoT → bifunctional (p)ppGpp synthase/hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,050,8720GA 9.0% 120.8 / 2.5 45P578P (CCG→CCA) spoTbifunctional (p)ppGpp synthase/hydrolase
Reads supporting (aligned to +/- strand):  ref base G (22/18);  new base A (2/2);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

GGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAATGTCG  >  NC_002947/6050732‑6051020
                                                                                                                                            |                                                                                                                                                       
ggCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGAc                                                                                                                                                   <  7:232198/147‑1 (MQ=255)
 gCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAATGCCCACTGGCGATCCGCGGCCCCGAAGGCCTGGTGAGAAGCTACGCCAAGTGCTGCACGCCCATACCGGGGGcccc                                                                                                                                                 >  5:91389/1‑148 (MQ=255)
                       cgcgcgCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTcggccg                                                                                                                          <  8:276591/149‑1 (MQ=255)
                        gcgcgcCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTcggccgg                                                                                                                         >  1:165048/1‑149 (MQ=255)
                          gcgcCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTcggccggc                                                                                                                        >  7:51256/1‑148 (MQ=255)
                                ccTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCa                                                                                                                 >  3:235492/1‑149 (MQ=255)
                                   gctgTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGCGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATgg                                                                                                              >  3:160459/1‑149 (MQ=255)
                                       tCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTc                                                                                                            <  3:274973/147‑1 (MQ=255)
                                       tCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGt                                                                                                          <  4:235492/149‑1 (MQ=255)
                                         ggCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTg                                                                                                         <  4:132694/148‑1 (MQ=255)
                                         ggCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGc                                                                                                        <  7:272913/149‑1 (MQ=255)
                                           ccGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCAc                                                                                                      <  6:65142/149‑1 (MQ=255)
                                                 gCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTg                                                                                                   >  1:164898/1‑146 (MQ=255)
                                                    aaCAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAAc                                                                                          >  5:42382/1‑149 (MQ=255)
                                                         cTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCg                                                                                     >  1:116799/1‑149 (MQ=255)
                                                           gCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGc                                                                                    >  2:198381/1‑148 (MQ=255)
                                                             cGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGCGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAAc                                                                                 <  4:160459/149‑1 (MQ=255)
                                                               gcgcCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACACGATTGTCGGCCACCTGTCGGCAGGCAAGGGCATGGTCGTGCCCCTTGA‑‑‑GAACTGCGGCAACAt                                                                               >  8:127591/1‑149 (MQ=255)
                                                                   cAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTgg                                                                                                  >  3:169849/1‑129 (MQ=255)
                                                                   cAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTgg                                                                                                  <  4:169849/129‑1 (MQ=255)
                                                                   cAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGt                                                                           <  1:218125/149‑1 (MQ=255)
                                                                           ccACTGGCGATCCGCGGCACCGAAGGCCTGGCGCGCAGCTACGCCTAGTGCTGCACTCCGATCCCGGGGGACCCGATTGTCGGCCAGCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCg                                                                   <  8:293439/149‑1 (MQ=255)
                                                                                ggCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCcac                                                               >  7:122237/1‑148 (MQ=255)
                                                                                 gCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACaa                                                             >  6:283836/1‑149 (MQ=255)
                                                                                 gCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACATGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACaa                                                             >  4:201969/1‑149 (MQ=255)
                                                                                             aCCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAgtg                                                 >  7:253749/1‑149 (MQ=255)
                                                                                                gAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGATGGGTACAGCCGCCACCTGAGAGATCTCCGCCACACCCCCGAagagag                                                 >  1:110564/1‑147 (MQ=255)
                                                                                                 aaGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGAGAAATCCTCCACAACCCCGAAAAgtgcgtt                                             >  8:60215/1‑148 (MQ=255)
                                                                                                    gCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCaa                                          <  7:60215/149‑1 (MQ=255)
                                                                                                         gTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACtctc                                     <  8:253749/149‑1 (MQ=255)
                                                                                                               aGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGc                               <  2:164898/149‑1 (MQ=255)
                                                                                                                 cTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGcc                              >  8:177140/1‑148 (MQ=255)
                                                                                                                    cGCCAAGTGCTGCACGCCGATCCCTGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAg                           <  7:127591/148‑1 (MQ=255)
                                                                                                                    cGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCaa                            <  8:122237/147‑1 (MQ=255)
                                                                                                                      ccAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGa                         >  4:143774/1‑148 (MQ=255)
                                                                                                                             ctgcACGCCGATCCCAGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTg                 >  3:306059/1‑149 (MQ=255)
                                                                                                                             ctgcACGCCGATCCCAGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTg                 >  3:301324/1‑149 (MQ=255)
                                                                                                                                cACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGc               <  6:42382/148‑1 (MQ=255)
                                                                                                                                 aCGCCGATCCCAGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGa             <  4:306059/149‑1 (MQ=255)
                                                                                                                                 aCGCCGATCCCAGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGa             <  4:301324/149‑1 (MQ=255)
                                                                                                                                  cGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGATATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAg            >  7:209645/1‑149 (MQ=255)
                                                                                                                                   gCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGt           >  5:102240/1‑149 (MQ=255)
                                                                                                                                        tCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAAt      >  2:11659/1‑149 (MQ=255)
                                                                                                                                          ccGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAATGt    <  8:51256/149‑1 (MQ=255)
                                                                                                                                           cGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAATGTc   <  2:165048/149‑1 (MQ=255)
                                                                                                                                            gggTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCTTGGTCGTGCACCTGGA‑‑‑GATCTGCCGCAACATGAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAAGTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAATGTCg  <  2:116799/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                       
GGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGA‑‑‑GAACTGCCGCAACATCAGTGAAATCCGCCACAACCCCGAAAAGTGCGTGCAACTCTCCTGGGCCAAGGACATCACTGGCGAGTTCAATGTCG  >  NC_002947/6050732‑6051020

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑