Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,897,050 +G coding (277/666 nt) yieH → phosphoenolpyruvate and 6‑phosphogluconate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,897,0461.G92.1% 122.7 / 3.5 38coding (273/666 nt)yieHphosphoenolpyruvate and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (14/21);  total (15/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTT  >  NC_000913/3896908‑3897156
                                                                                                                                           |                                                                                                               
tttCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                  <  2:108319/139‑1 (MQ=255)
          gtAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAg                                                                                                      >  2:67544/1‑139 (MQ=255)
                  gTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCt                                                                                               >  1:120142/1‑139 (MQ=255)
                    aCGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGt                                                                                             <  2:120142/139‑1 (MQ=255)
                           tATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGAt                                                                                      <  2:830684/139‑1 (MQ=255)
                             tCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCa                                                                                    >  2:756651/1‑139 (MQ=255)
                                            ttGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                     <  2:608856/139‑1 (MQ=255)
                                            tggAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTCCGATTCAGAACGGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                     <  2:631850/137‑1 (MQ=255)
                                                aCATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGt                                                                 <  2:89027/139‑1 (MQ=255)
                                                    ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTcc                                                                                                >  1:759581/1‑104 (MQ=255)
                                                    ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTcc                                                                                                <  2:759581/104‑1 (MQ=255)
                                                         tACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACgg                                                        >  2:437242/1‑139 (MQ=255)
                                                             ttAGCGAAAACAGAAGCTGAACACGTTTACCGTGCATAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCTTCGAAGGGCGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTATCGGCCCa                                                    <  2:392285/139‑1 (MQ=255)
                                                                gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCg                                                                                        <  2:550305/100‑1 (MQ=255)
                                                                gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCg                                                                                        >  1:550305/1‑100 (MQ=255)
                                                                     aaCAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACaa                                            >  1:474690/1‑139 (MQ=255)
                                                                         gAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATg                                        <  1:756651/139‑1 (MQ=255)
                                                                             cTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcagc                                    <  1:230904/139‑1 (MQ=255)
                                                                                aaCACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTc                                <  1:67544/139‑1 (MQ=255)
                                                                                        tACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                         <  2:429666/139‑1 (MQ=255)
                                                                                        tACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                         <  2:689407/139‑1 (MQ=255)
                                                                                        tACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                         >  1:689407/1‑139 (MQ=255)
                                                                                             tGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATGACAAAATGCAGCATTCTATGGGc                         >  2:723052/1‑134 (MQ=255)
                                                                                             tGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATGACAAAATGCAGCATTCTATGGGc                         <  1:723052/134‑1 (MQ=255)
                                                                                               cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATTTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGa                  <  2:530602/139‑1 (MQ=255)
                                                                                                  aaGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCa                                                                            >  2:725236/1‑78 (MQ=255)
                                                                                                  aaGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCa                                                                            <  1:725236/78‑1 (MQ=255)
                                                                                                       gCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAgcg                                                                                                    >  1:538920/1‑49 (MQ=255)
                                                                                                       gCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAgcg                                                                                                    <  2:538920/49‑1 (MQ=255)
                                                                                                        cTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACa                                             <  1:117875/103‑1 (MQ=255)
                                                                                                        cTCGGCTGTTCGATTCAGAACTGCAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACa                                             >  2:117875/1‑103 (MQ=255)
                                                                                                               gTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACtt  <  2:474690/139‑1 (MQ=255)
                                                                                                               gTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACtt  <  2:300453/139‑1 (MQ=255)
                                                                                                                          aaCTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCa                                                                            >  1:396006/1‑54 (MQ=255)
                                                                                                                          aaCTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCa                                                                            <  2:396006/54‑1 (MQ=255)
                                                                                                                          aaCTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                     >  1:511367/1‑61 (MQ=255)
                                                                                                                          aaCTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                     <  2:511367/61‑1 (MQ=255)
                                                                                                                               gAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgt                                                                        >  2:229901/1‑53 (MQ=255)
                                                                                                                               gAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgt                                                                        <  1:229901/53‑1 (MQ=255)
                                                                                                                                           |                                                                                                               
TTTCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTT  >  NC_000913/3896908‑3897156

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: