Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,279,217	T→A	39.3%	G332G (GGT→GGA) ‡	yeeJ →	adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,279,217	0	T	A	39.3%	13.8 / 22.7	28	G332G (GGT→GGA) ‡	yeeJ	adhesin
Reads supporting (aligned to +/- strand):  ref base T (3/14);  new base A (9/2);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 1.48e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-04

CACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAA  >  minE/1279150‑1279285                                                                    |                                                                     cACGCCCTGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                                                     <  1:3195472/69‑1 (MQ=255) cACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                                                     <  1:2309720/69‑1 (MQ=255) cACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                                                     <  1:3079171/69‑1 (MQ=255) cACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                                                     <  1:3178636/69‑1 (MQ=255) cACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                                                     <  1:434437/69‑1 (MQ=255) cACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                                                     <  1:740351/69‑1 (MQ=255) cACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                                                     <  1:1786963/69‑1 (MQ=255) cACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                                                     <  1:1656299/69‑1 (MQ=255) cACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                                                     <  1:1044041/69‑1 (MQ=255) cACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                                                     <  1:1490463/69‑1 (MQ=255) cACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                                                     <  1:1408253/69‑1 (MQ=255) cACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCCTGGCGGTa                                                                     <  1:2420305/69‑1 (MQ=255) cACGCCCGGCCAATGGCTGGGATGTAAGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                                                     <  1:3160401/69‑1 (MQ=255)          ccAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCt                                                            >  1:78877/1‑69 (MQ=255)                        gTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGtattat                                              >  1:2401024/1‑69 (MQ=255)                        gTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGtatta                                               >  1:2509561/1‑68 (MQ=255)                                         tttACCTGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCg                                          <  1:2665779/54‑1 (MQ=255)                                           tACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTg                            <  1:1130252/68‑1 (MQ=255)                                               tgCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCg                                          <  1:1989603/49‑1 (MQ=255)                                                    ggCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                   >  1:50163/1‑68 (MQ=255)                                                    ggCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                   >  1:820397/1‑68 (MQ=255)                                                    ggCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                   >  1:3154212/1‑68 (MQ=255)                                                    ggCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                   >  1:3063262/1‑68 (MQ=255)                                                    ggCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                   >  1:3040986/1‑68 (MQ=255)                                                    ggCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                   >  1:3003956/1‑68 (MQ=255)                                                    ggCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                   >  1:2092435/1‑68 (MQ=255)                                                    ggCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                   >  1:1980576/1‑68 (MQ=255)                                                    ggCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                   >  1:1811717/1‑68 (MQ=255)                                                                    tAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTaa  <  1:3095696/69‑1 (MQ=255)                                                                    tAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTaa  <  1:3025185/69‑1 (MQ=255)                                                                    tAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTaa  <  1:2666087/69‑1 (MQ=255)                                                                    tAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTaa  <  1:403228/69‑1 (MQ=255)                                                                    tAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTaa  <  1:2542769/69‑1 (MQ=255)                                                                    tAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTaa  <  1:1583538/69‑1 (MQ=255)                                                                    |                                                                     CACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAA  >  minE/1279150‑1279285

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