Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,929,812 C→T 8.7% L120L (CTG→CTA) ‡ cheA ← chemotaxis histidine kinase CheA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,929,8120CT 8.7% 126.1 / 2.7 46L120L (CTG→CTA) ‡cheAchemotaxis histidine kinase CheA
Reads supporting (aligned to +/- strand):  ref base C (20/22);  new base T (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.36e-01

GCTCGAACTCACTGTCAGTGATGTCCTGGTGTTGCTCCTCGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTATCCAGTGCCTCGAGCACCACATCCATCAGTTCC‑GCATCGACCCGGCGCTCACCTTTGCGCAGGACGTCGAACACGTTC  >  NC_002947/4929670‑4929959
                                                                                                                                              |                                                                                                                                                    
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                                                                             <  3:97691/138‑1 (MQ=255)
                                                                              gTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTATCCa                                                                             >  4:97691/1‑138 (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)
                                                                                                                                              |                                                                                                                                                    
GCTCGAACTCACTGTCAGTGATGTCCTGGTGTTGCTCCTCGGCCTCGGCCGGCGCTGCAGCCACCACGTCGGCAACCGGTGCCGGTACGGGCGCTGCAACAGCTTCGGCTGCCCCAGGTTCGGCCAGGCGCGACAGCGCTGCCAGCAGCGCCGGGGTGGCCGGTGTCACTTCGGCACGCTCGCGCACCTGGCCGAACATGCTGTTGACCGTATCCAGTGCCTCGAGCACCACATCCATCAGTTCC‑GCATCGACCCGGCGCTCACCTTTGCGCAGGACGTCGAACACGTTC  >  NC_002947/4929670‑4929959

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