Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,080,941	A→G	50.6%	A212A (GCA→GCG)	exuT →	hexuronate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,080,941	0	A	G	50.6%	‑1.5 / 84.6	93	A212A (GCA→GCG)	exuT	hexuronate transporter
Reads supporting (aligned to +/- strand):  ref base A (38/8);  new base G (39/8);  total (77/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GCTGGTGGTATGGGCAATCGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCTATAAACATCCGCGCGACCAGAAGCATCTG  >  minE/2080881‑2081007                                                             |                                                                   gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:1051113/1‑69 (MQ=255) gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:3176098/1‑69 (MQ=255) gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:2846525/1‑69 (MQ=255) gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:3207102/1‑69 (MQ=255) gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:2837666/1‑69 (MQ=255) gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:1995766/1‑69 (MQ=255) gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:328353/1‑69 (MQ=255) gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:282854/1‑69 (MQ=255) gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:166267/1‑69 (MQ=255) gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:162148/1‑69 (MQ=255) gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:348555/1‑69 (MQ=255) gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:1546807/1‑69 (MQ=255) gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:1518754/1‑69 (MQ=255) gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:1449880/1‑69 (MQ=255) gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:5266/1‑69 (MQ=255) gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:1210730/1‑69 (MQ=255) gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:1268970/1‑69 (MQ=255) gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:2520287/1‑69 (MQ=255) gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:2792876/1‑69 (MQ=255) gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:2742136/1‑69 (MQ=255) gcTGGTGGTATGGGCAATCGTGATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCGTTGAGCtt                                                            >  1:680434/1‑69 (MQ=255)    ggtggtATGGGCAATCGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTAt                                                         >  1:2870623/1‑69 (MQ=255)    ggtggtATGGGCAATCGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTAt                                                         >  1:2484880/1‑69 (MQ=255)    ggtggtATGGGCAATCGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTAt                                                         >  1:1708862/1‑69 (MQ=255)    ggtggtATGGGCAATCGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTAt                                                         >  1:1624467/1‑69 (MQ=255)    ggtggtATGGGCAATCGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTAt                                                         >  1:1441466/1‑69 (MQ=255)    ggtggtATGGGCAATCGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTAt                                                         >  1:508276/1‑69 (MQ=255)    ggtggtATGCGCAATCGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTAt                                                         >  1:645586/1‑69 (MQ=255)             ggCAATCGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAg                                                               >  1:733316/1‑54 (MQ=255)             ggCAATCGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTAt                                                         >  1:329110/1‑60 (MQ=255)             ggCAATCGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGat                                                >  1:2306010/1‑69 (MQ=255)             ggCAATCGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGat                                                >  1:2664561/1‑69 (MQ=255)             ggCAATCGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGat                                                >  1:2143152/1‑69 (MQ=255)             ggCAATCGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGat                                                >  1:2029536/1‑69 (MQ=255)             ggCAATCGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGat                                                >  1:1287096/1‑69 (MQ=255)                   cGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGatggcatg                                          >  1:3164235/1‑69 (MQ=255)                   cGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGatggcatg                                          >  1:2189441/1‑69 (MQ=255)                   cGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGatggcatg                                          >  1:236367/1‑69 (MQ=255)                   cGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGatggcatg                                          >  1:1879581/1‑69 (MQ=255)                   cGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGatggcatg                                          >  1:646260/1‑69 (MQ=255)                   cGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGatggcatg                                          >  1:797075/1‑69 (MQ=255)                   cGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGatggcatg                                          >  1:2729861/1‑69 (MQ=255)                   cGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGatggcat                                           >  1:2048337/1‑68 (MQ=255)                   cGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGatggcat                                           >  1:3363947/1‑68 (MQ=255)                   cGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGatggca                                            >  1:2044476/1‑67 (MQ=255)                       aTGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTg                                      <  1:2861867/69‑1 (MQ=255)                       aTGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTg                                      <  1:936660/69‑1 (MQ=255)                       aTGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTg                                      <  1:2554946/69‑1 (MQ=255)                       aTGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTg                                      <  1:2567277/69‑1 (MQ=255)                       aTGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTg                                      <  1:2448803/69‑1 (MQ=255)                       aTGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTg                                      <  1:1603326/69‑1 (MQ=255)                          cacaGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTGAtt                                   >  1:1285169/1‑69 (MQ=255)                          cacaGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTGAtt                                   >  1:2038891/1‑69 (MQ=255)                          cacaGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTGAtt                                   >  1:1501401/1‑69 (MQ=255)                          cacaGCTGGCAGATGGCATTTATCATATCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTGAt                                    >  1:703349/1‑68 (MQ=255)                                 ggCAGATGGCATTTATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            <  1:2150963/69‑1 (MQ=255)                                 ggCAGATGGCATTTATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            <  1:2314349/69‑1 (MQ=255)                                 ggCAGATGGCATTTATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            <  1:3283621/69‑1 (MQ=255)                                 ggCAGATGGCATTTATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            <  1:3306559/69‑1 (MQ=255)                                 ggCAGATGGCATTTATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            <  1:599852/69‑1 (MQ=255)                                 ggCAGATGGCATTTATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            <  1:716588/69‑1 (MQ=255)                                 ggCAGATGGCATTTATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            <  1:2642115/69‑1 (MQ=255)                                 ggCAGATGGCATTTATCATCTCCGGGGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            <  1:408553/69‑1 (MQ=255)                                         gCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCTATaaa                          >  1:492429/1‑63 (MQ=255)                                         gCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCTATAAACATcc                     >  1:1059280/1‑68 (MQ=255)                                         gCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCTATAAACATCcg                    >  1:2137267/1‑69 (MQ=255)                                         gCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCTATAAACATCcg                    >  1:2339628/1‑69 (MQ=255)                                         gCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCTATAAACATCcg                    >  1:1303971/1‑69 (MQ=255)                                         gCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCTATAAACATCcg                    >  1:2352720/1‑69 (MQ=255)                                         gCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCTATAAACATCcg                    >  1:1428735/1‑69 (MQ=255)                                         gCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCTATAAACATCcg                    >  1:2514745/1‑69 (MQ=255)                                             ttATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            >  1:2932613/1‑57 (MQ=255)                                             ttATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            >  1:2587972/1‑57 (MQ=255)                                             ttATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            >  1:828451/1‑57 (MQ=255)                                             ttATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            >  1:2603762/1‑57 (MQ=255)                                             ttATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            >  1:2724762/1‑57 (MQ=255)                                             ttATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            >  1:1314313/1‑57 (MQ=255)                                             ttATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            >  1:1429971/1‑57 (MQ=255)                                             ttATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            >  1:519012/1‑57 (MQ=255)                                             ttATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            >  1:1559244/1‑57 (MQ=255)                                             ttATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            >  1:2494867/1‑57 (MQ=255)                                             ttATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            >  1:453137/1‑57 (MQ=255)                                             ttATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            >  1:2830264/1‑57 (MQ=255)                                             ttATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            >  1:1830628/1‑57 (MQ=255)                                             ttATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            >  1:2845202/1‑57 (MQ=255)                                             ttATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            >  1:186476/1‑57 (MQ=255)                                             ttATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            >  1:1968673/1‑57 (MQ=255)                                             ttATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            >  1:2285767/1‑57 (MQ=255)                                             ttATCATCTCCGGTGCGTTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCtata                            >  1:3057777/1‑57 (MQ=255)                                                 catcTCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCTATAAACATCCGCGCGa               >  1:847150/1‑66 (MQ=255)                                                         gTGCATTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCTATAAACATCCGCGCGACCAGAAGc       <  1:516311/66‑1 (MQ=255)                                                         gTGCATTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCTATAAACATCCGCGCGACCAGAAGc       <  1:2685730/66‑1 (MQ=255)                                                           gCATTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCTATAAACATCCGCGCGACCAGAAGCATCTg  >  1:2526834/1‑69 (MQ=255)                                                             |                                                                   GCTGGTGGTATGGGCAATCGTAATGCACAGCTGGCAGATGGCATTTATCATCTCCGGTGCATTGAGCTTTATCTGGGCGATGGCATGGCTGATTTTCTATAAACATCCGCGCGACCAGAAGCATCTG  >  minE/2080881‑2081007

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