Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,976,012	C→T	34.8%	N245N (AAC→AAT)	galP →	D‑galactose transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,976,012	0	C	T	34.8%	28.1 / 21.2	43	N245N (AAC→AAT)	galP	D‑galactose transporter
Reads supporting (aligned to +/- strand):  ref base C (18/10);  new base T (5/10);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 6.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

GAAATCCGTGAAAGTTTGCAGGTTAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAA  >  minE/1975956‑1976080                                                         |                                                                     gAAATCCGTGAAAGTTTGCAGGTTAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAATAGCa                                                                  <  1:1122679/61‑1 (MQ=255) gAAATCCGTGAAAGTTTGCAGGTTAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAATAGCa                                                                  <  1:1184355/61‑1 (MQ=255) gAAATCCGTGAAAGTTTGCAGGTTAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAATAGCa                                                                  <  1:1282473/61‑1 (MQ=255) gAAATCCGTGAAAGTTTGCAGGTTAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAATAGCa                                                                  <  1:2776848/61‑1 (MQ=255) gAAATCCGTGAAAGTTTGCAGGTTAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAATAGCa                                                                  <  1:558507/61‑1 (MQ=255) gAAATCCGTGAAAGTTTGCAGGTTAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAATAGCa                                                                  <  1:1861575/61‑1 (MQ=255) gAAATCCGTGAAAGTTTGCAGGTTAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAATAGCa                                                                  <  1:161132/61‑1 (MQ=255) gAAATCCGTGAAAGTTTGAAGGTTAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAATAGCa                                                                  <  1:2930103/61‑1 (MQ=255)  aaaTCCGTGAAAGTTTGCAGGTTAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAATAGCa                                                                  <  1:969341/60‑1 (MQ=255)          gAAAGTTTGCAGGTTAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAATAGCa                                                                  <  1:976792/52‑1 (MQ=255)                       ttAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTAc                                    >  1:2286101/1‑69 (MQ=255)                       ttAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTAc                                    >  1:1709209/1‑69 (MQ=255)                        tAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACt                                   >  1:3029784/1‑69 (MQ=255)                        tAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACt                                   >  1:3283299/1‑69 (MQ=255)                        tAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACt                                   >  1:2923843/1‑69 (MQ=255)                        tAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACt                                   >  1:498694/1‑69 (MQ=255)                        tAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACt                                   >  1:956258/1‑69 (MQ=255)                        tAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACt                                   >  1:2520609/1‑69 (MQ=255)                        tAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACt                                   >  1:990303/1‑69 (MQ=255)                        tAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACt                                   >  1:1578123/1‑69 (MQ=255)                        tAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACt                                   >  1:139654/1‑69 (MQ=255)                        tAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACt                                   >  1:1983040/1‑69 (MQ=255)                        tAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTCGGCGTACt                                   >  1:1832247/1‑69 (MQ=255)                                   gCTGGGCGCTGTTTAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTaa                        <  1:490759/69‑1 (MQ=255)                                   gCTGGGCGCTGTTTAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTaa                        <  1:1860383/69‑1 (MQ=255)                                         cgcTGTTTAAAGAGAACAGCAACTTCCGCCGCTCGGTGTTCCTTGGCGTACTGTTGCAGGTAATgcagc                  <  1:2929930/69‑1 (MQ=255)                                         cgcTGTTTAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATgcagc                  <  1:2408807/69‑1 (MQ=255)                                         cgcTGTTTAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATgcagc                  <  1:330679/69‑1 (MQ=255)                                         cgcTGTTTAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATgcagc                  <  1:664490/69‑1 (MQ=255)                                              tttAAAGAGAATAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGt                                      >  1:2563306/1‑44 (MQ=255)                                              tttAAAGAGAATAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGt                                      >  1:753353/1‑44 (MQ=255)                                              tttAAAGAGAATAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGt                                      >  1:1183389/1‑44 (MQ=255)                                              tttAAAGAGAATAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGt                                      >  1:2221054/1‑44 (MQ=255)                                              tttAAAGAGAATAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGt                                      >  1:3435213/1‑44 (MQ=255)                                                tAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATgcagc                  <  1:1339565/62‑1 (MQ=255)                                                tAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATgcagc                  <  1:893765/62‑1 (MQ=255)                                                tAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATgcagc                  <  1:1470252/62‑1 (MQ=255)                                                tAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATgcagc                  <  1:3390415/62‑1 (MQ=255)                                                tAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATgcagc                  >  1:483493/1‑62 (MQ=255)                                                tAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATgcagc                  >  1:3398329/1‑62 (MQ=255)                                                tAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATgcagc                  >  1:3143181/1‑62 (MQ=255)                                                tAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATgcagc                  >  1:2035238/1‑62 (MQ=255)                                                tAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATgcagc                  >  1:2408496/1‑62 (MQ=255)                                                         cAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATg    >  1:1197233/1‑67 (MQ=255)                                                         cAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGaa  >  1:154673/1‑69 (MQ=255)                                                         cAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGaa  >  1:463500/1‑69 (MQ=255)                                                         cAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGaa  >  1:3100308/1‑69 (MQ=255)                                                         cAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGaa  >  1:1596191/1‑69 (MQ=255)                                                         |                                                                     GAAATCCGTGAAAGTTTGCAGGTTAAACAGAGTGGCTGGGCGCTGTTTAAAGAGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAA  >  minE/1975956‑1976080

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