Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,872,656	G→A	19.2%	A427T (GCA→ACA)	ygcE →	predicted kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,872,656	0	G	A	19.2%	143.1 / 21.3	93	A427T (GCA→ACA)	ygcE	predicted kinase
Reads supporting (aligned to +/- strand):  ref base G (21/54);  new base A (12/6);  total (33/60)
Fisher's exact test for biased strand distribution p-value = 4.81e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATG  >  minE/1872590‑1872716                                                                   |                                                             aaCGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCa                                                            >  1:347793/1‑69 (MQ=255)               tcTGGGGGCTGCGATCAATACGGCGGTAGGTCTGGGGCTATACCCGGATTACACAACGGCTGTCGATaa                                              <  1:2899630/68‑1 (MQ=255)               tcTGGGGGCTGCGATCAATACGGCGGTAGGTCTGGGGCTATACCCGGATTACACAACGGCTGTCGATaa                                              <  1:1199489/68‑1 (MQ=255)               tcTGGGGGCTGCGATCAATACGGCGGTAGGTCTGGGGCTATACCCGGATTACACAACGGCTGTCGATaa                                              <  1:2001243/68‑1 (MQ=255)               tcTGGGGGCTGCGATCAATACGGCGGTAGGTCTGGGGCTATACCCGGATTACACAACGGCTGTCGATaa                                              <  1:1900192/68‑1 (MQ=255)               tcTGGGGGCTGCGATCAATACGGCGGTAGGGCTGGGGCTATACCCGGATTACACAACGGCTGTCGATaa                                              <  1:2063511/68‑1 (MQ=255)               tcTGGGGGCTGCGACCAATACGGCGGTAGGTCTGGGGCTATACCCGGATTACACAACGGCTGTCGATaa                                              <  1:3303850/68‑1 (MQ=255)               ccTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATaa                                              <  1:2189897/69‑1 (MQ=255)                      gCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGAt                                                <  1:111536/60‑1 (MQ=255)                      gCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGAt                                                <  1:2413390/60‑1 (MQ=255)                      gCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGAt                                                <  1:1637856/60‑1 (MQ=255)                      gCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGAt                                                <  1:1410913/60‑1 (MQ=255)                                       ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAAGATGGTTCGCGTGAAAGATATCtt                      <  1:2001960/69‑1 (MQ=255)                                       ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:926472/69‑1 (MQ=255)                                       ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:271771/69‑1 (MQ=255)                                       ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:2708246/69‑1 (MQ=255)                                       ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:1014479/69‑1 (MQ=255)                                       ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:2377137/69‑1 (MQ=255)                                       ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:2920826/69‑1 (MQ=255)                                       ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:2260819/69‑1 (MQ=255)                                       ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:1180047/69‑1 (MQ=255)                                       ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:1254766/69‑1 (MQ=255)                                       ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:2026466/69‑1 (MQ=255)                                       ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:1981500/69‑1 (MQ=255)                                       ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:1402168/69‑1 (MQ=255)                                       ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:356783/69‑1 (MQ=255)                                       ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:1634574/69‑1 (MQ=255)                                       ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:1569386/69‑1 (MQ=255)                                       ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:1667164/69‑1 (MQ=255)                                       ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:1519356/69‑1 (MQ=255)                                       ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:470986/69‑1 (MQ=255)                                       ggTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGATCGCGTGAAAGATATCtt                      <  1:2061144/69‑1 (MQ=255)                                        gTAGGTCTGTGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:2347966/68‑1 (MQ=255)                                        gTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCtt                      <  1:2844952/68‑1 (MQ=255)                                           ggTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTtata                  >  1:1947710/1‑69 (MQ=255)                                           ggTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTtata                  >  1:908499/1‑69 (MQ=255)                                           ggTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTtata                  >  1:2231693/1‑69 (MQ=255)                                           ggTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTtata                  >  1:2147206/1‑69 (MQ=255)                                           ggTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTtata                  >  1:2515628/1‑69 (MQ=255)                                             tCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATAcc                <  1:1783625/69‑1 (MQ=255)                                               tGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGa              <  1:301596/69‑1 (MQ=255)                                               tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCg               >  1:230725/1‑68 (MQ=255)                                               tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCg               >  1:2233681/1‑68 (MQ=255)                                               tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa              >  1:3326747/1‑69 (MQ=255)                                               tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa              >  1:3238365/1‑69 (MQ=255)                                               tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa              >  1:2553891/1‑69 (MQ=255)                                               tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa              >  1:3371099/1‑69 (MQ=255)                                               tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa              >  1:1663494/1‑69 (MQ=255)                                               tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa              >  1:1335231/1‑69 (MQ=255)                                               tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa              >  1:62038/1‑69 (MQ=255)                                               tGGGGCTATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa              >  1:847866/1‑69 (MQ=255)                                               tGGGGCTATACCCGCATTACACAACGGCTGTCGATAATATGGTTCGCGTGAAAGATATCTTCATACCGa              >  1:2431865/1‑69 (MQ=255)                                               tGGGGCGATACCCGGATTACACAACGGCTGTCGATAAAATGGTTCGCGTGAAAGATATCTTCATACCGa              >  1:3020/1‑69 (MQ=255)                                                 gggCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGa              >  1:2714913/1‑67 (MQ=255)                                                      ataCCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTgaga        >  1:1582180/1‑68 (MQ=255)                                                       taCCCGGATTACGCAACGGCTGTCGATAACATGGTTCGGGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:2464983/69‑1 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:388596/67‑1 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:363335/67‑1 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:984664/67‑1 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:3327540/67‑1 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:2844417/67‑1 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:3371270/67‑1 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:726139/67‑1 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:1143153/67‑1 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:2655540/67‑1 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:1425078/67‑1 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:1512264/67‑1 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:2139977/67‑1 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:168967/67‑1 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:2761953/67‑1 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:3381351/1‑69 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:1360916/1‑69 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:1481573/1‑69 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:1524093/1‑69 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:1675949/1‑69 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:2826826/1‑69 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:2086497/1‑69 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:3274752/1‑69 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:268767/1‑69 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:221201/1‑69 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:205960/1‑69 (MQ=255)                                                         cccGGATTACGCAACGGCTGTCGATAACATGGATCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:2132856/1‑69 (MQ=255)                                                         ccaggATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCaa    >  1:2277287/4‑69 (MQ=255)                                                         cacGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGc      <  1:931381/65‑1 (MQ=255)                                                           cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg  <  1:434827/69‑1 (MQ=255)                                                           cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg  <  1:1271559/69‑1 (MQ=255)                                                           cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg  <  1:1395448/69‑1 (MQ=255)                                                           cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg  <  1:2583056/69‑1 (MQ=255)                                                           cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg  <  1:2248126/69‑1 (MQ=255)                                                           cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg  <  1:1616353/69‑1 (MQ=255)                                                           cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg  <  1:349137/69‑1 (MQ=255)                                                           cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg  <  1:1833466/69‑1 (MQ=255)                                                           cGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATg  <  1:200735/69‑1 (MQ=255)                                                                   |                                                             AACGGTTGTGCAAGCCTGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGATTGAGAGCAATG  >  minE/1872590‑1872716

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