Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	906,094	A→C	37.5%	I41L (ATC→CTC)	ycjR →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	906,094	0	A	C	37.5%	22.2 / 40.4	48	I41L (ATC→CTC)	ycjR	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/11);  new base C (10/8);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGGTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGC  >  minE/906055‑906156                                        |                                                               gATGGTTTTGAGATCGACGGCAAACTGCTGGTCAACAATCTCGAAGAAGTCAAAGCAGCAATCAAAGa                                    <  1:3426852/68‑1 (MQ=255) gATGGTTTTGAGATCGACGGCAAACTGCTGGTCAACAATCTCGAAGAAGTCAAAGCAGCAATCAAAGa                                    <  1:350953/68‑1 (MQ=255) gATGGTTTTGAGATCGACGGCAAACTGCTGGTCAACAATCTCGAAGAAGTCAAAGCAGCAATCAAAGa                                    <  1:441063/68‑1 (MQ=255) gATGGTTTTGAGATCGACGGCAAACTGCTGGTCAACAATCTCGAAGAAGTCAAAGCAGCAATCAAAGa                                    <  1:2381443/68‑1 (MQ=255) gATGGTTTTGAGATCGACGGCAAACTGCTGGTCAACAATCTCGAAGAAGTCAAAGCAGCAATCAAAGa                                    <  1:255387/68‑1 (MQ=255) gATGGTTTTGAGATCGACGGCAAACTGCTGGTCAACAATCTCGAAGAAGTCAAAGCAGCAATCAAAGa                                    <  1:149922/68‑1 (MQ=255) gATGGTTTTGAGATCGACGGCAAACTGCTGGTCAACAATCTCGAAGAAGTCAAAGCAGCAATCAAAGa                                    <  1:2984994/68‑1 (MQ=255) gATGGTTTTGAGATCGACGGCAAACTGCTGGTCAACAATCTCGAAGAAGTCAAAGAAGCAATCAAAGa                                    <  1:1985041/68‑1 (MQ=255)           agaTTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTa                          <  1:1920486/68‑1 (MQ=255)           agaTTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTa                          <  1:2756098/68‑1 (MQ=255)           agaTTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTa                          <  1:3051267/68‑1 (MQ=255)           agaTTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTa                          <  1:1435484/68‑1 (MQ=255)           agaTTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTa                          <  1:1041213/68‑1 (MQ=255)           agaTTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTa                          <  1:3324551/68‑1 (MQ=255)                gACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTg                    >  1:3113868/1‑69 (MQ=255)                gACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTg                    >  1:3208861/1‑69 (MQ=255)                gACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTg                    >  1:32312/1‑69 (MQ=255)                gACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTg                    >  1:2233325/1‑69 (MQ=255)                gACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTg                    >  1:3335790/1‑69 (MQ=255)                gACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTg                    >  1:1685191/1‑69 (MQ=255)                gACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTg                    >  1:758333/1‑69 (MQ=255)                         ctgctgGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGt                     <  1:1862229/59‑1 (MQ=255)                         ctgctgGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGt                     <  1:1841911/59‑1 (MQ=255)                         ctgctgGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGt                     >  1:1825397/1‑59 (MQ=255)                         ctgctgGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGt                     >  1:1164016/1‑59 (MQ=255)                         ctgctgGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGt                     >  1:614832/1‑59 (MQ=255)                         ctgctgGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGt                     <  1:1182187/59‑1 (MQ=255)                         ctgctgGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGt                     >  1:1305755/1‑59 (MQ=255)                         ctgctgGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGt                     <  1:3200988/59‑1 (MQ=255)                         ctgctgGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGt                     <  1:3123958/59‑1 (MQ=255)                         ctgctgGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGcc           >  1:269636/1‑69 (MQ=255)                                 tAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTg    >  1:2077967/1‑68 (MQ=255)                                 tAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTg    >  1:154501/1‑68 (MQ=255)                                 caacaatCTCGAAGAAGTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTgtggtg    >  1:729141/2‑68 (MQ=255)                                 caacaatCTCGAAGAAGTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTgtggtg    >  1:633348/2‑68 (MQ=255)                                 caacaatCTCGAAGAAGTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTgtggtg    >  1:3410224/2‑68 (MQ=255)                                 caacaatCTCGAAGAAGTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTgtggtg    >  1:933956/2‑68 (MQ=255)                                 caacaatCTCGAAGAAGTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTgtggtg    >  1:2418660/2‑68 (MQ=255)                                 caacaatCTCGAAGAAGTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTgtggtg    >  1:1602340/2‑68 (MQ=255)                                 caacaatCTCGAAGAAGTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTgtggtg    >  1:1562618/2‑68 (MQ=255)                                 caacaatCTCGAAGAAGTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTgtggtg    >  1:1474390/2‑68 (MQ=255)                                 caacaatCTCGAAGAAGTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTgtggtg    >  1:13578/2‑68 (MQ=255)                                 caacaatCTCGAAGAAGTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTgtggtg    >  1:1332247/2‑68 (MQ=255)                                  aacaacATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGc  >  1:2264096/1‑69 (MQ=255)                                  aacaacATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGc  >  1:1728899/1‑69 (MQ=255)                                  aacaacATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGc  >  1:412109/1‑69 (MQ=255)                                  aacaacATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGc  >  1:438526/1‑69 (MQ=255)                                  aacaacATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGc  >  1:1725614/1‑69 (MQ=255)                                        |                                                               GATGGTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGC  >  minE/906055‑906156

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