Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	906,111	G→A	48.6%	A46A (GCG→GCA)	ycjR →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	906,111	0	G	A	48.6%	1.4 / 62.0	72	A46A (GCG→GCA)	ycjR	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/15);  new base A (27/8);  total (49/23)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGGTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGAC  >  minE/906055‑906177                                                         |                                                                   gATGGTTTTGAGATCGACGGCAAACTGCTGGTCAACAATCTCGAAGAAGTCAAAGCAGCAATCAAAGa                                                         <  1:255387/68‑1 (MQ=255) gATGGTTTTGAGATCGACGGCAAACTGCTGGTCAACAATCTCGAAGAAGTCAAAGCAGCAATCAAAGa                                                         <  1:2381443/68‑1 (MQ=255) gATGGTTTTGAGATCGACGGCAAACTGCTGGTCAACAATCTCGAAGAAGTCAAAGCAGCAATCAAAGa                                                         <  1:2984994/68‑1 (MQ=255) gATGGTTTTGAGATCGACGGCAAACTGCTGGTCAACAATCTCGAAGAAGTCAAAGCAGCAATCAAAGa                                                         <  1:441063/68‑1 (MQ=255) gATGGTTTTGAGATCGACGGCAAACTGCTGGTCAACAATCTCGAAGAAGTCAAAGCAGCAATCAAAGa                                                         <  1:350953/68‑1 (MQ=255) gATGGTTTTGAGATCGACGGCAAACTGCTGGTCAACAATCTCGAAGAAGTCAAAGCAGCAATCAAAGa                                                         <  1:149922/68‑1 (MQ=255) gATGGTTTTGAGATCGACGGCAAACTGCTGGTCAACAATCTCGAAGAAGTCAAAGCAGCAATCAAAGa                                                         <  1:3426852/68‑1 (MQ=255) gATGGTTTTGAGATCGACGGCAAACTGCTGGTCAACAATCTCGAAGAAGTCAAAGAAGCAATCAAAGa                                                         <  1:1985041/68‑1 (MQ=255)           agaTTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTa                                               <  1:1041213/68‑1 (MQ=255)           agaTTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTa                                               <  1:2756098/68‑1 (MQ=255)           agaTTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTa                                               <  1:1920486/68‑1 (MQ=255)           agaTTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTa                                               <  1:3051267/68‑1 (MQ=255)           agaTTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTa                                               <  1:3324551/68‑1 (MQ=255)           agaTTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTa                                               <  1:1435484/68‑1 (MQ=255)                gACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTg                                         >  1:758333/1‑69 (MQ=255)                gACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTg                                         >  1:3335790/1‑69 (MQ=255)                gACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTg                                         >  1:32312/1‑69 (MQ=255)                gACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTg                                         >  1:1685191/1‑69 (MQ=255)                gACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTg                                         >  1:3208861/1‑69 (MQ=255)                gACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTg                                         >  1:2233325/1‑69 (MQ=255)                gACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTg                                         >  1:3113868/1‑69 (MQ=255)                         ctgctgGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGt                                          <  1:1841911/59‑1 (MQ=255)                         ctgctgGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGt                                          >  1:1164016/1‑59 (MQ=255)                         ctgctgGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGt                                          <  1:1182187/59‑1 (MQ=255)                         ctgctgGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGt                                          >  1:1305755/1‑59 (MQ=255)                         ctgctgGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGt                                          <  1:1862229/59‑1 (MQ=255)                         ctgctgGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGt                                          >  1:1825397/1‑59 (MQ=255)                         ctgctgGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGt                                          <  1:3123958/59‑1 (MQ=255)                         ctgctgGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGt                                          <  1:3200988/59‑1 (MQ=255)                         ctgctgGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGt                                          >  1:614832/1‑59 (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:3410224/2‑68 (MQ=255)                                 caacaatCTCGAAGAAGTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTgtggtg                         >  1:633348/2‑68 (MQ=255)                                 caacaatCTCGAAGAAGTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTgtggtg                         >  1:729141/2‑68 (MQ=255)                                 caacaatCTCGAAGAAGTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTgtggtg                         >  1:933956/2‑68 (MQ=255)                                 caacaatCTCGAAGAAGTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTgtggtg                         >  1:13578/2‑68 (MQ=255)                                 caacaatCTCGAAGAAGTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTgtggtg                         >  1:1474390/2‑68 (MQ=255)                                 caacaatCTCGAAGAAGTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTgtggtg                         >  1:1562618/2‑68 (MQ=255)                                 caacaatCTCGAAGAAGTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTgtggtg                         >  1:1602340/2‑68 (MQ=255)                                 caacaatCTCGAAGAAGTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTgtggtg                         >  1:1332247/2‑68 (MQ=255)                                 caacaatCTCGAAGAAGTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTgtggtg                         >  1:2418660/2‑68 (MQ=255)                                  aacaacATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGc                       >  1:2264096/1‑69 (MQ=255)                                  aacaacATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGc                       >  1:438526/1‑69 (MQ=255)                                  aacaacATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGc                       >  1:1728899/1‑69 (MQ=255)                                  aacaacATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGc                       >  1:1725614/1‑69 (MQ=255)                                  aacaacATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGc                       >  1:412109/1‑69 (MQ=255)                                         tCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACg                <  1:3134901/69‑1 (MQ=255)                                           gaagaaGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACgga              >  1:1369226/1‑69 (MQ=255)                                            aagaagTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGTGGTGGCTATGACggt              >  1:63942/1‑67 (MQ=255)                                            aagaagTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGTGGTGGCTATGACggt              >  1:63271/1‑67 (MQ=255)                                            aagaagTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGTGGTGGCTATGACggt              >  1:1141493/1‑67 (MQ=255)                                            aagaagTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGTGGTGGCTATGACggt              >  1:567594/1‑67 (MQ=255)                                            aagaagTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGTGGTGGCTATGACggt              >  1:1012056/1‑67 (MQ=255)                                            aagaagTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGTGGTGGCTATGACggt              >  1:431206/1‑67 (MQ=255)                                            aagaagTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGTGGTGGCTATGACggt              >  1:852669/1‑67 (MQ=255)                                            aagaagTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGTGGTGGCTATGACggt              >  1:2461161/1‑67 (MQ=255)                                            aagaagTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGTGGTGGCTATGACggt              >  1:1498466/1‑67 (MQ=255)                                            aagaagTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGTGGTGGCTATGACggt              >  1:1650093/1‑67 (MQ=255)                                            aagaagTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGTGGTGGCTATGACggt              >  1:183786/1‑67 (MQ=255)                                            aagaagTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGTGGTGGCTATGACggt              >  1:3078572/1‑67 (MQ=255)                                            aagaagTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGTGGTGGCTATGACggt              >  1:3035431/1‑67 (MQ=255)                                            aagaagTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGTGGTGGCTATGACggt              >  1:2270196/1‑67 (MQ=255)                                            aagaagTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGTGGTGGCTATGACggt              >  1:2737768/1‑67 (MQ=255)                                            aagaagTCAAAGCAGCAATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGTGGTGGCTATGACgg               >  1:553497/1‑67 (MQ=255)                                            aagaagTCAAAGCAG‑AATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGTGGTGGCTATGACggttg            >  1:1754184/1‑69 (MQ=255)                                                    aaaGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGc                       <  1:3430496/51‑1 (MQ=255)                                                    aaaGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGc                       >  1:1884145/1‑51 (MQ=255)                                                    aaaGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGc                       <  1:2781110/51‑1 (MQ=255)                                                    aaaGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGc                       <  1:1185117/51‑1 (MQ=255)                                                    aaaGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGc                       >  1:2512938/1‑51 (MQ=255)                                                        cggcgATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGAc  >  1:1232850/1‑68 (MQ=255)                                                        cggcgATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGAc  >  1:2439304/1‑68 (MQ=255)                                                         |                                                                   GATGGTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGAC  >  minE/906055‑906177

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