Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,049,910 C→G 9.3% intergenic (+43/‑58) 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,9100CG 9.3% 147.9 / 6.8 54intergenic (+43/‑58)quiA/oprB‑IIIquinate dehydrogenase/carbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base C (22/27);  new base G (2/3);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

CTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAG  >  NC_002947/4049762‑4050056
                                                                                                                                                    |                                                                                                                                                  
cTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGc                                                                                                                                                    >  4:121966/1‑149 (MQ=255)
    gTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCt                                                                                                                                                >  7:280403/1‑149 (MQ=255)
     tctcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCGGCCAGGCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTg                                                                                                                                               <  8:339593/149‑1 (MQ=255)
            gAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGc                                                                                                                                        >  1:59133/1‑149 (MQ=255)
             aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCTTACAAGGCCCCTGTAGGAGCg                                                                                                                                       <  4:10039/149‑1 (MQ=255)
             aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCg                                                                                                                                       <  8:326920/149‑1 (MQ=255)
              aGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCgg                                                                                                                                      <  3:121966/149‑1 (MQ=255)
                aCCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCGGGc                                                                                                                                    <  8:338578/149‑1 (MQ=255)
                  cGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCtt                                                                                                                                  <  8:267283/149‑1 (MQ=255)
                    gaaaGCAGTACGTGGTGATCACTGCTGGCGGCGCCGGCCAGTCGTCTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGGGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCgg                                                                                                                                      <  3:58748/141‑1 (MQ=38)
                    gCAAGCAGTACGTGGTGATCACTGCTGGCTGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGAATTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCACCTACACCTGCCTTGTAGGAGCgt                                                                                                                                      >  4:58748/1‑142 (MQ=255)
                               gtggtgATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGcc                                                                                                                     >  6:169168/1‑149 (MQ=255)
                                      tCACTGCTGGGGGCGCCCGCGAGTCGACTGCCCGTGGCGACTATGTGATTTCTTCCTCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCaaa                                                                                                              <  4:105132/149‑1 (MQ=38)
                                             tgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCaaa                                                                                                              >  4:50871/1‑142 (MQ=255)
                                             tgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCaaa                                                                                                              <  3:50871/142‑1 (MQ=255)
                                               gcggcgCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGGGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAAccc                                                                                                      <  8:49225/148‑1 (MQ=255)
                                                cggcgCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCga                                                                                                    <  2:59133/149‑1 (MQ=255)
                                                  tcgcccgccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCgaga                                                                                                  <  5:84010/148‑1 (MQ=255)
                                                      ccgccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCt                                                                                              <  3:92975/149‑1 (MQ=255)
                                                      ccgccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCt                                                                                              <  5:169168/149‑1 (MQ=255)
                                                      ccgccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAACCCCCCTACAAAAAAGCCCCAGTAGGAGCCGGCGTGCCCGCGAAAGCCCCTCAAAGACAACACCAGAttat                                                                                              >  4:247946/1‑147 (MQ=255)
                                                          cAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTaa                                                                                          <  8:280403/149‑1 (MQ=255)
                                                              ctcctccccGTGGCCCATATGTTATTCCTTCCGCTCTCCCTAAAACCCGCGGCCCCTGCTGCGCGGCAACCCCCGCCCCAACACGCCCCCTGTCGGCGGGGGCTGCCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGcc                                                                                      <  6:251662/142‑1 (MQ=255)
                                                                aCTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGCCAACCCGAGATTCTGTAATGCCAt                                                                                    >  5:273667/1‑149 (MQ=255)
                                                                        tGGCGACTATGTGATTTCTTACCCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAAtt                                                                            <  2:62785/149‑1 (MQ=255)
                                                                            gACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGcc                                                                        >  5:167138/1‑148 (MQ=255)
                                                                              cTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATc                                                                      <  7:317743/149‑1 (MQ=255)
                                                                                   tgATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTACAACAACGCCCCTGTAGGAGCCGGCTTGCCCGCGAAAGCCCCGCAAACACCACCCGAGATTATGTAATGCACTCCGCAATTAGCATCAcccc                                                                 >  6:146398/1‑149 (MQ=255)
                                                                                    gATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCGACTACCAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTTATGCCATGCGCAATTGGCATCAccccc                                                                >  2:107353/1‑149 (MQ=39)
                                                                                      tttCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCGACAAGGCCCCGGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACCAACCGAGATTATGTAATGCCATCCGCAATTCGCATACCCCCCAc                                                              >  8:218636/1‑149 (MQ=38)
                                                                                         cTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCAcgc                                                            <  6:367378/148‑3 (MQ=255)
                                                                                         cTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTACTACAATGCCCCTGTAGGAGCGGGCTTGCCCGCGAAAGCCCCTCAAAGACAACCCCAGATTCTGTAATGCCATCCGCAATTGGCATCACCCCCACCct                                                           >  5:106818/1‑149 (MQ=39)
                                                                                          ttACGCTCTGCCTTAAAACCGCGGCGCCTGCTTCGCGGGAAAGCCCGCTCCTACACGGCCCGTGTAGTAGGGGGGTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCctg                                                          <  4:2074/149‑1 (MQ=255)
                                                                                          ttACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTACTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTATGTAATGCCCTCCGCAATTCGCATCACACccccccag                                                          >  7:150941/1‑147 (MQ=38)
                                                                                             cGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCctgct                                                        <  5:22811/148‑1 (MQ=255)
                                                                                                tctGCCTTAAAACCGCGGCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTATGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCt                                                                                              <  6:208822/107‑1 (MQ=255)
                                                                                                tctGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCCCTCCTACAAGGCCCCTGTAAGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCt                                                                                              >  5:208822/1‑107 (MQ=255)
                                                                                                  tGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCctgctgctg                                                    <  7:63640/147‑1 (MQ=255)
                                                                                                                  cGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACt                                  <  7:218636/149‑1 (MQ=255)
                                                                                                                     cTGCTTCGCGGGCAAGCCCCCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCGCAAAGCCAACCCGAGATTCTGTAATGGCAACCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCAGGCCAGCACCCCGGcc                               >  6:282647/1‑149 (MQ=38)
                                                                                                                        cTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTTTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCtt                            <  6:106818/149‑1 (MQ=255)
                                                                                                                        cTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGccc                              <  1:107353/147‑1 (MQ=255)
                                                                                                                            gcgGGCAAGCCCGCTCCTACAAGCCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCa                        <  5:217044/149‑1 (MQ=255)
                                                                                                                              gggCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGcac                                      <  1:61103/133‑1 (MQ=255)
                                                                                                                              gggCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGcac                                      >  2:61103/1‑133 (MQ=255)
                                                                                                                              gggCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGc                      >  3:277476/1‑149 (MQ=255)
                                                                                                                                     cccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtc               >  7:304933/1‑149 (MQ=255)
                                                                                                                                      ccGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct              >  1:23872/1‑149 (MQ=255)
                                                                                                                                       cGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGc                                                                                                                      <  8:274607/44‑1 (MQ=255)
                                                                                                                                       cGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGc                                                                                                                      >  7:274607/1‑44 (MQ=255)
                                                                                                                                             tACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccctg                                              <  8:308256/110‑1 (MQ=255)
                                                                                                                                             tACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccctg                                              >  7:308256/1‑110 (MQ=255)
                                                                                                                                               cAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCg     >  3:269166/1‑149 (MQ=255)
                                                                                                                                               cAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAGCTCATGACCCg     >  8:343164/1‑149 (MQ=255)
                                                                                                                                                  ggCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgcag  >  6:104322/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAG  >  NC_002947/4049762‑4050056

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