Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,783,786	T→C	53.0%	R286R (CGT→CGC)	yjbB →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,783,786	0	T	C	53.0%	10.1 / 71.3	85	R286R (CGT→CGC)	yjbB	predicted transporter
Reads supporting (aligned to +/- strand):  ref base T (20/20);  new base C (19/26);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TGGGGAAGTTGTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGCCATTTGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTCGC  >  minE/2783719‑2783846                                                                    |                                                             tGGGGAAGTTGTCATTGCCAAAAGCGGAACTTGTGATCTATTTCCACGTCTTCTACAACCTTGTACGtt                                                             <  1:1900771/69‑1 (MQ=255) tGGGGAAGTTGTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGtt                                                             <  1:3421663/69‑1 (MQ=255) tGGGGAAGTTGTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGtt                                                             <  1:573685/69‑1 (MQ=255) tGGGGAAGTTGTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGtt                                                             <  1:1715347/69‑1 (MQ=255) tGGGGAAGTTGTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGtt                                                             <  1:3407163/69‑1 (MQ=255) tGGGGAAGTTGTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGtt                                                             <  1:2210137/69‑1 (MQ=255) tGGGGAAGTTGTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGtt                                                             <  1:781740/69‑1 (MQ=255) tGGGGAAGTTGTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGtt                                                             <  1:2344423/69‑1 (MQ=255) tGGGGAAGTTGTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGtt                                                             <  1:1896952/69‑1 (MQ=255) tGGGGAAGTTGTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGtt                                                             <  1:1977188/69‑1 (MQ=255) tGGGGAAGTTGTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGtt                                                             <  1:2056047/69‑1 (MQ=255)  ggggAAGTTGTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGtt                                                             <  1:3276661/68‑1 (MQ=255)      aaGTTGTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTg                                                        >  1:1554307/1‑69 (MQ=255)      aaGTTGCCATTGCCAAAAGCGGAATTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGCTGCCTg                                                        >  1:545168/1‑69 (MQ=255)      aaGTTGCCATTGCCAAAAGCGGAATTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGCTGCCTg                                                        >  1:568400/1‑69 (MQ=255)      aaGTTGCCATTGCCAAAAGCGGAATTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGCTGCCTg                                                        >  1:1908790/1‑69 (MQ=255)      aaGTTGCCATTGCCAAAAGCGGAATTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGCTGCCTg                                                        >  1:473680/1‑69 (MQ=255)      aaGTTGCCATTGCCAAAAGCGGAATTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGCTGCCTg                                                        >  1:2171647/1‑69 (MQ=255)      aaGTTGCCATTGCCAAAAGCGGAATTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGCTGCCTg                                                        >  1:1718973/1‑69 (MQ=255)      aaGTTGCCATTGCCAAAAGCGGAATTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGCTGCCTg                                                        >  1:3432911/1‑69 (MQ=255)      aaGTTGCCATTGCCAAAAGCGGAATTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGCTGCCTg                                                        >  1:157972/1‑69 (MQ=255)      aaGTTGCCATTGCCAAAAGCGGAATTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGCTGCCTg                                                        >  1:680068/1‑69 (MQ=255)      aaGTTGCCATTGCCAAAAGCGGAATTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGCTGCCTg                                                        >  1:3356569/1‑69 (MQ=255)      aaGTTGCCATTGCCAAAAGCGGAATTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGCTGCCTg                                                        >  1:1561453/1‑69 (MQ=255)      aaGTTGCCATTGCCAAAAGCGGAATTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGCTGCCTg                                                        >  1:1403815/1‑69 (MQ=255)      aaGTTGCCATTGCCAAAAGCGGAATTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGCTGCCTg                                                        >  1:850069/1‑69 (MQ=255)      aaGTTGCCATTGCCAAAAGCGGAATTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGCTGCCTg                                                        >  1:3207991/1‑69 (MQ=255)      aaGTTGCCATTGCCAAAAGCGGAATTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGCTGCCTg                                                        >  1:2636249/1‑69 (MQ=255)      aaGTTGCCATTGCCAAAAGCGGAATTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGCTGCCTg                                                        >  1:1105495/1‑69 (MQ=255)      aaGTTGCCATTGCCAAAAGCGGAATTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGCTGCCTg                                                        >  1:2901487/1‑69 (MQ=255)      aaGTTGCCATTGCCAAAAGCGGAATTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGCTGCCTg                                                        >  1:2905406/1‑69 (MQ=255)      aaGTTGCCATTGCCAAAAGCGGAATTGGTGAGCTATTTCCACGTCTTCTACAACCTTGTACGCTGc                                                           >  1:2395394/1‑66 (MQ=255)           gTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACTTTGTACGTTGCCTGGTCAt                                                   <  1:333765/69‑1 (MQ=255)           gTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCAt                                                   <  1:866273/69‑1 (MQ=255)           gTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCAt                                                   <  1:2774177/69‑1 (MQ=255)           gTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCAt                                                   <  1:2042011/69‑1 (MQ=255)           gTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCAt                                                   >  1:1782949/1‑69 (MQ=255)           gTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCAt                                                   <  1:729801/69‑1 (MQ=255)           gTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCAt                                                   >  1:860912/1‑69 (MQ=255)           gTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCAt                                                   <  1:794147/69‑1 (MQ=255)           gTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCAt                                                   <  1:1237050/69‑1 (MQ=255)            tCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCAt                                                   <  1:1564723/68‑1 (MQ=255)                tGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCAt                                                   <  1:1536150/64‑1 (MQ=255)                 gCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGcc                                             <  1:2912701/69‑1 (MQ=255)                 gCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGcc                                             <  1:2938372/69‑1 (MQ=255)                 gCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGcc                                             <  1:2943878/69‑1 (MQ=255)                 gCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGcc                                             <  1:2721429/69‑1 (MQ=255)                 gCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGcc                                             <  1:1204150/69‑1 (MQ=255)                 gCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGcc                                             >  1:297141/1‑69 (MQ=255)                 gCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGcc                                             <  1:1251158/69‑1 (MQ=255)                 gCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGcc                                             >  1:3283767/1‑69 (MQ=255)                 gCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGcc                                             <  1:1325813/69‑1 (MQ=255)                             aCTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCAtgctgc                                              <  1:1116923/56‑1 (MQ=255)                                   gATCTATTTCCACGTCTTCTGCAACCTTGTACGTTGCCTGGTCATGCTGCCATTTGTTGACCCGATGGc                           >  1:183597/1‑69 (MQ=255)                                   gATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGCCATTTGTTGACCCGATGGc                           >  1:3155082/1‑69 (MQ=255)                                   gATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGCCATTTGTTGACCCGATGGc                           >  1:973785/1‑69 (MQ=255)                                   gATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGCCATTTGTTGACCCGATGGc                           >  1:1075693/1‑69 (MQ=255)                                   gATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGCCATTTGTTGACCCGATGGc                           >  1:1588183/1‑69 (MQ=255)                                   gATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGCCATTTGTTGACCCGATGGc                           >  1:1883676/1‑69 (MQ=255)                                   gATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGCCATTTGTTGACCCGATGGc                           >  1:2048222/1‑69 (MQ=255)                                   gATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGCCATTTGTTGACCCGATGGc                           >  1:2377191/1‑69 (MQ=255)                                   gATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGCCATTTGTTGACCCGATGGc                           >  1:2421407/1‑69 (MQ=255)                                    aTCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGCCATTTGTTGACCCGATGGCa                          >  1:2813989/1‑69 (MQ=255)                                    aTCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGCCATTTGTTGACCCGATGGCa                          >  1:3366610/1‑69 (MQ=255)                                           tCCACGTCTTCTTCAACCTTGTACGTTGCCTGGTCATGCTGCCATTTGTTGACCCGATGGCACGGttt                    >  1:1537066/1‑68 (MQ=255)                                           tCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGCCATTTGTTGACCCGATGGCACGGttt                    >  1:262419/1‑68 (MQ=255)                                                   ttctACAACCTTGTACGTTGCCTGGTCATGCTGCCATTTGTTGACCCGATGGCACg                        <  1:327405/56‑1 (MQ=255)                                                   ttctACAACCTTGTACGTTGCCTGGTCATGCTGCCATTTGTTGACCCGATGGCACg                        <  1:128748/56‑1 (MQ=255)                                                   ttctACAACCTTGTACGTTGCCTGGTCATGCTGCCATTTGTTGACCCGATGGCACg                        <  1:2272882/56‑1 (MQ=255)                                                            cTTGTACGTTGCCTGGTCATGCTGCCATTTGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  >  1:1033949/1‑69 (MQ=255)                                                            cTTGTACGTTGCCTGGTCATGCTGCCATTTGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  >  1:2899724/1‑69 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:74636/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:563411/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:1864332/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:1657551/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:585600/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:1598326/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:2961196/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:538402/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:1255526/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:9424/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:1031768/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:54557/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:3214133/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:502259/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:1926274/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:2037651/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:3045806/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:3419246/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:3367756/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:3052957/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:2296177/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:2348319/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:2435893/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:2411229/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:3229891/69‑1 (MQ=255)                                                            cTTGTACGCTGCCTGGTCATGTTGCCATTCGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTcgc  <  1:3216734/69‑1 (MQ=255)                                                                    |                                                             TGGGGAAGTTGTCATTGCCAAAAGCGGAACTGGTGATCTATTTCCACGTCTTCTACAACCTTGTACGTTGCCTGGTCATGCTGCCATTTGTTGACCCGATGGCACGGTTTTGCAAAACGATTATTCGC  >  minE/2783719‑2783846

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