Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,049,907 A→G 20.2% intergenic (+40/‑61) 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,9070AG20.2% 68.1 / 14.7 42intergenic (+40/‑61)quiA/oprB‑IIIquinate dehydrogenase/carbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base A (19/12);  new base G (4/4);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.94e-01

CTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCA  >  NC_002947/4049762‑4050055
                                                                                                                                                  |                                                                                                                                                     
cTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑Gc                                                                                                                                                   >  5:71231/1‑149 (MQ=255)
 tACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑Gcc                                                                                                                                                  >  7:14890/1‑149 (MQ=255)
    gTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑GCCCCt                                                                                                                                               >  8:173159/1‑149 (MQ=255)
    gTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑GCCCCt                                                                                                                                               >  7:267142/1‑149 (MQ=255)
       tcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑CCCGCtcc                                                                                                                                             >  7:290270/1‑146 (MQ=255)
         gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAGG‑GGCCTTGTAg                                                                                                                                           <  6:3996/148‑1 (MQ=255)
         gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAGG‑GGCCTTGTAg                                                                                                                                           <  3:287609/148‑1 (MQ=255)
         gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑GCCCCTGTa                                                                                                                                            <  7:173159/147‑1 (MQ=255)
           cGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCACCTACAAG‑GCCCCTGTAGGAg                                                                                                                                        >  8:68847/1‑149 (MQ=255)
                        gCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑GCCCCTGTAGGAGCGGGCTTGCCcgc                                                                                                                           <  8:267142/149‑1 (MQ=255)
                                  gtgATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑GCCCCTGTAGGTGCGGGCTTGCCCGCGAATGCCCCt                                                                                                                 <  8:14890/149‑1 (MQ=255)
                                     aTCACTGCTGGCGGCGCCCGCCAGTCGACGGCCCGTGGCGCCTATGTGATTTCGTACGCTCTGCCGTAAACCCGCGTCGCCTGCGTCGC‑GGGCCCGCCCGCGCCGACCAG‑GCCCCGGGAGGCGCGGGCTGGCCCGCGCATGCCCCGCa                                                                                                               >  5:269138/1‑148 (MQ=255)
                                         ctgctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAGG‑GGCCTTGTAGGAGc                                                                                                                                       >  1:154806/1‑120 (MQ=255)
                                         ctgctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAGG‑GGCCTTGTAGGAGc                                                                                                                                       <  2:154806/120‑1 (MQ=255)
                                            ctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAGCGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGAc                                                                                                          <  1:242691/147‑1 (MQ=255)
                                             tgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGGGGGGCCCCCCCCAACAAG‑GCCCCGGTGGGGGGGGGGGGGGGCGCGATGGCCCCGCAAAaacaac                                                                                                       >  1:345848/1‑149 (MQ=255)
                                                  atcaccgccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCTTAAAACCGCGGCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCgaga                                                                                                 <  7:119180/145‑1 (MQ=255)
                                                     cccgccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑GCCCCGGTAGGAGCGGGCTTGCCCCCGAATGCCCCTCAAAGACAACCCgacaata                                                                                              >  6:89096/1‑144 (MQ=255)
                                                        gccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCCAGCCCCCGGCGTCAAC‑GCGCCTGGAGGGGCCGGCGTGCCCGCGCCTGGCCCGCCAACACCAGCCGCGATTCTgg                                                                                           >  3:219817/1‑148 (MQ=255)
                                                           aGTCGCCTGCCCGGGGCGCCTCGGGGATTTCGTACGCGCTGCCGTACACCCGCGGCGCCGGCTTCGC‑GGGCAAGCCCGCTCCTACCCG‑GCCCCGTTCGGCGCGGGCGTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAAt                                                                                        <  6:269138/149‑1 (MQ=255)
                                                           aGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGAAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAAt                                                                                        <  7:68847/149‑1 (MQ=255)
                                                            gTCGACTGACCGTGGCGAGTTTGTGTTTTCTTCCGCTGTGCCGTAAAACCGCGGCGCCGGCTTCGG‑GGGCAAGCCCGCTCCTACAAG‑GCCCCTGTTGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATg                                                                                       <  8:301851/149‑1 (MQ=38)
                                                            gTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACACG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATg                                                                                       <  4:23309/149‑1 (MQ=255)
                                                                aCTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCAt                                                                                   <  5:89096/149‑1 (MQ=255)
                                                                    aCCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCaaa                                                                                                             <  3:319357/119‑1 (MQ=255)
                                                                                     aTTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCa                                                              <  6:71231/149‑1 (MQ=255)
                                                                                      tttCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCGCAAAGCCCACCCGAGCTTCTGTAATGCCATCCGGCATTCCCATCCCCCCCAc                                                             >  5:276566/1‑149 (MQ=37)
                                                                                          ttACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCCCTCCAACAAG‑GCCCCAGTAGGAGCGGGCTTGCCCGCGAATGCCCCACAAAGACAACCCCAGATTATGTAATGCCa                                                                                    >  6:330192/1‑122 (MQ=38)
                                                                                             cGCTCTGCCGTAAAACCGCGTCGCCTGCTGCGC‑GGGCAAGCCCGCTCCTACAAG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCctgctg                                                      <  4:219817/149‑1 (MQ=255)
                                                                                                           aCCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATcc                                                                                 <  2:60365/108‑1 (MQ=255)
                                                                                                           aCCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATcc                                                                                 >  1:60365/1‑108 (MQ=255)
                                                                                                                    ccTGCTTCGC‑GGGCAAGCCCGCTCCTACACG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCa                                         <  3:239061/139‑1 (MQ=255)
                                                                                                                         ttCGC‑GGGCAAGCCCGCTCCTACAAG‑GCCCCTGTAGGAGCGGGCTTGGCCGCGAATGGCGCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTg                          >  5:330472/1‑149 (MQ=255)
                                                                                                                                       ccGCTCCTACAAG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCa                                         >  4:275338/1‑121 (MQ=255)
                                                                                                                                       ccGCTCCTACAAG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtc              >  5:356975/1‑148 (MQ=255)
                                                                                                                                        cGCTCCTACAAG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct             >  8:213853/1‑148 (MQ=255)
                                                                                                                                         gCTCCTACAGG‑GGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa           >  3:183179/1‑149 (MQ=255)
                                                                                                                                          cTCCTACAGG‑GGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCtt                           >  7:300569/1‑132 (MQ=255)
                                                                                                                                          cTCCTACAGG‑GGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCtt                           <  8:300569/132‑1 (MQ=255)
                                                                                                                                          cTCCTACAGG‑GGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCtt                           >  3:206776/1‑132 (MQ=255)
                                                                                                                                               aCAAG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAccc     >  5:328545/1‑149 (MQ=255)
                                                                                                                                               aCAAG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAccc     >  3:73504/1‑149 (MQ=255)
                                                                                                                                                  aG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgca  <  6:276566/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCA  >  NC_002947/4049762‑4050055

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