Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,984,575	A→T	59.1%	V105V (GTA→GTT)	yggW →	predicted oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,984,575	0	A	T	59.1%	17.6 / 41.0	44	V105V (GTA→GTT)	yggW	predicted oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (13/5);  new base T (17/9);  total (30/14)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

GTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGT  >  minE/1984510‑1984641                                                                  |                                                                   atacgtgcgCGTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGa                                                                  <  1:1684878/65‑1 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGc                                                        >  1:1348131/1‑68 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGc                                                        >  1:982328/1‑68 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                                       >  1:1058121/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:2714020/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                                       >  1:3356916/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                                       >  1:2402524/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                                       >  1:2363537/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                                       >  1:2328701/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                                       >  1:380003/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                                       >  1:2197560/1‑69 (MQ=255)           gTTTGCCGCTGGCAGCGGATGCAGAAATTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCt                                                       >  1:434735/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)                                     aTTACTATGGAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGtgtg                             >  1:1792503/1‑69 (MQ=255)                                     aTTACTATGGAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGtgtg                             >  1:2688569/1‑69 (MQ=255)                                     aTTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCTTTGTTGATTATCAGCGTGCTGGtgtg                             <  1:584979/69‑1 (MQ=255)                                     aTTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCTTTGTTGATTATCAGCGTGCTGGtgtg                             <  1:3361234/69‑1 (MQ=255)                                     aTTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCTTTGTTGATTATCAGCGTGCTGGtgtg                             <  1:654522/69‑1 (MQ=255)                                     aTTACGATGGAAGCGAATCCTGGTACGGTTGAAGCCGATCGCTTTGTTGATTATCAGCGTGCTGGtgtg                             <  1:1121181/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:2198648/69‑1 (MQ=255)                                           aTGGAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGc                       <  1:2439192/69‑1 (MQ=255)                                           aTGGAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGc                       <  1:2903676/69‑1 (MQ=255)                                           aTGGAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTCAACCGc                       <  1:437877/69‑1 (MQ=255)                                              gAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCAtc                    >  1:3148886/1‑69 (MQ=255)                                              gAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCAtc                    >  1:2923588/1‑69 (MQ=255)                                              gAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCAtc                    >  1:2957307/1‑69 (MQ=255)                                              gAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCAtc                    >  1:2927444/1‑69 (MQ=255)                                              gAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCAtc                    >  1:1657898/1‑69 (MQ=255)                                                                gTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGt  >  1:3205014/1‑69 (MQ=255)                                                                gTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGt  >  1:2631793/1‑69 (MQ=255)                                                                 tAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTTCAGAGt  >  1:471210/1‑68 (MQ=255)                                                                 tAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGt  >  1:349809/1‑68 (MQ=255)                                                                 tAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGt  >  1:2011556/1‑68 (MQ=255)                                                                 tAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGt  >  1:1635008/1‑68 (MQ=255)                                                                  |                                                                   GTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCGAACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGT  >  minE/1984510‑1984641

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