Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,097,224	G→A	Q117* (CAG→TAG)	glpF ←	glycerol facilitator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,097,224	0	G	A	97.1%	124.4 / ‑3.0	34	Q117* (CAG→TAG)	glpF	glycerol facilitator
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (20/13);  total (21/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01

GATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCGCAGAAAGCGCCGGCAACT  >  REL606/4097152‑4097294                                                                         |                                                                       gATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTACt                                                                      <  2:1463995/75‑1 (MQ=255)  aTTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTACTc                                                                     >  1:1312977/1‑75 (MQ=255)  aTTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTACTc                                                                     >  1:1562010/1‑75 (MQ=255)    tAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTACTCGa                                                                   <  2:775326/75‑1 (MQ=255)    tAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTACTCGa                                                                   >  1:1359229/1‑75 (MQ=255)    tAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTACTCGa                                                                   <  1:9342/75‑1 (MQ=255)      gggTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTACTCGAAg                                                                 <  1:1365488/75‑1 (MQ=255)       ggTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTACTCGAAGt                                                                <  2:1128250/75‑1 (MQ=255)            gTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTACTCGAAGTCGaaa                                                           >  1:1004218/1‑75 (MQ=255)                 gaAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTACTCGAAGTCGAAaaataa                                                      >  2:1052518/1‑75 (MQ=255)                 gaAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTACTCGAAGTCGAAaaataa                                                      >  2:385252/1‑75 (MQ=255)                  aaaGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTACTCGAAGTCGAAaaataaa                                                     <  2:487644/75‑1 (MQ=255)                   aaGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTACTCGAAGTCGAAaaataaat                                                    >  1:626465/1‑75 (MQ=255)                         cagccagATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTACTCGAAGTCGAAAAATAAATTgtagt                                              >  1:145267/1‑75 (MQ=255)                           gccagATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTACTCGAAGTCGAAAAATAAATTGTAGTaa                                            >  1:1064987/1‑75 (MQ=255)                            ccagATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTACTCGAAGTCGAAAAATAAATTGTAGTaaa                                           <  2:1043639/75‑1 (MQ=255)                             cagATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTACTCGAAGTCGAAAAATAAATTGTAGTAAAg                                          >  1:900773/1‑75 (MQ=255)                                 tCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTACTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCg                                      >  2:329522/1‑75 (MQ=255)                                        tttCAACGCTGCCGCGAACAATGTGATGAGTCTACTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTa                               <  2:590112/75‑1 (MQ=255)                                          tCAACGCTGCCGCGAACAATGTGATGAGTCTACTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTaaa                             >  2:1342718/1‑75 (MQ=255)                                             aCGCTGCCGCGAACAATGTGATGAGTCTACTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGcc                          <  2:884618/75‑1 (MQ=255)                                             aCGCTGCCGCGAACAATGTGATGAGTCTACTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGcc                          <  2:1014966/75‑1 (MQ=255)                                              cGCTGCCGCGAACAATGTGATGAGTCTACTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCcg                         >  2:350195/1‑75 (MQ=255)                                                cTGCCGCGAACAATGTGATGAGTCTACTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCcgca                       >  1:1142924/1‑75 (MQ=255)                                                   ccGCGAACAATGTGATGAGTCTACTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCcgcagcg                    >  1:333410/1‑75 (MQ=255)                                                      cgAACAATGTGATGAGTCTACTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCcgcagcgcag                 >  1:1439620/1‑75 (MQ=255)                                                       gAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCGCAGa                >  1:777622/1‑75 (MQ=255)                                                                gatgaGTCTACTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCGCAGAAAGCGCCgg       >  2:248213/1‑75 (MQ=255)                                                                gatgaGTCTACTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCGCAGAAAGCGCCgg       >  1:228935/1‑75 (MQ=255)                                                                  tgaGTCTACTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCGCAGAAAGCGCCGGCa     <  2:1286599/75‑1 (MQ=255)                                                                  tgaGTCTACTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCGCAGAAAGCGCCGGCa     <  2:1272275/75‑1 (MQ=255)                                                                  tgaGTCTACTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCGCAGAAAGCGCCGGCa     >  1:881254/1‑75 (MQ=255)                                                                   gaGTCTACTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCGCAGAAAGCGCCGGCaa    >  1:10597/1‑75 (MQ=255)                                                                     gTCTACTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCGCAGAAAGCGCCGGCAACt  <  2:1573954/75‑1 (MQ=255)                                                                         |                                                                       GATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCGCAGAAAGCGCCGGCAACT  >  REL606/4097152‑4097294

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