Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,622,431 G→A 9.0% G449S (GGC→AGC)  PP_1421 → sensor histidine kinase TctE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,622,4310GA 9.0% 121.1 / 2.1 44G449S (GGC→AGC) PP_1421sensor histidine kinase TctE
Reads supporting (aligned to +/- strand):  ref base G (19/21);  new base A (2/2);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGT  >  NC_002947/1622292‑1622574
                                                                                                                                           |                                                                                                                                                    
cgGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAgg                                                                                                                                             <  3:147901/149‑1 (MQ=255)
         tATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTtgc                                                                                                                                    <  6:40526/149‑1 (MQ=255)
         tATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTtgc                                                                                                                                    <  6:45818/149‑1 (MQ=255)
             cccGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTtgcgggt                                                                                                                                >  3:211436/1‑149 (MQ=255)
             cccGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTtgcgggt                                                                                                                                >  3:173832/1‑149 (MQ=255)
             cccGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTtgcggg                                                                                                                                 <  4:50559/148‑1 (MQ=255)
              ccGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTtgcgggtg                                                                                                                               >  8:119156/1‑149 (MQ=255)
               cGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTtgcgggtgc                                                                                                                              >  5:257912/1‑149 (MQ=255)
                   gCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTgcgggtgcgggt                                                                                                                          <  1:269446/149‑1 (MQ=255)
                        gcgcGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCt                                                                                                                     >  2:205615/1‑149 (MQ=255)
                                gggTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATcgcg                                                                                                             >  7:32365/1‑149 (MQ=255)
                                 ggTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATcgcg                                                                                                             >  7:192621/1‑148 (MQ=255)
                                  gtgtTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGa                                                                                                           <  8:192621/149‑1 (MQ=255)
                                    gttggcgCGGTTCTACCGGTGTAGGGGGCAGGGCTGTGGGTTATGGCTGGCGTTGGTTGGCGATTTCTGTCGGGCGCACCTGGTGTAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGGGAGCTTTATCGCGGAtt                                                                                                         <  8:53055/143‑1 (MQ=255)
                                        gAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATg                                                                                                     >  1:320149/1‑149 (MQ=255)
                                                 tACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCtg                                                                                            >  3:271229/1‑149 (MQ=255)
                                                   ccGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCt                                                                                             <  4:7629/146‑1 (MQ=255)
                                                   ccGGCGAAGTGCGCAGGGCAGTGGGTTAGGACTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCtgtg                                                                                          >  1:208564/1‑149 (MQ=255)
                                                        gAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGcc                                                                                     <  8:32365/149‑1 (MQ=255)
                                                              cgcAGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTc                                                                          <  2:208564/149‑1 (MQ=255)
                                                                 aGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATg                                                                       >  7:77843/1‑149 (MQ=255)
                                                                   ggCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATggg                                                                     <  6:257912/149‑1 (MQ=255)
                                                                     cAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCt                                                                   <  4:271229/149‑1 (MQ=255)
                                                                          ggTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCTCAccc                                                              >  6:265743/1‑149 (MQ=255)
                                                                           gTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCTCACCcg                                                             >  1:172727/1‑149 (MQ=255)
                                                                            ttAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCTCACCcgc                                                            >  7:296214/1‑149 (MQ=255)
                                                                            ttAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCTCACCcgc                                                            <  2:320149/149‑1 (MQ=255)
                                                                                    ggCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCTCACCCGCGAAAgggg                                                    <  1:205615/149‑1 (MQ=255)
                                                                                           gTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGc                                                                    <  4:135006/126‑1 (MQ=255)
                                                                                           gTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGc                                                                    >  3:135006/1‑126 (MQ=255)
                                                                                                 gagaTCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTgg                                       <  5:265743/149‑1 (MQ=255)
                                                                                                    aTCTGTCGGGCGCACGTGGGGCAAGTGACCCTGCACGTTGGCGAGAAGGGCGGGTTGGGGGTGCGGGTGAGGTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTa                                     <  2:172727/148‑1 (MQ=255)
                                                                                                         tCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCTCACCCGCGAAAgg                                                      >  2:161311/1‑126 (MQ=255)
                                                                                                         tCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCTCACCCGCGAAAgg                                                      <  1:161311/126‑1 (MQ=255)
                                                                                                         tCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGc                                 <  7:119156/147‑1 (MQ=255)
                                                                                                                  ccTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACCGCTCAGCGCAACa                      >  5:15972/1‑149 (MQ=255)
                                                                                                                  ccTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACa                      >  3:108008/1‑149 (MQ=255)
                                                                                                                              gACCCTGCACGATAGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCt          >  2:61464/1‑149 (MQ=255)
                                                                                                                              gACCCTGCACGATAGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCt          >  2:56187/1‑149 (MQ=255)
                                                                                                                               aCCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCt          <  4:173832/148‑1 (MQ=255)
                                                                                                                               aCCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTc         <  4:211436/149‑1 (MQ=255)
                                                                                                                                 ccTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCGGGGGGTTTAATGGGCTAACACGCAAAATGGGCGAAAAAGGGTGGTACAGCTCATCGCAACATGCTCcgccc            >  1:112090/1‑147 (MQ=255)
                                                                                                                                      aCGATAGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGt  <  1:61464/149‑1 (MQ=255)
                                                                                                                                      aCGATAGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGt  <  1:56187/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
CGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCC‑‑‑‑‑CTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGT  >  NC_002947/1622292‑1622574

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