Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,622,188 G→A 12.0% E368K (GAG→AAG) ‡ 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,1880GA12.0% 120.8 / 6.4 51E368K (GAG→AAG) ‡PP_1421sensor histidine kinase TctE
Reads supporting (aligned to +/- strand):  ref base G (17/26);  new base A (3/3);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

GGCGGTGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGC  >  NC_002947/1622044‑1622335
                                                                                                                                                |                                                                                                                                                   
ggCGGTGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAgtt                                                                                                                                                 <  6:344045/149‑1 (MQ=255)
 gCGGTGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAgttg                                                                                                                                                <  1:51076/149‑1 (MQ=255)
   ggTGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAgttgtt                                                                                                                                              <  7:173111/149‑1 (MQ=255)
    gTGCGCAGCGGCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGttgttg                                                                                                                                             <  2:351884/149‑1 (MQ=255)
     tGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGa                                                                                                                                            <  6:355578/149‑1 (MQ=255)
        gcAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAgttg                                                                                                                                                <  1:112887/142‑1 (MQ=255)
        gcAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAgttg                                                                                                                                                >  2:112887/1‑142 (MQ=255)
        gcAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACAACTCGTTGAGCa                                                                                                                                         <  2:164858/149‑1 (MQ=255)
            cGTCTGGCCTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGACCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACAACTCGTTGAGCaatgt                                                                                                                                     <  4:236175/149‑4 (MQ=255)
            cGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCGTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCt                                                                                                                                     <  4:264300/149‑1 (MQ=255)
             gTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCtg                                                                                                                                    >  4:86439/1‑149 (MQ=255)
                   aCTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCATTAGCGCGGGGTGGCGCTGCCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGAc                                                                                                                              <  1:115315/149‑1 (MQ=255)
                   aCTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGAc                                                                                                                              <  5:285235/149‑1 (MQ=255)
                       gAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAAt                                                                                                                           <  1:135309/148‑1 (MQ=255)
                             gTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTAGATCAAGGTGGACCCCAAATAAATAAAAAAAATATTGAGCAACCTGGTGGACAATGCCTTgg                                                                                                                    >  4:167687/1‑149 (MQ=255)
                               tGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGgc                                                                                                                   <  3:129419/148‑1 (MQ=255)
                                 gCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCa                                                                                                                >  2:278560/1‑149 (MQ=255)
                                   ccGCGAGCTGGGCAGGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATa                                                                                                              <  6:224861/149‑1 (MQ=255)
                                        aGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCa                                                                                                         <  3:167687/149‑1 (MQ=255)
                                          cTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGc                                                                                                       <  1:278560/149‑1 (MQ=255)
                                             ggCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGTATGCGCCaaggca                                                                                                    <  5:344998/149‑6 (MQ=255)
                                              gCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTAGATCAAGGTGGAGCCGAAATACATAAAATAAATGTTGAGCAACCTGGTGGACAATGCCTTGGAGCATACGCCAGAGGGt                                                                                                   >  6:158985/1‑149 (MQ=255)
                                               catggccatggcGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTg                                                                                                  <  3:86439/149‑1 (MQ=255)
                                                    ccatggcGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAAt                                                                                             >  5:63103/1‑149 (MQ=255)
                                                         gcgcCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGa                                                                                         <  4:133546/148‑1 (MQ=255)
                                                                aggggcACAAGGGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCgg                                                                                 <  3:301058/145‑1 (MQ=255)
                                                                    gcACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTaa                                                                             >  7:275367/1‑149 (MQ=255)
                                                                     cacaAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAAt                                                                            <  7:247305/149‑1 (MQ=255)
                                                                     cacaAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAAt                                                                            <  7:294329/149‑1 (MQ=255)
                                                                                     cgctggcgctggAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtg                                                            <  6:63103/149‑1 (MQ=255)
                                                                                      gctggcgctggAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtg                                                            >  6:247256/1‑148 (MQ=255)
                                                                                         ggcgctggAGGCCGAGGCGCCGGTGTGGATCAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaag                                                        >  6:173304/1‑149 (MQ=255)
                                                                                          gcgctggAGGCCGAGGCGCCGGTGTGGTTTAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagt                                                       <  6:81943/149‑1 (MQ=255)
                                                                                            gctggAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagtg                                                      >  3:84288/1‑148 (MQ=255)
                                                                                                     cGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGgc                                                                                                                   <  5:61579/78‑1 (MQ=255)
                                                                                                     cGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGgc                                                                                                                   >  6:61579/1‑78 (MQ=255)
                                                                                                                  tgGTTGAAGGGGGAGCCGACATTGCTCAACTAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTTATTTTGCGGGTAATGGCGCCGGCAGTGCTTGAAGTGGAAGACGACGGGCCAGGTATTc                                <  6:328492/148‑1 (MQ=255)
                                                                                                                    gTTGAAGGGGAAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCg                             >  2:361388/1‑149 (MQ=255)
                                                                                                                        aaGGGGAATCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTgg                          <  1:361388/148‑1 (MQ=255)
                                                                                                                                 accaggttGCTCAACAACTCGTTGAGCAATCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGcgcg                 >  5:33271/7‑148 (MQ=255)
                                                                                                                                   gACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAgc                    <  4:243600/143‑1 (MQ=255)
                                                                                                                                   gACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAgc                    >  3:243600/1‑143 (MQ=255)
                                                                                                                                       ttGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGgt          <  5:13360/149‑1 (MQ=255)
                                                                                                                                       ttGCTCAACAACTCGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagt                                                       <  1:64662/104‑1 (MQ=255)
                                                                                                                                       ttGCTCAACAACTCGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagt                                                       >  2:64662/1‑104 (MQ=255)
                                                                                                                                         gCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGgtgt        >  5:341743/1‑149 (MQ=255)
                                                                                                                                             aaCGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTc      >  7:107329/1‑147 (MQ=255)
                                                                                                                                             aaCGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTc      >  3:241855/1‑147 (MQ=255)
                                                                                                                                             aaCGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGAGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTc      >  7:215659/1‑147 (MQ=255)
                                                                                                                                               cGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGc  >  6:152044/1‑149 (MQ=255)
                                                                                                                                               cGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGc  >  6:146442/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GGCGGTGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGC  >  NC_002947/1622044‑1622335

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