Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	798,949	C→T	34.1%	M214I (ATG→ATA)	treA ←	periplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	798,949	0	C	T	34.1%	47.9 / 40.4	44	M214I (ATG→ATA)	treA	periplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base C (27/2);  new base T (9/6);  total (36/8)
Fisher's exact test for biased strand distribution p-value = 1.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

CAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGAATA  >  minE/798898‑799012                                                    |                                                                cAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCg                                              >  1:1822738/1‑71 (MQ=255) cAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCg                                              >  1:2270361/1‑71 (MQ=255) cAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCg                                              >  1:2326717/1‑71 (MQ=255) cAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCg                                              >  1:2396192/1‑71 (MQ=255) cAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCg                                              >  1:2059266/1‑71 (MQ=255) cAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCg                                              >  1:1952472/1‑71 (MQ=255) cAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCg                                              >  1:2763324/1‑71 (MQ=255) cAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCg                                              >  1:2892460/1‑71 (MQ=255) cAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCg                                              >  1:57955/1‑71 (MQ=255) cAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCg                                              >  1:1665767/1‑71 (MQ=255) cAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCg                                              >  1:588487/1‑71 (MQ=255) cAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCg                                              >  1:71342/1‑71 (MQ=255) cAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCg                                              >  1:824292/1‑71 (MQ=255) cAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCg                                              >  1:1064111/1‑71 (MQ=255) cAGGTAATGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCg                                              >  1:1235756/1‑71 (MQ=255)   ggTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCgg                                             >  1:1798398/1‑70 (MQ=255)   ggTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGt                                            >  1:1131596/1‑71 (MQ=255)   ggTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGt                                            >  1:2814384/1‑71 (MQ=255)   ggTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGt                                            >  1:2742122/1‑71 (MQ=255)   ggTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGt                                            >  1:2596344/1‑71 (MQ=255)   ggTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGt                                            >  1:1369977/1‑71 (MQ=255)   ggTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGt                                            >  1:1687701/1‑71 (MQ=255)   ggTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGt                                            >  1:2127689/1‑71 (MQ=255)   ggTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGt                                            >  1:2127009/1‑71 (MQ=255)                  cgGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTaaa                             <  1:460721/71‑1 (MQ=255)                   ggCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTaaa                             <  1:53436/70‑1 (MQ=255)                      aTCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGt                         <  1:1885623/71‑1 (MQ=255)                      aTCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGt                         <  1:849988/71‑1 (MQ=255)                      aTCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGt                         <  1:1561607/71‑1 (MQ=255)                      aTCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGt                         <  1:336733/71‑1 (MQ=255)                      aTCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGt                         <  1:2548928/71‑1 (MQ=255)                       tCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGt                         <  1:1730259/70‑1 (MQ=255)                                             aCTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAATTAGTAACTGCGGTTGCCGTTGGGAata  >  1:1452130/1‑71 (MQ=255)                                             aCTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGaa    >  1:2699786/1‑69 (MQ=255)                                             aCTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGAata  >  1:2740852/1‑71 (MQ=255)                                             aCTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGAata  >  1:1887991/1‑71 (MQ=255)                                             aCTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGAata  >  1:2220873/1‑71 (MQ=255)                                             aCTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGAata  >  1:2873297/1‑71 (MQ=255)                                             aCTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGAata  >  1:2495980/1‑71 (MQ=255)                                             aCTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGAata  >  1:2421410/1‑71 (MQ=255)                                             aCTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGAata  >  1:2413436/1‑71 (MQ=255)                                             aCTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGAata  >  1:738162/1‑71 (MQ=255)                                             aCTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGAata  >  1:214428/1‑71 (MQ=255)                                             aCTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGAata  >  1:2362474/1‑71 (MQ=255)                                                    |                                                                CAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGAATA  >  minE/798898‑799012

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