Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,049,926 C→A 8.7% intergenic (+15/‑86) 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,9260CA 8.7% 136.1 / 3.4 46intergenic (+15/‑86)quiA/oprB‑IIIquinate dehydrogenase/carbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base C (9/33);  new base A (2/2);  total (11/35)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

CGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACC  >  NC_002947/4049779‑4050065
                                                                                                                                                   |                                                                                                                                           
cGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCtt                                                                                                                                           <  1:702066/150‑1 (MQ=255)
   gggCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTTGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTcgc                                                                                                                                        <  1:796549/150‑1 (MQ=255)
   gggCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTcgc                                                                                                                                        <  1:105562/150‑1 (MQ=255)
    ggCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTcgcg                                                                                                                                       <  1:1403653/150‑1 (MQ=255)
     gCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCgg                                                                                                                                      <  1:360912/150‑1 (MQ=255)
     gCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGACCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCgg                                                                                                                                      <  1:232865/150‑1 (MQ=255)
         ccAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGNCaa                                                                                                                                  <  2:1190431/150‑1 (MQ=255)
               cggcCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGc                                                                                                                            <  1:1680660/150‑1 (MQ=255)
                tgcCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCt                                                                                                                           <  2:1487580/149‑1 (MQ=255)
                ggcCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGACGATGGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCt                                                                                                                           <  2:1485001/150‑1 (MQ=255)
                 gcCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGTCCGTGGCGACTATTTTATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTc                                                                                                                          <  1:1179352/150‑1 (MQ=255)
                   ccGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCgg                                                                                                                                      <  2:1568989/136‑1 (MQ=255)
                   ccGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCgg                                                                                                                                      >  1:1568989/1‑136 (MQ=255)
                    cGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTa                                                                                                                       <  1:578635/150‑1 (MQ=255)
                       tGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACaa                                                                                                                    <  1:447438/150‑1 (MQ=255)
                            tACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACaggggcc                                                                                                               <  1:710545/150‑7 (MQ=255)
                             aCGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGgccc                                                                                                              <  2:222090/150‑1 (MQ=255)
                             aCGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGACAGCCCGCGCCTACCCGGcccc                                                                                                              >  2:1544098/1‑150 (MQ=255)
                                  tcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACaa                                                                                                                    <  2:1478319/139‑1 (MQ=255)
                                  tcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACaa                                                                                                                    >  1:1478319/1‑139 (MQ=255)
                                  tcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAg                                                                                                                   >  2:1497809/1‑140 (MQ=255)
                                  tcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAg                                                                                                                   <  1:1497809/140‑1 (MQ=255)
                                    gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGa                                                                                                       <  2:1549927/150‑1 (MQ=255)
                                    gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGa                                                                                                       <  2:1468228/150‑1 (MQ=255)
                                     ccGAAGACCGCCGAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGGCGGTCGACTGACCGGGGCGGCTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAg                                                                                                      <  2:534783/150‑1 (MQ=38)
                                      cGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGc                                                                                                     <  2:1749990/150‑1 (MQ=255)
                                        aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCgg                                                                                                   <  2:1227034/150‑1 (MQ=255)
                                        aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCgg                                                                                                   <  2:1255889/150‑1 (MQ=255)
                                        aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGAATATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCgg                                                                                                   <  2:1591194/150‑1 (MQ=255)
                                        aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCCGCAGTCTACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCgg                                                                                                   <  1:409109/150‑1 (MQ=255)
                                         aGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCggg                                                                                                  <  2:581610/150‑1 (MQ=255)
                                         aGACCGGCAAGCAGTACGTGGTGATCAATGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGCCCCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCggg                                                                                                  <  2:1011196/150‑1 (MQ=255)
                                          gACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGc                                                                                                 >  1:1202510/1‑150 (MQ=255)
                                                                aTCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCg                                                                                                                             >  2:1677681/1‑100 (MQ=255)
                                                                aTCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCg                                                                                                                             <  1:1677681/100‑1 (MQ=255)
                                                                            ggcgCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAg                                                                                                      >  1:248760/1‑111 (MQ=38)
                                                                            ggcgCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAg                                                                                                      <  2:248760/111‑1 (MQ=38)
                                                                                       gTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGc                                                                                                                 <  1:788723/89‑1 (MQ=255)
                                                                                       gTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGc                                                                                                                 >  2:788723/1‑89 (MQ=255)
                                                                                                                aTTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGATTcgcg                                                                                                                                       <  1:502705/42‑1 (MQ=255)
                                                                                                                aTTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGATTcgcg                                                                                                                                       <  1:502661/42‑1 (MQ=255)
                                                                                                                aTTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGATTcgcg                                                                                                                                       >  2:502661/1‑42 (MQ=255)
                                                                                                                aTTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGATTcgcg                                                                                                                                       >  2:502705/1‑42 (MQ=255)
                                                                                                                      tACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCctgct                     <  2:624594/150‑1 (MQ=255)
                                                                                                                        cGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCGCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCctgctgc                   <  1:1544098/150‑1 (MQ=255)
                                                                                                                                         gcgTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGcacc  >  1:1344283/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
CGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACC  >  NC_002947/4049779‑4050065

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