Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	912,048	A→T	24.9%	S54S (TCA→TCT)	ompG →	outer membrane porin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	912,048	0	A	T	24.9%	52.3 / 14.0	36	S54S (TCA→TCT)	ompG	outer membrane porin
Reads supporting (aligned to +/- strand):  ref base A (9/18);  new base T (7/2);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 4.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGCTCTGGCCTATTATCAGGAAGGGCCGG  >  minE/911985‑912112                                                                |                                                                 gATGTACGAAATAGAAAACGTCGAAGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCTGTCTa                                                             >  1:505241/1‑69 (MQ=255) gATGTACGAAATAGAAAACGTCGAAGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCTGTCTa                                                             >  1:1746220/1‑69 (MQ=255) gATGTACGAAATAGAAAACGTCGAAGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCTGTCTa                                                             >  1:1812920/1‑69 (MQ=255) gATGTACGAAATAGAAAACGTCGAAGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCTGTCTa                                                             >  1:1865563/1‑69 (MQ=255) gATGTACGAAATAGAAAACGTCGAAGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCTGTCTa                                                             >  1:2997741/1‑69 (MQ=255) gATGTACGAAATAGAAAACGTCGAAGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCTGTCTa                                                             >  1:2959482/1‑69 (MQ=255) gATGTACGAAATAGAAAACGTCGAAGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCTGTCTa                                                             >  1:2686459/1‑69 (MQ=255)     tACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTAtttt                                                         <  1:403411/69‑1 (MQ=255)     tACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTAtttt                                                         <  1:3369997/69‑1 (MQ=255)     tACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTAtttt                                                         <  1:2000355/69‑1 (MQ=255)     tACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTAtttt                                                         <  1:2551381/69‑1 (MQ=255)      aCGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTAtttt                                                         <  1:3080308/68‑1 (MQ=255)                       gAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCg                                       <  1:948967/69‑1 (MQ=255)                       gAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCg                                       <  1:2233417/69‑1 (MQ=255)                       gAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCg                                       <  1:198504/69‑1 (MQ=255)                                 gCGAAGATATGGATGGGCTGGCGGAGCCTTCTGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAAttt                             <  1:2447264/69‑2 (MQ=255)                                 gCGAAGATATGGATGGGCTGGCGGAGCCTTCTGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAAttt                             <  1:500949/69‑2 (MQ=255)                                 gCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGgaga                                 >  1:1393580/1‑65 (MQ=255)                                 gCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAAtt                              >  1:3064157/1‑68 (MQ=255)                                 gCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTg                             <  1:2368542/69‑1 (MQ=255)                                 gCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTg                             <  1:2214874/69‑1 (MQ=255)                                 gCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTg                             <  1:2148837/69‑1 (MQ=255)                                 gCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTg                             <  1:1044535/69‑1 (MQ=255)                                  cGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTg                             <  1:505498/68‑1 (MQ=255)                                        tatGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGct                           >  1:390063/1‑64 (MQ=255)                                        tatGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGCTCTgg                       >  1:2301681/1‑68 (MQ=255)                                        tatGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGCTCTgg                       >  1:2762866/1‑68 (MQ=255)                                        tatGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGCTCTgg                       >  1:2758790/1‑68 (MQ=255)                                         atggatggGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTg                                     >  1:1044721/1‑53 (MQ=255)                                         atggatggGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGCTCTGGcc                     >  1:3228746/1‑69 (MQ=255)                                         atggatggGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGCTCTGGcc                     >  1:1029177/1‑69 (MQ=255)                                                            cTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGCTCTGGCCTATTATCAGGAAGGGCCgg  <  1:1990643/69‑1 (MQ=255)                                                            cTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGCTCTGGCCTATTATCAGGAAGGGCCgg  <  1:456506/69‑1 (MQ=255)                                                            cTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGCTCTGGCCTATTATCAGGAAGGGCCgg  <  1:2095950/69‑1 (MQ=255)                                                            cTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGCTCTGGCCTATTATCAGGAAGGGCCgg  <  1:602630/69‑1 (MQ=255)                                                            cTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGCTCTGGCCTATTATCAGGAAGGGCCgg  <  1:823524/69‑1 (MQ=255)                                                                |                                                                 GATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGCTCTGGCCTATTATCAGGAAGGGCCGG  >  minE/911985‑912112

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