Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,096,697	(T)8→7	coding (8/1509 nt)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,096,690	0	T	.	94.1%	143.9 / 3.6	34	coding (15/1509 nt)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base . (15/17);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTT  >  REL606/4096616‑4096757                                                                           |                                                                    cTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAtt                                                                    <  1:33008/76‑1 (MQ=255) cTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAtt                                                                    >  2:1559768/1‑76 (MQ=255)  tAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAttt                                                                   >  2:842331/1‑76 (MQ=255)   aaTGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAtttt                                                                  <  2:870290/76‑1 (MQ=255)           tGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTc                                                         <  1:1270221/76‑1 (MQ=255)            ggCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCa                                                        <  1:804129/76‑1 (MQ=255)                    gATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAAt                                                >  2:1114872/1‑76 (MQ=255)                     aTCCATTACGACCGCGCGGGAGCGGGATGTGCCCTGGTCGAGCGCTACGATGGCTTTTTTTTCCCTCCTGTTTAAt                                                >  2:653298/1‑76 (MQ=255)                       ccATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGt                                             >  1:390012/1‑76 (MQ=255)                       ccATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGt                                             <  1:67259/76‑1 (MQ=255)                        cATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGt                                             >  1:1116870/1‑76 (MQ=255)                        cATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTc                                            >  2:1251630/1‑76 (MQ=255)                            aCGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCTGt                                        >  2:197518/1‑76 (MQ=255)                            aCGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCga                                        <  2:1411931/76‑2 (MQ=255)                            aCGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGt                                        <  2:976665/76‑1 (MQ=255)                                ccGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTc                                    <  2:1286140/76‑1 (MQ=255)                                   cgcgGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCata                                 <  2:451498/76‑1 (MQ=255)                                     cgGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATAtt                               >  2:583223/1‑76 (MQ=255)                                     cgGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATAtt                               <  1:415679/76‑1 (MQ=255)                                     cgGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATAtt                               >  2:379219/1‑76 (MQ=255)                                       ggAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTAc                             <  2:79674/76‑1 (MQ=255)                                            tggtggTGCCCTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGa                        >  1:1449957/1‑76 (MQ=255)                                               tggtgCCCTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGc                     <  2:927227/76‑1 (MQ=255)                                                   gCCCTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCtttt                 >  1:1523982/1‑76 (MQ=255)                                                    cccTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGGGAAGCtttta                <  2:1248764/76‑2 (MQ=255)                                                    cccTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCttttt                >  1:985678/1‑76 (MQ=255)                                                     ccTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTg               >  1:1471236/1‑76 (MQ=255)                                                     ccTGGTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTg               <  2:853036/76‑1 (MQ=255)                                                         gTCGAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCt           <  1:422962/76‑1 (MQ=255)                                                            gAGCGCAACGATATATTTTTTTTCAGTCATGGTTAATTGTCCCGTAGTCATATTACAGCGAAGGTTTTTGTTCTGa         <  2:72236/76‑1 (MQ=255)                                                            gAGCGCAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGaa        >  1:444102/1‑76 (MQ=255)                                                              gcgcAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAgg      <  2:736794/76‑1 (MQ=255)                                                              gcgcAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAgg      >  2:921547/1‑76 (MQ=255)                                                              gcgcAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAgg      <  2:985678/76‑1 (MQ=255)                                                              gcgcAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCAGATTACAGCGAAGCTTTTTGTTCTGAAgg      <  2:348933/76‑1 (MQ=255)                                                                 cAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGCAGGAGt   >  2:1408358/1‑76 (MQ=255)                                                                 cAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGt   <  1:655670/76‑1 (MQ=255)                                                                 cAACGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGt   >  2:1015825/1‑76 (MQ=255)                                                                  aaCGATATA‑TTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGtt  >  2:1557942/1‑76 (MQ=255)                                                                   aCGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGtt  >  2:665197/1‑76 (MQ=255)                                                                           |                                                                    CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTT  >  REL606/4096616‑4096757

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