Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,870,835	G→A	55.2%	S254S (TCG→TCA)	yqcE →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,870,835	0	G	A	55.2%	8.8 / 53.9	67	S254S (TCG→TCA)	yqcE	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (15/15);  new base A (2/35);  total (17/50)
Fisher's exact test for biased strand distribution p-value = 3.70e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCTACGCCATTCTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTT  >  minE/1870778‑1870901                                                          |                                                                   cATCTACGCCATTCTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTggcgtc                                                         <  1:1772354/69‑3 (MQ=255) cATCTACGCCATTCTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTggcggc                                                         <  1:950208/69‑1 (MQ=255) cATCTACGCCATTCTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTggcggc                                                         <  1:3370305/69‑1 (MQ=255) cATCTACGCCATTCTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTggcggc                                                         <  1:1221600/69‑1 (MQ=255) cATCTACGCCATTCTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTggcggc                                                         <  1:1320652/69‑1 (MQ=255) cATCTACGCCATTCTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTggcggc                                                         <  1:379829/69‑1 (MQ=255) cATCTACGCCATTCTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTggcggc                                                         <  1:470148/69‑1 (MQ=255) cATCTACGCCATTCTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTggcggc                                                         <  1:792983/69‑1 (MQ=255) cATCTACGCCATTCTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTggcggc                                                         <  1:1543892/69‑1 (MQ=255) cATCTACGCCATTCTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTggcggc                                                         <  1:1664581/69‑1 (MQ=255) cATCTACGCCATTCTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTggcg                                                           >  1:2516861/1‑67 (MQ=255) cATCTACGCCATTCTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTggcg                                                           >  1:623171/1‑67 (MQ=255) cATCTACGCCATTCTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTggcg                                                           >  1:1584213/1‑67 (MQ=255)     tACGCCATTCTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTGGCGGCGAg                                                      >  1:273467/1‑68 (MQ=255)          cATTCTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTGGCGGCGAGCTACAt                                                <  1:1432861/69‑1 (MQ=255)              cTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCg                                                                    >  1:688472/1‑45 (MQ=255)              cTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTGGCGGCGAGCTACATgggg                                            >  1:1571295/1‑69 (MQ=255)              cTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTGGCGGCGAGCTACATgggg                                            >  1:918071/1‑69 (MQ=255)              cTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTGGCGGCGAGCTACATgggg                                            >  1:1451218/1‑69 (MQ=255)                    tACTCCACCAACTATCTGACCTAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:745439/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGTGATtgtg                                      <  1:1032236/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:695013/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:623311/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:695443/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:81234/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:840155/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:3369128/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:326726/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:3154626/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:3126951/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:290895/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:2790042/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:2744738/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:963823/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:2472430/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:2320530/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:1125154/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:2295750/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:137997/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:1527199/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:2233496/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:2000796/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:2620922/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:2036699/69‑1 (MQ=255)                    tACTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:2127520/69‑1 (MQ=255)                     aCTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTGGCGGCGAGCTACATGGGGATtgtg                                      >  1:538294/1‑68 (MQ=255)                     aCTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTGGCGGCGAGCTACATGGGGATtgtg                                      >  1:787979/1‑68 (MQ=255)                     aCTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTGGCGGCGAGCTACATGGGGATtgtg                                      >  1:1691663/1‑68 (MQ=255)                     aCTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTGGCGGCGAGCTACATGGGGATtgtg                                      >  1:2138770/1‑68 (MQ=255)                     aCTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTGGCGGCGAGCTACATGGGGATtgtg                                      >  1:2937072/1‑68 (MQ=255)                     aCTCCACCAACTATCTGACCGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:2035122/68‑1 (MQ=255)                                        cGAAATGTACTGCATGTCACTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATCTTCcgc                    <  1:193491/67‑1 (MQ=255)                                        cGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATtgtg                                      <  1:861203/49‑1 (MQ=255)                                        cGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATCTTCcgc                    <  1:1099031/67‑1 (MQ=255)                                        cGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATCTTCcgc                    <  1:3433590/67‑1 (MQ=255)                                        cGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATCTTCcgc                    <  1:2795324/67‑1 (MQ=255)                                        cGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATCTTCcgc                    <  1:2279593/67‑1 (MQ=255)                                        cGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATCTTCcgc                    <  1:1861809/67‑1 (MQ=255)                                        cGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATCTTCcgc                    <  1:1880078/67‑1 (MQ=255)                                         gAAATGTATGGCATGTCGCTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATATTCcgcgcg                 <  1:1633831/69‑1 (MQ=255)                                         gAAATGTATGGCATGTCGCTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATATTCcgcgcg                 <  1:3127912/69‑1 (MQ=255)                                         gAAATGTATGGCATGTCGCTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATATTCcgcgcg                 <  1:120563/69‑1 (MQ=255)                                            aTGTATGGCATGTCGCTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCtg              <  1:132556/69‑1 (MQ=255)                                                acggCATGTCACTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATCTTCCGCGCGTTATGc           >  1:3067380/3‑68 (MQ=255)                                                acggCATGTCACTGGTGGCGGCGAGCTACATGAGGATTGTGATCAACAAAATCTTCCGCGCGTTATGc           >  1:71114/3‑68 (MQ=255)                                                   gCATGTCGCTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGcc       >  1:843240/1‑69 (MQ=255)                                                   gCATGTCGCTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGcc       >  1:2215134/1‑69 (MQ=255)                                                        tCGCTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACtt  >  1:321762/1‑69 (MQ=255)                                                          |                                                                   CATCTACGCCATTCTGAGTTACTCCACCAACTATCTGACCGAAATGTATGGCATGTCGCTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTT  >  minE/1870778‑1870901

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