Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,622,190 G→C 13.4% E368D (GAG→GAC) ‡ 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,1900GC13.4% 123.4 / 5.7 46E368D (GAG→GAC) ‡PP_1421sensor histidine kinase TctE
Reads supporting (aligned to +/- strand):  ref base G (16/23);  new base C (2/3);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

GGCGGTGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGT  >  NC_002947/1622044‑1622338
                                                                                                                                                  |                                                                                                                                                    
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                                                          >  2:64662/1‑104 (MQ=255)
                                                                                                                                       ttGCTCAACAACTCGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagt                                                          <  1:64662/104‑1 (MQ=255)
                                                                                                                                         gCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGgtgt           >  5:341743/1‑149 (MQ=255)
                                                                                                                                             aaCGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTc         >  3:241855/1‑147 (MQ=255)
                                                                                                                                             aaCGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTc         >  7:107329/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)
                                                                                                                                                  gttgttGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGgtgt           >  6:23795/1‑140 (MQ=255)
                                                                                                                                                  gttgttGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGgtgt           <  5:23795/140‑1 (MQ=255)
                                                                                                                                                  gttgttGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCAAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGt  >  5:278441/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGCGGTGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGT  >  NC_002947/1622044‑1622338

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