Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,090,490	C→T	17.6%	A234V (GCG→GTG)	yggJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,090,490	0	C	T	17.6%	199.5 / 33.2	85	A234V (GCG→GTG)	yggJ	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (15/55);  new base T (7/8);  total (22/63)
Fisher's exact test for biased strand distribution p-value = 5.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

AATTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGT  >  W3110S.gb/3090421‑3090553                                                                      |                                                                aaTTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACcgcg                                                                >  1:1717657/1‑71 (MQ=255) aaTTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACcgcg                                                                >  1:2236983/1‑71 (MQ=255) aaTTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACcgcg                                                                >  1:1923128/1‑71 (MQ=255) aaTTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACcgcg                                                                >  1:3120494/1‑71 (MQ=255) aaTTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACcgcg                                                                >  1:864163/1‑71 (MQ=255) aaTTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACcgcg                                                                >  1:1562524/1‑71 (MQ=255) aaTTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACcgcg                                                                >  1:1538267/1‑71 (MQ=255) aaTTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACcgcg                                                                >  1:1454231/1‑71 (MQ=255) aaTTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACcgcg                                                                >  1:2905989/1‑71 (MQ=255)     tACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTa                                                             >  1:3128457/1‑70 (MQ=255)     tACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTa                                                             >  1:1560148/1‑70 (MQ=255)     tACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTa                                                             >  1:548299/1‑70 (MQ=255)     tACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTa                                                             >  1:372966/1‑70 (MQ=255)     tACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTa                                                             >  1:1113163/1‑70 (MQ=255)     tACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTa                                                             >  1:1881054/1‑70 (MQ=255)     aaCTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTa                                                             >  1:94220/2‑70 (MQ=255)             tCCTGTTGGGACCTCGCGTTTTGCGTACATAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcg                                                     <  1:1613687/70‑1 (MQ=255)             tCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcga                                                    <  1:3030376/71‑1 (MQ=255)             tCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcga                                                    <  1:985687/71‑1 (MQ=255)             tCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcga                                                    <  1:2718217/71‑1 (MQ=255)             tCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcg                                                     <  1:635825/70‑1 (MQ=255)             tCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcg                                                     <  1:370100/70‑1 (MQ=255)             tCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcg                                                     <  1:2392472/70‑1 (MQ=255)             tCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcg                                                     <  1:1847026/70‑1 (MQ=255)             tCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcg                                                     <  1:3099856/70‑1 (MQ=255)             tCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcg                                                     <  1:2765628/70‑1 (MQ=255)             tCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcg                                                     <  1:2706302/70‑1 (MQ=255)             tCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcg                                                     <  1:129589/70‑1 (MQ=255)             tCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcg                                                     <  1:2982480/70‑1 (MQ=255)             tCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGATACAAGTAcg                                                     <  1:3103570/70‑1 (MQ=255)              ccTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcg                                                     <  1:1234669/69‑1 (MQ=255)                        ccTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcgatttggcgatt                                          >  1:3307129/1‑70 (MQ=255)                        ccTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcgatttggcgatt                                          >  1:601510/1‑70 (MQ=255)                        ccTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcgatttggcgatt                                          >  1:2825562/1‑70 (MQ=255)                        ccTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcgatttggcgatt                                          >  1:202932/1‑70 (MQ=255)                        ccTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcgatttggcgatt                                          >  1:1385212/1‑70 (MQ=255)                        ccTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcgatttggcgatt                                          >  1:1856447/1‑70 (MQ=255)                        ccTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcgatttggcgat                                           >  1:1544366/1‑69 (MQ=255)                        ccTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcgatttggagatt                                          >  1:1979006/1‑70 (MQ=255)                                        cAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGagaag                          <  1:2874177/70‑1 (MQ=255)                                        cAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGagaag                          <  1:606422/70‑1 (MQ=255)                                        cAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGagaag                          <  1:2359859/70‑1 (MQ=255)                                        cAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGagaag                          <  1:1972743/70‑1 (MQ=255)                                        cAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGagaag                          <  1:263181/70‑1 (MQ=255)                                        cAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGagaag                          <  1:1912364/70‑1 (MQ=255)                                                aCTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGgaga                            <  1:2236227/60‑1 (MQ=255)                                                aCTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGgaga                            <  1:980110/60‑1 (MQ=255)                                                aCTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGgaga                            <  1:1205814/60‑1 (MQ=255)                                                aCTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGgaga                            <  1:2659478/60‑1 (MQ=255)                                                aCTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGgaga                            <  1:1686265/60‑1 (MQ=255)                                                aCTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGaaga                            <  1:3289784/60‑1 (MQ=255)                                                              cATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGgaga                            <  1:1843985/46‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGgaga                            <  1:793764/45‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAGTAATGATCAAGCTCGGCATCGt  <  1:985867/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:3123685/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:3150034/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:3161362/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:959475/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:2966112/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:2477918/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:525041/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:471997/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:2807641/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:1154943/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:1167078/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:1399381/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:1436621/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:2203036/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:2213972/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:2225137/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:234782/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:2357541/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:2817887/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:2762656/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:2495473/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:2739329/71‑1 (MQ=255)                                                               aTTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:2602536/71‑1 (MQ=255)                                                                ttACCGTGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCAt     <  1:2603741/67‑1 (MQ=255)                                                                ttACCGTGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCAt     <  1:1158263/67‑1 (MQ=255)                                                                ttACCGTGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCAt     <  1:1740955/67‑1 (MQ=255)                                                                ttACCGTGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCAt     <  1:2455764/67‑1 (MQ=255)                                                                ttACCGTGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCAt     <  1:2280050/67‑1 (MQ=255)                                                                ttACCGTGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCAt     <  1:2471100/67‑1 (MQ=255)                                                                ttACCGTGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCAt     <  1:3158843/67‑1 (MQ=255)                                                                ttACCGTGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCAt     <  1:2452067/67‑1 (MQ=255)                                                                ttACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:2517257/70‑1 (MQ=255)                                                                ttACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:1950512/70‑1 (MQ=255)                                                                  aCCGTGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCg   >  1:382091/1‑67 (MQ=255)                                                                  aCCGTGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCg   >  1:2259790/1‑67 (MQ=255)                                                                  aCCGTGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCg   >  1:2090077/1‑67 (MQ=255)                                                                  aCCGTGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCg   >  1:722518/1‑67 (MQ=255)                                                                  aCCGTGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCg   >  1:2327520/1‑67 (MQ=255)                                                                  aCCGTGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCg   >  1:2780543/1‑67 (MQ=255)                                                                  aCCGTGCTACAAGTACGATTTGGCGATTTGGGATAACGGAGAAGAATAATGATCAAGCTCGGCATCg   >  1:2175205/1‑67 (MQ=255)                                                                      cgcTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGt  <  1:2192746/64‑1 (MQ=255)                                                                      |                                                                AATTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGT  >  W3110S.gb/3090421‑3090553

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: ‑