Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,314,904	A→G	7.6%	Q74Q (CAA→CAG)	glpK →	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,314,904	0	A	G	7.6%	192.9 / 5.5	66	Q74Q (CAA→CAG)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base A (38/23);  new base G (2/3);  total (40/26)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAAATCTGGGCCACCCAAAGCTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTATTGTCTGGGAAAAAGAAACCGGCA  >  Exported/2314835‑2314968                                                                      |                                                                 ggAAATCTGGGCCACCCAAAGCTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCaaa                                                                 >  1:1532599/1‑71 (MQ=255)  gAAATCTGGGCCACCCAAAGCTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAAt                                                                <  1:3021537/71‑1 (MQ=255)  gAAATCTGGGCCACCCAAAGCTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAAt                                                                <  1:3311009/71‑1 (MQ=255)  gAAATCTGGGCCACCCAAAGCTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAAt                                                                >  1:2871081/1‑71 (MQ=255)       cTGGGCCACCCAAAGCTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCa                                                            >  1:1475165/1‑70 (MQ=255)       cTGGGCCACCCAAAGCTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAg                                                           >  1:632097/1‑71 (MQ=255)               cccAAAGCTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCgg                                                    >  1:2602632/1‑70 (MQ=255)               cccAAAGCTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCgg                                                    >  1:4365802/1‑70 (MQ=255)               cccAAAGCTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCgg                                                    >  1:1980097/1‑70 (MQ=255)                ccAAAGCTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTa                                                  <  1:2870517/71‑1 (MQ=255)                 cAAAGCTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTAt                                                 >  1:3073277/1‑71 (MQ=255)                 cAAAGCTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTAt                                                 >  1:4403366/1‑71 (MQ=255)                 cAAAGCTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTAt                                                 >  1:4540442/1‑71 (MQ=255)                 cAAAGCTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTAt                                                 >  1:1404914/1‑71 (MQ=255)                 cAAAGCTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTAt                                                 >  1:627029/1‑71 (MQ=255)                     gCTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCaaa                                                                 >  1:4407894/1‑51 (MQ=255)                     gCTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCaaa                                                                 >  1:291783/1‑51 (MQ=255)                     gCTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCaaa                                                                 >  1:2723971/1‑51 (MQ=255)                      cTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACg                                             >  1:4960928/1‑70 (MQ=255)                      cTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACg                                             >  1:2616450/1‑70 (MQ=255)                      cTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACg                                             >  1:2708002/1‑70 (MQ=255)                        ccACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGaa                                           >  1:3081111/1‑70 (MQ=255)                        ccACGCTGGTAGAAGTGCTGGCGAAAGCAGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAAc                                          >  1:3606351/1‑71 (MQ=255)                            gCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGc                                      <  1:2017858/71‑1 (MQ=255)                            gCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGc                                      <  1:2180413/71‑1 (MQ=255)                            gCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGc                                      <  1:763977/71‑1 (MQ=255)                            gCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGc                                      <  1:4978095/71‑1 (MQ=255)                                gTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAGATTGCAGCTATCGGTATTACGAACCAGCGTg                                   >  1:1741403/1‑70 (MQ=255)                                gTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAGATTGCAGCTATCGGTATTACGAACCAGCGTg                                   >  1:3605023/1‑70 (MQ=255)                                  aGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCaaa                                                                 >  1:3961951/1‑38 (MQ=255)                                  aGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCaaa                                                                 >  1:4056216/1‑38 (MQ=255)                                  aGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCaaa                                                                 >  1:4347793/1‑38 (MQ=255)                                  aGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCaaa                                                                 >  1:4651310/1‑38 (MQ=255)                                  aGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTAtt                                                >  1:1219666/1‑55 (MQ=255)                                  aGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGaaa                                >  1:1608924/1‑71 (MQ=255)                                  aGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGaaa                                >  1:503332/1‑71 (MQ=255)                                  aGAAGTGCAGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGaaa                                >  1:1621597/1‑71 (MQ=255)                                   gAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGAAAc                               >  1:1746993/1‑71 (MQ=255)                                   gAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGAAAc                               >  1:4004965/1‑71 (MQ=255)                                         cTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTa                                               >  1:1744748/1‑49 (MQ=255)                                         cTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTa                                               >  1:323070/1‑49 (MQ=255)                                         cTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTa                                               >  1:430193/1‑49 (MQ=255)                                         cTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTa                                               >  1:4089615/1‑49 (MQ=255)                                         cTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTAt                         >  1:2212998/1‑71 (MQ=255)                                            gCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCgg                                                    <  1:1703230/41‑1 (MQ=255)                                            gCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCgg                                                    <  1:4852794/41‑1 (MQ=255)                                            gCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCgg                                                    <  1:3969037/41‑1 (MQ=255)                                            gCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCgg                                                    <  1:3914603/41‑1 (MQ=255)                                            gCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCgg                                                    <  1:1452250/41‑1 (MQ=255)                                            gCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCgg                                                    <  1:5143910/41‑1 (MQ=255)                                            gCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCgg                                                    <  1:2649296/41‑1 (MQ=255)                                            gCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTAtt                        <  1:2645188/69‑1 (MQ=255)                                            gCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTAtt                        <  1:4505351/69‑1 (MQ=255)                                              gAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTa                                                  >  1:362879/1‑41 (MQ=255)                                              gAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTa                                                  >  1:2910051/1‑41 (MQ=255)                                              gAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTa                                                  >  1:3136104/1‑41 (MQ=255)                                              gAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTa                                                  >  1:4519238/1‑41 (MQ=255)                                               aaaGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTATTGTCTg                   >  1:4182040/1‑71 (MQ=255)                                               aaaGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTATTGTCTg                   >  1:2954530/1‑71 (MQ=255)                                               aaaGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTATTGTCTg                   >  1:326818/1‑71 (MQ=255)                                                  gCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTATTGTCTGGGa                <  1:3429526/71‑1 (MQ=255)                                                  gCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTATTGTCTGGGa                <  1:322937/71‑1 (MQ=255)                                                    cGATATCAGTTCCTATCAGATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTATTGTCTGGGaaa              <  1:2876508/71‑1 (MQ=255)                                                    cGATATCAGTTCCGATCAGATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTATTGTCTGGGaaa              <  1:1082730/71‑1 (MQ=255)                                                    cGATATCAGTTCCGATCAGATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTATTGTCTGGGaaa              <  1:4240411/71‑1 (MQ=255)                                                    cGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTATTGTCTGTGaa               >  1:3659800/1‑70 (MQ=255)                                                    cGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTATTGTCTGGGaaa              <  1:4387806/71‑1 (MQ=255)                                                    cGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTATTGTCTGGGaaa              <  1:264151/71‑1 (MQ=255)                                                    cGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTATTGTCTGGGaaa              <  1:1571992/71‑1 (MQ=255)                                                    cGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTATTGTCTGGGaa               >  1:4086468/1‑70 (MQ=255)                                                               ccGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTATTGTCTGGGAAAAAGAAACCgg    <  1:4972685/70‑1 (MQ=255)                                                               ccGATCAAATTGAAGCTATCGGTATTACGAACCAGCGTGAAACCACTATTGTCTGGGAAAAAGAAACCgg    <  1:2987286/70‑1 (MQ=255)                                                                    cAAATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTATTGTCTGGGAAAAAGAAACCGGCa  >  1:3562395/1‑67 (MQ=255)                                                                    cAAATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTATTGTCTGGGAAAAAGAAACCGGCa  >  1:143611/1‑67 (MQ=255)                                                                      |                                                                 GGAAATCTGGGCCACCCAAAGCTCCACGCTGGTAGAAGTGCTGGCGAAAGCCGATATCAGTTCCGATCAAATTGCAGCTATCGGTATTACGAACCAGCGTGAAACCACTATTGTCTGGGAAAAAGAAACCGGCA  >  Exported/2314835‑2314968

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