Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,223,190	G→A	22.2%	V186V (GTG→GTA)	ebgA →	evolved beta‑D‑galactosidase, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,223,190	0	G	A	22.2%	33.4 / 4.8	18	V186V (GTG→GTA)	ebgA	evolved beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (4/10);  new base A (2/2);  total (6/12)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTATTTTGAAGTCTATGTTAACGGTCAGTATGTGGGTTTCAGCAAGGGCAGTCGCCTGACCGCAGAGTTTGACATCAGCGCGATGGTTAAAACCGGCGACAACCTGTTGTGTGTGCGCGTGATGCAGTGGGCGGACTCTACCTACGTGGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGCGATGCCACGCTTTCCTGCGAAGTGGTGCTGG  >  NC_000913/3223043‑3223338                                                                                                                                                    |                                                                                                                                                     cTATTTTGAAGTCTATGTTAACGGTCAGTATGTGGGTTTCAGCAAGGGCAGTCGCCTGACCGCAGAGTTTGACATCAGCGCGATGGTTAAAACCGGCGACAACCTGTTGTGTGTGCGCGTGATGCAGTGGGCGGACTCTACCTACGTGGa                                                                                                                                                    <  2:225767/150‑1 (MQ=255)                                                          gACCGCAGAGTTTGACATCAGCGCGATGGTTAAAACCGGCGACAACCTGTTGTGTGTGCGCGTGATGCAGTGGGCGGACTCTACCTACGTGGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGaaa                                                                                           <  2:157199/150‑1 (MQ=255)                                                             cGCAGAGTTTGACATCAGCGCGATGGTTAAAACCGGCGACAACCTGTTGTGTGTGCGCGTGATGCAGTGGGCGGACTCTACCTACGTGGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAaca                                                                                        >  2:167490/1‑150 (MQ=255)                                                               cAGAGTTTGACATCAGCGCGATGGTTAAAACCGGCGACAACCTGTTGTGTGTGCGCGTGATGCAGTGGGCGGACTCTACCTACGTGGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACAcc                                                                                      <  1:244953/150‑1 (MQ=255)                                                                  agTTTGACATCAGCGCGATGGTTAAAACCGGCGACAACCTGTTGTGTGTGCGCGTGATGCAGTGGGCGGACTCTACCTACGTGGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTaa                                                                                   >  2:138050/1‑150 (MQ=255)                                                                                gcgATGGTTAAAACCGGCGACAACCTGTTGTGTGTGCGCGTGATGCAGTGGGCGGACTCTACCTACGTAGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGAt                                                                     >  2:103440/1‑150 (MQ=255)                                                                                gcgATGGTTAAAACCGGCGACAACCTGTTGTGTGTGCGCGTGATGCAGTGGGCGGACTCTACCTACGTAGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGAt                                                                     >  2:113390/1‑150 (MQ=255)                                                                                    tGGTTAAAACCGGCGACAACCTGTTGTGTGTGCGCGTGATGCAGTGGGCGGACTCTACCTACGTGGAAGACCAGGATATGCGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCa                                                                 <  1:60660/150‑1 (MQ=255)                                                                                          aaaCCGGCGACAACCTGTTGTGTGTGCGCGTGATGCAGTGGGCGGACTCTACCTACGTGGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGc                                                           >  2:57822/1‑150 (MQ=255)                                                                                             ccGGCGACAACCTGTTGTGTGTGCGCGTGATGCAGTGGGCGGACTCTACCTACGTAGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTa                                                        <  1:113390/150‑1 (MQ=255)                                                                                             ccGGCGACAACCTGTTGTGTGTGCGCGTGATGCAGTGGGCGGACTCTACCTACGTAGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTa                                                        <  1:103440/150‑1 (MQ=255)                                                                                                    cAACCTGTTGTGTGTGCGCGTGATGCAGTGGGCGGACTCTACCTCCGTGGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACtt                                                 <  2:32639/150‑1 (MQ=255)                                                                                                        cTGTTGTGTGTGCGCGTGATGCAGTGGGCGGACTCTACCTACGTGGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGAc                                             <  2:170799/150‑1 (MQ=255)                                                                                                        cTGTTGTGTGTGCGCGTGATGCAGTGGGCGGACTCTACCTACGTGGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGAc                                             <  2:226990/150‑1 (MQ=255)                                                                                                                      cgTGATGCAGTGGGCGGACTCTACCTACGTGGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGCGa                               <  2:68751/150‑1 (MQ=255)                                                                                                                      cgTGATGCAGTGGGCGGACACTACCTACGTGGAAGACCAGGATATGTGGTGGTCAGCGGGGAACTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGCGa                               <  2:68760/150‑1 (MQ=255)                                                                                                                                           tacctacGTGGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTACCGCGATGTTTATCTGGTAGGAAAACACCTAACGCATATTATCGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGCGATGCCACGCTTTCCTGCGACgt          >  2:103077/1‑150 (MQ=255)                                                                                                                                                   tGGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGCGATGCCACGCTTTCCTGCGAAGTGGTGCTgg  <  1:184481/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CTATTTTGAAGTCTATGTTAACGGTCAGTATGTGGGTTTCAGCAAGGGCAGTCGCCTGACCGCAGAGTTTGACATCAGCGCGATGGTTAAAACCGGCGACAACCTGTTGTGTGTGCGCGTGATGCAGTGGGCGGACTCTACCTACGTGGAAGACCAGGATATGTGGTGGTCAGCGGGGATCTTCCGCGATGTTTATCTGGTCGGAAAACACCTAACGCATATTAACGATTTCACTGTGCGTACCGACTTTGACGAAGCCTATTGCGATGCCACGCTTTCCTGCGAAGTGGTGCTGG  >  NC_000913/3223043‑3223338

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