Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,565,671	G→T	45.2%	V125V (GTC→GTA)	glk ←	glucokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,565,671	0	G	T	45.2%	26.8 / 46.1	55	V125V (GTC→GTA)	glk	glucokinase
Reads supporting (aligned to +/- strand):  ref base G (23/7);  new base T (13/12);  total (36/19)
Fisher's exact test for biased strand distribution p-value = 8.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCAGATGCGCAACCCCAAGCCCCGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCTTTTTTCAGCATCGGGATCGCCATC  >  minE/1565611‑1565731                                                             |                                                             gACCAGATGCGCAACCCCAAGCCCCGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCTGTTc                                                      >  1:518060/1‑69 (MQ=255) gACCAGATGCGCAACCCCAAGCCCCGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGtt                                                       >  1:2896575/1‑68 (MQ=255) gACCAGATGCGCAACCCCAAGCCCCGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTc                                                      >  1:328179/1‑69 (MQ=255)         gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                                                           >  1:3195015/1‑56 (MQ=255)         gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                                                           >  1:21665/1‑56 (MQ=255)         gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                                                           >  1:2352263/1‑56 (MQ=255)         gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                                                           >  1:2407332/1‑56 (MQ=255)         gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                                                           >  1:196746/1‑56 (MQ=255)         gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                                                           >  1:2463883/1‑56 (MQ=255)         gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                                                           >  1:1796469/1‑56 (MQ=255)         gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                                                           >  1:619540/1‑56 (MQ=255)         gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                                                           >  1:292950/1‑56 (MQ=255)         gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                                                           >  1:1382287/1‑56 (MQ=255)         gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                                                           >  1:1271615/1‑56 (MQ=255)         gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                                                           >  1:3271471/1‑56 (MQ=255)         gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTAACTTCTAcc                                                           >  1:2905941/1‑56 (MQ=255)                         cGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCATATg                              <  1:3059752/69‑1 (MQ=255)                         cGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATg                              <  1:1621031/69‑1 (MQ=255)                         cGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATg                              <  1:1504294/69‑1 (MQ=255)                         cGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATg                              <  1:1866578/69‑1 (MQ=255)                         cGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATg                              <  1:2345259/69‑1 (MQ=255)                         cGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATg                              <  1:23268/69‑1 (MQ=255)                         cGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTGCTGCGCCACCAAACTGAATCAGATg                              <  1:1263423/69‑1 (MQ=255)                          gTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGCCCGGTTCTGCGCCACCAAACTGAATCAGATGc                             >  1:2181488/1‑69 (MQ=255)                          gTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATGc                             >  1:2507147/1‑69 (MQ=255)                          gTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATGc                             >  1:2127235/1‑69 (MQ=255)                          gTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATGc                             >  1:1241503/1‑69 (MQ=255)                          gTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATGc                             >  1:1167929/1‑69 (MQ=255)                             ccGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGcca                                               >  1:2305095/1‑48 (MQ=255)                             ccGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATGctct                          >  1:145453/1‑69 (MQ=255)                             ccGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATGctct                          >  1:3017696/1‑69 (MQ=255)                                  aCCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCtttt                       >  1:1892998/1‑67 (MQ=255)                                  aCCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCtttt                       >  1:648223/1‑67 (MQ=255)                                  aCCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCtttt                       >  1:771282/1‑67 (MQ=255)                                  aCCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCttt                        >  1:369734/1‑66 (MQ=255)                                  aCCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCttt                        >  1:393124/1‑66 (MQ=255)                                               aTAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCTTTTTTCAGCATCGGGATc        >  1:2339798/1‑69 (MQ=255)                                               aTAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCTTTTTTCAGCATCGGGATc        >  1:1657158/1‑69 (MQ=255)                                               aTAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCTTTTTTCAGCATCGGGATc        >  1:355376/1‑69 (MQ=255)                                                  ggCTTACCTTCTACCGGTTCTGCGCCACCAAACTGAATCAGATTCTCTTTTTTCAGCATCGGGATCGcc     <  1:1240620/69‑1 (MQ=255)                                                  ggCTTACCTTCTACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCTTTTTTCAGCATCGGGATCGcc     <  1:1096667/69‑1 (MQ=255)                                                  ggCTTACCTTCTACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCTTTTTTCAGCATCGGGATCGcc     <  1:55990/69‑1 (MQ=255)                                                  ggCTTACCTTCTACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCTTTTTTCAGCATCGGGATCGcc     <  1:3010712/69‑1 (MQ=255)                                                  ggCTTACCTTCTACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCTTTTTTCAGCATCGGGATCGcc     <  1:2954451/69‑1 (MQ=255)                                                  ggCTTACCTTCTACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCTTTTTTCAGCATCGGGATCGcc     <  1:2380365/69‑1 (MQ=255)                                                  ggCTTACCTTCTACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCTTTTTTCAGCATCGGGATCGcc     <  1:1977758/69‑1 (MQ=255)                                                  ggCTTACCTTCTACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCTTTTTTCAGCATCGGGATCGcc     <  1:1742523/69‑1 (MQ=255)                                                  ggCTTACCTTCTACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCTTTTTTCAGCATCGGGATCGcc     <  1:1681521/69‑1 (MQ=255)                                                  ggCTTACCTTCTACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCTTTTTTCAGCATCGGGATCGcc     <  1:1432048/69‑1 (MQ=255)                                                  ggCTTACCTTCTACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCTTTTTTCAGCATCGGGATCGcc     <  1:102836/69‑1 (MQ=255)                                                  ggCTTACCTTCTACCGGTTCTGCGCCACCAAACAGAATCAGATGCTCTTTTTTCAGCATCGGGATCGcc     <  1:3322835/69‑1 (MQ=255)                                                     ttACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCTTTTTTCAGCATCGGGATCGCCAt   >  1:1702603/1‑68 (MQ=255)                                                     ttACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCTTTTTTCAGCATCGGGATCGCCATc  >  1:3008437/1‑69 (MQ=255)                                                     ttACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCTTTTTTCAGCATCGGGATCGCCATc  >  1:532373/1‑69 (MQ=255)                                                     ttACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCTTTTTTCAGCATCGGGATCGCCATc  >  1:88690/1‑69 (MQ=255)                                                             |                                                             GACCAGATGCGCAACCCCAAGCCCCGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATGCTCTTTTTTCAGCATCGGGATCGCCATC  >  minE/1565611‑1565731

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑