Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	799,144	C→G	39.3%	P149P (CCG→CCC)	treA ←	periplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	799,144	0	C	G	39.3%	32.4 / 51.7	56	P149P (CCG→CCC)	treA	periplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base C (32/2);  new base G (7/15);  total (39/17)
Fisher's exact test for biased strand distribution p-value = 1.00e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAAC  >  minE/799080‑799208                                                                 |                                                                 tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATGCCTCGCGAAAACGTCCGCCGGGc                                                               >  1:81341/1‑68 (MQ=255) tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGc                                                               >  1:1369795/1‑68 (MQ=255) tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGc                                                               >  1:223241/1‑68 (MQ=255) tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGc                                                               >  1:3381263/1‑68 (MQ=255) tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGc                                                               >  1:1477765/1‑68 (MQ=255) tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGc                                                               >  1:2396153/1‑68 (MQ=255) tttCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGc                                                               >  1:2501869/1‑68 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGt                                                   <  1:2721777/69‑1 (MQ=255)            ccTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGt                                                   <  1:204758/69‑1 (MQ=255)                     gTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGc                                          >  1:1163460/1‑69 (MQ=255)                           tAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAg                                                     >  1:1653604/1‑52 (MQ=255)                           tAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGt                                       >  1:2637191/1‑66 (MQ=255)                           tAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGt                                       >  1:2887917/1‑66 (MQ=255)                               tGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACaga                                 >  1:2449044/1‑68 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:2368546/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:2615958/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:2791089/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:2263411/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:222910/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:1981825/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:1972956/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:3005558/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:1814340/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:1745785/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:1313035/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:615884/69‑1 (MQ=255)                                gTCCCAGTAATATACCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAgaga                               <  1:80524/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:1991619/1‑69 (MQ=255)                                       tAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAt                        >  1:1070586/1‑69 (MQ=255)                                       tAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAt                        >  1:185862/1‑69 (MQ=255)                                       tAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAt                        >  1:1265937/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:1013861/1‑69 (MQ=255)                                         atatACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAttt                      >  1:2291147/1‑69 (MQ=255)                                         atatACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAttt                      >  1:773300/1‑69 (MQ=255)                                         atatACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCAttt                      >  1:367176/1‑69 (MQ=255)                                              ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg                 >  1:2527826/1‑69 (MQ=255)                                              ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg                 >  1:416574/1‑69 (MQ=255)                                              ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg                 >  1:1318125/1‑69 (MQ=255)                                              ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg                 >  1:342748/1‑69 (MQ=255)                                              ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg                 >  1:1556474/1‑69 (MQ=255)                                              ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg                 >  1:1863319/1‑69 (MQ=255)                                              ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg                 >  1:2992244/1‑69 (MQ=255)                                              ccTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCgg                 >  1:2921012/1‑69 (MQ=255)                                                             cGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAAc  >  1:1582394/1‑69 (MQ=255)                                                             cGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAAc  >  1:1112481/1‑69 (MQ=255)                                                             cGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAAc  >  1:2531539/1‑69 (MQ=255)                                                                 |                                                                 TTTCGGCAAGTCCTAACATGGTGAAGTAACTGTCCCAGTAATATACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCCATTTTTCGGTGTTTTCGGTAGAAC  >  minE/799080‑799208

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