Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,422,263	G→A	G186S (GGC→AGC)	yfcH →	epimerase family protein YfcH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,422,263	0	G	A	92.0%	71.0 / ‑6.0	26	G186S (GGC→AGC)	yfcH	epimerase family protein YfcH
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (13/10);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 2.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATCTTCAGTACTGGTATTAGG  >  NC_000913/2422042‑2422507                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     ccTCGTCAACGGGCTATTGGGGGGACTTAGGTGAAGTGGTGGTCACCGAAGAGGAACCGCCGCATAACGAATTGACCCATAACCTCTGCGCGCGCTGGGAAGAAATTGCCGGCGGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTTGCTCTGGTCGGCAGTAt                                                                                                                                                                                                                         <  2:290294/250‑1 (MQ=255)   ccGGCAACGGGCTATGATGGTGACTTAGGTGAAGTGGTGGTTCCCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCt                                                                                                                                                                                                                       <  1:319154/250‑1 (MQ=255)     ggaaaagggatattatGGTGACTTAGGTGAAGTGGTGGTTAACGAAGGGGAACCGCCGGATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTACCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGATGCGTACCGGTGTAGAGCTGGCACCGGATGGAGGTATTCTCGGTAAAATGCTGACGCCGTTTAGTCTTGGCATGCGCGGGCCGATTGGCTCTGGTCGGCAGTATCTgg                                                                                                                                                                                                                     >  2:13866/11‑251 (MQ=255)                              ggTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGt                                                                                                                                                                                            >  2:75972/1‑251 (MQ=255)                               gtgaatgggtggtTACCGAAGAGAAACCGCCGCATAACGAATGTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTc                                                                                                                                                                                           <  2:26350/244‑1 (MQ=255)                                               cGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCGGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGGTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGgc                                                                                                                                                                           <  1:215767/251‑1 (MQ=255)                                                               cATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCt                                                                                                                                                           >  2:186627/1‑251 (MQ=255)                                                                           aCCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTgg                                                                                                                                                                            >  2:82746/1‑222 (MQ=255)                                                                           aCCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTgg                                                                                                                                                                            <  1:82746/222‑1 (MQ=255)                                                                                 aatctctGCGCTCTCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCGTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTTTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGGATGCTCTGGCGGCTGGATAACAAGCTGCGGGGGCCATTTAatat                                                                                                                                         <  1:62037/248‑1 (MQ=255)                                                                                                     agaaATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCTTACCGGGGTAGTGCTGGCACCGGATGGTGTTATGCTCGGTAAAATGCGTCCGCCGTTTCGTCTTGTCCTTGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTTGCCTGGATTCATATGGGTGATATGGTCAACGGTATTCTCGGGCTTCTGGAGAACGTGCTCCGCGTGCCATTTAATATGGTTTCGCCCTACCCGGTAc                                                                                                                     <  1:341749/251‑1 (MQ=255)                                                                                                                gccGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGAGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACa                                                                                                           >  1:28655/1‑250 (MQ=255)                                                                                                                      cgcAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTACAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGACCTACCCGGTACGCAATGAACAATTTGcc                                                                                                    >  2:256226/1‑251 (MQ=255)                                                                                                                                          gtgggCCTGCTGCGTCCCGGGGGAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAACTGCTGCCGCCGTTTCGTCTTGGCCTGAGCGTGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATGCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCt                                                                                <  1:244888/247‑1 (MQ=255)                                                                                                                                                     gcggaCCGGTGTAGTGCTGGCACCGGAGGGCGGTATTCTCGGTAAAATTCTGCCTCCGGTTCGGCTTGGCCTGAGCGGGCCGATGGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAAGATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCAGGCGCTCGGTCAGGCGCTGCATCGCCCgg                                                                     <  2:114641/247‑1 (MQ=255)                                                                                                                                                     gcGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGATCTGGTCGGCAGTATCTGGCCTGGAGTCATATTGATGCTATGGTCAACGGTATTCTCTGGCTGCTGGCTAACAAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCACTGAACAATTTACCCATGCGCTCGCTCATGCGCTGCATCGCCCgg                                                                     >  2:128173/1‑251 (MQ=255)                                                                                                                                                        tACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGcc                                                                   >  2:98126/1‑250 (MQ=255)                                                                                                                                                            ggTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCAtttt                                                              >  2:127281/1‑251 (MQ=255)                                                                                                                                                                  gTGCTGGCACCGGATGGCGGTATTCTCGGTAAACTGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTAACAGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAAAAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCAGGCCATTTTGCGCGt                                                        >  1:148768/1‑251 (MQ=255)                                                                                                                                                                      tGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTcccc                                                    >  2:278889/1‑251 (MQ=255)                                                                                                                                                                                                   aTGCTGCCGCCGTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATGGGCg                        >  2:157991/1‑250 (MQ=255)                                                                                                                                                                                                              gTTTCGTCTTGGCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATCTTCAGTAc            >  1:197183/1‑251 (MQ=255)                                                                                                                                                                                                                   gTCTTGGCCTTAGCTGGCCGATTGGCTCTGGGCGTGAGTGTCTGGCCTGGATTCATATTGATGATATGGTCAACGTTATGCTCTGTGTGCTTGATAACTTGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGTTACGCAATGAACAATTGTCCCATGCTCTCGGTCATGCGCTGCAGCGCGCGGCCATTGTGCTCGGCCCCGCGACCGTTATTCTACTGTTAATGGGCGAATCTTCAGTACTGGTa       <  1:173614/251‑1 (MQ=255)                                                                                                                                                                                                                       gggCCTGGGCGTGCCGATTGGCTCTGGTCGGCAGTACCGGGCCGGGAGTCATATGGAGGATATGGTCAACGGTATTCTCTGGCTGCTGGAAAACGAGCGGCGCGGGCCATTTAATAGGGTTTCGCCCTACCCGGTCCGCAAGGAACAATTTGCCCATGCGCTCGTTCATGCGCTGCATCGCCCGGCAATTTTGCGCGGCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATCTTCAGTACTGGTATTAg   <  2:248491/250‑1 (MQ=255)                                                                                                                                                                                                                        ggCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTACTGGGCGAATCTTCAGTACTGGTATTAgg  >  2:287664/1‑251 (MQ=255)                                                                                                                                                                                                                        ggCCTGAGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATCTTCAGTACTGATATTAg   <  1:98126/250‑1 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     GCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTACCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATCTTCAGTACTGGTATTAGG  >  NC_000913/2422042‑2422507

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