Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	908,257	C→T	63.8%	S145L (TCA→TTA)	ycjT →	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	908,257	0	C	T	63.8%	33.0 / 33.6	58	S145L (TCA→TTA)	ycjT	predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base C (13/8);  new base T (18/19);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TTTCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTCAGTGCTGATTTCAACAGGCATTGACGCCACGCAAACCAATCACGGTCGCCAACATCTCG  >  minE/908191‑908316                                                                   |                                                            tttCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTc                                                             <  1:2434562/67‑1 (MQ=255) tttCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTc                                                             <  1:2063434/67‑1 (MQ=255) tttCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTc                                                             <  1:729377/67‑1 (MQ=255) tttCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTc                                                             <  1:604092/67‑1 (MQ=255) tttCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTc                                                             <  1:2319744/67‑1 (MQ=255) tttCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTc                                                             <  1:1854875/67‑1 (MQ=255) tttCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTc                                                             <  1:588132/67‑1 (MQ=255) tttCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTc                                                             <  1:2941650/67‑1 (MQ=255) tttCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTc                                                             <  1:1433600/67‑1 (MQ=255) tttCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTc                                                             <  1:1023805/67‑1 (MQ=255) tttCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTc                                                             <  1:348254/67‑1 (MQ=255)   tCGGCAGCCCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGCCGCGTTAGTg                                                         >  1:1220538/1‑69 (MQ=255)   tCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTg                                                         >  1:2492906/1‑69 (MQ=255)   tCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTg                                                         >  1:2793160/1‑69 (MQ=255)   tCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTg                                                         >  1:19864/1‑69 (MQ=255)   tCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTg                                                         >  1:258215/1‑69 (MQ=255)   tCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTg                                                         >  1:2322540/1‑69 (MQ=255)   tCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTg                                                         >  1:2595138/1‑69 (MQ=255)   tCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTg                                                         >  1:2606558/1‑69 (MQ=255)   tCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTg                                                         >  1:2127615/1‑69 (MQ=255)   tCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTg                                                         >  1:64134/1‑69 (MQ=255)   tCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTg                                                         >  1:749408/1‑69 (MQ=255)   tCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTg                                                         >  1:1955616/1‑69 (MQ=255)   tCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTg                                                         >  1:1761054/1‑69 (MQ=255)   tCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTg                                                         >  1:45359/1‑69 (MQ=255)   tCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTg                                                         >  1:844435/1‑69 (MQ=255)   tCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTg                                                         >  1:873415/1‑69 (MQ=255)   tCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTg                                                         >  1:137924/1‑69 (MQ=255)   tCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTg                                                         >  1:3226791/1‑69 (MQ=255)       cAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTCAGTGCTg                                                      <  1:783254/68‑1 (MQ=255)            cAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTCAGTGCTGATTTc                                                 >  1:2916352/1‑68 (MQ=255)              aCTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTCAGTGCTGATTTCAAc                                              >  1:406712/1‑69 (MQ=255)              aCTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTCAGTGCTGATTTCAAc                                              >  1:543639/1‑69 (MQ=255)              aCTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTCAGTGCTGATTTCAAc                                              >  1:457372/1‑69 (MQ=255)              aCTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTCAGTGCTGATTTCAAc                                              >  1:3423896/1‑69 (MQ=255)              aCTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTCAGTGCTGATTTCAAc                                              >  1:2767515/1‑69 (MQ=255)              aCTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTCAGTGCTGATTTCAAc                                              >  1:194561/1‑69 (MQ=255)              aCTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTCAGTGCTGATTTCAAc                                              >  1:822806/1‑69 (MQ=255)              aCTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTCAGTGCTGATTTCAAc                                              >  1:1484994/1‑69 (MQ=255)              aCTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTCAGTGCTGATTTCAAc                                              >  1:936860/1‑69 (MQ=255)                    gCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTCAGTGCTGATTTCAACAGGc                                          <  1:755986/67‑1 (MQ=255)                    gCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTCAGTGCTGATTTCAACAGGc                                          <  1:1112075/67‑1 (MQ=255)                        aTTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTCAGTGCTGATTTCAACAGGCATTGa                                     <  1:2439250/68‑1 (MQ=255)                          tGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTCAGTGCTGATTTCAACAGGCATTGAc                                    <  1:2896513/67‑1 (MQ=255)                          tGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTCAGTGCTGATTTCAACAGGCATTGAc                                    <  1:1121487/67‑1 (MQ=255)                               tGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTGCTGATTTCAACAGGCATTGACGCCACGc                             <  1:587932/69‑1 (MQ=255)                               tGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTGCTGATTTCAACAGGCATTGACGCCACGc                             <  1:783898/69‑1 (MQ=255)                               tGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTGCTGATTTCAACAGGCATTGACGCCACGc                             <  1:447593/69‑1 (MQ=255)                               tGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTGCTGATTTCAACAGGCATTGACGCCACGc                             <  1:2685511/69‑1 (MQ=255)                               tGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTGCTGATTTCAACAGGCATTGACGCCACGc                             <  1:392347/69‑1 (MQ=255)                               tGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTGCTGATTTCAACAGGCATTGACGCCACGc                             <  1:1415177/69‑1 (MQ=255)                               tGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTGCTGATTTCAACAGGCATTGACGCCACGc                             <  1:2689349/69‑1 (MQ=255)                               tGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTGCTGATTTCAACAGGCATTGACGCCACGc                             <  1:2267223/69‑1 (MQ=255)                               tGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTGCTGATTTCAACAGGCATTGACGCCACGc                             <  1:2103187/69‑1 (MQ=255)                               tGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTGCTGATTTCAACAGGCATTGACGCCACGc                             <  1:2066589/69‑1 (MQ=255)                               tGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTGCTGATTTCAACAGGCATTGACGCCACGc                             <  1:2010090/69‑1 (MQ=255)                               tGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTGCTGATTTCAACAGGCATTGACGCCACGc                             <  1:1937791/69‑1 (MQ=255)                               tGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTGCTGATTTCAACAGGCATTGACGCCACGc                             <  1:1450250/69‑1 (MQ=255)                               tGGAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTGCTGATTTCAACAGGCATTGACGCCACGc                             <  1:1604206/69‑1 (MQ=255)                                ggAAATCACTATTACGCCACTGGGCGCCGACGCGTTAGTGCTGATTTCAACAGGCATTGACGCCACGc                             <  1:234867/68‑1 (MQ=255)                                ggAAATCACTATTACGCCACTGGACGCCGACGCGTTAGTGCTGATTTCAACAGGCATTGACGCCACGc                             <  1:1784171/68‑1 (MQ=255)                                      cACTATTACGCCACTGGACGCCGACGCGTTAGTGCTGATTTCAACAGGCATTGACGCCACGc                             <  1:3180760/62‑1 (MQ=255)                                             aCGCCACTGGACGCCGACGCGTCAGTGCTGATTTCAACAGGCATTGACGCCACGCAAACCAATCACgg                >  1:910629/1‑68 (MQ=255)                                             aCGCCACTGGACGCCGACGCGTCAGTGCTGATTTCAACAGGCATTGACGCCACGCAAACCAATCACGGt               >  1:753099/1‑69 (MQ=255)                                                 cACTGGACGCCGACGCGTCAGTGCTGATTTCAACAGGCATTGACGCCACGCAAACCAATCACGGTCGcc           <  1:3290432/69‑1 (MQ=255)                                                     ggACGCCGACGCGTTAGTGCTGATTTCAACAGGCATTGACGCCACGc                             <  1:1732812/47‑1 (MQ=255)                                                     ggACGCCGACGCGTCAGTGCTGATTTCAACAGGCATTGAc                                    <  1:3197994/40‑1 (MQ=255)                                                          ccGACGCGTCAGTGCTGATTTCAACAGGCATTGACGCCACGCAAACCATTCACGGTCGCCAACATCTcg  >  1:1037292/1‑69 (MQ=255)                                                              cgcgTTAGTGCTGATTTCAACAGGCATTGACGCCACGc                             <  1:2923070/38‑1 (MQ=255)                                                                   |                                                            TTTCGGCAGACCAACTGCCGCTCATTGCGCTGGAAATCACTATTACGCCACTGGACGCCGACGCGTCAGTGCTGATTTCAACAGGCATTGACGCCACGCAAACCAATCACGGTCGCCAACATCTCG  >  minE/908191‑908316

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