Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,487,535	G→A	31.2%	P329S (CCT→TCT)	valS ←	valyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,487,535	0	G	A	31.2%	67.0 / 38.1	48	P329S (CCT→TCT)	valS	valyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (15/18);  new base A (11/4);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGGCGCTTTCACG  >  W3110S.gb/4487466‑4487598                                                                      |                                                                cAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAgg                                                                 >  1:1249551/1‑70 (MQ=255) cAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAgg                                                                 >  1:1473346/1‑70 (MQ=255) cAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAgg                                                                 >  1:4043490/1‑70 (MQ=255)  aaGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGa                                                               >  1:4801934/1‑71 (MQ=255)  aaGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGa                                                               <  1:3054709/71‑1 (MQ=255)     cgcgTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGAttt                                                            >  1:464834/1‑71 (MQ=255)     cgcgTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGAttt                                                            >  1:2760022/1‑71 (MQ=255)        gTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGAGAt                                                              <  1:70261/66‑1 (MQ=255)          cAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTg                                                       >  1:2702261/1‑71 (MQ=255)           aaCTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTgg                                                      <  1:3007193/71‑1 (MQ=255)           aaCTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTg                                                       >  1:533317/1‑70 (MQ=255)             cTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTc                                                           <  1:996433/64‑1 (MQ=255)             cTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTc                                                           <  1:4669880/64‑1 (MQ=255)               gCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGc                                                         >  1:2701600/1‑64 (MQ=255)                cGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAAt                                                   <  1:1096857/69‑1 (MQ=255)                cGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTAGAACTCTGCAGGGATTTCGCTGGAATag                                                 <  1:3217378/71‑2 (MQ=255)                 ggCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGAGATTTCGCTg                                                       <  1:3232393/64‑1 (MQ=255)                   cAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGAGATTTCGCTGGa                                                     >  1:4584715/1‑64 (MQ=255)                   cAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGAGATTTCGCTGGa                                                     >  1:1093327/1‑64 (MQ=255)                   cAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGAGATTTCGCTGGa                                                     >  1:2145555/1‑64 (MQ=255)                      cGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCa                                           <  1:1682725/71‑1 (MQ=255)                       gACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGt                                             >  1:1542561/1‑68 (MQ=255)                       gACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGt                                             >  1:3090748/1‑68 (MQ=255)                       gACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGAGATTTCGCTGGAATAAACGt                                             >  1:3927870/1‑68 (MQ=255)                       gACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGAGATTTCGCTGGAATAAACGt                                             >  1:1546655/1‑68 (MQ=255)                       gACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGAGATTTCGCTGGAATAAACGTCAg                                          <  1:716535/71‑1 (MQ=255)                       gACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGAGATTTCGCTGGAATAAACGTCAg                                          <  1:2808773/71‑1 (MQ=255)                            cTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCg                                     >  1:4836699/1‑71 (MQ=255)                               tACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATag                                                 <  1:3436659/56‑2 (MQ=255)                                     cagcaAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCttt                             <  1:475206/70‑1 (MQ=255)                                     cagcaAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCttt                             <  1:4574832/70‑1 (MQ=255)                                     cagcaAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCttt                             <  1:4580768/70‑1 (MQ=255)                                         aaaaCGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGAtt                                       >  1:591419/1‑56 (MQ=255)                                          aaaCGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCttt                             <  1:2596800/65‑1 (MQ=255)                                          aaaCGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCttt                             <  1:1652156/65‑1 (MQ=255)                                          aaaCGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATc                       >  1:130215/1‑71 (MQ=255)                                          aaaCGCTCCAGTTGCTGGAACTCTGCAGAGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATc                       >  1:1491662/1‑71 (MQ=255)                                              gCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAac                   <  1:2538235/71‑1 (MQ=255)                                                tCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTc                                            <  1:3584136/44‑1 (MQ=255)                                                 ccAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCg                                     >  1:946262/1‑50 (MQ=255)                                                  cAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCATATTCGTTACCTTTGGTATCGAacac                 <  1:1740994/69‑1 (MQ=255)                                                  cAGTTTCTGGAACTCTGCAGAGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCt               >  1:4152233/1‑71 (MQ=255)                                                  cAGTTTCTGGAACTCTGCAGAGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCt               >  1:922695/1‑71 (MQ=255)                                                      ttCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGgc           >  1:1511275/1‑71 (MQ=255)                                                          ggAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAac                   >  1:3879471/1‑59 (MQ=255)                                                           gAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGGCGCt        <  1:2870868/69‑1 (MQ=255)                                                                 tGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGGCGCTTTCACg  <  1:1443783/69‑1 (MQ=255)                                                                  gCAGGGATTTCGCTGGAATAAACGTCAGATTCGtt                                   >  1:439594/1‑35 (MQ=255)                                                                   cAGAGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGaa                    >  1:4900065/1‑49 (MQ=255)                                                                   cAGAGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGaa                    >  1:1342705/1‑49 (MQ=255)                                                                   cAGAGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGaa                    >  1:227296/1‑49 (MQ=255)                                                                      |                                                                CAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGGCGCTTTCACG  >  W3110S.gb/4487466‑4487598

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