Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,187,440	T→C	32.5%	K139K (AAA→AAG)	nanT ←	sialic acid transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,187,440	0	T	C	32.5%	31.3 / 20.3	43	K139K (AAA→AAG)	nanT	sialic acid transporter
Reads supporting (aligned to +/- strand):  ref base T (15/14);  new base C (5/9);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-01

GACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATG  >  minE/2187383‑2187506                                                          |                                                                   gacgGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCt                                                         >  1:737438/1‑69 (MQ=255) gacgGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCt                                                         >  1:397084/1‑69 (MQ=255) gacgGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCt                                                         >  1:2635097/1‑69 (MQ=255) gacgGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCt                                                         >  1:2302203/1‑69 (MQ=255) gacgGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCCTTGTTACGCAGATGTTTTGGCCAGc                                                          >  1:2364701/1‑68 (MQ=255)  acgGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCt                                                         >  1:951724/1‑68 (MQ=255)  acgGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCt                                                         >  1:2595220/1‑68 (MQ=255)  acgGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCt                                                         >  1:2090065/1‑68 (MQ=255)  acgGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCt                                                         >  1:1920292/1‑68 (MQ=255)       ccccACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAAt                                                   <  1:1145151/69‑1 (MQ=255)         ccACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGa                                                 <  1:3052035/69‑1 (MQ=255)            ctgagaAGCCGGAAATCAGAAAACCACTGGCTTTGTTACTCAGATGCTTTGGCCAGCTTTCAATGACAt                                              >  1:1812555/3‑69 (MQ=255)            ctgagaAGCCGGAAATCAGAAAACCACTGGCTTTGTTACGCAGATGCTTTGGCCAGCTTTCAATGACAt                                              >  1:3005024/3‑69 (MQ=255)            ctgagaAGCCGGAAATCAGAAAACCACTGGCTTTGTTACGCAGATGCTTTGGCCAGCTTTCAATGACAt                                              >  1:2767981/3‑69 (MQ=255)            ctgagaAGCCGGAAATCAGAAAACCACTGGCTTTGTTACGCAGATGCTTTGGCCAGCTTTCAATGACAt                                              >  1:1453889/3‑69 (MQ=255)            ctgagaAGCCGGAAATCAGAAAACCACTGGCTTTGTTACGCAGATGCTTTGGCCAGCTTTCAATGACAt                                              >  1:1432890/3‑69 (MQ=255)              gagaAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCaa                                                    >  1:2688798/1‑61 (MQ=255)              gagaAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAg                                            >  1:3334648/1‑69 (MQ=255)                      ggAAATCAGAAAACCACTGGCTTTGTTACGCAGATGCTTTGGCCAGCTTTCAATGACATATGTGgcg                                      <  1:3401559/66‑1 (MQ=255)                      ggAAATCAGAAAACCACTGGCTTTGTTACGCAGATGCTTTGGCCAGCTTTCAATGACATAGGTGgcg                                      <  1:1138184/66‑1 (MQ=255)                      ggAAATCAGAAAACCACTGGCTTTGTTACGCAGATGCTTTGGCCAGCTTTCAATGACATAGGTGgcg                                      <  1:1191660/66‑1 (MQ=255)                      ggAAATCAGAAAACCACTGGCTTTGTTACGCAGATGCTTTGGCCAGCTTTCAATGACATAGGTGgcg                                      <  1:2931959/66‑1 (MQ=255)                      ggAAATCAGAAAACCACTGGCTTTGTTACGCAGATGCTTTGGCCAGCTTTCAATGACATAGGTGgcg                                      <  1:1947866/66‑1 (MQ=255)                      ggAAATCAGAAAACCACTGGCTTTGTTACGCAGATGCTTTGGCCAGCTTTCAATGACATAGGTGgcg                                      <  1:2373797/66‑1 (MQ=255)                      ggAAATCAGAAAACCACTGGCTTTGTTACGCAGATGCTTTGGCCAGCTTTCAATGACATAGGTGgcg                                      <  1:2029365/66‑1 (MQ=255)                      ggAAATCAGAAAACCACTGGCTTTGTTACGCAGATGCTTTGGCCAGCTTTCAATGACATAGGTGgcg                                      <  1:2467223/66‑1 (MQ=255)                      ggAAATCAGAAAACCACTGGCTTTGTTACGCAGATGCTTTGGCCAGCTTTCAATGACATAGGTGgcg                                      <  1:1713951/66‑1 (MQ=255)                               aaaaCCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGt                            <  1:2046515/68‑1 (MQ=255)                               aaaaCCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGt                            <  1:700850/68‑1 (MQ=255)                               aaaaCCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGt                            <  1:360404/68‑1 (MQ=255)                               aaaaCCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGt                            <  1:354379/68‑1 (MQ=255)                               aaaaCCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGt                            <  1:299829/68‑1 (MQ=255)                                   ccACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCa                       <  1:2118092/69‑1 (MQ=255)                                   ccACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCa                       <  1:2451950/69‑1 (MQ=255)                                           tttGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCa                <  1:251201/68‑1 (MQ=255)                                                  cGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATg        >  1:1402530/1‑69 (MQ=255)                                                         gTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGt                            <  1:2313218/42‑1 (MQ=255)                                                         gTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGt                            <  1:520109/42‑1 (MQ=255)                                                         gTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGt                            <  1:565104/42‑1 (MQ=255)                                                         gTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGt                            <  1:687913/42‑1 (MQ=255)                                                         gTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATg  >  1:1568235/1‑68 (MQ=255)                                                         gTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATg  >  1:1975604/1‑68 (MQ=255)                                                         gTTTTGCCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATg  >  1:1264022/1‑68 (MQ=255)                                                          |                                                                   GACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATG  >  minE/2187383‑2187506

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