Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	484,668	C→T	40.3%	G158G (GGC→GGT)	tolB →	periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	484,668	0	C	T	40.3%	21.7 / 32.9	52	G158G (GGC→GGT)	tolB	periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (24/7);  new base T (6/15);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 6.56e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAAGCAGTGGCTGCGTTATGCTGGTCATACCGCCAGTGATGAAGTGTTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATG  >  minE/484603‑484735                                                                  |                                                                    aaCAAGCAGTGGCTGCGTTATGCTGGTCATACCGCCAGTGACGAAGTATTTGAAAAACTGACCGGTAt                                                                   >  1:101264/1‑68 (MQ=255) aaCAAGCAGTGGCTGCGTTATGCTGGTCATACCGCCAGTGACGAAGTATTTGAAAAACTGACCGGTAt                                                                   >  1:2328061/1‑68 (MQ=255) aaCAAGCAGTGGCTGCGTTATGCTGGTCATACCGCCAGTGACGAAGTATTTGAAAAACTGACCGGTAt                                                                   >  1:2187822/1‑68 (MQ=255) aaCAAGCAGTGGCTGCGTTATGCTGGTCATACCGCCAGTGACGAAGTATTTGAAAAACTGACCGGTAt                                                                   >  1:2968138/1‑68 (MQ=255) aaCAAGCAGTGGCTGCGTTATGCTGGTCATACCGCCAGTGACGAAGTATTTGAAAAACTGACCGGTAt                                                                   >  1:1783469/1‑68 (MQ=255) aaCAAGCAGTGGCTGCGTTATGCTGGTCATACCGCCAGTGACGAAGTATTTGAAAAACTGACCGGTAt                                                                   >  1:2941738/1‑68 (MQ=255)                  ttATGCTGGTCATACCGCCAGTGACGAAGTATTTGAAAAACTGACCGGTATTAAAGGTGCGTTCCGTa                                                  <  1:2429957/68‑1 (MQ=255)                  ttATGCTGGTCATACCGCCAGTGACGAAGTATTTGAAAAACTGACCGGTATTAAAGGTGCGTTCCGTa                                                  <  1:252342/68‑1 (MQ=255)                  ttATGCTGGTCATACCGCCAGTGACGAAGTATTTGAAAAACTGACCGGTATTAAAGGTGCGTTCCGTa                                                  <  1:2549774/68‑1 (MQ=255)                  ttATGCTGGTCATACCGCCAGTGACGAAGTATTTGAAAAACTGACCGGTATTAAAGGTGCGTTCCGTa                                                  <  1:2202629/68‑1 (MQ=255)                  ttATGCTGGTCATACCGCCAGTGACGAAGTATTTGAAAAACTGACCGGTATTAAAGGTGCGTTCCGTa                                                  <  1:2151363/68‑1 (MQ=255)                  ttATGCTGGTCATACCGCCAGTGACGAAGTATTTGAAAAACTGACCGGTATTAAAGGTGCGTTCCGTa                                                  <  1:2556490/68‑1 (MQ=255)                  ttATGCTGGTCATACCGCCAGTGACGAAGTATTTGAAAAACTGACCGGTATTAAAGGTGCGTTCCGTa                                                  <  1:2602504/68‑1 (MQ=255)                  ttATGCTGGTCATACCGCCAGTGACGAAGTATTTGAAAAACTGACCGGTATTAAAGGTGCGTTCCGTa                                                  <  1:1842731/68‑1 (MQ=255)                  ttATGCTGGTCATACCGCCAGTGACGAAGTATTTGAAAAACTGACCGGTATTAAAGGTGCGTTCCGTa                                                  <  1:2959420/68‑1 (MQ=255)                  ttATGCTGGTCATACCGCCAGTGACGAAGTATTTGAAAAACTGACCGGTATTAAAGGTGCGTTCCGTa                                                  <  1:2992004/68‑1 (MQ=255)                  ttATGCTGGTCATACCGCCAGTGACGAAGTATTTGAAAAACTGACCGGTATTAAAGGTGCGTTCCGTa                                                  <  1:426881/68‑1 (MQ=255)                  ttATGCTGGTCATACCGCCAGTGACGAAGTATTTGAAAAACTGACCGGTATTAAAGGTGCGGTCCGTa                                                  <  1:1888/68‑1 (MQ=255)                  ttATGCTGGTCATACCGCCAGTGACGAAGTATTTGAAAAACTGACCGGTATTAAAGGTGCGGTCCGTa                                                  <  1:715168/68‑1 (MQ=255)                   tATGCTGGTCATACCGCCAGTGACGAAGTATTTGAAAAACTGACCGGTATTAAAGGTGCGTTCCGTa                                                  <  1:1675098/67‑1 (MQ=255)                       cTGGTCATACCGCCAGTGATGAAGTGTTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTa                                            >  1:2977244/1‑69 (MQ=255)                       cTGGTCATACCGCCAGTGATGAAGTGTTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTa                                            >  1:1407376/1‑69 (MQ=255)                       cTGGTCATACCGCCAGTGATGAAGTGTTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTa                                            >  1:2295828/1‑69 (MQ=255)                       cTGGTCATACCGCCAGTGATGAAGTGTTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTa                                            >  1:1636327/1‑69 (MQ=255)                       cTGGTCATACCGCCAGTGATGAAGTGTTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTa                                            >  1:1702821/1‑69 (MQ=255)                       cTGGTCATACCGCCAGTGATGAAGTGTTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTa                                            >  1:1706926/1‑69 (MQ=255)                       cTGGTCATACCGCCAGTGATGAAGTGTTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTAGTa                                            >  1:926583/1‑69 (MQ=255)                               aCCGCCAGTGATGAAGTGTTTGAAAAACTGACCTGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTAc                                    >  1:2619374/1‑69 (MQ=255)                               aCCGCCAGTGATGAAGTGTTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTAc                                    >  1:686047/1‑69 (MQ=255)                               aCCGCCAGTGATGAAGTGTTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTAc                                    >  1:701161/1‑69 (MQ=255)                               aCCGCCAGTGATGAAGTGTTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTAc                                    >  1:3303067/1‑69 (MQ=255)                               aCCGCCAGTGATGAAGTGTTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTAc                                    >  1:827444/1‑69 (MQ=255)                               aCCGCCAGTGATGAAGTGTTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTAc                                    >  1:3045102/1‑69 (MQ=255)                               aCCGCCAGTGATGAAGTGTTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTAc                                    >  1:1989411/1‑69 (MQ=255)                                        gatgaAGTGTTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAg                           >  1:2586897/1‑69 (MQ=255)                                           gaAGTGTTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGAcc                        >  1:2019118/1‑69 (MQ=255)                                              gtgtTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAAc                     >  1:1831973/1‑69 (MQ=255)                                              gtgtTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAAc                     >  1:1893786/1‑69 (MQ=255)                                              gtgtTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAAc                     >  1:342162/1‑69 (MQ=255)                                              gtgtTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAAc                     >  1:1415487/1‑69 (MQ=255)                                              gtgtTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAAc                     >  1:588524/1‑69 (MQ=255)                                                  ttGAAAAACTGACCGGTATTAAAGGTGCGTTCCGTa                                                  <  1:954780/36‑1 (MQ=255)                                                   tGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGt                                             >  1:3210512/1‑40 (MQ=255)                                                   tGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGt                                             >  1:2220510/1‑40 (MQ=255)                                                   tGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGt                                             >  1:534374/1‑40 (MQ=255)                                                     aaaaaCTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTT‑AGACCAACGGCGGTc              <  1:2486023/68‑1 (MQ=255)                                                      aaaaCTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCTGTc              <  1:3054094/68‑1 (MQ=255)                                                      aaaaCTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTc              <  1:684885/68‑1 (MQ=255)                                                      aaaaCTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTc              <  1:1513210/68‑1 (MQ=255)                                                      aaaaCTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTc              <  1:2508122/68‑1 (MQ=255)                                                      aaaaCTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTc              <  1:2447630/68‑1 (MQ=255)                                                                 gCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATg  <  1:519900/69‑1 (MQ=255)                                                                  |                                                                    AACAAGCAGTGGCTGCGTTATGCTGGTCATACCGCCAGTGATGAAGTGTTTGAAAAACTGACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATG  >  minE/484603‑484735

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