Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,973,263 T→C 8.7% D26G (GAT→GGT)  cheW ← purine‑binding chemotaxis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,973,2630TC 8.7% 134.4 / 3.9 46D26G (GAT→GGT) cheWpurine‑binding chemotaxis protein
Reads supporting (aligned to +/- strand):  ref base T (21/21);  new base C (2/2);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAG  >  NC_000913/1973116‑1973408
                                                                                                                                                   |                                                                                                                                                 
aaCTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTtcat                                                                                                                                                   <  2:92952/148‑1 (MQ=255)
 aCTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGCGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTtcatc                                                                                                                                                  <  3:123983/148‑1 (MQ=255)
    tAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTTCTCTTCATCAc                                                                                                                                                >  2:181029/1‑147 (MQ=255)
       ttCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAg                                                                                                                                            <  4:253036/148‑1 (MQ=255)
        tCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTtcatc                                                                                                                                                  >  7:63883/1‑141 (MQ=255)
        tCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTtcatc                                                                                                                                                  <  8:63883/141‑1 (MQ=255)
             aGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTaa                                                                                                                                       <  3:90044/147‑1 (MQ=255)
             aGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTATCAGGATATCAATACCGTACTCTTCATCACCAAGGGTaa                                                                                                                                       <  8:96375/147‑1 (MQ=255)
                  acaaTCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCa                                                                                                                                 <  5:152533/148‑1 (MQ=255)
                        ggCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTtcatc                                                                                                                                                  >  2:325442/1‑125 (MQ=255)
                        ggCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTtcatc                                                                                                                                                  <  1:325442/125‑1 (MQ=255)
                        ggCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTtcatc                                                                                                                                                  <  3:124617/125‑1 (MQ=255)
                        ggCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTtcatc                                                                                                                                                  >  4:124617/1‑125 (MQ=255)
                        ggCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAAtt                                                                                                                           <  2:195032/148‑1 (MQ=255)
                        ggCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAAtt                                                                                                                           <  4:93209/148‑1 (MQ=255)
                          cacaATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTcc                                                                                                                         <  6:351360/148‑1 (MQ=255)
                             aaTAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTgg                                                                                                                      >  3:95672/1‑148 (MQ=255)
                                            aTTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTg                                                                                                       >  7:47280/1‑148 (MQ=255)
                                            aTTCGTGACGCCCTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCAcc                                                                                                                                               >  1:93690/1‑108 (MQ=255)
                                            aTTCGTGACGCCCTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCAcc                                                                                                                                               <  2:93690/108‑1 (MQ=255)
                                                         ttGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCt                                                                                                                   <  7:50864/123‑1 (MQ=255)
                                                         ttGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCt                                                                                                                   >  8:50864/1‑123 (MQ=255)
                                                         ttGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGtt                                                                                           >  2:251281/1‑147 (MQ=255)
                                                         ttGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTa                                                                                          >  4:193944/1‑148 (MQ=255)
                                                           gATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTAc                                                                                         <  8:47280/147‑1 (MQ=255)
                                                                 cGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTc                                                                                  <  1:251281/148‑1 (MQ=255)
                                                                 cGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTc                                                                                  <  3:193944/148‑1 (MQ=255)
                                                                          gtTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTa                                                                                          >  4:137881/1‑131 (MQ=255)
                                                                          gtTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTa                                                                                          <  3:137881/131‑1 (MQ=255)
                                                                          gtTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGt                                                                          >  7:374531/1‑147 (MQ=255)
                                                                             cGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCat                                                                       <  4:1691/147‑1 (MQ=255)
                                                                                  tCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGt                                                                  >  4:286771/1‑147 (MQ=255)
                                                                                   ccGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTa                                                                >  6:150242/1‑148 (MQ=255)
                                                                                      tgtTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGTTCATATTGTTACCt                                                             >  5:88704/1‑148 (MQ=255)
                                                                                        ttACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCtt                                                            <  3:286771/147‑1 (MQ=255)
                                                                                            cTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCACCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCttttt                                                         >  5:66066/1‑146 (MQ=255)
                                                                                            cTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCACCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCttttt                                                         <  6:66066/146‑1 (MQ=255)
                                                                                              gATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTg                                                                                                                  <  6:147787/87‑1 (MQ=255)
                                                                                              gATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTg                                                                                                                  >  5:147787/1‑87 (MQ=255)
                                                                                                tCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCAt                                                   >  1:109548/1‑148 (MQ=255)
                                                                                                         cGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAggcggc                                          <  8:374531/148‑1 (MQ=255)
                                                                                                                    cACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGcc                                <  4:80527/147‑1 (MQ=255)
                                                                                                                     aCTTTCAGGATATCAATACCGTACTCTTCACCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCa                               >  1:189346/1‑147 (MQ=255)
                                                                                                                       tttCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCGGGCGGCGGTGTTCGCCATAc                            >  4:241567/1‑148 (MQ=255)
                                                                                                                       tttCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATAc                            >  7:19875/1‑148 (MQ=255)
                                                                                                                        ttCAGGATATCAATACCGTACTCTTCATCACAAAGGGTAAATACCATAAATTCCTGTCCTGAAGGCTCGCTGGCAAGCTTTGTTACATTCGTCATACAGGTCATATTGTTACCTTTTTACTCATTCATCCGGCGgtg                                      >  3:268970/1‑137 (MQ=255)
                                                                                                                           aGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACgttg                        >  7:244838/1‑148 (MQ=255)
                                                                                                                                 tCAATACCGTACTCGTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTcgcg                  <  5:150242/148‑1 (MQ=255)
                                                                                                                                   aaTTCCGTACTCTTCATCACGATGGGTAAATACCAGAAATTCCTTGCCTTACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATTTTGTTACCTTTTTACTCATTCAGGCGGCGgt                                       <  3:142146/125‑1 (MQ=255)
                                                                                                                                   aaTTCCGTACTCTTCATCACCAAGGTTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCgg                 <  1:181029/147‑1 (MQ=255)
                                                                                                                                      aCCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTa             >  4:8084/1‑148 (MQ=255)
                                                                                                                                       ccGTACTCTTCACCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGGTCGCGGTTTAt            <  2:189346/148‑1 (MQ=255)
                                                                                                                                         gTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCg          >  5:275327/1‑148 (MQ=255)
                                                                                                                                           aCTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGcc        >  5:99993/1‑148 (MQ=255)
                                                                                                                                                tcatcaCCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTc                                                 <  4:294952/102‑1 (MQ=255)
                                                                                                                                                tcatcaCCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTc                                                 >  3:294952/1‑102 (MQ=255)
                                                                                                                                                 catcaCCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGc     <  6:88704/145‑1 (MQ=255)
                                                                                                                                                 catcaCCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAg  >  3:313530/1‑148 (MQ=255)
                                                                                                                                                 catcaCCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAg  >  2:260729/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAG  >  NC_000913/1973116‑1973408

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: