Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,460,224	T→C	12.2%	S83P (TCA→CCA)	paaG →	acyl‑CoA hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,460,224	0	T	C	12.2%	97.1 / 5.9	41	S83P (TCA→CCA)	paaG	acyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base T (14/22);  new base C (3/2);  total (17/24)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGT  >  W3110S.gb/1460176‑1460294                                                 |                                                                       aaCGATCGTAACGTTGATCCCAAAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATaa                                                     <  1:2749600/68‑1 (MQ=255)  aCGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCtata                                                      <  1:2520179/66‑1 (MQ=255)  aCGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATaa                                                     <  1:276685/67‑1 (MQ=255)  aCGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCa                                                 >  1:3285277/1‑71 (MQ=255)       cGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTgg                                             >  1:735753/1‑70 (MQ=255)       cGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTgg                                             >  1:2839941/1‑70 (MQ=255)       cGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTgg                                             <  1:2008430/70‑1 (MQ=255)        gTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGCCAGTTGAACGTTTCTATAACCCACTGGTa                                           <  1:2706725/71‑1 (MQ=255)        gTAACGTTGATCCCACAGGCCCCGCACAGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTa                                           >  1:2239263/1‑71 (MQ=255)           aCGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCa                                                 >  1:1490980/1‑62 (MQ=255)           aCGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCa                                                 >  1:3978572/1‑62 (MQ=255)              ttGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTg                                                                  <  1:46845/42‑1 (MQ=255)              ttGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTg                                                                  <  1:502898/42‑1 (MQ=255)              ttGATCCCACAGGCCCCGCACCGGATTTAGGGATGCCAGTTGAACGTTTCTATAACCCACTGGTACGTCGc                                     >  1:1717218/1‑71 (MQ=255)              ttGATCCCACAGGCCACGCACCGGATTTAGGGATGCCAGTTGAACGTTTCTATAACCCACTGGTACGTCGc                                     >  1:2656838/1‑71 (MQ=255)                    ccACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCt                                                         <  1:3106432/45‑1 (MQ=255)                    ccACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCt                                                         <  1:556486/45‑1 (MQ=255)                      acaGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCaaa                             <  1:2066645/71‑1 (MQ=255)                      acaGGCCCCGCACCGGATTTAGGGATGCCAGTTGAACGTTTCTATAAcc                                                   <  1:865451/49‑1 (MQ=255)                        aGGCCCCGCACCGGATTTAGGGATGCCAGTTGAACGTTTCTATAACCCa                                                 >  1:1133483/1‑49 (MQ=255)                               gCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCaaa                    <  1:1227010/71‑1 (MQ=255)                               gCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCaa                     >  1:3636329/1‑70 (MQ=255)                                  ccGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAAcc                 <  1:2207420/71‑1 (MQ=255)                                      aTTTAGGGATTTCAGTTGAACGTTTCTATAACCCAc                                                <  1:3008403/36‑1 (MQ=255)                                      aTTTAGGGATGTCAGTTGAACGTTTCTATAACCCAc                                                <  1:1681839/36‑1 (MQ=255)                                      aTTTAGGGATGTCAGTTGAACGTTTCTATAACCCAc                                                <  1:1829112/36‑1 (MQ=255)                                      aTTTAGGGATGTCAGTTGAACGTTTCTATAACCCAc                                                <  1:712111/36‑1 (MQ=255)                                       tttAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTAc                                          <  1:2171238/41‑1 (MQ=255)                                       tttAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTAc                                          <  1:824911/41‑1 (MQ=255)                                       tttAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTAc                                          <  1:733480/41‑1 (MQ=255)                                       tttAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTAc                                          <  1:543804/41‑1 (MQ=255)                                       tttAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTAc                                          <  1:506494/41‑1 (MQ=255)                                       tttAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTAc                                          <  1:1952970/41‑1 (MQ=255)                                       tttAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTAc                                          <  1:204102/41‑1 (MQ=255)                                        ttAGGGATGTCAGTTGAACGTTTCTATAACCCACTg                                              >  1:3747352/1‑36 (MQ=255)                                        ttAGGGATGTCAGTTGAACGTTTCTATAACCCACTg                                              >  1:79828/1‑36 (MQ=255)                                        ttAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAAc                           >  1:3349055/1‑55 (MQ=255)                                        ttAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAAc                           >  1:217892/1‑55 (MQ=255)                                        ttAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAAc                           >  1:3150639/1‑55 (MQ=255)                                         tAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATc          >  1:3707164/1‑71 (MQ=255)                                         tAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATc          >  1:3373615/1‑71 (MQ=255)                                                 tCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGt  >  1:4020371/1‑71 (MQ=255)                                                 tCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGt  >  1:2770267/1‑71 (MQ=255)                                                 tCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGt  >  1:3612629/1‑71 (MQ=255)                                                 tCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGt  >  1:1933366/1‑71 (MQ=255)                                                 tCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGt  >  1:2288170/1‑71 (MQ=255)                                                 |                                                                       AACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGT  >  W3110S.gb/1460176‑1460294

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