Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,246,943	T→G	100%	S122S (TCT→TCG)	znuC →	zinc transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,246,943	0	T	G	90.1%	65.1 / ‑0.9	30	S122S (TCT→TCG)	znuC	zinc transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base G (16/11);  total (16/14)
Fisher's exact test for biased strand distribution p-value = 8.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGTGTCCAGGCCGGGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGT  >  minE/1246896‑1247005                                                |                                                               cGTGTCCAGGCCGGGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTa                                           <  1:1375634/69‑1 (MQ=255)  gtgtCCAGGCCGGGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTa                                           <  1:2285406/68‑1 (MQ=255)                  tcTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:2940353/68‑1 (MQ=255)                  tcTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:3049573/68‑1 (MQ=255)                  tcTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:400016/68‑1 (MQ=255)                  tcTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:647382/68‑1 (MQ=255)                  tcTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:2066186/68‑1 (MQ=255)                  tcTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:2062162/68‑1 (MQ=255)                  tcTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:199684/68‑1 (MQ=255)                  tcTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:1946544/68‑1 (MQ=255)                  tcTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:1684465/68‑1 (MQ=255)                  tcTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:738896/68‑1 (MQ=255)                   cTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:676046/68‑1 (MQ=255)                              aCCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCgac               >  1:3243932/1‑66 (MQ=255)                              aCCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCgac               >  1:3055903/1‑66 (MQ=255)                              aCCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCgac               >  1:894/1‑66 (MQ=255)                              aCCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCgac               >  1:2502318/1‑66 (MQ=255)                              aCCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCgac               >  1:1277028/1‑66 (MQ=255)                              aCCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCgac               >  1:2300454/1‑66 (MQ=255)                              aCCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCgac               >  1:2181258/1‑66 (MQ=255)                              aCCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCgac               >  1:1571005/1‑66 (MQ=255)                                    gCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCa           <  1:2303009/66‑1 (MQ=255)                                          gCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAgcg                                  >  1:2008672/1‑37 (MQ=37)                                          gCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCGACCGCAATTATTAGt  >  1:2439770/1‑69 (MQ=255)                                          gCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCGACCGCAATTATTAGt  >  1:3075790/1‑69 (MQ=255)                                          gCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCGACCGCAATTATTAGt  >  1:316605/1‑69 (MQ=255)                                          gCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCGACCGCAATTATTAGt  >  1:676705/1‑69 (MQ=255)                                          gCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCGACCGCAATTATTAGt  >  1:1428572/1‑69 (MQ=255)                                          gCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCGACCGCAATTATTAGt  >  1:995268/1‑69 (MQ=255)                                          gCTCTCGGGAGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCGACCGCAATTATTAGt  >  1:700543/1‑69 (MQ=255)                                                |                                                               CGTGTCCAGGCCGGGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGT  >  minE/1246896‑1247005

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