Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	62,820	G→C	58.7%	A213A (GCC→GCG)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	62,820	0	G	C	58.7%	20.9 / 42.6	46	A213A (GCC→GCG)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base G (15/4);  new base C (19/8);  total (34/12)
Fisher's exact test for biased strand distribution p-value = 7.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

TTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAG  >  minE/62754‑62888                                                                   |                                                                     tttttGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTgg                                                                     >  1:1281197/1‑68 (MQ=255) tttttGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTgg                                                                     >  1:2904841/1‑68 (MQ=255) tttttGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTgg                                                                     >  1:163195/1‑68 (MQ=255)         gcaTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAg                                                             >  1:550040/1‑68 (MQ=255)          caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                                           >  1:385853/1‑69 (MQ=255)          caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                                           >  1:2133030/1‑69 (MQ=255)          caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                                           >  1:2135923/1‑69 (MQ=255)          caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                                           >  1:2173530/1‑69 (MQ=255)          caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                                           >  1:2291475/1‑69 (MQ=255)          caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                                           >  1:2765880/1‑69 (MQ=255)          caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                                           >  1:3132254/1‑69 (MQ=255)          caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                                           >  1:2657457/1‑69 (MQ=255)          caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                                           >  1:53269/1‑69 (MQ=255)          caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                                           >  1:714925/1‑69 (MQ=255)          caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                                           >  1:732424/1‑69 (MQ=255)          caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                                           >  1:955319/1‑69 (MQ=255)          caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                                           >  1:981789/1‑69 (MQ=255)          caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                                           >  1:2092325/1‑69 (MQ=255)          caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                                           >  1:1129254/1‑69 (MQ=255)          caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                                           >  1:1156732/1‑69 (MQ=255)          caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                                           >  1:1917570/1‑69 (MQ=255)          caTTCGCAGATCGAACTGCACCACGTTCCAACAGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                                           >  1:82231/1‑69 (MQ=255)          caTTCGCAGATAGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCATGAGTCGCAAACCAGGc                                                           >  1:2148703/1‑69 (MQ=255)                    tCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGtttt                                                 >  1:1044803/1‑69 (MQ=255)                    tCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGtttt                                                 >  1:1055302/1‑69 (MQ=255)                     cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGTCAAACCAGGCGTTGAGtttt                                                 >  1:2711010/1‑68 (MQ=255)                     cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGtttt                                                 >  1:2062702/1‑68 (MQ=255)                     cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGtttt                                                 >  1:624519/1‑68 (MQ=255)                     cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGtttt                                                 >  1:2774838/1‑68 (MQ=255)                     cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGtttt                                                 >  1:1028482/1‑68 (MQ=255)                     cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGtttt                                                 >  1:2269854/1‑68 (MQ=255)                     cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGtttt                                                 >  1:2344173/1‑68 (MQ=255)                     cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGtttt                                                 >  1:2497394/1‑68 (MQ=255)                     cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGttt                                                  >  1:2578415/1‑67 (MQ=255)                          tGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGc                                           <  1:877215/69‑1 (MQ=255)                          tGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGc                                           <  1:737671/69‑1 (MQ=255)                                acGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTgc                                     >  1:2013876/1‑69 (MQ=255)                                acGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTgc                                     >  1:353229/1‑69 (MQ=255)                                            atgatCACATCAGGATCATGAGTCGCAAACCAGGCGTTAAGTTTTTCCAGTAGCAGCGGGCGGCTGGCg                         <  1:1490062/69‑1 (MQ=255)                                            atgatCACATCAGGATCATGAGTCGCAAACCAGGCGTTAAGTTTTTCCAGTAGCAGCGGGCGGCTGGCg                         <  1:2060673/69‑1 (MQ=255)                                            atgatCACATCAGGATCATGAGTCGCAAACCAGGCGTTAAGTTTTTCCAGTAGCAGCGGGCGGCTGGCg                         <  1:901390/69‑1 (MQ=255)                                            atgatCACATCAGGATCATGAGTCGCAAACCAGGCGTTAAGTTTTTCCAGTAGCAGCGGGCGGCTGGCg                         <  1:2671661/69‑1 (MQ=255)                                            atgatCACATCAGGATCATGAGTCGCAAACCAGGCGTTAAGTTTTTCCAGTAGCAGCGGGCGGCTGGCg                         <  1:3284486/69‑1 (MQ=255)                                            atgatCACATCAGGATCATGAGTCGCAAACCAGGCGTTAAGTTTTTCCAGTAGCAGCGGGCGGCTGGCg                         <  1:1516117/69‑1 (MQ=255)                                            atgatCACATCAGGATCATGAGTCGCAAACCAGGCGTTAAGTTTTTCCAGTAGCAGCGGGCGGCTGGCg                         <  1:3021665/69‑1 (MQ=255)                                            atgatCACATCAGGATCATGAGTCGCAAACCAGGCGTTAAGTTTTTCCAGTAGCAGCGGGCGGCTGGCg                         <  1:28236/69‑1 (MQ=255)                                               atCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACg                      <  1:3272295/69‑1 (MQ=255)                                               atCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACg                      <  1:3270849/69‑1 (MQ=255)                                                 cacaTCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTgcggg                                  >  1:324677/1‑55 (MQ=255)                                                                   ggCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCaa   >  1:1184982/1‑68 (MQ=255)                                                                   ggCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCaa   >  1:265199/1‑68 (MQ=255)                                                                   ggCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAg  <  1:1429901/69‑1 (MQ=255)                                                                   ggCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAg  <  1:561185/69‑1 (MQ=255)                                                                   ggCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAg  <  1:1935552/69‑1 (MQ=255)                                                                   |                                                                     TTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAG  >  minE/62754‑62888

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