Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,588,185 A→G 10.2% E11G (GAG→GGG)  benE‑II → benzoate transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,588,1850AG10.2% 97.7 / 3.5 40E11G (GAG→GGG) benE‑IIbenzoate transport protein
Reads supporting (aligned to +/- strand):  ref base A (17/18);  new base G (2/2);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-01

GCAGCGCAGCGGGCGGCCGCGAGCTTTGCAGGAGGCCTGATCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCGGCAGCATTGCCACGACGATCTCCTACGCCGGCCCGTTGGTGATCATTTTCCACGCCGCCG  >  NC_002947/3588040‑3588331
                                                                                                                                                 |                                                                                                                                                  
gcagcgcagcGGGCGGCCGCGAGCTTTGCAGGAGGCCTGATCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGCAGCAGCCAACACCCAAACAAGAGcc                                                                                                                                                 >  1:131709/1‑149 (MQ=255)
    cgcagcgGCCGGCCGCGAGCTTTGCAGGAGGCCTGATCCACCGAGGGGGCTGTTTCGCGGGCAGGCCCGCTCCCACAGGGGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGcc                                                                                                                                             <  4:58337/149‑1 (MQ=255)
        gcgggcggCCGCGAGCTTTGCAGGAGGCCTGATCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGCAGCAGCCAACACCCAAACAAGAGCCCGCcaca                                                                                                                                          >  7:206647/1‑148 (MQ=255)
          gggcggCCGCGAGCTTTGCAGGAGGCCTGATCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAg                                                                                                                                       <  7:137685/149‑1 (MQ=255)
               gccggtaGCTTTGCTGGTGGCCTTATCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCTGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAggccgg                                                                                                                                  <  2:131709/143‑1 (MQ=255)
               gcCGCGAGGTTTGCAGGTGGCCTGATCCACCCAGGGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCATTCGTTGTTGGAGCGGGCATGCCCGCGAAGCATCCAAAACCCAAACAAGAGCCCGCCACAAAggccgg                                                                                                                                  <  8:177618/149‑1 (MQ=255)
                        tttGCAGGAGGCCTGATCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTACCACAGGTGAAGCCCGTGCAGTCGTTGTTGGAGCCGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAACGCCCGCCCCAGcca                                                                                                                         >  8:45148/1‑148 (MQ=255)
                          tGCAGGAGGCCTGATCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGCAGCAGCCCACCCCCAAACAAGCGCCCCCCCCAAAGGCCCGCCCCGGccggg                                                                                                                       >  8:265593/1‑146 (MQ=255)
                           gcatgcggCCTGATCCACCGCGTGGGCTGTTTCGCGGGCATGCCCGCACCCACAGGTTAAGCCTGTGCAGTCGTTGTTGGAGCGGGAATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTg                                                                                                                       <  8:206647/142‑1 (MQ=255)
                           gCAGGAGGCCTGATCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTg                                                                                                                       <  2:64283/148‑1 (MQ=255)
                           gCAGGAGGCCTGATCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTg                                                                                                                       <  7:260429/148‑1 (MQ=255)
                            cAGGAGGCCTGATCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGc                                                                                                                      <  1:43019/148‑1 (MQ=255)
                               gaggCCTGATCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCgaga                                                                                                                  <  7:265593/149‑1 (MQ=255)
                                   ccTGATCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCt                                                                                                              >  5:293980/1‑149 (MQ=255)
                                     tGATCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCtt                                                                                                             >  8:3727/1‑148 (MQ=255)
                                       aTCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCCGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTc                                                                                                            <  7:45148/147‑1 (MQ=255)
                                                      tGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAACCGCTGATc                                                                                             >  1:43010/1‑147 (MQ=255)
                                                            gcgGGCATGCCCGATCCCACAGGTGAAGACTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCACCCCCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGAGAGAAGCTTCATGAAAACGCTGAACAAGGACTg                                                                                     >  7:130517/1‑149 (MQ=255)
                                                                gCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGctct                                                                                 <  6:293980/149‑1 (MQ=255)
                                                                       cGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCgg                                                                          >  2:175656/1‑149 (MQ=255)
                                                                                    gAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCgg                                                             >  5:281326/1‑149 (MQ=255)
                                                                                               aGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGGGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGc                                                                <  1:11056/135‑1 (MQ=255)
                                                                                               aGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGGGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGc                                                                <  1:15837/135‑1 (MQ=255)
                                                                                               aGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGGGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGc                                                                >  2:11056/1‑135 (MQ=255)
                                                                                               aGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGGGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGc                                                                >  2:15837/1‑135 (MQ=255)
                                                                                                 tCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCAtcgtc                                                                  <  4:295848/131‑1 (MQ=255)
                                                                                                 tCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCAtcgtc                                                                  >  3:295848/1‑131 (MQ=255)
                                                                                                   gttgttGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCGGCAGCATTGCCacgac                                              >  8:84352/1‑149 (MQ=255)
                                                                                                    ttgttgGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAACCGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCGGCAGCATTGCCacgacg                                             <  2:43010/149‑1 (MQ=255)
                                                                                                      gttgGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAg                                                                                          <  2:43706/102‑1 (MQ=255)
                                                                                                      gttgGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAg                                                                                          >  1:43706/1‑102 (MQ=255)
                                                                                                              gggCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCtgtg                                                                              >  3:255892/1‑106 (MQ=255)
                                                                                                              gggCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCtgtg                                                                              <  4:255892/106‑1 (MQ=255)
                                                                                                                     cccGCTAAGCAGCCAACACCCAAACAATAGCCCTACACAAAGGCCTGCCTCTGCCTGGCTTGAATCTTCATGAAAACGCTTCCCAAGGACTGCTCTGTGCCGTATTTGGTCGCCGGCAGCATGGCCACGACGATCTCCTACGCCGGccc                            <  8:130517/149‑1 (MQ=255)
                                                                                                                         cgAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCGGCAGCATTGCCACGACGATCTCCTACGCCGGCCCGTTg                        <  7:3727/149‑1 (MQ=255)
                                                                                                                          gAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCGGCAGCATTGCCACGACGATCTCCTACGCCGGCCCGTTgg                       >  2:17585/1‑149 (MQ=255)
                                                                                                                                      acCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCGGCAGCATTGCCACGACGATCTCCTACGCCGGCCCGTTGGTGa                    <  6:287849/140‑1 (MQ=255)
                                                                                                                                      acCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCGGCAGCATTGCCACGACGATCTCCTACGCCGGCCCGTTGGTGa                    >  5:287849/1‑140 (MQ=255)
                                                                                                                                        ccAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCGGCAGCATTGCCACGACGAtc                                          >  7:276418/1‑116 (MQ=255)
                                                                                                                                         cAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCGGCAGCATTGCCACGACGATCTCCTACGCCGGCCCGTTGGTGATCATTTTCCAcg        >  3:142321/1‑149 (MQ=255)
                                                                                                                                                gagCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCGGCAGCATTGCCACGACGATCTCCTACGCCGGCCCGTTGGTGATCATTTTCCACGccgccg  >  7:293006/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCAGCGCAGCGGGCGGCCGCGAGCTTTGCAGGAGGCCTGATCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCGGCAGCATTGCCACGACGATCTCCTACGCCGGCCCGTTGGTGATCATTTTCCACGCCGCCG  >  NC_002947/3588040‑3588331

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