Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,698,856:1 +C 22.6% coding (562/1014 nt) 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,8561.C22.6% 72.5 / 26.8 31coding (562/1014 nt)cheBCchemotaxis response regulator protein‑glutamate methylesterase
Reads supporting (aligned to +/- strand):  ref base . (9/15);  new base C (4/3);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACAGAATG  >  NC_002947/1698714‑1699003
                                                                                                                                               |                                                                                                                                                   
ggTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCa                                                                                                                                              >  2:331975/1‑150 (MQ=255)
         aCCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGc                                                                                                                                     <  1:267092/150‑1 (MQ=255)
              ggcgcgcGGCTCCCCAGCGACGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGTGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAg                                                                                                                                <  2:236976/149‑1 (MQ=255)
                     ggCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGc                                                                                                                         <  1:429220/150‑1 (MQ=255)
                          cccAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAAt                                                                                                                    <  2:501539/150‑1 (MQ=255)
                                      ggTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAgcagcg                                                                                                        >  2:530646/1‑150 (MQ=255)
                                               cAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGcc                                                                                               <  2:296079/150‑1 (MQ=255)
                                                aaTCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCt                                                                                              >  2:327263/1‑150 (MQ=255)
                                                aaTCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCACATGCTGGACCa                                                                                                                                              >  2:456316/1‑102 (MQ=255)
                                                aaTCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCACATGCTGGACCa                                                                                                                                              <  1:456316/102‑1 (MQ=255)
                                                    gggTCTTCGGCGGGTGGCCCGGCAGCAGTCGAAGTGCTGCTCAAGCGTTTGCCGGCAGCGTTTACGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCTCAGCCATGGCCTAATGGGTCACCAGCGCAGCAGGCCTGCCg                                                                                          >  1:361494/1‑150 (MQ=255)
                                                      gTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCACATGCTGGACCAgc                                                                                                                                            >  1:492098/1‑97 (MQ=255)
                                                      gTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCACATGCTGGACCAgc                                                                                                                                            <  2:492098/98‑2 (MQ=255)
                                                                ggTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTTCGCCttcac                                                                              >  2:396414/1‑146 (MQ=255)
                                                                     cccGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGa                                                                         >  1:353824/1‑150 (MQ=255)
                                                                        ggCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAgg                                                                      <  1:331975/150‑1 (MQ=255)
                                                                         gcagcaCTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGc                                                                     <  1:579810/150‑1 (MQ=255)
                                                                                             cAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCa                                                 <  1:530646/150‑1 (MQ=255)
                                                                                              aaGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAg                                                >  1:519376/1‑150 (MQ=255)
                                                                                              aaGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAg                                                >  2:116758/1‑150 (MQ=255)
                                                                                              aaGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAg                                                <  1:441625/150‑1 (MQ=255)
                                                                                                         ggCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGc                                     <  1:116758/150‑1 (MQ=255)
                                                                                                               gTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTac                               <  1:327263/150‑1 (MQ=255)
                                                                                                                        ggcCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAAccacca                      <  2:502799/150‑1 (MQ=255)
                                                                                                                               caccTGGTCCACATGCTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTac                               <  2:334753/131‑1 (MQ=255)
                                                                                                                               caccTGGTCCACATGCTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTac                               >  1:334753/4‑134 (MQ=255)
                                                                                                                                 ccTGGTCCACATGCTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCt             >  2:119463/2‑150 (MQ=255)
                                                                                                                                        cAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAaccac                        >  2:301477/1‑132 (MQ=255)
                                                                                                                                        cAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAaccac                        <  1:301477/132‑1 (MQ=255)
                                                                                                                                            aTGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACAGAATg  <  2:519376/150‑1 (MQ=255)
                                                                                                                                            aTGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACCCCTGGCCGGTACAAACAACCAGATTCGCCTGCTACAGAATg  <  2:361494/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACAGAATG  >  NC_002947/1698714‑1699003

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