Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,493,440 T→C V324A (GTG→GCG)  trg → methyl‑accepting chemotaxis protein Trg

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,493,4400TC98.3% 204.8 / ‑4.4 59V324A (GTG→GCG) trgmethyl‑accepting chemotaxis protein Trg
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (28/30);  total (28/31)
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.72e-01

CGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGTGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGA  >  NC_000913/1493304‑1493571
                                                                                                                                        |                                                                                                                                   
cGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGaa                                                                                                                                  >  2:1002971/1‑140 (MQ=255)
  gcAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAAc                                                                                                                                <  2:250165/139‑1 (MQ=255)
  aCAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAAc                                                                                                                                >  1:1155155/1‑140 (MQ=255)
   cAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACa                                                                                                                               >  2:36950/1‑140 (MQ=255)
        cgcgGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAAt                                                                                                                           <  1:605319/139‑1 (MQ=255)
        cgcgGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAAt                                                                                                                           >  2:605319/1‑139 (MQ=255)
         gcgGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGc                                                                                                                         <  1:626841/140‑1 (MQ=255)
          cgGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCg                                                                                                                        <  1:1002971/140‑1 (MQ=255)
           ggAAGAGATTTATCGTGGCACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCgg                                                                                                                       <  2:1212071/140‑1 (MQ=255)
             aaGAGATTTATCGTGGCACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGAt                                                                                                                     <  1:485860/140‑1 (MQ=255)
                agaTTTATCGTGGCACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAAc                                                                                                                  >  1:238344/1‑140 (MQ=255)
                 gaTTTATCGTGGCACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAAcg                                                                                                                 <  1:884700/140‑1 (MQ=255)
                  aTTTATCGTGGCACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCTAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAAcgc                                                                                                                >  2:228260/1‑140 (MQ=255)
                     tATCGTGGCACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGca                                                                                                             <  1:564359/140‑1 (MQ=255)
                       tCGTGGCACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGACTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGcatc                                                                                                           >  2:76571/1‑140 (MQ=255)
                          tGGCACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGcat                                                                                                            >  2:54282/1‑136 (MQ=255)
                          tGGCACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGcat                                                                                                            <  1:54282/136‑1 (MQ=255)
                            gCACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGcc                                                                                                      >  2:1129594/1‑140 (MQ=255)
                             cACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCa                                                                                                     <  2:1155155/140‑1 (MQ=255)
                             cACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCa                                                                                                     <  1:1039277/140‑1 (MQ=255)
                                cAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCa                                                                                                  >  1:829327/1‑140 (MQ=255)
                                 aGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCaa                                                                                                 >  2:572238/1‑140 (MQ=255)
                                  gCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCaaa                                                                                                >  1:801759/1‑140 (MQ=255)
                                     aaaTTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCTAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTg                                                                                             >  2:228107/1‑140 (MQ=255)
                                         ttCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGgcgc                                                                                         <  1:325870/140‑1 (MQ=255)
                                           cAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCaa                                                                                       <  1:441032/140‑1 (MQ=255)
                                           cAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCaa                                                                                       <  1:36950/140‑1 (MQ=255)
                                           cAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCaa                                                                                       <  1:778818/140‑1 (MQ=255)
                                            aGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCTAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAag                                                                                      <  1:228260/140‑1 (MQ=255)
                                                ggCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGc                                                                                  <  2:721815/140‑1 (MQ=255)
                                                 gCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCt                                                                                 <  1:1129594/140‑1 (MQ=255)
                                                           ccTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCtt                                                                                >  2:1019046/1‑131 (MQ=255)
                                                           ccTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCtt                                                                                <  1:1019046/131‑1 (MQ=255)
                                                           ccTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAg                                                                       >  2:336739/1‑140 (MQ=255)
                                                               tCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTAtt                                                                           <  2:1080304/132‑1 (MQ=255)
                                                               tCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTAtt                                                                           >  1:1080304/1‑132 (MQ=255)
                                                               tCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAg                                                                   >  2:34957/1‑140 (MQ=255)
                                                                    tcGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATg                                                              >  2:163257/1‑140 (MQ=255)
                                                                      gCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGc                                                            >  1:722416/1‑140 (MQ=255)
                                                                         ccGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCggg                                                         >  1:229978/1‑140 (MQ=255)
                                                                         ccGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCggg                                                         >  1:876829/1‑140 (MQ=255)
                                                                            aagaaCAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAg                                                      >  1:869258/1‑140 (MQ=255)
                                                                            aagaaCAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAg                                                      <  1:412333/140‑1 (MQ=255)
                                                                            aagaaCAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAg                                                      <  1:34957/140‑1 (MQ=255)
                                                                             agaaCAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAGa                                                     <  1:165283/140‑1 (MQ=255)
                                                                              gaaCAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGTGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAGAc                                                    <  1:224409/140‑1 (MQ=255)
                                                                                 cAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAGACGGt                                                 <  2:876829/140‑1 (MQ=255)
                                                                                    gcggcgGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAGACGGTTTc                                              >  2:337129/1‑140 (MQ=255)
                                                                                    gcggcgGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAGACGGTTTc                                              <  2:801759/140‑1 (MQ=255)
                                                                                      ggcggcTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAGACGGTTTCCg                                            <  2:233528/140‑1 (MQ=255)
                                                                                         ggcTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAGACg                                                   <  1:904096/130‑1 (MQ=255)
                                                                                         ggcTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAGACg                                                   >  2:904096/1‑130 (MQ=255)
                                                                                                    aaaCTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACg                              <  1:337129/140‑1 (MQ=255)
                                                                                                     aaCTGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGa                             >  1:998917/1‑140 (MQ=255)
                                                                                                        tGCCGCCAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATgg                          <  2:229978/140‑1 (MQ=255)
                                                                                                          ccgccAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGgc                        >  1:719319/1‑140 (MQ=255)
                                                                                                           cgccAGCATGGAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGgcg                       <  1:163257/140‑1 (MQ=255)
                                                                                                                    ggAGCAACTCACTGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCa              >  1:28297/1‑140 (MQ=255)
                                                                                                                                tGCGACGGCGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCgaaga  >  1:828649/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                   
CGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGCTGGCAATGCGGACCTGTCATCTCGCACCGAAGAACAAGCGGCGGCTATCGAACAAACTGCCGCCAGCATGGAGCAACTCACTGCGACGGTGAAACAGAATGCGGATAACGCGCATCATGCCAGCAAACTGGCGCAAGAGGCTTCTATTAAAGCCAGCGATGGCGGGCAGACGGTTTCCGGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGA  >  NC_000913/1493304‑1493571

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