Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,096,926	A→G	L216P (CTG→CCG)	glpF ←	glycerol facilitator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,096,926	0	A	G	97.1%	125.9 / ‑3.4	34	L216P (CTG→CCG)	glpF	glycerol facilitator
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (18/15);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.58e-01

GCCAAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAGAT  >  REL606/4096856‑4096997                                                                       |                                                                        gcaaaaaaGCGCCGCGAGGAAGTAAGGAATGTCTCGGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGccggccag                                                                    <  1:364038/73‑1 (MQ=255)  caaaaaaGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCGGCCAgg                                                                   <  2:814171/74‑1 (MQ=255)  caaaaaaGCGGCACCAGGAAGTAAGGAATGTCTCGGCCGCCGGTAAAGGCGACATTGCCCCAGCACGCCGGCCAgg                                                                   <  1:344274/74‑1 (MQ=255)   cAAAAAGCGGCACCAGAAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCGAGCCCGCCAGCCAGGc                                                                  <  2:1011493/76‑1 (MQ=255)    aaaaaGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCGGCCAGGCa                                                                 >  2:463591/1‑76 (MQ=255)        aggggCACCAGGAATAAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCGGCCAGGCAAAGa                                                             <  1:168592/73‑1 (MQ=255)                 aGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCGGCCAGGCAAAGGCTTTCGGAc                                                    <  2:549987/76‑1 (MQ=255)                       tAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGt                                              >  2:217734/1‑76 (MQ=255)                            aaTGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCAcgc                                         <  1:252908/76‑1 (MQ=255)                               gTCTCTGNCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCt                                      >  2:980490/1‑76 (MQ=255)                                tcncCGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTTACGCCNTg                                     <  2:993565/76‑1 (MQ=255)                                 ntctGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTgg                                    >  1:352197/2‑76 (MQ=255)                                 ctctGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTgg                                    >  2:303889/1‑76 (MQ=255)                                 ctctGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTgg                                    <  2:1061241/76‑1 (MQ=255)                                   ctGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGt                                  >  2:498483/1‑76 (MQ=255)                                   ctGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGt                                  >  2:340335/1‑76 (MQ=255)                                   ctGCCGCCGCTAAAGGCGACATTGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGt                                  >  2:212433/1‑76 (MQ=255)                                     gccgccGGTAAAGGCGACATTGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTc                                >  2:121964/1‑76 (MQ=255)                                     gccgccGGTAAAGGCGACATTGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTc                                >  2:431510/1‑76 (MQ=255)                                        gccgGTAAAGGCGACATTGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATg                             >  2:870614/1‑76 (MQ=255)                                         ccgGTAAAGGCGACATTGCCCCAGCNCTCCGGCCAGGCAAAGACTTANGGAGCGAAGTCACGCGCTGGGTTCATgg                            <  1:175082/76‑1 (MQ=255)                                          cgGTAAAGGCGACATTGCCCCAGCCCGCCGGCCGGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGc                           >  2:666982/1‑76 (MQ=255)                                           ggTAAAGGCGACATTGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCa                          >  1:53335/1‑76 (MQ=255)                                              aaaGGCGACATTGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCaaaa                       >  2:442567/1‑76 (MQ=255)                                                  gCGACATTGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCg                   <  2:493108/76‑1 (MQ=255)                                                    gACATTGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTc                 >  1:559629/1‑76 (MQ=255)                                                    gACATTGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTc                 >  2:747617/1‑76 (MQ=255)                                                     aCATTGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCa                <  1:748347/76‑1 (MQ=255)                                                     aCATTGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCa                <  1:870614/76‑1 (MQ=255)                                                        ttGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATg             >  1:485101/1‑76 (MQ=255)                                                         tGCCCCAGCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATgg            >  2:403854/1‑76 (MQ=255)                                                            cccACCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGcc         <  2:630718/76‑1 (MQ=255)                                                                gCCCGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATa     <  2:625983/76‑1 (MQ=255)                                                                 cccGCCGGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAg    <  1:618865/76‑1 (MQ=255)                                                                  ccgccgGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAACGGGCCCATAGa   >  1:332568/1‑76 (MQ=255)                                                                   cGCCGGCCAGGCAAAGACTTTCGGACCGAGGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAGAt  <  2:390770/76‑1 (MQ=255)                                                                       |                                                                        GCCAAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAGAT  >  REL606/4096856‑4096997

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