Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,984,569	C→T	68.5%	G103G (GGC→GGT)	yggW →	predicted oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,984,569	0	C	T	68.5%	28.2 / 16.1	38	G103G (GGC→GGT)	yggW	predicted oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (7/5);  new base T (17/9);  total (24/14)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATC  >  minE/1984510‑1984623                                                            |                                                       atacgtgcgCGTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGa                                                <  1:1684878/65‑1 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGc                                      >  1:1348131/1‑68 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGc                                      >  1:982328/1‑68 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                     >  1:2328701/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                     >  1:2986658/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                     >  1:3340718/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                     >  1:3348558/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                     >  1:3356916/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                     >  1:380003/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                     >  1:2714020/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                     >  1:434735/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                     >  1:2402524/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                     >  1:2363537/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                     >  1:2197560/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                     >  1:911447/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                     >  1:1433944/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                     >  1:1123495/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                     >  1:1058121/1‑69 (MQ=255)                                     aTTACTATGGAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGtgtg           >  1:2688569/1‑69 (MQ=255)                                     aTTACTATGGAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGtgtg           >  1:1792503/1‑69 (MQ=255)                                     aTTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCTTTGTTGATTATCAGCGTGCTGGtgtg           <  1:1121181/69‑1 (MQ=255)                                     aTTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCTTTGTTGATTATCAGCGTGCTGGtgtg           <  1:584979/69‑1 (MQ=255)                                     aTTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCTTTGTTGATTATCAGCGTGCTGGtgtg           <  1:654522/69‑1 (MQ=255)                                     aTTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCTTTGTTGATTATCAGCGTGCTGGtgtg           <  1:3361234/69‑1 (MQ=255)                                     aTTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCTTTGTTGATTATCAGCGTGCTGGtgtg           <  1:2937181/69‑1 (MQ=255)                                     aTTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCTTTGTTGATTATCAGCGTGCTGGtgtg           <  1:281897/69‑1 (MQ=255)                                      ttACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCTTTGTTGATTATCAGCGTGCTGGtgtg           <  1:1782563/68‑1 (MQ=255)                                       tACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCTTTGTTGATTATCAGCGTGCTGGtgtg           <  1:1998473/67‑1 (MQ=255)                                           aTGGAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGc     <  1:2954730/69‑1 (MQ=255)                                           aTGGAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGc     <  1:2903676/69‑1 (MQ=255)                                           aTGGAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGc     <  1:2439192/69‑1 (MQ=255)                                           aTGGAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGc     <  1:2198648/69‑1 (MQ=255)                                           aTGGAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTCAACCGc     <  1:437877/69‑1 (MQ=255)                                              gAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCAtc  >  1:2957307/1‑69 (MQ=255)                                              gAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCAtc  >  1:3148886/1‑69 (MQ=255)                                              gAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCAtc  >  1:2927444/1‑69 (MQ=255)                                              gAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCAtc  >  1:2923588/1‑69 (MQ=255)                                              gAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCAtc  >  1:1657898/1‑69 (MQ=255)                                                            |                                                       GTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATC  >  minE/1984510‑1984623

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