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.C92.3% 123.4 / 4.0 40coding (896/1050 nt)cheBfused chemotaxis regulator: protein‑glutamate methylesterase in two‑component regulatory system with CheA
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base C (19/17);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGCCGATGACGGT  >  NC_000913/1967458‑1967741
                                                                                                                                                      |                                                                                                                                        
aCGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGTGTCCATGc                                                                                                                                         <  1:908385/150‑1 (MQ=255)
aCGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGcc                                                                                                                                        >  2:259531/1‑151 (MQ=255)
aCGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGcc                                                                                                                                        >  1:870162/1‑151 (MQ=255)
    aTCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCGGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGcccccc                                                                                                                                   <  1:902585/151‑2 (MQ=255)
     tCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGcccccc                                                                                                                                   >  1:429623/1‑149 (MQ=255)
        ccTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGccccccgcct                                                                                                                               >  1:607887/1‑146 (MQ=255)
           gTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGa                                                                                                                              <  1:239495/150‑1 (MQ=255)
           gTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGAc                                                                                                                             >  1:478225/1‑151 (MQ=255)
              cGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGc                                                                                                                           >  1:450662/1‑150 (MQ=255)
                gCACTAATCTTTGCCAACATTTGCTGGCTTACCTGGCTAAGTTCGATCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCCAACCCCACGCACCCTTCTTCGCTTTGCGCACGGGTCCATGCCCCCCGCCTGCCGCa                                                                                                                          >  2:387518/1‑149 (MQ=255)
                  aCTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGc                                                                                                                      >  1:106855/1‑151 (MQ=255)
                     aaTTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGcc                                                                                                                     >  1:510603/1‑149 (MQ=255)
                            gCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTccc                                                                                                            <  1:394245/151‑1 (MQ=255)
                            gCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTcc                                                                                                             <  1:56127/150‑1 (MQ=255)
                              cAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCccgc                                                                                                          <  1:544650/151‑1 (MQ=255)
                                          gCTAACCTGTCTAGGATCGAACACTCCGAAGAGCACACCCATATTATAGCCCTCGCCCGGCAGGCCGAACACCACGCAGTGGGCTGGGTTTGGCGCAAGGGCCCAAGCCCCCCGACGGGCGGAACGCCAACATTCCCGCCGCGCcgtc                                                                                                 <  2:449892/148‑1 (MQ=255)
                                                     tAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCgg                                                                                   >  1:607183/1‑151 (MQ=255)
                                                     tAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCgg                                                                                   >  1:189004/1‑151 (MQ=255)
                                                     tAAGATCGACCACTCCGCAGACCCCCCCCATATTGATGGCCTCGCGCGCCATGCCGACCCCCACCCAACTTGCTTCGTTTTGCCCAGGGTTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATa                                                                                       <  2:38645/147‑1 (MQ=255)
                                                         aTCGACCACTTCGCAGACACCCCCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCACCTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAg                                                                               <  1:671857/151‑1 (MQ=255)
                                                         aTCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTc                                                                                 >  2:715226/1‑149 (MQ=255)
                                                            gACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGAt                                                                            >  1:22321/1‑151 (MQ=255)
                                                              ccccTTCGCAGCCACCCCCCATATTGATGGCCTCGCGCGGCATGCCGAACCCCACGCACCTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGa                                                                             <  2:152605/145‑1 (MQ=255)
                                                                 cTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCAccc                                                                       >  2:265320/1‑151 (MQ=255)
                                                                    cGCAGCCCCCACCCATATTGATGGCCTCGCGCGGCATGCCGAACCCCACGCAACTTTCTTCTTTTTGCGCAGGGGTCCATGCCCCCCGCCTGACGCATCGCCACCATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCAccc                                                                       <  2:408990/148‑1 (MQ=255)
                                                                      cAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCCCGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCAAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCAccc                                                                       <  1:100452/146‑1 (MQ=255)
                                                                        gcccccccccATATGGATGCCCTCGCGCGGCATGCCAAACCCCACGAACCTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCCCCAACATCCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAAcc                                                                  <  2:523237/142‑1 (MQ=255)
                                                                            ccccccATATTGATGGCCTCGCGCGGCATGCCGAACCCCACGCAACTTGCTTCGTTTTGCGCAGGGGCCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCAtt                                                             <  2:382077/147‑1 (MQ=255)
                                                                              accCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGTCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTAc                                                           >  2:318815/1‑150 (MQ=255)
                                                                               cccATATTGATGGCCTAGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCTTCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCCGGATCACCCCAACCGCATTACGc                                                         >  2:551727/1‑151 (MQ=255)
                                                                                 cATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGccc                                                       >  2:266836/1‑151 (MQ=255)
                                                                                     ttGATGCCCTCCCGCGGCATGCCGACCACCCCGCAACTTGCTTCGTTTGGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCg                                                     <  2:624896/149‑1 (MQ=255)
                                                                                       gATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTg                                                 >  1:111736/1‑151 (MQ=255)
                                                                                                gcgcgGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGAc                                        >  1:590366/1‑151 (MQ=255)
                                                                                                             acccccccGCACCTTCCTCCGTTTTGCGCAGGGGCCCATGCCCCCCGCCTGACGCACCCCCAACTTCCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGaacaa                            <  2:30772/143‑1 (MQ=255)
                                                                                                                 ccacGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACAt                       >  2:125122/1‑151 (MQ=255)
                                                                                                                  cccccAATTTGCTCCTTTTTGCCCAGGGGCCCATCCCCCCCCCCTGCCCCTCCGCCAACATCACCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATc                      <  2:303520/147‑1 (MQ=255)
                                                                                                                       aaCTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTa                    <  1:380117/148‑1 (MQ=255)
                                                                                                                       aaCTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCg                 >  1:254646/1‑151 (MQ=255)
                                                                                                                         cgtgctgcGTTTTGCGCAAGGGCCCATGCCCCCCGCCTGACCCATGGCCACCATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTAc                   <  2:839153/145‑1 (MQ=255)
                                                                                                                           tGCTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTGCGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAgg             >  2:236482/1‑151 (MQ=255)
                                                                                                                             cTTCGTTTTGCGCAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGcc           >  1:885532/1‑151 (MQ=255)
                                                                                                                                  ttttGCGCAGGGGCCCATCCCCCCCGCTTGCCGCACCCCCAACATCCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGCCg          <  2:188979/147‑1 (MQ=255)
                                                                                                                                  ttttGCGCACGGGCCCAGCCCCCGCGCCTCCCGCAGCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATAACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAg              <  2:370654/143‑1 (MQ=255)
                                                                                                                                  tttgccgcAGGGGCCCAGGCCCCCCCCCTGACGCACCCCAACCATCCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGCCGAt        <  2:694995/144‑1 (MQ=255)
                                                                                                                                      gcgcAAGGGTCCATGCCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGCCGATGACGGt  >  1:479021/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                        
ACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCATGCCCCCGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGCCGATGACGGT  >  NC_000913/1967458‑1967741

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: