Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,901,362 A→G intergenic (‑62/+26) eno ← / ← pyrG phosphopyruvate hydratase/CTP synthase (glutamine hydrolyzing)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,901,3620AG96.8% 104.1 / ‑3.3 31intergenic (‑62/+26)eno/pyrGphosphopyruvate hydratase/CTP synthase (glutamine hydrolyzing)
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (14/16);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

TACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCA  >  NZ_CP009273/2901226‑2901497
                                                                                                                                        |                                                                                                                                       
tACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTGcg                                                                                                                                       <  2:537961/139‑1 (MQ=255)
  cTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTGcgcg                                                                                                                                     >  2:138597/1‑139 (MQ=255)
       gCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTGCGCGTTGcc                                                                                                                                >  1:814509/1‑139 (MQ=255)
                                          ttCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTc                                                                                             >  2:108317/1‑139 (MQ=255)
                                                    atgatTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGc                                                                                   >  2:255877/1‑139 (MQ=255)
                                                    atgatTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGc                                                                                   >  2:533581/1‑139 (MQ=255)
                                                                    tGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCGGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACaa                                                                    <  2:547741/138‑1 (MQ=255)
                                                                         aTTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTg                                                              <  2:653487/139‑1 (MQ=255)
                                                                           tAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGc                                                             >  1:652670/1‑138 (MQ=255)
                                                                            aGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCa                                                            <  1:533581/138‑1 (MQ=255)
                                                                             ggTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCaaa                                                          >  1:30439/1‑139 (MQ=255)
                                                                              gTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAAc                                                         <  2:404466/139‑1 (MQ=255)
                                                                                 ttCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAg                                                       <  2:42256/138‑1 (MQ=255)
                                                                                    cTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCgg                                                    <  2:652670/138‑1 (MQ=255)
                                                                                       aaGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGt                                                  <  1:255877/137‑1 (MQ=255)
                                                                                              tAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCAc                                         >  1:681367/1‑139 (MQ=255)
                                                                                                gTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTGCGCGTTGCCTCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACg                                        <  1:580366/138‑1 (MQ=255)
                                                                                                gTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACg                                        <  2:814509/138‑1 (MQ=255)
                                                                                                    aaCTGAAACTCCAGACAAACAACGCGTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAg                                   <  2:346041/139‑1 (MQ=255)
                                                                                                         aaaCTCCAGACAAACAACGCGTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAg                                >  1:442526/1‑137 (MQ=255)
                                                                                                                aGACAAACAACGCGTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCa                                             >  2:112699/1‑117 (MQ=255)
                                                                                                                aGACAAACAACGCGTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCa                                             <  1:112699/117‑1 (MQ=255)
                                                                                                                aGACAAACAACGCGTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACt                       <  2:285294/139‑1 (MQ=255)
                                                                                                                     aacaacGCGTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGa                  <  1:759438/139‑1 (MQ=255)
                                                                                                                       caacGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATg                >  2:452194/1‑139 (MQ=255)
                                                                                                                        aacGCGTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATg                >  2:300379/1‑138 (MQ=255)
                                                                                                                            cgTACCCAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACt           >  1:261386/1‑139 (MQ=255)
                                                                                                                                cccAGGGTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGc        >  2:213292/1‑138 (MQ=255)
                                                                                                                                     ggTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCa  <  2:442526/139‑1 (MQ=255)
                                                                                                                                     ggTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCa  >  1:548389/1‑139 (MQ=255)
                                                                                                                                     ggTGCGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCa  <  2:786096/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
TACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCA  >  NZ_CP009273/2901226‑2901497

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