Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,397,181 G→C 10.5% P826A (CCG→GCG)  PP_2100 ← two‑component system sensor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,397,1810GC10.5% 123.3 / 5.3 47P826A (CCG→GCG) PP_2100two‑component system sensor protein
Reads supporting (aligned to +/- strand):  ref base G (24/18);  new base C (2/3);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

TCGCGCACCTCCAGCTCCAGTACGGCTGGTGTGCTGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGACCATCGAAAATGCGCGCCACCGATTCGACCAGCGCCGCCAGGTCGACCACCTCCGGGGCCAGGGAC  >  NC_002947/2397042‑2397320
                                                                                                                                           |                                                                                                                                           
tCGCGCACCTCCAGCTCCAGTACGGCTGGTGTGCTGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCa                                                                                                                                     <  7:307536/148‑1 (MQ=255)
tCGCGCACCTCCAGCTCCAGTACGGCTGGTGTGCTGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAg                                                                                                                                    <  3:290300/149‑1 (MQ=255)
 cgcgCACCTCCAGCTCCAGTACGGCTGGTGTGCGGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGACCCGCTGCGCTTCAAGCAgg                                                                                                                                   <  8:211025/149‑1 (MQ=255)
  gcgcACCTCCAGCTCCAGTACGGCTGGTGTGCTGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGACCCGCTGCGCTTCAAGCAgg                                                                                                                                   <  6:170312/148‑1 (MQ=255)
  gcgcACCTCCAGCTCCAGTACGGCTGGTGTGCTGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGa                                                                                                                                  <  7:129130/149‑1 (MQ=255)
        ctccagctccagTACGGCTGGTGTGCTGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACAAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATgg                                                                                                                             >  3:176197/1‑148 (MQ=255)
         tccagctccagTACGGCTGGTGTGCTGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATgg                                                                                                                             >  1:303736/1‑147 (MQ=255)
         tccagctccagTACGGCTGGTGTGCTGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATgg                                                                                                                             >  3:159943/1‑147 (MQ=255)
             gctccagTACGGCTGGTGTGCTGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGt                                                                                                                       >  8:92100/1‑149 (MQ=255)
                 cagTACGGCTGGTGTGCTGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCgg                                                                                                                   >  1:10164/1‑149 (MQ=255)
                    tACGGCTGGTGTGCTGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCg                                                                                                                >  7:278404/1‑149 (MQ=255)
                      cGGCTGGTGTGCTGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCggg                                                                                                              <  1:223759/149‑1 (MQ=255)
                      cGGCTGGTGTGCTGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCgg                                                                                                               <  4:159943/148‑1 (MQ=255)
                            gtgtGCTGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGAt                                                                                                         <  6:159077/148‑1 (MQ=255)
                             tgtgCTGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCt                                                                                                       >  1:99216/1‑149 (MQ=255)
                               tgCTGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGc                                                                                                     <  7:92100/149‑1 (MQ=255)
                                 cgggTTCCGTCATCCAGCAGCCTGACGCTTATGCGCACCTGACCGTGCTCGTTGAACTTGATGGCATTGCTCACCAGGGTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCa                                                                                                    <  4:176843/146‑1 (MQ=255)
                                  tGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCAcc                                                                                                  >  8:230848/1‑149 (MQ=255)
                                  tGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCAcc                                                                                                  >  5:306405/1‑149 (MQ=255)
                                    gTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACATGATGGCATTGCTCCCCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGAAGTACATGGCAGCGTGCGGCCGGGGCGATCTGAACCTc                                                                                                >  6:345766/1‑149 (MQ=255)
                                          tCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGc                                                                                          >  1:86176/1‑149 (MQ=255)
                                            aTCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGACCCGCTGCGCTTCAAGCAg                                                                                                                                    <  1:19084/105‑1 (MQ=38)
                                            aTCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGACCCGCTGCGCTTCAAGCAg                                                                                                                                    >  2:19084/1‑105 (MQ=38)
                                             tCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAg                                                                                       <  4:176197/149‑1 (MQ=255)
                                                       tGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGc                                                                             >  5:19179/1‑149 (MQ=255)
                                                                 gcgcACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGAc                                                                   >  8:27207/1‑149 (MQ=255)
                                                                  cgcACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGAcc                                                                  <  2:33160/149‑1 (MQ=255)
                                                                         gACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAggg                                                                                           >  6:71450/1‑117 (MQ=255)
                                                                         gACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAggg                                                                                           <  5:71450/117‑1 (MQ=255)
                                                                                   cgTGAACATGATGGAATTGATCACCAGGTTTGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCCGCGTGCGGCCGGGCCGATCTTCACCTCAAGGGCCAGGCCTTTCTGCCGCGAATGACCATCGAAAATGCGAGCCa                                                 >  8:216321/2‑149 (MQ=255)
                                                                                     tGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGACCATCGAAAATGCGCGCCAcc                                               >  3:163793/1‑149 (MQ=255)
                                                                                        aCTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGACCATCGAAAATGCGCGCCACCGAt                                            >  4:147941/1‑149 (MQ=255)
                                                                                         cTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGAcc                                                                  >  5:63234/1‑126 (MQ=255)
                                                                                         cTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGAcc                                                                  <  6:63234/126‑1 (MQ=255)
                                                                                                         aCCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGACCATCg                                                              >  8:304889/1‑114 (MQ=255)
                                                                                                               ttGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCATGCCTTTATGCCGCGCCTGACCATAGAAAATGCGCGCCAACGATTCGACCAGCGCCTCCAGGTCGac                     <  2:99216/149‑1 (MQ=255)
                                                                                                               ttGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGACCATCGAAAATGCGCGCCACCGATTCGACCAGCGCCGCCAGGTCGac                     <  7:27207/149‑1 (MQ=255)
                                                                                                                 ggACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGACCATCGAAAATGCGCGCCACCGATTCGACCAGCGCCGCCAGGTCGacc                    <  2:303736/148‑1 (MQ=255)
                                                                                                                   aCAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGACCATCGAAAATGCGCGCCACCGATTCGACCAGCGCCGCCAGGTCGaccac                  >  3:74767/1‑148 (MQ=255)
                                                                                                                           tGCTTGACCCGCTGCGCTTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGACCATCGAAAATGCGCGCCACCGATTCGACCAGCGCCGCCAGGTCGACCACCTCCGGGGc         >  3:236545/1‑149 (MQ=255)
                                                                                                                              ttGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGACCATCGAAAATGCGCGCCACCGATTCGACCAGCGCCGCCAGGTCGACCACCTCCGGGGCCAg      >  1:63759/1‑149 (MQ=255)
                                                                                                                                gAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGACCATCGAAAATGCGCGCCACCGATTCGACCAGCGCCGCCAGGTCGACCACCTCCGGGGCCAgg     <  2:63759/148‑1 (MQ=255)
                                                                                                                                gAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGACCATCGAAAATGCGCGCCACCGATTCGACCAGCGCCGCCAGGTCGACCACCTACGGGGCCAgg     >  1:319326/1‑148 (MQ=255)
                                                                                                                                  aGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGACCATCGAAAATGCGCGCCACCGATTCGACCAGCGCCGCCAGGTCGACCACCTCCGGGGCCAGGGAc  <  2:319326/149‑1 (MQ=255)
                                                                                                                                   gcgcAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGTCGCGCCTGACCATCGAAAATGCGCGACACCGATTCGACCAGCGCCGCCAGGTCGACCACCTCCGGGGCCAGGGAc  >  4:264991/1‑148 (MQ=255)
                                                                                                                                      cAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGACCATCGAAAATGCGCGCCACCGATTCGACCAGCGCCGCCAGGTCGACCACCTCCGGGGCCAg      <  6:198959/141‑1 (MQ=255)
                                                                                                                                      cAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGACCATCGAAAATGCGCGCCACCGATTCGACCAGCGCCGCCAGGTCGACCACCTCCGGGGCCAg      >  5:198959/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                           
TCGCGCACCTCCAGCTCCAGTACGGCTGGTGTGCTGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGACCATCGAAAATGCGCGCCACCGATTCGACCAGCGCCGCCAGGTCGACCACCTCCGGGGCCAGGGAC  >  NC_002947/2397042‑2397320

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