Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,328,372	G→A	50.0%	F131F (TTC→TTT)	wzxC ←	colanic acid exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,328,372	0	G	A	50.0%	‑3.6 / 67.9	46	F131F (TTC→TTT)	wzxC	colanic acid exporter
Reads supporting (aligned to +/- strand):  ref base G (15/8);  new base A (14/9);  total (29/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

CCGCCAGCACCGCGCTGGTTTCGATCATGCCGATTTTGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATGTTTTAATCAACGGTGCCAGGTCCG  >  minE/1328302‑1328439                                                                       |                                                                    ccGCCAGCACCGCGCTGGTTTCGATCATGCCGATTTTGTTGAACTCCAGCTCTTTTTGCATCAACGCGCgg                                                                     <  1:3734936/71‑1 (MQ=255)     cAGCACCGCGCTGGTTTCGATCATGCCGATTTTGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGAAACt                                                                 <  1:3021435/71‑1 (MQ=255)              gcTGGTTTCGATCATGCCGATTTTGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGAAACTGTTGCCCGt                                                        >  1:1815732/1‑71 (MQ=255)              gcTGGTTTCGATCATGCCGATTTTGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGAAACTGTTGCCCGt                                                        >  1:2137085/1‑71 (MQ=255)                  gTTTCGATCATGCCGATTTTGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGAAACTGTTGCCCGTgggg                                                    <  1:546801/71‑1 (MQ=255)                    ttCGATCATGCCGATTTTGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGAAACTGTTGCCCGTgggg                                                    <  1:976099/69‑1 (MQ=255)                    ttCGATCATGCCGATTTTGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGAAACTGTTGCCCGTGGGGGa                                                  >  1:2339980/1‑71 (MQ=255)                     tCGATCATGCCGATTTTGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGAAACTGTTGCCCGTGGGGGAt                                                 >  1:3481066/1‑71 (MQ=255)                      cGATCATGCCGATTTTGTTGAACTCCAGCTCTTTTTGCATCAACGCTCGAAACTGTTGCCCGTggggg                                                   >  1:68311/1‑68 (MQ=255)                      cGATCATGCCGATTTTGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGAAACTGTTGCCCGTGGGGGa                                                  >  1:147973/1‑69 (MQ=255)                      cGATCATGCCGATTTTGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGAAACTGTTGCCCGTGGGGGa                                                  >  1:3308075/1‑69 (MQ=255)                      cGATCATGCCGATTTTGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGAAACTGTTGCCCGTGGGGGa                                                  >  1:2068190/1‑69 (MQ=255)                        aTCATGCCGATTTTGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGAAACt                                                                 <  1:2472050/52‑1 (MQ=255)                          cATGCCGATTTTGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGGAACTGTTGCCCGTGGGGGATTACca                                            >  1:2044198/1‑71 (MQ=255)                          cATGCCGATTTTGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGAAACTGt                                                               >  1:133907/1‑52 (MQ=255)                          cATGCCGATTTTGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGAAACTGt                                                               >  1:593467/1‑52 (MQ=255)                          cATGCCGATTTTGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGAAACTGt                                                               >  1:1784908/1‑52 (MQ=255)                          aaTGCCGATTTTGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGGAACTGTTGCCCGTGGGGGATTACca                                            >  1:199710/2‑71 (MQ=255)                                   tttGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGc                                     >  1:4107526/1‑69 (MQ=255)                                   tttGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGc                                     >  1:3716758/1‑69 (MQ=255)                                   tttGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGc                                     >  1:1381727/1‑69 (MQ=255)                                        tgAACTCCAGCTCTTTTTGCATCAACGCGCGGAACTGTTGc                                                            <  1:1698528/41‑1 (MQ=255)                                            cTCCAGCTCTTTTTGCATCAACGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATg                          <  1:3589955/71‑1 (MQ=255)                                            cTCCAGCTCTTTTTGCATCAACGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATg                          <  1:3318022/71‑1 (MQ=255)                                            cTCCAGCTCTTTTTGCATCAACGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATg                          <  1:1765665/71‑1 (MQ=255)                                                  ctctTTTTGCATCAACGCGCGGAACTGTTGTCCGTGGGGGATTACCACAAACGCCAGCGATAATGtttt                      <  1:1909330/69‑1 (MQ=255)                                                  ctctTTTTGCATCAACGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATGtttt                      <  1:1871662/69‑1 (MQ=255)                                                  ctctTTTTGCATCAACGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATGtttt                      <  1:1607301/69‑1 (MQ=255)                                                  ctctTTTTGCATCAACGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATGtttt                      <  1:784064/69‑1 (MQ=255)                                                     tttttGCATCAACGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATGTTTTAATc                  >  1:377321/1‑70 (MQ=255)                                                       tttGCATCAACGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATGTTTTAATCAAc               >  1:2616915/1‑71 (MQ=255)                                                       tttGCATCAACGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATGTTTTAATCAAc               >  1:3187817/1‑71 (MQ=255)                                                       tttGCATCAACGCGCGAAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATGTTTTAATCAAc               >  1:3421281/1‑71 (MQ=255)                                                          gCATCAACGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATGTTTTAATCAACGGt            >  1:1817016/1‑71 (MQ=255)                                                          gCATCAACGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATGTTTTAATCAACGGt            >  1:4041566/1‑71 (MQ=255)                                                          gCATCAACGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATGTTTTAATCAACGGt            >  1:4130760/1‑71 (MQ=255)                                                               aaCGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATg                          >  1:3199911/1‑52 (MQ=255)                                                               aaCGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATg                          >  1:3067128/1‑52 (MQ=255)                                                               aaCGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATg                          >  1:2629664/1‑52 (MQ=255)                                                               aaCGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATg                          >  1:1310036/1‑52 (MQ=255)                                                               aaCGCGCGAAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATGTTTTAATCAACGGTGCCAg       >  1:1581698/1‑71 (MQ=255)                                                               aaCGCGCGAAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGAGATAATGTTTTAATc                  >  1:2936129/1‑60 (MQ=255)                                                                    gcgAAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATGTTTTAATCAACGGTGCCAGGTCCg  <  1:2117255/71‑1 (MQ=255)                                                                    gcgAAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATGTTTTAATCAACGGTGCCAGGTCCg  <  1:3262695/71‑1 (MQ=255)                                                                    gcgAAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATGTTTTAATCAACGGTGCCAGGTCCg  <  1:2640323/71‑1 (MQ=255)                                                                    gcgAAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATGTTTTAATCAACGGTGCCAGGTCCg  <  1:593157/71‑1 (MQ=255)                                                                    gcgAAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATGTTTTAATCAACGGTGCCAGGTCCg  <  1:224505/71‑1 (MQ=255)                                                                       gAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATGt                         <  1:411747/45‑1 (MQ=255)                                                                       gAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATGt                         <  1:480805/45‑1 (MQ=255)                                                                       |                                                                    CCGCCAGCACCGCGCTGGTTTCGATCATGCCGATTTTGTTGAACTCCAGCTCTTTTTGCATCAACGCGCGGAACTGTTGCCCGTGGGGGATTACCACAAACGCCAGCGATAATGTTTTAATCAACGGTGCCAGGTCCG  >  minE/1328302‑1328439

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