Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,049,913 C→T 20.6% intergenic (+46/‑55) quiA → / → oprB‑III quinate dehydrogenase/carbohydrate‑selective porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,049,9130CT20.6% 63.6 / 11.2 34intergenic (+46/‑55)quiA/oprB‑IIIquinate dehydrogenase/carbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base C (11/16);  new base T (4/3);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.74e-01

TCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAG  >  NC_002947/4049767‑4050056
                                                                                                                                                  |                                                                                                                                                
tctcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTg                                                                                                                                                >  5:266696/1‑149 (MQ=255)
    gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTa                                                                                                                                              <  1:132163/147‑1 (MQ=255)
       gAAGACCGTCAAGCAGTACGTGGTGATCACTGCTGGCGGCTCCCGCAAGTCGACTGACCGTGGCGACTATGTGATTTGTTACGCTCTGCCGTATAACCTCGTCGCCTGATTCTCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGc                                                                                                                                         <  6:355149/149‑1 (MQ=37)
       gAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGc                                                                                                                                         <  6:266696/149‑1 (MQ=255)
          gACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCggg                                                                                                                                      <  2:316134/149‑1 (MQ=255)
                         cGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGGCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGc                                                                                                                       <  1:342008/149‑1 (MQ=255)
                          gtggtgATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCCCGCCGACAAGGGCCCAGTAGGAGCGGGCGTGGCCGCGaaaggc                                                                                                                      >  2:119318/1‑145 (MQ=255)
                           gggtgATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGGGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGccc                                                                                                                     <  5:291346/148‑1 (MQ=255)
                                 tCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGATTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCaaa                                                                                                               <  1:119318/149‑1 (MQ=255)
                                 tCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCAGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCaaa                                                                                                               <  4:221869/149‑1 (MQ=255)
                                       ctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑c                                                                                                        <  5:362627/149‑1 (MQ=255)
                                                           aCTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGc                                                                                                                                         <  6:146950/97‑1 (MQ=255)
                                                           aCTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGc                                                                                                                                         >  5:146950/1‑97 (MQ=255)
                                                               aCCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCgg                                                                                                                                       >  5:199090/1‑95 (MQ=255)
                                                               aCCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCgg                                                                                                                                       <  6:199090/95‑1 (MQ=255)
                                                                       gACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTg                                                                                             <  5:179934/128‑1 (MQ=255)
                                                                       gACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTg                                                                                             >  6:179934/1‑128 (MQ=255)
                                                                               gATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGTGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGGAATGCCATCCGCAATTCGCATCAccccc                                                                >  6:184474/1‑149 (MQ=255)
                                                                               gATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCCGAGATTCTGTAATGCCATCCGCAATTCGCATCAccc                                                                  <  1:59759/148‑1 (MQ=255)
                                                                                               ccGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggcc                                                 <  5:165836/148‑1 (MQ=255)
                                                                                                                   cTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGTGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGccc                              <  5:184474/147‑1 (MQ=255)
                                                                                                                   cTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCCGCTGGCCAGGGCCAgcagcagggtgg                              <  3:71324/147‑7 (MQ=38)
                                                                                                                     acgcgGGCAAGCCCGCTCCTACAAGGCCCCTTGAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc                          <  7:280605/148‑1 (MQ=255)
                                                                                                                             aaGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcgac                   >  1:149698/1‑148 (MQ=255)
                                                                                                                                 ccGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct              >  7:76395/1‑149 (MQ=255)
                                                                                                                                 ccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCCCCACTGCCCTTGCCAGCGACGAtct              >  2:251748/1‑149 (MQ=255)
                                                                                                                                 ccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct              >  1:284867/1‑149 (MQ=255)
                                                                                                                                  cGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCCCCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct              >  8:251769/1‑148 (MQ=255)
                                                                                                                                  cGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc             <  8:76395/149‑1 (MQ=255)
                                                                                                                                      ccTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGa         >  4:296783/1‑149 (MQ=255)
                                                                                                                                           aaGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc                          >  5:87852/1‑127 (MQ=255)
                                                                                                                                           aaGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc                          <  6:87852/127‑1 (MQ=255)
                                                                                                                                           aaGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCAGCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCCCTGCCCTTGCCAGCGCCGATCTCAGGACCCgc    >  4:247806/1‑149 (MQ=255)
                                                                                                                                              gCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgcag  >  5:9219/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAA‑CCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAG  >  NC_002947/4049767‑4050056

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