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

ACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑G‑CCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGCCGATGACGGT  >  NC_000913/1967458‑1967741
                                                                                                                                                      |                                                                                                                                        
aCGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGTGTCCAT‑G‑c                                                                                                                                         <  1:908385/150‑1 (MQ=255)
aCGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑G‑cc                                                                                                                                        >  2:259531/1‑151 (MQ=255)
aCGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑G‑cc                                                                                                                                        >  1:870162/1‑151 (MQ=255)
    aTCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCGGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑G‑cccc‑cc                                                                                                                                   <  1:902585/151‑2 (MQ=255)
     tCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑G‑cccc‑cc                                                                                                                                   >  1:429623/1‑149 (MQ=255)
        ccTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑G‑cccc‑ccgcct                                                                                                                               >  1:607887/1‑146 (MQ=255)
           gTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGa                                                                                                                              <  1:239495/150‑1 (MQ=255)
           gTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGAc                                                                                                                             >  1:478225/1‑151 (MQ=255)
              cGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGc                                                                                                                           >  1:450662/1‑150 (MQ=255)
                gCACTAATCTTTGCCAACATTTGCTGGCTTACCTGGCTAAGTTCGATCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCCAACCCCACGCACCCTTCTTCGCTTTGCGCACGGGTCCAT‑GCCCCC‑CGCCTGCCGCa                                                                                                                          >  2:387518/1‑149 (MQ=255)
                  aCTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGc                                                                                                                      >  1:106855/1‑151 (MQ=255)
                     aaTTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGcc                                                                                                                     >  1:510603/1‑149 (MQ=255)
                            gCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTccc                                                                                                            <  1:394245/151‑1 (MQ=255)
                            gCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTcc                                                                                                             <  1:56127/150‑1 (MQ=255)
                              cAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTCccgc                                                                                                          <  1:544650/151‑1 (MQ=255)
                                          gCTAACCTGTCTAGGATCGAACACTCCGAAGAGCACACCCATATTATAGCCCTCGCCCGGCAGGCCGAACACCACGCAGTGGGCTGGGTTTGGCGCAAGGGCCCAA‑GCCCCC‑CGACGGGCGGAACGCCAACATTCCCGCCGCGCcgtc                                                                                                 <  2:449892/148‑1 (MQ=255)
                                                     tAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCgg                                                                                   >  1:607183/1‑151 (MQ=255)
                                                     tAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCgg                                                                                   >  1:189004/1‑151 (MQ=255)
                                                     tAAGATCGACCACTCCGCAGACCCCCCCCATATTGATGGCCTCGCGCGCCATGCCGACCCCCACCCAACTTGCTTCGTTTTGCCCAGGGTTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATa                                                                                       <  2:38645/147‑1 (MQ=255)
                                                         aTCGACCACTTCGCAGACACCCCCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCACCTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAg                                                                               <  1:671857/151‑1 (MQ=255)
                                                         aTCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTc                                                                                 >  2:715226/1‑149 (MQ=255)
                                                            gACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGAt                                                                            >  1:22321/1‑151 (MQ=255)
                                                              ccccTTCGCAGCCACCCCCCATATTGATGGCCTCGCGCGGCATGCCGAACCCCACGCACCTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGa                                                                             <  2:152605/145‑1 (MQ=255)
                                                                 cTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCAccc                                                                       >  2:265320/1‑151 (MQ=255)
                                                                    cGCAGCCCCCACCCATATTGATGGCCTCGCGCGGCATGCCGAACCCCACGCAACTTTCTTCTTTTTGCGCAGGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCACCATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCAccc                                                                       <  2:408990/148‑1 (MQ=255)
                                                                      cAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCCCGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCAAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCAccc                                                                       <  1:100452/146‑1 (MQ=255)
                                                                        gcccccccccATATGGATGCCCTCGCGCGGCATGCCAAACCCCACGAACCTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCCCCAACATCCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAAcc                                                                  <  2:523237/142‑1 (MQ=255)
                                                                            ccccccATATTGATGGCCTCGCGCGGCATGCCGAACCCCACGCAACTTGCTTCGTTTTGCGCAGGGGCCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCAtt                                                             <  2:382077/147‑1 (MQ=255)
                                                                              accCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑GTCCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTAc                                                           >  2:318815/1‑150 (MQ=255)
                                                                               cccATATTGATGGCCTAGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCTTCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCCGGATCACCCCAACCGCATTACGc                                                         >  2:551727/1‑151 (MQ=255)
                                                                                 cATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGccc                                                       >  2:266836/1‑151 (MQ=255)
                                                                                     ttGATGCCCTCCCGCGGCATGCCGACCACCCCGCAACTTGCTTCGTTTGGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCg                                                     <  2:624896/149‑1 (MQ=255)
                                                                                       gATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTg                                                 >  1:111736/1‑151 (MQ=255)
                                                                                                gcgcgGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGAc                                        >  1:590366/1‑151 (MQ=255)
                                                                                                             acccccccGCACCTTCCTCCGTTTTGCGCAGGGGCCCAT‑GCCCCC‑CGCCTGACGCACCCCCAACTTCCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGaacaa                            <  2:30772/143‑1 (MQ=255)
                                                                                                                 ccacGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACAt                       >  2:125122/1‑151 (MQ=255)
                                                                                                                  cccccAATTTGCTCCTTTTTGCCCAGGGGCCCATCC‑CCCC‑CCCCTGCCCCTCCGCCAACATCACCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATc                      <  2:303520/147‑1 (MQ=255)
                                                                                                                       aaCTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTa                    <  1:380117/148‑1 (MQ=255)
                                                                                                                       aaCTTGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCg                 >  1:254646/1‑151 (MQ=255)
                                                                                                                         cgtgctgcGTTTTGCGCAAGGGCCCAT‑GCCCCC‑CGCCTGACCCATGGCCACCATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTAc                   <  2:839153/145‑1 (MQ=255)
                                                                                                                           tGCTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTGCGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAgg             >  2:236482/1‑151 (MQ=255)
                                                                                                                             cTTCGTTTTGCGCAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGcc           >  1:885532/1‑151 (MQ=255)
                                                                                                                                  ttttGCGCAGGGGCCCATCC‑CCCC‑CGCTTGCCGCACCCCCAACATCCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGCCg          <  2:188979/147‑1 (MQ=255)
                                                                                                                                  ttttGCGCACGGGCCCAG‑C‑CCCCGCGCCTCCCGCAGCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATAACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAg              <  2:370654/143‑1 (MQ=255)
                                                                                                                                  tttgccgcAGGGGCCCAG‑GCCCCC‑CCCCTGACGCACCCCAACCATCCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGCCGAt        <  2:694995/144‑1 (MQ=255)
                                                                                                                                      gcgcAAGGGTCCAT‑GCCCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGCCGATGACGGt  >  1:479021/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                        
ACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACACCACGCAACTTGCTTCGTTTTGCGCAAGGGTCCAT‑G‑CCCC‑CGCCTGACGCATCGCCAACATTCCCGCCGCGCCGTCGTTGCCCATACCGGTCAGGATCACCCCAACCGCATTACGCCCCGCCTGTTTGGCGACAGAATGGAACAACACATCTACCGAAGGCCGATGACGGT  >  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: ‑