Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,565,578	G→T	29.7%	G156G (GGC→GGA)	glk ←	glucokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,565,578	0	G	T	29.7%	115.5 / 46.5	85	G156G (GGC→GGA)	glk	glucokinase
Reads supporting (aligned to +/- strand):  ref base G (16/43);  new base T (5/20);  total (21/64)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

TTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCAACCCCAAGCCCCGTTCCGGCAC  >  minE/1565520‑1565645                                                           |                                                                    tttCGAGGATAATGGCCTCTTCTTCCCTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAg                                                            >  1:2708333/1‑68 (MQ=255) tttCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAg                                                            >  1:2347506/1‑68 (MQ=255) tttCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAg                                                            >  1:382933/1‑68 (MQ=255) tttCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAg                                                            >  1:2399546/1‑68 (MQ=255) tttCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAg                                                            >  1:1365057/1‑68 (MQ=255) tttCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAg                                                            >  1:1686136/1‑68 (MQ=255)       ggATAATGCCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTAccc                                                     <  1:1121353/69‑1 (MQ=255)       ggATAATGCCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTAccc                                                     <  1:3417746/69‑1 (MQ=255)       ggATAATGCCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTAccc                                                     <  1:1280265/69‑1 (MQ=255)       ggATAATGCCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTAccc                                                     <  1:817197/69‑1 (MQ=255)          tAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAAc                                                  >  1:2145787/1‑69 (MQ=255)          tAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAAc                                                  >  1:2402732/1‑69 (MQ=255)           aaTGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACg                                                 >  1:1681805/1‑69 (MQ=255)           aaTGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACg                                                 >  1:1488242/1‑69 (MQ=255)           aaTGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACg                                                 >  1:2725399/1‑69 (MQ=255)           aaTGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACg                                                 >  1:351892/1‑69 (MQ=255)           aaTGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACg                                                 >  1:2393294/1‑69 (MQ=255)             tGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCt                                               >  1:3362642/1‑69 (MQ=255)                 ctctTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATc                                           >  1:1290630/1‑69 (MQ=255)                   cttcttCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGa                                         <  1:3001253/69‑1 (MQ=255)                   cttcttCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGa                                         <  1:1586009/69‑1 (MQ=255)                   cttcttCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGa                                         <  1:1399704/69‑1 (MQ=255)                   cttcttCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGa                                         <  1:1713041/69‑1 (MQ=255)                   cttcttCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGa                                         <  1:820196/69‑1 (MQ=255)                   cttcttCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGa                                         <  1:3056539/69‑1 (MQ=255)                   cttcttCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGa                                         <  1:596853/69‑1 (MQ=255)                        tcACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                     <  1:2507476/68‑1 (MQ=255)                        tcACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                     <  1:2883524/68‑1 (MQ=255)                        tcACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                     <  1:2673303/68‑1 (MQ=255)                        tcACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                     <  1:2689727/68‑1 (MQ=255)                        tcACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                     <  1:274636/68‑1 (MQ=255)                        tcACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                     <  1:2151774/68‑1 (MQ=255)                        tcACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                     <  1:3165762/68‑1 (MQ=255)                        tcACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                     <  1:3273396/68‑1 (MQ=255)                        tcACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                     <  1:783320/68‑1 (MQ=255)                        tcACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                     <  1:102287/68‑1 (MQ=255)                        tcACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                     <  1:1153941/68‑1 (MQ=255)                        tcACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                     <  1:1547865/68‑1 (MQ=255)                        tcACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                     <  1:78777/68‑1 (MQ=255)                        tcACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                     <  1:740655/68‑1 (MQ=255)                        tcACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                     <  1:464794/68‑1 (MQ=255)                        tcACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                     <  1:2216052/68‑1 (MQ=255)                        tcACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                     <  1:451561/68‑1 (MQ=255)                         cACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                     <  1:1079670/67‑1 (MQ=255)                           cTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATGGAc                                  <  1:1731413/68‑1 (MQ=255)                             aTTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATGGAcct                                <  1:1791562/68‑2 (MQ=255)                             aTTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCa                                <  1:585423/68‑1 (MQ=255)                             aTTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCa                                <  1:300854/68‑1 (MQ=255)                             aTTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCa                                <  1:2311838/68‑1 (MQ=255)                             aTTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCa                                <  1:921503/68‑1 (MQ=255)                             aTTCGGCGCAAAATCAACGTGACCGCCTTCACCTGGCAAGCTTACCCAACGCTTATCGACATGGACCa                                <  1:2902152/68‑1 (MQ=255)                              ttCGGCGCAAAATCAACGTGACCGCCTTCGCCTTGCAAGCTTACCCAACGCTTATCGAc                                        <  1:3281707/59‑1 (MQ=255)                              ttCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGAc                                        <  1:547490/59‑1 (MQ=255)                              ttCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGAc                                        <  1:1056818/59‑1 (MQ=255)                              ttCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGAc                                        <  1:395545/59‑1 (MQ=255)                              ttCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGAc                                        <  1:411107/59‑1 (MQ=255)                              ttCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGAc                                        <  1:309047/59‑1 (MQ=255)                              ttCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGAc                                        <  1:2722719/59‑1 (MQ=255)                              ttCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGAc                                        <  1:1691974/59‑1 (MQ=255)                              ttCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGAc                                        <  1:1695964/59‑1 (MQ=255)                              ttCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGAc                                        <  1:1732245/59‑1 (MQ=255)                                cGGCGCAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCaa                                                   <  1:991205/46‑1 (MQ=255)                                cGGCGCAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCaa                                                   <  1:3092653/46‑1 (MQ=255)                                cGGCGCAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCaa                                                   <  1:872286/46‑1 (MQ=255)                                cGGCGCAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCaa                                                   <  1:1221936/46‑1 (MQ=255)                                cGGCGCAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCaa                                                   <  1:766758/46‑1 (MQ=255)                                cGGCGCAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCaa                                                   <  1:1786962/46‑1 (MQ=255)                                cGGCGCAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCaa                                                   <  1:676738/46‑1 (MQ=255)                                cGGCGCAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCaa                                                   <  1:1885774/46‑1 (MQ=255)                                cGGCGCAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCaa                                                   <  1:589903/46‑1 (MQ=255)                                cGGCGCAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCaa                                                   <  1:588800/46‑1 (MQ=255)                                cGGCGCAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCaa                                                   <  1:1965318/46‑1 (MQ=255)                                cGGCGCAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCaa                                                   <  1:230755/46‑1 (MQ=255)                                cGGCGCAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCaa                                                   <  1:2316829/46‑1 (MQ=255)                                cGGCGCAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCaa                                                   <  1:2841620/46‑1 (MQ=255)                                cGGCGCAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCaa                                                   <  1:3404492/46‑1 (MQ=255)                                cGGCGCAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCaa                                                   <  1:3403242/46‑1 (MQ=255)                                    gcAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCa                        >  1:1925955/1‑69 (MQ=255)                                    gcAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCa                        >  1:632015/1‑69 (MQ=255)                                    gcAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCa                        >  1:791056/1‑69 (MQ=255)                                    gcAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCa                        >  1:2662761/1‑69 (MQ=255)                                    gcAAAATCAACGTGACCGCCTTCTCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCa                        >  1:1139821/1‑69 (MQ=255)                                           cAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATg                                     <  1:1957723/49‑1 (MQ=255)                                                   ccGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGAc                                        <  1:88249/38‑1 (MQ=255)                                                          cGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCAACCCCAAGCCCCGTTCCGGCAc  >  1:2381690/1‑69 (MQ=255)                                                           |                                                                    TTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCGCAAAATCAACGTGACCGCCTTCGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCAACCCCAAGCCCCGTTCCGGCAC  >  minE/1565520‑1565645

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