Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	799,135	G→A	57.6%	R152R (CGC→CGT)	treA ←	periplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	799,135	0	G	A	57.6%	19.0 / 54.9	73	R152R (CGC→CGT)	treA	periplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base G (29/2);  new base A (27/15);  total (56/17)
Fisher's exact test for biased strand distribution p-value = 4.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGG  >  minE/799080‑799193                                                        |                                                           tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATGCCTCGCGAAAACGTCCGCCGGGc                                                >  1:81341/1‑68 (MQ=255) tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGc                                                >  1:223241/1‑68 (MQ=255) tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGc                                                >  1:3381263/1‑68 (MQ=255) tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGc                                                >  1:2501869/1‑68 (MQ=255) tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGc                                                >  1:2396153/1‑68 (MQ=255) tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGc                                                >  1:1477765/1‑68 (MQ=255) tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGc                                                >  1:1369795/1‑68 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGt                                    <  1:2721777/69‑1 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGt                                    <  1:204758/69‑1 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACg                                                          >  1:86074/1‑47 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACg                                                          >  1:2985967/1‑47 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACg                                                          >  1:2872834/1‑47 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACg                                                          >  1:3126697/1‑47 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACg                                                          >  1:2830943/1‑47 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACg                                                          >  1:2807457/1‑47 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACg                                                          >  1:540540/1‑47 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACg                                                          >  1:2519152/1‑47 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACg                                                          >  1:2485043/1‑47 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACg                                                          >  1:2476858/1‑47 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACg                                                          >  1:2404915/1‑47 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACg                                                          >  1:919279/1‑47 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACg                                                          >  1:606123/1‑47 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACg                                                          >  1:891713/1‑47 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACg                                                          >  1:1254772/1‑47 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACg                                                          >  1:1528611/1‑47 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACg                                                          >  1:847041/1‑47 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACg                                                          >  1:3004736/1‑47 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACg                                                          >  1:1821250/1‑47 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACg                                                          >  1:638334/1‑47 (MQ=255)                     gTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGc                           >  1:1163460/1‑69 (MQ=255)                           tAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAg                                      >  1:1653604/1‑52 (MQ=255)                           tAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGt                        >  1:2887917/1‑66 (MQ=255)                           tAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGt                        >  1:2637191/1‑66 (MQ=255)                               tGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACaga                  >  1:2449044/1‑68 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                <  1:1814340/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                <  1:3005558/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                <  1:1313035/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                <  1:2791089/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                <  1:2615958/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                <  1:1745785/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                <  1:80524/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                <  1:1972956/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                <  1:2368546/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                <  1:1981825/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                <  1:615884/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                <  1:2263411/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                <  1:222910/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAcaga                <  1:3110854/69‑1 (MQ=255)                                      gTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCa          >  1:358886/1‑69 (MQ=255)                                      gTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCa          >  1:119942/1‑69 (MQ=255)                                       tAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCa          >  1:3063869/1‑68 (MQ=255)                                       tAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAt         >  1:347214/1‑69 (MQ=255)                                       tAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAt         >  1:957497/1‑69 (MQ=255)                                       tAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAt         >  1:1070586/1‑69 (MQ=255)                                       tAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAt         >  1:1265937/1‑69 (MQ=255)                                       tAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAt         >  1:185862/1‑69 (MQ=255)                                       tAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAt         >  1:1991619/1‑69 (MQ=255)                                       tAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAAAAGAGAATCCCAt         >  1:1185827/1‑69 (MQ=255)                                        aaTATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAtt        >  1:2965852/1‑69 (MQ=255)                                        aaTATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                <  1:2281198/61‑1 (MQ=255)                                         atatACCTCGCGAAAGCGTCCGCCCGGCACGACATaa                                       >  1:2859876/1‑37 (MQ=255)                                         atatACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAttt       >  1:2291147/1‑69 (MQ=255)                                         atatACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAttt       >  1:367176/1‑69 (MQ=255)                                         atatACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAttt       >  1:773300/1‑69 (MQ=255)                                         atatACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAttt       >  1:1013861/1‑69 (MQ=255)                                              ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg  >  1:416574/1‑69 (MQ=255)                                              ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg  >  1:2527826/1‑69 (MQ=255)                                              ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg  >  1:2992244/1‑69 (MQ=255)                                              ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg  >  1:1863319/1‑69 (MQ=255)                                              ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg  >  1:1556474/1‑69 (MQ=255)                                              ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg  >  1:1318125/1‑69 (MQ=255)                                              ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg  >  1:2921012/1‑69 (MQ=255)                                              ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg  >  1:342748/1‑69 (MQ=255)                                                        |                                                           TTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGG  >  minE/799080‑799193

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