Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,049,913 C→T 19.5% 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,9130CT19.5% 79.3 / 12.8 41intergenic (+46/‑55)quiA/oprB‑IIIquinate dehydrogenase/carbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base C (17/16);  new base T (4/4);  total (21/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 = 4.92e-01

GTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCA  >  NC_002947/4049766‑4050055
                                                                                                                                                     |                                                                                                                                              
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                                                                                                                                           <  3:287609/148‑1 (MQ=255)
     gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAGG‑GGCCTTGTAg                                                                                                                                           <  6:3996/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                           >  3:206776/1‑132 (MQ=255)
                                                                                                                                      cTCCTACAGG‑GGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCtt                           >  7:300569/1‑132 (MQ=255)
                                                                                                                                      cTCCTACAGG‑GGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCtt                           <  8:300569/132‑1 (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)
                                                                                                                                                     |                                                                                                                                              
GTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGC‑GGGCAAGCCCGCTCCTACAAG‑GCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCA  >  NC_002947/4049766‑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: ‑