Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,393,647 C→G 7.7% A254P (GCC→CCC)  PP_2997 ← protein CnmA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,393,6470CG 7.7% 160.2 / 4.4 52A254P (GCC→CCC) PP_2997protein CnmA
Reads supporting (aligned to +/- strand):  ref base C (25/23);  new base G (2/2);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

GCTCACTGGAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCTTCCAGCAGCTTCTGCCC  >  NC_002947/3393505‑3393794
                                                                                                                                              |                                                                                                                                                   
gctCACTGGAACAATCCGCGAACGACCTCCCTCAACCATCTATTTCCATTGTCCTTATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAgt                                                                                                                                                 <  5:180308/147‑1 (MQ=255)
gctCACTGGAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAgt                                                                                                                                                 <  6:353234/147‑1 (MQ=255)
 ctCACTGGAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCCTGCATATCTGcacgggtc                                                                                                                                              >  6:272896/1‑142 (MQ=255)
    aCTGGAACAATCCGCGAACGTCCTCCCTCAACCCTCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCt                                                                                                                                             <  4:5372/147‑1 (MQ=255)
    aCTGGAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGc                                                                                                                                           >  6:233085/1‑149 (MQ=255)
        gAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGa                                                                                                                                       >  7:372021/1‑149 (MQ=255)
        gAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGCCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGa                                                                                                                                       <  8:343436/149‑1 (MQ=255)
             aTCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGc                                                                                                                                  <  8:372021/149‑1 (MQ=255)
             aTCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGc                                                                                                                                  <  5:390933/149‑1 (MQ=255)
                cgcgAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCcgc                                                                                                                               >  1:65588/1‑149 (MQ=255)
                 gcgAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCcgcg                                                                                                                              <  3:232498/149‑1 (MQ=255)
                   gAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACATCGACCGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGc                                                                                                                            <  5:272896/149‑1 (MQ=255)
                      cGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAg                                                                                                                         <  8:60232/149‑1 (MQ=255)
                      cGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAg                                                                                                                         <  5:233085/149‑1 (MQ=255)
                        aCCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTc                                                                                                                                              >  7:297060/1‑126 (MQ=255)
                        aCCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTc                                                                                                                                              <  8:297060/126‑1 (MQ=255)
                             cctcAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGcc                                                                                                                    >  5:155666/1‑147 (MQ=255)
                                    cATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTc                                                                                                           >  5:158560/1‑149 (MQ=255)
                                    cATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTc                                                                                                           >  8:88463/1‑149 (MQ=255)
                                            tCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTcc                                                                                                   >  8:86503/1‑149 (MQ=255)
                                              cATTGTCCTGATGATAGCGGGCGTGCCAATGGTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCTTATCGGCCAGTGTCTGCCCGATATGCCGCGGCTAGGTGCCCAGTa                                                                                                               >  1:181700/1‑135 (MQ=255)
                                              cATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTa                                                                                                               >  1:175752/1‑135 (MQ=255)
                                              cATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTa                                                                                                               <  2:175752/135‑1 (MQ=255)
                                              cATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTccc                                                                                                  >  1:363580/1‑148 (MQ=255)
                                                   tCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGGCAGACACTGGCCGatat                                                                                                                                    >  5:142233/1‑109 (MQ=39)
                                                   tCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGGCAGACACTGGCCGatat                                                                                                                                    <  6:142233/109‑1 (MQ=39)
                                                       gatgatAGCGGGTGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGGCAGACACTGGCCGATATGc                                                                                                                                  <  8:32232/107‑1 (MQ=38)
                                                       gatgatAGCGGGTGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGGCAGACACTGGCCGATATGc                                                                                                                                  >  7:32232/1‑107 (MQ=38)
                                                        atgatAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTc                                                                                        <  6:155666/148‑1 (MQ=255)
                                                                           tGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTgg                                                                    <  2:65588/149‑1 (MQ=255)
                                                                               gTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGc                                                                >  1:391681/1‑149 (MQ=255)
                                                                                               cTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATccgcc                                                <  2:343805/149‑1 (MQ=255)
                                                                                               cTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGAGTCTGCCCGATATGCCGCGGCAAGGTGCCAAGTAGGTCGGTAGTCACGATAAACTTCCCCAGCCCCAGAAACGCAGGCAGCTCGAGCATGATccgcc                                                >  7:325179/1‑149 (MQ=255)
                                                                                                tCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATccgccg                                               >  7:70792/1‑149 (MQ=255)
                                                                                                  cACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGtt                                             >  5:164871/1‑149 (MQ=255)
                                                                                                   aCTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTc                                            <  7:86503/149‑1 (MQ=255)
                                                                                                     ttCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGAtt                                                    <  3:253633/139‑2 (MQ=255)
                                                                                                     ttCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGAtt                                                    >  4:253633/1‑138 (MQ=255)
                                                                                                     ttCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGa                                          >  2:268120/1‑149 (MQ=255)
                                                                                                     ttCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGa                                          >  7:85166/1‑149 (MQ=255)
                                                                                                       cGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGAt                                         <  2:363580/148‑1 (MQ=255)
                                                                                                        gAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATcc                                       >  5:122234/1‑149 (MQ=255)
                                                                                                        gAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATcc                                       >  4:154913/1‑149 (MQ=255)
                                                                                                        gAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATc                                        <  7:88463/148‑1 (MQ=255)
                                                                                                              gCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGccc                                 <  8:85166/149‑1 (MQ=255)
                                                                                                                     tGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGGATGATCCGCCGTTCGATCCCGGCCCGCGCCaa                          >  7:162336/1‑149 (MQ=255)
                                                                                                                      ggACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTTTCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAAt                         >  2:383864/1‑149 (MQ=255)
                                                                                                                        aCCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATcc                       >  8:115315/1‑149 (MQ=255)
                                                                                                                            gTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCt                   >  3:259721/1‑149 (MQ=255)
                                                                                                                            ctAGCCCATGCATATCGGCCAGTGTCTGCCCGCTATGCCGCGGCAAGGTGCCCACTAGGTCGGTAGTCCCGCTTATCTTCCCCAGCCCGAGAAACGCTCCCAGCTCGACCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCt                   <  2:358311/148‑1 (MQ=255)
                                                                                                                             tAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGc                    <  8:70792/147‑1 (MQ=255)
                                                                                                                                    tGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCTTCcagcag           <  6:122234/149‑1 (MQ=255)
                                                                                                                                             gCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCTTCCAGCAGCTTCTGccc  >  5:110132/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GCTCACTGGAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCTTCCAGCAGCTTCTGCCC  >  NC_002947/3393505‑3393794

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