Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,629,040	G→T	49.3%	R2R (CGC→CGA)	ypfG ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,629,040	0	G	T	49.3%	9.6 / 83.0	75	R2R (CGC→CGA)	ypfG	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/6);  new base T (24/13);  total (56/19)
Fisher's exact test for biased strand distribution p-value = 6.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCAC  >  minE/1628975‑1629081                                                                  |                                          gcTGGTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAgcg                                         <  1:13044/68‑1 (MQ=255)   tGGTGCCGCCCACACCAAAGACGTCGGCAACAAAGCG‑AAAAAAAGAGAAGGAAAATGCGATAGCGCATg                                     <  1:596254/69‑1 (MQ=255)    ggTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATg                                     <  1:1288560/69‑1 (MQ=255)    ggTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATg                                     <  1:3250872/69‑1 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACgc                           >  1:449204/1‑68 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:29625/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:1907371/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:3102749/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:873835/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:1983485/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:892829/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:3035912/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:3102940/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:2882371/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:2832538/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:2505191/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:2725065/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:2717076/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:3197360/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:1109099/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:1142622/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:470750/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:985774/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:1382878/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:833812/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:3379347/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:1605384/1‑69 (MQ=255)              acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:3227228/1‑69 (MQ=255)                   aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:2974952/1‑68 (MQ=255)                   aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:3386579/1‑68 (MQ=255)                   aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:2954250/1‑68 (MQ=255)                   aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:2800407/1‑68 (MQ=255)                   aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:509999/1‑68 (MQ=255)                   aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:2796437/1‑68 (MQ=255)                   aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:2722529/1‑68 (MQ=255)                   aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:715590/1‑68 (MQ=255)                   aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:2667286/1‑68 (MQ=255)                   aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:1532940/1‑68 (MQ=255)                   aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:2207634/1‑68 (MQ=255)                   aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:2161193/1‑68 (MQ=255)                   aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:1622843/1‑68 (MQ=255)                   aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGA‑GGTACGAcccc                       >  1:1952535/1‑67 (MQ=255)                            ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:1391063/1‑68 (MQ=255)                            ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:1323694/1‑68 (MQ=255)                            ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:35977/1‑68 (MQ=255)                            ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:1297515/1‑68 (MQ=255)                            ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:1604147/1‑68 (MQ=255)                            ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:3347799/1‑68 (MQ=255)                            ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:3196388/1‑68 (MQ=255)                            ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:1787434/1‑68 (MQ=255)                            ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:1935759/1‑68 (MQ=255)                            ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:719441/1‑68 (MQ=255)                            ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:2248409/1‑68 (MQ=255)                            ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:2288851/1‑68 (MQ=255)                            ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:2398763/1‑68 (MQ=255)                            ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:2465133/1‑68 (MQ=255)                               aacaaAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTcaaaca          <  1:1673946/69‑1 (MQ=255)                               aacaaAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTcaaaca          <  1:1285093/69‑1 (MQ=255)                               aacaaAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTcaaaca          <  1:2445267/69‑1 (MQ=255)                                  aaaGCGAAAAAAAAGAGAAGGAAAATGCGATAgcg                                         >  1:3188361/1‑35 (MQ=255)                                  aaaGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCaaacaaac       >  1:66064/1‑69 (MQ=255)                                  aaaGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCaaacaaac       >  1:2608136/1‑69 (MQ=255)                                  aaaGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCaaacaaac       >  1:3050649/1‑69 (MQ=255)                                  aaaGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCaaacaaac       >  1:2581472/1‑69 (MQ=255)                                             aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:805772/63‑1 (MQ=255)                                             aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:898566/63‑1 (MQ=255)                                             aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:895447/63‑1 (MQ=255)                                             aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:994605/63‑1 (MQ=255)                                             aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:655544/63‑1 (MQ=255)                                             aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:419012/63‑1 (MQ=255)                                             aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:366493/63‑1 (MQ=255)                                             aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:3080040/63‑1 (MQ=255)                                             aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:1671455/63‑1 (MQ=255)                                             aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:1537532/63‑1 (MQ=255)                                             aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:1358291/63‑1 (MQ=255)                                             aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:1113967/63‑1 (MQ=255)                                             aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACAACCCCAGAAGTTCAAACAAACAGcac  <  1:2055067/63‑1 (MQ=255)                                                                  |                                          GCTGGTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCAC  >  minE/1628975‑1629081

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