Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	454,293	C→A	64.8%	R264R (CGG→CGT)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	454,293	0	C	A	64.8%	62.7 / 75.5	91	R264R (CGG→CGT)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (21/11);  new base A (35/24);  total (56/35)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTC  >  minE/454230‑454356                                                                |                                                                ttgtGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTg                                                            >  1:58892/1‑69 (MQ=255) ttgtGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTg                                                            >  1:1747834/1‑69 (MQ=255)             ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGtttt                                                 >  1:2305692/1‑68 (MQ=255)             ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGtttt                                                 >  1:1790977/1‑68 (MQ=255)             ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc                                                >  1:11234/1‑69 (MQ=255)             ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc                                                >  1:2110715/1‑69 (MQ=255)             ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc                                                >  1:2266002/1‑69 (MQ=255)             ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc                                                >  1:2382367/1‑69 (MQ=255)             ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc                                                >  1:2927921/1‑69 (MQ=255)             ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc                                                >  1:1940833/1‑69 (MQ=255)             ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc                                                >  1:1872910/1‑69 (MQ=255)             ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc                                                >  1:311570/1‑69 (MQ=255)             ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc                                                >  1:395850/1‑69 (MQ=255)             ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc                                                >  1:1662344/1‑69 (MQ=255)             ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc                                                >  1:969270/1‑69 (MQ=255)             ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc                                                >  1:772244/1‑69 (MQ=255)             ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc                                                >  1:1169824/1‑69 (MQ=255)             ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTc                                                >  1:911619/1‑69 (MQ=255)             ccAGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTAAGCTGTTTTTGGtttt                                                 >  1:633953/1‑68 (MQ=255)               aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGtt                                                   >  1:321731/1‑64 (MQ=255)               aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg                                              >  1:971492/1‑69 (MQ=255)               aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg                                              >  1:3383827/1‑69 (MQ=255)               aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg                                              >  1:2116057/1‑69 (MQ=255)               aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg                                              >  1:2211205/1‑69 (MQ=255)               aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg                                              >  1:2632299/1‑69 (MQ=255)               aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg                                              >  1:2652927/1‑69 (MQ=255)               aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg                                              >  1:89835/1‑69 (MQ=255)               aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg                                              >  1:2989567/1‑69 (MQ=255)               aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg                                              >  1:3020108/1‑69 (MQ=255)               aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg                                              >  1:3140766/1‑69 (MQ=255)               aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg                                              >  1:354016/1‑69 (MQ=255)               aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg                                              >  1:1947926/1‑69 (MQ=255)               aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg                                              >  1:2022828/1‑69 (MQ=255)               aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg                                              >  1:1724399/1‑69 (MQ=255)               aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg                                              >  1:1511849/1‑69 (MQ=255)               aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg                                              >  1:2050359/1‑69 (MQ=255)               aGAACAACATGCTGAAGAAGGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTg                                              >  1:166451/1‑69 (MQ=255)                  acaacaTACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCtt                                           <  1:277675/69‑1 (MQ=255)                  acaacaTACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCtt                                           <  1:1179325/69‑1 (MQ=255)                          cTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGAGATaa                                   >  1:574652/1‑69 (MQ=255)                            gaagaaGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGGTTCTGGTTTTCTGCTTTGCGATAAAt                                 >  1:2918933/1‑69 (MQ=255)                                aagGTGAGGGTAATAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:903632/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTAGGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:3430357/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:1658804/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:599329/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:1559521/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:795652/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:3288789/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:3252072/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:176228/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:1771575/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:1771753/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:1798430/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:2088730/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:934281/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:1078636/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:2954967/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:2815624/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:2778205/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:2753931/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:2736166/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:2128877/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:2697907/69‑1 (MQ=255)                                aagGTGAGGGTAACAATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:2526161/69‑1 (MQ=255)                                  ggTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAAttttt                             >  1:2545399/1‑67 (MQ=255)                                     gAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAg                        >  1:2870460/1‑69 (MQ=255)                                     gAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAg                        >  1:848124/1‑69 (MQ=255)                                     gAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAg                        >  1:665121/1‑69 (MQ=255)                                     gAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAg                        >  1:1966800/1‑69 (MQ=255)                                     gAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAg                        >  1:1725012/1‑69 (MQ=255)                                     gAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAg                        >  1:1166749/1‑69 (MQ=255)                                      aGCGTCACAATCAGCCCCAGCTCGTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCa                         >  1:920826/1‑67 (MQ=255)                                      aGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAg                        >  1:213550/1‑68 (MQ=255)                                      aGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAg                        >  1:2674125/1‑68 (MQ=255)                                      aGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAg                        >  1:2142653/1‑68 (MQ=255)                                      aGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAg                        >  1:240434/1‑68 (MQ=255)                                      aGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAg                        >  1:2345067/1‑68 (MQ=255)                                      aGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAg                        >  1:128802/1‑68 (MQ=255)                                      aGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAg                        >  1:2210798/1‑68 (MQ=255)                                         gtaacaATCAGCCCCAGCTCCTTACGCTGTTTTTGGTTTTCTGCTTTGCGATAAAttttt                             <  1:2029195/57‑1 (MQ=255)                                         gTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTaca                    <  1:2916414/69‑1 (MQ=255)                                              aaTCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAg               <  1:508749/69‑1 (MQ=255)                                              aaTCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAg               <  1:2671097/69‑1 (MQ=255)                                              aaTCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAg               <  1:2708345/69‑1 (MQ=255)                                              aaTCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAg               <  1:2712241/69‑1 (MQ=255)                                                    ccccAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCaat         <  1:208672/69‑1 (MQ=255)                                                    ccccAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCaat         <  1:2049142/69‑1 (MQ=255)                                                    caccAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCaat         <  1:2121961/67‑1 (MQ=255)                                                      ccAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCaataa       <  1:1749483/69‑1 (MQ=255)                                                           tCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTc  >  1:749914/1‑69 (MQ=255)                                                           tCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTc  >  1:709734/1‑69 (MQ=255)                                                                |                                                                TTGTGCGAAGGCCCAGAACAACATACTGAAGAAGGTGAGCGTCACAATCAGCCCCAGCTCCTTCCGCTGTTTCTGGTTTTCTGCTTTGCGATAAATTTTTGCCAGTACACCCAGGCCAATAATAGTC  >  minE/454230‑454356

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