Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	859,070	G→A	Y138Y (TAC→TAT)	ybiU ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	859,070	0	G	A	78.1%	62.6 / 14.4	32	Y138Y (TAC→TAT)	ybiU	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (7/0);  new base A (20/5);  total (27/5)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

TTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC  >  W3110S.gb/859014‑859117                                                         |                                                ttCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGATAAATCTCGGGACg                                   >  1:1111032/1‑71 (MQ=255) ttCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGATAAATCTCGGGACg                                   >  1:827861/1‑71 (MQ=255) ttCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGATAAATCTCGGGACg                                   >  1:616191/1‑71 (MQ=255) ttCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGATAAATCTCGGGACg                                   >  1:463963/1‑71 (MQ=255) ttCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGATAAATCTCGGGACg                                   >  1:3684426/1‑71 (MQ=255) ttCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGATAAATCTCGGGACg                                   >  1:2955765/1‑71 (MQ=255)   cgcgATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGATAAATCTCGGGACGt                                  >  1:980722/1‑70 (MQ=255)        tttCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGATAAATCTCGGGACGTGAAgcg                            <  1:549674/71‑1 (MQ=255)        tttCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGATAAATCTCGGGACGTGAAgcg                            <  1:3504525/71‑1 (MQ=255)          tcttcACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGATAAAt                                            >  1:1138419/1‑53 (MQ=255)          tcttcACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGATAAAt                                            >  1:2828274/1‑53 (MQ=255)          tcttcACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGATAAAt                                            >  1:1953177/1‑53 (MQ=255)              cACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGATAAATCTCGGGACGTGAAGCGCTGAGc                      >  1:3777658/1‑71 (MQ=255)                  ctGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGATAAATCTCGGGACGTGAAGCGCTGAGCGTcc                  >  1:844121/1‑71 (MQ=255)                  ctGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGATAAATCTCGGGACGTGAAGCGCTGAGCGTcc                  >  1:2388077/1‑71 (MQ=255)                  ctGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGATAAATCTCGGGACGTGAAGCGCTGAGCGTcc                  >  1:3222629/1‑71 (MQ=255)                  ctGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGATAAATCTCGGGACGTGAAGCGCTGAGCGTcc                  >  1:4068690/1‑71 (MQ=255)                  ctGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGATAAATCTCGGGACGTGAAGCGCTGAGCGTcc                  >  1:1874177/1‑71 (MQ=255)                  ctGGCGGGCCTGCATTTCCGCCTGCGACCAGTAGATGGGATAAATCTCGGGACGTGAAGCGCTGAGCGTcc                  >  1:137803/1‑71 (MQ=255)                                  tgcgcctgcgccCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAAt         >  1:3283495/1‑64 (MQ=255)                                  tGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAAt         >  1:2848394/1‑64 (MQ=255)                                  tGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAAt         >  1:960410/1‑64 (MQ=255)                                  tGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAAt         >  1:690319/1‑64 (MQ=255)                                  tGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAAt         >  1:43802/1‑64 (MQ=255)                                  tGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAAt         >  1:189546/1‑64 (MQ=255)                                  tGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAAt         >  1:2884009/1‑64 (MQ=255)                                  tGCGCCTGCGACCAGTAGATGGGATAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGcc  >  1:504542/1‑71 (MQ=255)                                  tGCGCCTGCGACCAGTAGATGGGATAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGcc  >  1:1521657/1‑71 (MQ=255)                                  tGCGCCTGCGACCAGTAGATGGGATAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGcc  >  1:1379622/1‑71 (MQ=255)                                    cgcCTGCGACCAGTAGATGGGATAAATCTCGGGACGTg                                 <  1:759320/38‑1 (MQ=255)                                    cgcCTGCGACCAGTAGATGGGATAAATCTCGGGACGTg                                 <  1:2452601/38‑1 (MQ=255)                                    cgcCTGCGACCAGTAGATGGGATAAATCTCGGGACGTg                                 <  1:1341507/38‑1 (MQ=255)                                                         |                                                TTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC  >  W3110S.gb/859014‑859117

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