Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,835,870	C→T	22.9%	A14V (GCT→GTT)	ydjY →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,835,870	0	C	T	22.9%	77.2 / 19.6	40	A14V (GCT→GTT)	ydjY	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/10);  new base T (2/7);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 2.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

CCAAAAAATGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGAATGCGGC  >  W3110S.gb/1835808‑1835930                                                               |                                                             ccAAAAAATGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCAct                                                        >  1:563441/1‑69 (MQ=255) ccAAAAAATGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCAct                                                        >  1:1759821/1‑69 (MQ=255) ccAAAAAATGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCAct                                                        >  1:50475/1‑69 (MQ=255)    aaaaaTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGTTGCACTCTg                                                     <  1:3460044/69‑1 (MQ=255)    aaaaaTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGTTGCACTCTg                                                     <  1:617609/69‑1 (MQ=255)      aaaTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGTTGCACTCTGtt                                                   <  1:2846243/69‑1 (MQ=255)      aaaTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGtt                                                   <  1:2089916/69‑1 (MQ=255)      aaaTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGtt                                                   <  1:1655772/69‑1 (MQ=255)      aaaTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGtt                                                   <  1:2085487/69‑1 (MQ=255)           aaGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTg                                              >  1:1907957/1‑69 (MQ=255)           aaGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTg                                              >  1:2446962/1‑69 (MQ=255)                 aGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGTTGCACTCTGTTTACTGGCTGtt                                        <  1:793801/69‑1 (MQ=255)                 aGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGt                                         <  1:1762221/68‑1 (MQ=255)                  gAACGATGTTGGAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGttt                                                  >  1:2280197/1‑58 (MQ=255)                  gAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGTTTGCACTCTGt                                                    <  1:760134/56‑1 (MQ=255)                  gAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGTTGCACTCTGt                                                    >  1:891974/1‑56 (MQ=255)                  gAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTg                                       <  1:1371672/69‑1 (MQ=255)                  gAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTg                                       <  1:3716015/69‑1 (MQ=255)                  gAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTg                                       >  1:3437560/1‑69 (MQ=255)                  gAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTg                                       <  1:2289362/69‑1 (MQ=255)                  gAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGATGCACTCTGt                                                    >  1:2494184/1‑56 (MQ=255)                         gTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGTTGCACTCTGTTTACTGGCTGTTGCTGGGCt                                <  1:2142694/69‑1 (MQ=255)                            gCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAg                             >  1:741927/1‑69 (MQ=255)                             cAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTa                              <  1:1125312/67‑1 (MQ=255)                                cATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCtt                               >  1:2958065/1‑63 (MQ=255)                                cATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCtt                               >  1:305971/1‑63 (MQ=255)                                  ttattCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTg                                    >  1:2289246/1‑56 (MQ=255)                                     ttCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCg                     >  1:1006046/1‑68 (MQ=255)                                       cAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATc                  <  1:3183810/69‑1 (MQ=255)                                       cAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATc                  <  1:1501078/69‑1 (MQ=255)                                            tCATGGAAAAAAGGACTGGTTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATcaacaa             <  1:2143309/69‑1 (MQ=255)                                            tCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATcaac               >  1:867100/1‑67 (MQ=255)                                              aTGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAag            >  1:468074/1‑68 (MQ=255)                                              aTGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAag            >  1:152692/1‑68 (MQ=255)                                                 gAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATcaacaa             >  1:195735/1‑64 (MQ=255)                                                 gAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGaa        >  1:2616100/1‑69 (MQ=255)                                                     aaaGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAg                             >  1:1085718/1‑44 (MQ=255)                                                     aaaGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTTCGATCAACAAGAGAATgcg    >  1:353675/1‑69 (MQ=255)                                                     aaaGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGAATgcg    >  1:2800203/1‑69 (MQ=255)                                                       aGGACTGGTTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGAATgcggc  >  1:1259033/1‑69 (MQ=255)                                                               |                                                             CCAAAAAATGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGAATGCGGC  >  W3110S.gb/1835808‑1835930

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