Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,128,907	C→A	A278A (GCC→GCA)	yhfZ →	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,128,907	0	C	A	81.0%	72.0 / 22.3	43	A278A (GCC→GCA)	yhfZ	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base A (34/0);  total (39/4)
Fisher's exact test for biased strand distribution p-value = 6.25e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

GTCGATGATTACCCAATGCAACAACTGCTGCGTGCCGTTGTAGATAAACACGCCCTGCTCGCCCATCAACAACGGGTAGTGAGTGGGGAACAGGA  >  W3110S.gb/4128872‑4128966                                    |                                                            gTCGATGATTACCCAATGCAACAACTGCTGCGTGCGGTTg                                                         >  1:2599417/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCCGTTGTAGATAAACACGCCCTGCTCg                                    >  1:1595410/1‑61 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCCGTTGTAGATAAACACGCCCTGCTCg                                    >  1:1825548/1‑61 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCCGTTGTAGATAAACACGCCCTGCTCg                                    >  1:19452/1‑61 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGtt                                                          >  1:852638/1‑39 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGtt                                                          >  1:422353/1‑39 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:907587/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:2632541/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:2797766/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:2816103/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:2907532/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:294687/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:3070388/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:379550/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:386987/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:405209/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:895540/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:557216/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:571328/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:619440/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:717526/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:803131/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:1518478/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:1091086/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:1208951/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:1222780/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:1224813/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:1231187/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:1233148/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:1489479/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:2666539/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:1843178/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:1864423/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:195026/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:1967281/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:912529/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:2192570/1‑40 (MQ=255) gTCGATGATTACCCAATGCAACAACTGCTGCGTGCAGTTg                                                         >  1:2597952/1‑40 (MQ=255)                              gcGTGCCGTTGTAGATAAACACGCCCTGCTCGCCCATCAACAACGGGTAGTGAGTGGGGAAc      >  1:592786/1‑62 (MQ=255)                                  gCCGTTGTAGATAAACACGCCCTGCTCGCCCATCAACAACGGGTAGTGAGTGGGGAACAGGa  <  1:2668650/62‑1 (MQ=255)                                  gCCGTTGTAGATAAACACGCCCTGCTCGCCCATCAACAACGGGTAGTGAGTGGGGAACAGGa  <  1:535165/62‑1 (MQ=255)                                  gCCGTTGTAGATAAACACGCCCTGCTCGCCCATCAACAACGGGTAGTGAGTGGGGAACAGGa  <  1:1078995/62‑1 (MQ=255)                                  gCCGTTGTAGATAAACACGCCCTGCTCGCCCATCAACAACGGGTAGTGAGTGGGGAACAGGa  <  1:2169240/62‑1 (MQ=255)                                   ccGTTGTAGATAAACACGCCCTGCTCGCCCATCAACAACGGGTAGTGAGTGGGGAACAGGa  <  1:928155/61‑1 (MQ=255)                                    |                                                            GTCGATGATTACCCAATGCAACAACTGCTGCGTGCCGTTGTAGATAAACACGCCCTGCTCGCCCATCAACAACGGGTAGTGAGTGGGGAACAGGA  >  W3110S.gb/4128872‑4128966

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