Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,872,133 T→C 7.5% Q83R (CAA→CGA)  PP_3420 ← sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,872,1330TC 7.5% 151.4 / 2.7 53Q83R (CAA→CGA) PP_3420sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (27/22);  new base C (2/2);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

CATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCAACAGGCTGATCA  >  NC_002947/3871987‑3872280
                                                                                                                                                   |                                                                                                                                                   
cATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGc                                                                                                                                                   <  3:148125/149‑1 (MQ=255)
 aTTGGCCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGc                                                                                                                                                   >  6:21554/1‑148 (MQ=255)
                   cGGCAGGTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCt                                                                                                                                >  3:3824/1‑149 (MQ=255)
                    ggCAGGTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGTCCGGTGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTg                                                                                                                               <  4:223264/149‑1 (MQ=255)
                     gCAGGTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTgg                                                                                                                              <  8:53023/149‑1 (MQ=255)
                       aGGTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGc                                                                                                                            <  2:28122/149‑1 (MQ=255)
                         gTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCa                                                                                                                          >  2:136271/1‑149 (MQ=255)
                             tCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGc                                                                                                                      >  5:223022/1‑149 (MQ=255)
                               gACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTgg                                                                                                                              >  4:243827/1‑139 (MQ=255)
                               gACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTgg                                                                                                                              <  3:243827/139‑1 (MQ=255)
                                      tCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGcc                                                                                                                           <  6:89208/135‑1 (MQ=255)
                                      tCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGcc                                                                                                                           >  5:89208/1‑135 (MQ=255)
                                         aGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTg                                                                                                          >  7:191982/1‑149 (MQ=255)
                                             aCTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCATTCGCAATTGCCGc                                                                                                                                          <  4:46728/113‑1 (MQ=255)
                                             aCTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCATTCGCAATTGCCGc                                                                                                                                          >  3:46728/1‑113 (MQ=255)
                                                 gATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGt                                                                                                    >  8:235222/1‑147 (MQ=255)
                                                 gATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTaa                                                                                                  <  1:136271/149‑1 (MQ=255)
                                                  aTCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGCACTGCGCGGTAAg                                                                                                 <  4:3824/149‑1 (MQ=255)
                                                      gggTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATg                                                                                             >  2:40364/1‑149 (MQ=255)
                                                      gggTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATg                                                                                             <  7:235222/149‑1 (MQ=255)
                                                      gggTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATg                                                                                             <  1:43887/149‑1 (MQ=255)
                                                                    aGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCAc                                                                               >  3:55915/1‑149 (MQ=255)
                                                                                aCCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGca                                                                    <  5:21554/148‑1 (MQ=255)
                                                                                  cTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGcacca                                                                 <  6:223022/149‑1 (MQ=255)
                                                                                  cTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGcacca                                                                 >  5:101440/1‑149 (MQ=255)
                                                                                  cTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGcacca                                                                 >  2:243033/1‑149 (MQ=255)
                                                                                    gCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTcgc                                                                       <  6:104731/141‑1 (MQ=255)
                                                                                    gCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTcgc                                                                       >  5:104731/1‑141 (MQ=255)
                                                                                     cAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGcaccac                                                                <  4:55915/147‑1 (MQ=255)
                                                                                              cAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTc                                                      <  1:243033/148‑1 (MQ=255)
                                                                                                ggTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGa                                                    >  3:80553/1‑148 (MQ=255)
                                                                                                 gTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGa                                                    >  1:44361/1‑148 (MQ=255)
                                                                                                      gCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTaa                                             >  1:129129/1‑149 (MQ=255)
                                                                                                       ccGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAAc                                            <  2:129129/149‑1 (MQ=255)
                                                                                                       ccGCGGGTCCAGAAACTCGCCCATGGTGGTGCCAGGCGGCAATTTAGAATGCAGCACTTAGGCCTGGGCCATTGCATGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCAAGAAGCGCGTCGAAACGTAAc                                            >  2:241287/1‑149 (MQ=255)
                                                                                                         gcgGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGACCTGCGCGGTAAGCATGGCATCGCGCTCCCCGCTGTCGCGCACCACGACGCGCGTCGAAACCTCAcc                                           >  5:195819/1‑148 (MQ=255)
                                                                                                                       tCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCAc                            >  8:75291/1‑149 (MQ=255)
                                                                                                                         gCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATg                                                                                             <  5:34854/82‑1 (MQ=255)
                                                                                                                         gCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATg                                                                                             >  6:34854/1‑82 (MQ=255)
                                                                                                                           ccATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAg                        >  2:218541/1‑149 (MQ=255)
                                                                                                                            cATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAgg                       >  3:20216/1‑149 (MQ=255)
                                                                                                                            cATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAgg                       >  6:99768/1‑149 (MQ=255)
                                                                                                                               ggtggtGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTaa                    >  7:47064/1‑149 (MQ=255)
                                                                                                                                 tggtgCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAAt                   >  8:221884/1‑148 (MQ=255)
                                                                                                                                 tggtgCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATg                  <  4:20216/149‑1 (MQ=255)
                                                                                                                                     tgCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATc               <  7:75291/149‑1 (MQ=255)
                                                                                                                                     tgCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATc               <  2:44361/149‑1 (MQ=255)
                                                                                                                                       ccAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCaa             >  2:177664/1‑149 (MQ=255)
                                                                                                                                       aagtgcggcATTCGCAATTGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGt                                                       <  5:199786/103‑1 (MQ=255)
                                                                                                                                       aagtgcggcATTCGCAATTGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGt                                                       >  6:199786/5‑107 (MQ=255)
                                                                                                                                            cggcAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTg                                                                           <  6:252320/82‑1 (MQ=255)
                                                                                                                                            cggcAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTg                                                                           >  5:252320/1‑82 (MQ=255)
                                                                                                                                             ggcAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCAACAGGCt       <  5:99768/149‑1 (MQ=255)
                                                                                                                                                  ttGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCAACAGGCTGATCa  >  6:227350/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCAACAGGCTGATCA  >  NC_002947/3871987‑3872280

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: