Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,931,947 C→T 9.7% G310D (GGC→GAC)  cheBA ← chemotaxis response regulator protein‑glutamate methylesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,931,9470CT 9.7% 121.9 / 2.6 41G310D (GGC→GAC) cheBAchemotaxis response regulator protein‑glutamate methylesterase
Reads supporting (aligned to +/- strand):  ref base C (27/10);  new base T (2/2);  total (29/12)
Fisher's exact test for biased strand distribution p-value = 5.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCCAGGTTCGCCTTGACGATGGCCATGGGCATGCCATAAATCACGCAGCTGGCTTCATCCTGGGCCCACACCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTGTAGTTCAGGCGCTCGTCGCCGGGCAGGATCTTCACCGTGCCACGGCCA  >  NC_002947/4931819‑4932072
                                                                                                                                |                                                                                                                             
ccGCCAGGTTCGCCTTGACGATGGCCATGGGCATGCCATAAATCACGCAGCTGGCTTCATCCTGGGCCCACACCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAg                                                                                                                        >  1:243288/1‑136 (MQ=255)
   ccAGGTTCGCCTTGACGATGGCCATGGGCATGCCATAAATCACGCAGCTGGCTTCATCCTGGGCCCACCCCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGcac                                                                                                                     >  2:202432/1‑136 (MQ=255)
   ccAGGTTCGCCTTGACGATGGCCATGGGCATGCCATAAATCACGCAGCTGGCTTCATCCTGGGCCCACACCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGcac                                                                                                                     >  1:126833/1‑136 (MQ=255)
   ccAGGTTCGCCTTGACGATGGCCATGGGCATGCCATAAATCACGCAGCTGGCTTCATCCTGGGCCCACACCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGcac                                                                                                                     >  2:143920/1‑136 (MQ=255)
   ccAGGTTCGCCTTGACGATGGCCATGGGCATGCCATAAATCACGCAGCTGGCTTCATCCTGGGCCCACACCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGcac                                                                                                                     >  1:65817/1‑136 (MQ=255)
   ccAGGTTCGCCTTGACGATGGCCATGGGCATGCCATAAATCACGCAGCTGGCTTCATCCTGGGCCCACACCGAGCTGCCCCCCTGCGTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATCCCGGTGAGcac                                                                                                                     >  2:258405/1‑136 (MQ=255)
          cGCCTTGACGATGGCCATGGGCATGCCATAAATCACGCAGCTGGCTTCATCCTGGGCCCACACCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGac                                                                                                              >  1:148528/1‑136 (MQ=255)
          cGCCTTGACGATGGCCATGGGCATGCCATAAATCACGCAGCTGGCTTCATCCTGGGCCCACACCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGac                                                                                                              >  1:387006/1‑136 (MQ=255)
              ttGACGATGGCCATGGGCATGCCATAAATCACGCAGCTGGCTTCATCCTGGGCCCACACCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGacaacc                                                                                                          >  2:294297/1‑135 (MQ=255)
                    aTGGCCATGGGCATGCCATAAATCACGCAGCTGGCTTCATCCTGGGCCCACACCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGt                                                                                                    >  1:83954/1‑136 (MQ=255)
                         cATGGGCATGCCATAAATCACGCAGCTGGCTTCATCCTGGGCCCACACCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCg                                                                                               >  1:386676/1‑136 (MQ=255)
                                 tGCCATAAATCACGCAGCTGGCTTCATCCTGGGCCCACACCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACtt                                                                                       <  1:256918/136‑1 (MQ=255)
                                     aTAAATCACGCAGCTGGCTTCATCCTGGGCCCACACCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGTCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGgcc                                                                                   >  2:383188/1‑136 (MQ=255)
                                     aTAAATCACGCAGCTGGCTTCATCCTGGGCCCACACCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGTCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGgcc                                                                                   >  2:384760/1‑136 (MQ=255)
                                                 gCTGGCTTCATCCTGGGCCCACACCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGTCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAg                                                                       <  1:383188/136‑1 (MQ=255)
                                                 gCTGGCTTCATCCTGGGCCCACACCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGTCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAg                                                                       <  1:384760/136‑1 (MQ=255)
                                                 gCTGGCTTCATCCTGGGCCCACACCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAg                                                                       <  2:368339/136‑1 (MQ=255)
                                                            ccTGGGCCCACACCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCcac                                                            <  2:66291/136‑1 (MQ=255)
                                                              tGGGCCCACACCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCTGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCcacgc                                                          <  1:8006/136‑1 (MQ=255)
                                                                 gCCCACACCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCcacgcacg                                                       <  1:43989/136‑1 (MQ=255)
                                                                        ccGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTGTa                                                >  1:200991/1‑136 (MQ=255)
                                                                        ccGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTGTa                                                >  2:394919/1‑136 (MQ=255)
                                                                        ccGTGCTGCCCCCCTGCATGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACATAGTCACCGTACGCCTTGGCCGCCGAACCGAACGTGATGTCCACGCACCGCTTGTa                                                >  1:66525/1‑136 (MQ=255)
                                                                             ctgcCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTGTAGTTCa                                           >  1:370904/1‑136 (MQ=255)
                                                                             ctgcCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTGTAGTTCa                                           >  1:295297/1‑136 (MQ=255)
                                                                              tgcCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGACCCGAAGGTGATGTCCACGCACGGCTTGTAGTTCAg                                          >  1:136805/1‑136 (MQ=255)
                                                                              tgcCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTGTAGTTCAg                                          >  1:111261/1‑136 (MQ=255)
                                                                              tgcCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTGTAGTTCAg                                          >  1:13369/1‑136 (MQ=255)
                                                                              tgcCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTGTAGTTCAg                                          >  1:293091/1‑136 (MQ=255)
                                                                                        cTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTGTAGTTCAGGCGCtcgtcg                                >  1:126901/1‑136 (MQ=255)
                                                                                          tGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCCCCGACAACCCCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTGTAGTTCAGGCGCTCGTCGcc                              >  1:179717/1‑136 (MQ=255)
                                                                                           gAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTGTAGTTCAGGCGCTCGTCGCCg                             <  2:1338/136‑1 (MQ=255)
                                                                                           gAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTGTAGTTCAGGCGCTCGTCGCCg                             <  2:93430/136‑1 (MQ=255)
                                                                                             gCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTGTAGTTCAGGCGCTCGTCGCCggg                           <  2:165472/136‑1 (MQ=255)
                                                                                               aGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTGTAGTTCAGGCGCTCGTCGCCGGGCa                         >  2:350237/1‑136 (MQ=255)
                                                                                               aGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTGTAGTTCAGGCGCTCGTCGCCGGGCa                         >  1:95846/1‑136 (MQ=255)
                                                                                                    cgcgcCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTGTAGTTCAGGCGCTCGTCGCCGGGCAGGATc                    >  2:99697/1‑136 (MQ=255)
                                                                                                     gcgcCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTGTAGTTCAGGCGCTCGTCGCCGGGCAGGATCt                   >  2:303727/1‑136 (MQ=255)
                                                                                                                    tCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTGTAGTTCAGGCGCTCGTCGCCGGGCAGGATCTTCACCGTGCCACGGc    <  2:29458/136‑1 (MQ=255)
                                                                                                                    tCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTATAGTTCAGGCGCTCGTCGCCGGGCAGGATCTTCACCGTGCCACGGc    <  2:387006/136‑1 (MQ=255)
                                                                                                                      ggCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTGTAGTTCAGGCGCTCGTCGCCGGGCAGGATCTTCACCGTGCCACGGCCa  >  1:43987/1‑136 (MQ=255)
                                                                                                                                |                                                                                                                             
CCGCCAGGTTCGCCTTGACGATGGCCATGGGCATGCCATAAATCACGCAGCTGGCTTCATCCTGGGCCCACACCGTGCTGCCCCCCTGCTTGAGCAGGCGCGCGCCCTCGCGGCCATCGGCGCCCATGCCGGTGAGCACCACCGACAACACCTTGTCACCGTACGACTTGGCCGCCGAACCGAAGGTGATGTCCACGCACGGCTTGTAGTTCAGGCGCTCGTCGCCGGGCAGGATCTTCACCGTGCCACGGCCA  >  NC_002947/4931819‑4932072

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: ‑