Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,049,913 C→T 9.4% 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,9130CT 9.4% 138.6 / 5.1 54intergenic (+46/‑55)quiA/oprB‑IIIquinate dehydrogenase/carbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base C (22/26);  new base T (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

GTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCAC  >  NC_002947/4049766‑4050059
                                                                                                                                                   |                                                                                                                                                  
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)
                                                                                                                                                 cccTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCAc  >  8:101667/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCAC  >  NC_002947/4049766‑4050059

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