Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,698,641 T→C 10.0% V130A (GTC→GCC)  cheBC → chemotaxis response regulator protein‑glutamate methylesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,698,6410TC10.0% 100.0 / 3.4 41V130A (GTC→GCC) cheBCchemotaxis response regulator protein‑glutamate methylesterase
Reads supporting (aligned to +/- strand):  ref base T (18/18);  new base C (2/2);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

CCGTAAACAGAACGTCCACCGGGTATTCGAGGCCATGGGCCATGGCGCCCTCGATGTGGTCGACACCCCTGCGTTGGGCGCTGGCGATGCCCGCGAAGCGGCGGCGCCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTT  >  NC_002947/1698498‑1698785
                                                                                                                                               |                                                                                                                                                 
ccGTAAACAGAACGTCCACCGGGTATTCGAGGCCATGGGCCATGGCGCCCTCGATGTGGTCGACACCCCTGCGTTGGGCGCTGGCGATGCCCGCGAAGCGGCGGCGCCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGcc                                                                                                                                              <  7:287335/149‑1 (MQ=255)
                         ttCGAGGCCATGGGCCATGGCGCCCTCGATGTGGTCGACACCCCTGCGTTGGGCGCTGGCGATGCCCGCGAAGCGGCGGCGCCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGt                                                                                                                     <  5:175877/148‑1 (MQ=255)
                                 catgggccatgtCGCCCTCGATGTGGTCGACACCCCTGCGTTGGGCGCTGGCGATGCCCGCGAAGCGGCGGCGCCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTg                                                                                                            <  8:20678/149‑1 (MQ=255)
                                 catgggccatggCGCCCTCGATGTGGTCGACACCCCTGCGTTGGGCGCTGGCGATGCCCGCGAAGCGGCGGCGCCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCCGACcagccagccga                                                                                                                                       >  4:237884/1‑119 (MQ=255)
                                 catgggccatggCGCCCTCGATGTGGTCGACACCCCTGCGTTGGGCGCTGGCGATGCCCGCGAAGCGGCGGCGCCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCCGACcagccagccga                                                                                                                                       <  3:237884/122‑4 (MQ=255)
                                  atgggccatggCGCCCTCGATGTGGTCGACACCCCTGCGTTGGGCGCTGGCGATGCCCGCGAAGCGGCGGCGCCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTg                                                                                                            >  4:47454/1‑148 (MQ=255)
                                  atgggccatggCGCCCTCGATGTGGTCGACACCCCTGCGTTGGGCGCTGGCGATGCCCGCGAAGCGGCGGCGCCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGc                                                                                                           >  3:194005/1‑149 (MQ=255)
                                              gcCCTCGATGTGGTCGACACCCCTGCGTTGGGCGCTGGCGATGCCCGCGAAGCGGCGGCGCCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTgcgc                                                                                                 >  2:14097/1‑147 (MQ=255)
                                                      tgtgGTCGACACCCCTGCGTTGGGCGCTGGCGATGCCCGCGAAGCGGCGGCGCCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTcc                                                                                        >  8:39356/1‑148 (MQ=255)
                                                      tgtgGTCGACACACCTGCGATGGGCGCTGGCGATGCCCGCGAAGCGCCGGCCCCCTTCCAGCCCAAGATCCTCAATATCGGCTCGCTGGTCGGCCAGCAGCGTGCAACTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCATcc                                                                                        >  4:239625/1‑148 (MQ=255)
                                                        tgGTCGACACCCCTGCGTTGGGCGCTGGCGATGCCCGCGAAGCGGCGGCGCCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTg                                                                                                                                    <  1:277596/102‑1 (MQ=255)
                                                        tgGTCGACACCCCTGCGTTGGGCGCTGGCGATGCCCGCGAAGCGGCGGCGCCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTg                                                                                                                                    >  2:277596/1‑102 (MQ=255)
                                                                 ccccTGCGTTGGGCGCTGGCGATGCCCGCGAAGCGGCGGCGCCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCggg                                                                            <  7:219787/149‑1 (MQ=255)
                                                                  cccTGCGTTGGGCGCTGGCGATGCCCGCGAAGCGGCGGCGCCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGt                                                                           >  6:118820/1‑149 (MQ=255)
                                                                  cccTGCGTTGGGCGCTGGCGATGCCCGCGAAGCGGCGGCGCCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGt                                                                           >  2:270582/1‑149 (MQ=255)
                                                                  cccTGCGTTGGGCGCTGGCGATGCCCGCGAAGCGGCGGCGCCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGt                                                                           >  5:105863/1‑149 (MQ=255)
                                                                        gTTGGGCGCTGGCGATGCCCGCGAAGCGGCGGCGCCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTtggtgg                                                                     <  1:14097/149‑1 (MQ=255)
                                                                                          ccGCGAAGCGGCGGCGCCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTcggcgg                                                   >  3:213734/1‑149 (MQ=255)
                                                                                           cgcgAAGCGGCGGCGCCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGTAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCggg                                                  >  2:99117/1‑149 (MQ=255)
                                                                                                        cgcCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGgcagca                                     <  5:118820/149‑1 (MQ=255)
                                                                                                               gctgcGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGa                                <  3:47454/147‑1 (MQ=255)
                                                                                                                   cgcAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGtgc                           <  1:153788/148‑1 (MQ=255)
                                                                                                                    gcAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGGTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGtgcgg                         >  4:239484/1‑147 (MQ=255)
                                                                                                                      aaGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGTAAGGTCGGTGGCTGCACCGCTCCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTTCGAGtgctgcc                       <  1:99117/149‑2 (MQ=255)
                                                                                                                      aaGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTTCGAGtgctgcc                       <  8:146612/149‑2 (MQ=255)
                                                                                                                          tCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGgca                                        >  6:316383/1‑128 (MQ=255)
                                                                                                                          tCCTCAATATCGGCTGGCTGGGCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGGCTTCGGCGGGTGGCCCGgca                                        <  5:316383/128‑1 (MQ=255)
                                                                                                                           ccTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAgg                  >  7:230323/1‑149 (MQ=255)
                                                                                                                            cTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGGAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTTCTCAAgg                  <  3:239625/148‑1 (MQ=255)
                                                                                                                            cTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACGGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGTGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAgg                  <  7:39356/148‑1 (MQ=255)
                                                                                                                            cTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAggg                 >  3:218549/1‑149 (MQ=255)
                                                                                                                             tCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGt                >  8:220178/1‑149 (MQ=255)
                                                                                                                               aaTATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGttt              >  6:99866/1‑149 (MQ=255)
                                                                                                                               aaTATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGttt              <  6:105863/149‑1 (MQ=255)
                                                                                                                               aaTATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGttt              <  4:218549/149‑1 (MQ=255)
                                                                                                                                 tatCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGc            <  4:213734/149‑1 (MQ=255)
                                                                                                                                   tCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCg          <  7:220178/149‑1 (MQ=255)
                                                                                                                                     ggctggctGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGc        <  8:230323/149‑1 (MQ=255)
                                                                                                                                     acgctgctggCCGACCAGCCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGaa                               >  4:159905/6‑127 (MQ=255)
                                                                                                                                     acgctgctggCCGACCAGCCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGaa                               <  3:159905/122‑1 (MQ=255)
                                                                                                                                           ctggTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGtt  >  4:253849/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CCGTAAACAGAACGTCCACCGGGTATTCGAGGCCATGGGCCATGGCGCCCTCGATGTGGTCGACACCCCTGCGTTGGGCGCTGGCGATGCCCGCGAAGCGGCGGCGCCCTTGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTT  >  NC_002947/1698498‑1698785

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