Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,967,611 (C)5→6 coding (891/1050 nt) cheB ← fused chemotaxis regulator: protein‑glutamate methylesterase in two‑component regulatory system with CheA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,967,6061.C81.8% 125.7 / 22.9 44coding (896/1050 nt)cheBfused chemotaxis regulator: protein‑glutamate methylesterase in two‑component regulatory system with CheA
Reads supporting (aligned to +/- strand):  ref base . (1/7);  new base C (17/19);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGCCGATGACGGTTAACCGCC  >  NC_000913/1967459‑1967749
                                                                                                                                                      |                                                                                                                                               
cGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGccc                                                                                                                                              <  2:370272/151‑1 (MQ=255)
cGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGccc                                                                                                                                              >  2:768492/1‑151 (MQ=255)
  tATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGccccc                                                                                                                                            >  2:55894/1‑151 (MQ=255)
   aTCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGcccccc                                                                                                                                           <  2:429871/151‑2 (MQ=255)
     cGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGccccccg                                                                                                                                          <  1:509277/150‑3 (MQ=255)
       ccTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGccccccgcct                                                                                                                                       >  1:787104/1‑146 (MQ=255)
       ccTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGccccccgcct                                                                                                                                       <  2:411791/151‑6 (MQ=255)
          gTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTg                                                                                                                                       <  1:123778/149‑1 (MQ=255)
          gTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGAc                                                                                                                                     <  2:678828/151‑1 (MQ=255)
            ccGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGc                                                                                                                                   >  2:616270/1‑151 (MQ=255)
               gCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATc                                                                                                                                >  2:592510/1‑151 (MQ=255)
                         ttGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTc                                                                                                                      >  1:792087/1‑151 (MQ=255)
                                cATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCcgccgc                                                                                                               >  1:713610/1‑151 (MQ=255)
                                 aTTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCcgccgcg                                                                                                              >  2:44601/1‑151 (MQ=255)
                                     gCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCCCCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCcgt                                                                                                          <  2:171512/151‑1 (MQ=255)
                                      cTGGCTTCCCTGGCTAAGATCGACCACTTCGCAGCCACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCGCGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCcgccgcgccg                                                                                                           <  2:341191/149‑1 (MQ=255)
                                         gCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGtt                                                                                                      >  1:174029/1‑151 (MQ=255)
                                           ttACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGcc                                                                                                   >  2:178996/1‑151 (MQ=255)
                                              ccTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCACGCGCGGCATGACGAACACCACGCAACTTGCTTCGTTTTGCGCAAGCGGCCATGCCCCCCGCCTGACGCATCGCCATCATTCCCGCCGCGCCGTCGTTGccc                                                                                                  >  2:432413/1‑150 (MQ=255)
                                                tGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCAt                                                                                                >  1:304907/1‑150 (MQ=255)
                                                 ggCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATAc                                                                                              >  2:278698/1‑151 (MQ=255)
                                                  gCTAAGATCGACCACTTCGCAGACACCCCCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATa                                                                                               <  1:634781/149‑1 (MQ=255)
                                                   caaaGATCGACCACCTCGCAGACACCCCCCATATGGATGGCCACGCGCGCATTGCCGAACACCTCTCACCTTGCTTCTTTTGGCGCTAGGTCCGCTACCCCCCGCCGTCCGCACCGCCAAAATTCCCGCCGCGCCGTCGTTGCCCATAcc                                                                                             <  1:768435/148‑1 (MQ=255)
                                                     aaGATCGACCACTTCGCAGACACCCCCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCg                                                                                            <  2:180870/149‑1 (MQ=255)
                                                      aGATCGCCCACTTCGCAGACCCCCCCCATATTGATGGCCTCGCGCGGCATGCCGAACCCCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTc                                                                                         <  1:393872/151‑1 (MQ=255)
                                                          cGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGa                                                                                     <  2:583380/151‑1 (MQ=255)
                                                             ccACTTCGCAGACCCCACCCATATTGATGGCCTCGCGCGGCATGCCGAACCCCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCa                                                                                  <  2:595053/151‑1 (MQ=255)
                                                                        acacCCCCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCAAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAAc                                                                           <  1:727284/147‑1 (MQ=255)
                                                                         ccccccccATATTGAGGGCCTCGCGCGGCATGCCGAACCCCACGCAACTTGCTTCGTTTTGCGCAGGGGTCCATGCCCCCCGCCTGACGCACCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGc                                                                        <  1:647783/144‑1 (MQ=255)
                                                                          accaccCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCAtt                                                                     >  2:25272/1‑151 (MQ=255)
                                                                           ccaccCATATTGATGGCCCCGCGCGGCAGGCCGAACACCACGCACCTTGTTTCGTTTTGCCCAGGGGTCCATGCCCCCCGCCTGACGCTTCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTa                                                                    <  2:44955/151‑1 (MQ=255)
                                                                              cccAAATTGGGGGCCTCGCGCGGCATGCCGAAACCCACGCAACTTCCTTTGTTTTCCGAAAGGGCCCAGGCCCCCCCCGGACGCACCCCCAACATCCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGcc                                                                <  1:849855/151‑1 (MQ=255)
                                                                                 atatTGATGGCCTCGCGCGGCATGCCGAACACCACCCACCTGGCTTCGTTTTGCGCAAGGTTCCATGCCCCCCGCCGGACGCACCGCCACCATACCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGccc                                                               <  1:262978/150‑1 (MQ=255)
                                                                                         ggCCCCGCGCGGCATGCCGAACCCCACCCACCTGGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATCCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGttt                                                      <  1:32412/151‑1 (MQ=255)
                                                                                           ccTCCCGCGCCATGCCGACCCCCCCGCACCTTGCTTCTTTTTGCGAAAGGTTCCATGCCCCCCCCCTACCGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGttt                                                      <  2:736881/149‑1 (MQ=255)
                                                                                            cTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGc                                                   >  2:277757/1‑151 (MQ=255)
                                                                                                  cgGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGAc                                                >  2:411818/1‑148 (MQ=255)
                                                                                                         ccGACCCCCCCGCAACTTGCTTCTTTTTGCGAAGGGGTCCATGCCCCCCGCCTGACCCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGaac                                      <  2:563382/151‑1 (MQ=255)
                                                                                                            aaCCCCCCGCAACTTGCTTCGTTTTGCCCAGGGGTCCATGCCCCCCGCCGGCCGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGaacaa                                    <  2:590846/150‑1 (MQ=255)
                                                                                                            aaCACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGaacaac                                   >  1:74160/1‑151 (MQ=255)
                                                                                                              caccccGAACCTTGCTTCGTTTTGCGCAAGGGCCCATCCCCCCCGCCTGCCGCATCGCCAACATCCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAacaaca                                  <  2:743294/150‑1 (MQ=255)
                                                                                                                ccccGCACCTTGCTTCGTTTTGCGCAGGGGTCCATGCCCCCCCCCTGACGCATCGCCACCATCCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAacac                                 <  2:222079/146‑1 (MQ=255)
                                                                                                                      aaCTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCg                         >  2:172683/1‑151 (MQ=255)
                                                                                                                      aaCTTCCTTCTTTTTGCGCAAGGCTCCATCCCCCCCGCCTGTCGCCCCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCACCCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTAc                           <  2:570043/149‑1 (MQ=255)
                                                                                                                            cTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGcc                   >  1:581671/1‑151 (MQ=255)
                                                                                                                             ttCGTTTTGCGCAAGGGTCCATCCCCCCCGCCTGACGCATCGCCAACATCCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGc                    <  2:246797/149‑1 (MQ=255)
                                                                                                                             ttCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGCCg                  >  1:658884/1‑151 (MQ=255)
                                                                                                                                 ttttGCGAAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGCCGAt                <  2:792143/149‑1 (MQ=255)
                                                                                                                                      cgcAGGGGTCCATCCCCCCCGCCGGCCGCACCGCCAACATCCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCACATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGCCGATGACg            <  1:194868/148‑1 (MQ=255)
                                                                                                                                           gggTCCATCCCCCCCCCCTGACGCACCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGCCGATGACGGTTAAccg    <  2:545537/151‑1 (MQ=255)
                                                                                                                                             gTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGCCGATGACGGTTAAccgcc  >  1:754001/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                               
CGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGCCGATGACGGTTAACCGCC  >  NC_000913/1967459‑1967749

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