Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,397,182 G→T 9.9% D825E (GAC→GAA)  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,1820GT 9.9% 136.6 / 4.8 47D825E (GAC→GAA) PP_2100two‑component system sensor protein
Reads supporting (aligned to +/- strand):  ref base G (24/18);  new base T (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.66e-01

TCGCGCACCTCCAGCTCCAGTACGGCTGGTGTGCTGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGACCATCGAAAATGCGCGCCACCGATTCGACCAGCGCCGCCAGGTCGACCACCTCCGGGGCCAGGGACAAGTGCCCG  >  NC_002947/2397042‑2397329
                                                                                                                                            |                                                                                                                                                   
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                                                                                                    <  5:71450/117‑1 (MQ=255)
                                                                         gACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAggg                                                                                                    >  6:71450/1‑117 (MQ=255)
                                                                                   cgTGAACATGATGGAATTGATCACCAGGTTTGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCCGCGTGCGGCCGGGCCGATCTTCACCTCAAGGGCCAGGCCTTTCTGCCGCGAATGACCATCGAAAATGCGAGCCa                                                          >  8:216321/2‑149 (MQ=255)
                                                                                     tGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGACCATCGAAAATGCGCGCCAcc                                                        >  3:163793/1‑149 (MQ=255)
                                                                                        aCTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGACCATCGAAAATGCGCGCCACCGAt                                                     >  4:147941/1‑149 (MQ=255)
                                                                                         cTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGAcc                                                                           <  6:63234/126‑1 (MQ=255)
                                                                                         cTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGAcc                                                                           >  5:63234/1‑126 (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)
                                                                                                                                            gTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGACCATCGAAAATGCGCGCCACCGATTCGACCAGCGCCGCCAGGTCGACCACCTCCGGGGCCAGGGACAAGTGCCCg  >  6:293513/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
TCGCGCACCTCCAGCTCCAGTACGGCTGGTGTGCTGGTTCCGTCATCCAGCAGCCTGACGCTGATGCGCACCTGACCGTGCTCGGTGAACTTGATGGCATTGCTCACCAGGTTGGACAGCACCTGCTTGAAGCGCAGCGGGTCAAGCAGGACATGGCAGCGTGCGGCCGGGGCGATCTGCACCTCCAGGGCCAGGCCTTTCTGCCGCGCCTGACCATCGAAAATGCGCGCCACCGATTCGACCAGCGCCGCCAGGTCGACCACCTCCGGGGCCAGGGACAAGTGCCCG  >  NC_002947/2397042‑2397329

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