Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	125,135	G→A	32.3%	T277T (ACC→ACT)	gcd ←	glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	125,135	0	G	A	32.3%	60.8 / 34.2	71	T277T (ACC→ACT)	gcd	glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (25/23);  new base A (18/5);  total (43/28)
Fisher's exact test for biased strand distribution p-value = 4.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

CAGTCGACCATCATTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTT  >  minE/125069‑125201                                                                   |                                                                   cAGTCGACCATCATTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCgct                                                                  <  1:2505827/69‑3 (MQ=255) cAGTCGACCATCATTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGt                                                                  <  1:3317262/69‑1 (MQ=255) cAGTCGACCATCATTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGt                                                                  <  1:1147323/69‑1 (MQ=255) cAGTCGACCATCATTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGt                                                                  <  1:2924934/69‑1 (MQ=255) cAGTCGACCATCATTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGt                                                                  <  1:2921500/69‑1 (MQ=255) cAGTCGACCATCATTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGt                                                                  <  1:282394/69‑1 (MQ=255) cAGTCGACCATCATTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGt                                                                  <  1:1582702/69‑1 (MQ=255) cAGTCGACCATCATTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGt                                                                  <  1:1556581/69‑1 (MQ=255) cAGTCGACCATCATTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGt                                                                  <  1:511879/69‑1 (MQ=255)  aGTCGACCATCATTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGtt                                                                 >  1:2452261/1‑69 (MQ=255)  aGTCGACCATCATTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGtt                                                                 >  1:1391652/1‑69 (MQ=255)  aGTCGACCATCATTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGt                                                                  <  1:2001195/68‑1 (MQ=255)     cGACCATCATTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCt                                                              >  1:1057412/1‑69 (MQ=255)     cGACCATCATTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCt                                                              >  1:52634/1‑69 (MQ=255)     cGACCATCATTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCt                                                              >  1:1271888/1‑69 (MQ=255)        ccATCATTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGc                                                            <  1:1588631/68‑1 (MQ=255)        ccATCATTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGc                                                            <  1:232767/68‑1 (MQ=255)         catcatTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGc                                                            <  1:2890381/67‑1 (MQ=255)          atcatTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCttt                                                         <  1:598109/69‑1 (MQ=255)            catTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTgg                                                       <  1:3151725/69‑1 (MQ=255)             atTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTgg                                                       >  1:2352222/1‑68 (MQ=255)             atTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTgg                                                       >  1:2163913/1‑68 (MQ=255)             atTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGc                                                      >  1:2865763/1‑69 (MQ=255)             atTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGc                                                      >  1:1327265/1‑69 (MQ=255)             atTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGc                                                      >  1:1961354/1‑69 (MQ=255)             atTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAACCGGTTTCTGCTTTGGc                                                      >  1:2895791/1‑69 (MQ=255)             atTGACCGGAAGAATGATACGACGCGGGAAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGc                                                      >  1:2310789/1‑69 (MQ=255)               tGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCtt                                                    <  1:1258192/69‑1 (MQ=255)                gACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCtt                                                    <  1:2098415/68‑1 (MQ=255)                  ccGGAAGAATGATACGACGCGGGCAGTCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCAt                                                 >  1:596704/1‑69 (MQ=255)                  ccGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCa                                                  >  1:1299427/1‑68 (MQ=255)                  ccGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCAt                                                 >  1:122068/1‑69 (MQ=255)                   cGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCAt                                                 <  1:1721623/68‑1 (MQ=255)                    ggAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCtt                                                    <  1:1201923/64‑1 (MQ=255)                        gaaTGATACGACGCGGTCAATCCGCCATCACTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAag                                           <  1:1096591/69‑1 (MQ=255)                        gaaTGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAag                                           <  1:2080651/69‑1 (MQ=255)                        gaaTGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAag                                           <  1:1251335/69‑1 (MQ=255)                        gaaTGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAag                                           <  1:1680048/69‑1 (MQ=255)                           tGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCgg                                                                   >  1:571293/1‑42 (MQ=255)                           tGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCgg                                                                   >  1:671205/1‑42 (MQ=255)                              tACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGt                                                                  <  1:2535550/40‑1 (MQ=255)                               acgacgCGGGCAATCCGCCATCACTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAag                                           <  1:3355430/62‑1 (MQ=255)                                  acgCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACg                                 <  1:3251571/69‑1 (MQ=255)                                    gcgGGAAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGc                               <  1:2645657/69‑1 (MQ=255)                                    gcgGGAAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGc                               <  1:3426939/69‑1 (MQ=255)                                                    cACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTa                         <  1:2257107/59‑1 (MQ=255)                                                     aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa              >  1:690617/1‑69 (MQ=255)                                                     aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa              >  1:3279672/1‑69 (MQ=255)                                                     aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa              >  1:1441309/1‑69 (MQ=255)                                                     aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa              >  1:2828240/1‑69 (MQ=255)                                                     aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa              >  1:1713968/1‑69 (MQ=255)                                                     aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa              >  1:602357/1‑69 (MQ=255)                                                     aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa              >  1:1923513/1‑69 (MQ=255)                                                     aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa              >  1:1994441/1‑69 (MQ=255)                                                     aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaa               >  1:2847457/1‑68 (MQ=255)                                                     aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag              >  1:337114/1‑68 (MQ=255)                                                     aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag              >  1:588566/1‑68 (MQ=255)                                                     aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag              >  1:3350012/1‑68 (MQ=255)                                                     aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag              >  1:774835/1‑68 (MQ=255)                                                     aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag              >  1:324447/1‑68 (MQ=255)                                                     aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag              >  1:3083329/1‑68 (MQ=255)                                                     aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag              >  1:2992560/1‑68 (MQ=255)                                                     aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag              >  1:2896890/1‑68 (MQ=255)                                                     aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag              >  1:2624873/1‑68 (MQ=255)                                                     aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag              >  1:2541583/1‑68 (MQ=255)                                                     aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag              >  1:2463468/1‑68 (MQ=255)                                                     aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag              >  1:2066344/1‑68 (MQ=255)                                                     aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag              >  1:1946691/1‑68 (MQ=255)                                                     aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag              >  1:1910709/1‑68 (MQ=255)                                                     aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag              >  1:134501/1‑68 (MQ=255)                                                     aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag              >  1:1263375/1‑68 (MQ=255)                                                     aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag              >  1:1058307/1‑68 (MQ=255)                                                     aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaa               >  1:2860879/1‑68 (MQ=255)                                                                 gCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGtt  >  1:3145991/1‑69 (MQ=255)                                                                   |                                                                   CAGTCGACCATCATTGACCGGAAGAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTT  >  minE/125069‑125201

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