Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,929,813 A→G 8.7% L120P (CTG→CCG) ‡ cheA ← chemotaxis histidine kinase CheA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,929,8130AG 8.7% 123.1 / 2.2 46L120P (CTG→CCG) ‡cheAchemotaxis histidine kinase CheA
Reads supporting (aligned to +/- strand):  ref base A (20/22);  new base G (2/2);  total (22/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 = 9.63e-01

GCTCGAACTCACTGTCAGTGATGTCCTGGTGTTGCTCCTCGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTATCCAGTGCCTCGAGCACCACATCCATCAGTTCC‑GCATCGACCCGGCGCTCACCTTTGCGCAGGACGTCGAACACGTTCTC  >  NC_002947/4929670‑4929961
                                                                                                                                               |                                                                                                                                                     
gCTCGAACTCACTGTCAGTGATGTCCTGGTGTTGCACCTCGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTACAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCagcagc                                                                                                                                                  >  2:184043/1‑149 (MQ=255)
    gAACTCACTGTCAGTGATGTCCTGGTGTTGCTCCTCGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACCGCGCTGCCAGCAGCGCCg                                                                                                                                              >  6:211443/1‑149 (MQ=255)
      aCTCACTGTCAGTGATGTCCTGGTGTTGCTCCTCGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCggg                                                                                                                                            <  5:183905/149‑1 (MQ=255)
           cTGTCAGTGATGTCCTGGTGTTGCTCCTCGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGc                                                                                                                                       <  2:287592/149‑1 (MQ=255)
               cAGTGATGTCCTGGTGTTGCTCCTCGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGgt                                                                                                                                   <  6:103521/149‑1 (MQ=255)
                aGTGATGTCCTGGTGTTGCTCCTCGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCcagc                                                                                                                                                     <  5:85808/130‑1 (MQ=255)
                aGTGATGTCCTGGTGTTGCTCCTCGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCcagc                                                                                                                                                     >  6:85808/1‑130 (MQ=255)
                aGTGATGTCCTGGTGTTGCTCCTCGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGcc                                                                                                                                               <  1:195019/136‑1 (MQ=255)
                aGTGATGTCCTGGTGTTGCTCCTCGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGcc                                                                                                                                               >  2:195019/1‑136 (MQ=255)
                      gTCCTGGTGTTGCTCCTCGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACtt                                                                                                                            >  5:195793/1‑149 (MQ=255)
                         cTGGTGTTGCTCCTCGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGcc                                                                                                                                               >  8:42906/1‑127 (MQ=255)
                         cTGGTGTTGCTCCTCGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGcc                                                                                                                                               <  7:42906/127‑1 (MQ=255)
                         cTGGTGTTGCTCCTCGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCgg                                                                                                                         >  4:204535/1‑149 (MQ=255)
                            gtgtTGCTCCTCGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGc                                                                                                                        <  5:16677/147‑1 (MQ=255)
                                 gCTCCTCGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACGGCGCTGCTGGCAgcgctgtcgc                                                                                                                                          <  4:190398/124‑7 (MQ=255)
                                 gCTCCTCGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACGGCGCTGCTGGCAgcgctgtcgc                                                                                                                                          >  3:190398/1‑118 (MQ=255)
                                   tcctcGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTcgcg                                                                                                               >  6:54576/1‑149 (MQ=255)
                                    cctcggcctcggCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGCTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTc                                                                                                                                <  4:150443/131‑1 (MQ=255)
                                     ctcggcctcggcCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCa                                                                                                             >  3:296729/1‑149 (MQ=255)
                                     ctcggcctcggcCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCa                                                                                                             >  6:256688/1‑149 (MQ=255)
                                        ggcctcggccGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCt                                                                                                          <  3:204535/149‑1 (MQ=255)
                                         gcctcggccGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTg                                                                                                         >  2:210057/1‑149 (MQ=255)
                                          cctcggccGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGACAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTgg                                                                                                        >  7:51339/1‑149 (MQ=255)
                                           ctcggccGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGc                                                                                                       <  8:51339/149‑1 (MQ=255)
                                           ctcggccGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGc                                                                                                       <  3:233602/149‑1 (MQ=255)
                                                   gcgcTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATg                                                                                               <  5:54576/149‑1 (MQ=255)
                                                   gcgcTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATg                                                                                               >  4:256817/1‑149 (MQ=255)
                                                         cAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGtt                                                                                          >  3:197958/1‑148 (MQ=255)
                                                           gCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGAc                                                                                       <  3:256817/149‑1 (MQ=255)
                                                            ccaccaCGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGAcc                                                                                      >  6:241369/1‑149 (MQ=255)
                                                            ccaccaCGTCGGCAACCGGTGCCGGTACGGGCGCTACAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGAcc                                                                                      <  1:184043/149‑1 (MQ=255)
                                                              accacGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCg                                                                                     >  4:276126/1‑148 (MQ=255)
                                                               ccacGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTcgc                                                                                                                <  7:246818/120‑1 (MQ=255)
                                                               ccacGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTcgc                                                                                                                >  8:246818/1‑120 (MQ=255)
                                                                cacGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTAt                                                                                  <  4:296729/149‑1 (MQ=255)
                                                                    tCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTATCCAg                                                                              <  8:263418/149‑1 (MQ=255)
                                                                              gTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTATCCa                                                                               >  4:97691/1‑138 (MQ=255)
                                                                              gTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTATCCa                                                                               <  3:97691/138‑1 (MQ=255)
                                                                                                   cAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTATCCAGTGCCTCGAGCACCACATCCATCAGTTCC‑GCa                                              >  2:183992/1‑149 (MQ=255)
                                                                                                               ccccAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTATCCAGTGCCTCGAGCACCACATCCATCAGTTCC‑GCATCGACCCGgcgc                                  >  5:337768/1‑149 (MQ=255)
                                                                                                                          gCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTATCCAGTGCCTCGAGCACCACATCGATCAGTTCC‑GCATCGACCCGGCGCTCACCTTTgcg                       >  6:228676/1‑149 (MQ=255)
                                                                                                                              ggCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTATCCAGTGCCTCGAGCACCACATCCATCAGTTCCAGCATCGACCCGGCGCTCACCTTTGCGCAg                    <  1:210057/149‑1 (MQ=255)
                                                                                                                              ggCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTATCCAGTGCCTCGAGCACCACATCCATCAGTTCC‑GCATCGACCCGGCGCTCACCTTTGCGCAgg                   <  5:241369/149‑1 (MQ=255)
                                                                                                                                 gcgACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTATCCAGTGCCTCGAGCACCACATCCATCAGTTCC‑GCATCGACCCGGCGCTCACCTTTGCGCAGGACg                <  6:195793/149‑1 (MQ=255)
                                                                                                                                 accccggcgcTGCTGGCAGCGCTGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTATCCAGTGCCTCGAGCACCACATCCATCAGTTCC‑GCATCGa                                          >  4:353206/7‑123 (MQ=255)
                                                                                                                                 accccggcgcTGCTGGCAGCGCTGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTATCCAGTGCCTCGAGCACCACATCCATCAGTTCC‑GCATCGa                                          <  3:353206/117‑1 (MQ=255)
                                                                                                                                   gACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTATCCAGTGCCTCGAGCACCACATCCATCAGTTCC‑GCATCGACCCGGCGCTCACCTTTGCGCAGGACGTc              <  5:256688/149‑1 (MQ=255)
                                                                                                                                   gACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTATCCAGTGCCTCGAGCACCACATCCATCAGTTCC‑GCATCGACCCGGCGCTCACCTTTGCGCAGGACGTc              <  1:183992/149‑1 (MQ=255)
                                                                                                                                            gacagcagCGCCGGGGTGGCCGGTGTCACTTAGTCACGCTCGCGCACCTGGCCGAACATGCTCTTGACCGTATCCAGTGCCTCGAGCACCCCATACATCAGTTCC‑GCATCGACCCCGCGATCACCTTTTCGCAGGACGTCGAACACGt      >  7:273587/3‑148 (MQ=255)
                                                                                                                                             ccAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTATCCAGTGCCTCGAGCACCACATCCATCAGTTCC‑GCATCGACCCGGCGCTCACCTTTGCGCAGGACGTCGAACACGTtc    >  8:272409/1‑149 (MQ=255)
                                                                                                                                               agcagcGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTATCCAGTGCCTCGAGCACCACATCCATCAGTTCC‑GCATCGACCCGGCGCTCACCTTTGCGCAGGATGTCGAACACGTtctc  <  5:211443/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GCTCGAACTCACTGTCAGTGATGTCCTGGTGTTGCTCCTCGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTATCCAGTGCCTCGAGCACCACATCCATCAGTTCC‑GCATCGACCCGGCGCTCACCTTTGCGCAGGACGTCGAACACGTTCTC  >  NC_002947/4929670‑4929961

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