Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,286,015	A→G	37.3%	G119G (GGT→GGC)	argB ←	acetylglutamate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,286,015	0	A	G	37.3%	31.8 / 40.4	67	G119G (GGT→GGC)	argB	acetylglutamate kinase
Reads supporting (aligned to +/- strand):  ref base A (35/7);  new base G (16/9);  total (51/16)
Fisher's exact test for biased strand distribution p-value = 8.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.96e-02

TAACCGTTCTCCAGCAAGGAGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCGAGAAACAAACCTACGGCCGC  >  minE/2285953‑2286083                                                               |                                                                     tAACCGTTCTCCAGCAAGGAGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGACCTAAc                                                                >  1:81346/1‑69 (MQ=255)  aaCCGTTCTCCAGCAAGGAGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGACCTAAct                                                               >  1:3424396/1‑69 (MQ=255)  aaCCGTTCTCCAGCAAGGAGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGACCTAAct                                                               >  1:1812376/1‑69 (MQ=255)  aaCCGTTCTCCAGCAAGGAGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGACCTAAct                                                               >  1:1778227/1‑69 (MQ=255)  aaCCGTTCTCCAGCAAGGAGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGACCTAAct                                                               >  1:3012168/1‑69 (MQ=255)  aaCCGTTCTCCAGCAAGGAGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGACCTAAct                                                               >  1:3029710/1‑69 (MQ=255)  aaCCGTTCTCCAGCAAGGAGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGACCTAAct                                                               >  1:3030863/1‑69 (MQ=255)  aaCCGTTCTCCAGCAAGGAGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGACCTAAct                                                               >  1:208930/1‑69 (MQ=255)  aaCCGTTCTCCAGCAAGGAGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGACCTAAct                                                               >  1:914015/1‑69 (MQ=255)  aaCCGTTCTCCAGCAAGGAGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGACCTAAct                                                               >  1:2148426/1‑69 (MQ=255)  aaCCGTTCTCCAGCAAGGAGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGACCTAAc                                                                >  1:858965/1‑68 (MQ=255)              gCAAGCTGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGGcc                                                                    >  1:2048249/1‑52 (MQ=255)              gCAAGCTGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGGCGTAACTCTTCATCAAGCt                                                   >  1:1180616/1‑69 (MQ=255)              gCAAGCTGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCt                                                   >  1:2221883/1‑69 (MQ=255)              gCAAGCTGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCt                                                   >  1:1688193/1‑69 (MQ=255)              gCAAGCTGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCt                                                   >  1:1006029/1‑69 (MQ=255)              gCAAGCTGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCt                                                   >  1:199293/1‑69 (MQ=255)              gCAAGCTGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCt                                                   >  1:26663/1‑69 (MQ=255)              gCAAGCTGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCt                                                   >  1:2702554/1‑69 (MQ=255)              gCAAGCTGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCt                                                   >  1:2941330/1‑69 (MQ=255)              gCAAGCTGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCt                                                   >  1:3020437/1‑69 (MQ=255)              gCAAGCTGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCt                                                   >  1:3181632/1‑69 (MQ=255)              gCAAGCTGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCt                                                   >  1:1464525/1‑69 (MQ=255)              gCAAGCTGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCt                                                   >  1:1358649/1‑69 (MQ=255)              gCAAGCTGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCt                                                   >  1:828445/1‑69 (MQ=255)              gCAAGCTGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCt                                                   >  1:1265692/1‑69 (MQ=255)              gCAAGCTGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCt                                                   >  1:998146/1‑69 (MQ=255)                        gATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGa                                                                      >  1:2564280/1‑40 (MQ=255)                        gATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACttt                                         >  1:1031148/1‑69 (MQ=255)                              cTTAGGCGAACCTGGCTGCGCCAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCt                                   >  1:3236377/1‑69 (MQ=255)                                tAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGt                                 <  1:276358/69‑1 (MQ=255)                                tAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGt                                 <  1:1921760/69‑1 (MQ=255)                                tAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGt                                 <  1:2128673/69‑1 (MQ=255)                                tAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGt                                 <  1:3073874/69‑1 (MQ=255)                                tAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGt                                 <  1:258831/69‑1 (MQ=255)                                tAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGt                                 <  1:2359161/69‑1 (MQ=255)                                tAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGt                                 <  1:3299974/69‑1 (MQ=255)                                tAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGt                                 <  1:3427310/69‑1 (MQ=255)                                tAGGCGAACCTGGCTGCGCCAGTCCAACATGGCCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGt                                 <  1:977938/69‑1 (MQ=255)                                           ggCTGCGCCAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCa                         >  1:2644769/1‑66 (MQ=255)                                              tGCGCCAGTCCGACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCgaga                   >  1:2002034/1‑69 (MQ=255)                                              tGCGCCAGTCCAACATGACCTAACTCTTCATCAAg                                                     >  1:3240709/1‑35 (MQ=255)                                              tGCGCCAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCgaga                   >  1:325292/1‑69 (MQ=255)                                              tGCGCCAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCgaga                   >  1:3214156/1‑69 (MQ=255)                                              tGCGCCAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCgaga                   >  1:1380621/1‑69 (MQ=255)                                              tGCGCCAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCgaga                   >  1:3029363/1‑69 (MQ=255)                                              tGCGCCAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCgaga                   >  1:2721081/1‑69 (MQ=255)                                              tGCGCCAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCgaga                   >  1:2352517/1‑69 (MQ=255)                                              tGCGCCAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCgaga                   >  1:2359912/1‑69 (MQ=255)                                              tGCGCCAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCgag                    >  1:1518907/1‑68 (MQ=255)                                               gcgcCAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCgaga                   >  1:2755562/1‑68 (MQ=255)                                               gcgcCAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCgaga                   >  1:2238962/1‑68 (MQ=255)                                               gcgcCAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCgaga                   >  1:996200/1‑68 (MQ=255)                                               gcgcCAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCgaga                   >  1:1449809/1‑68 (MQ=255)                                               gcgcCAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCgaga                   >  1:2021770/1‑68 (MQ=255)                                               gcgcCAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCgaga                   >  1:1932007/1‑68 (MQ=255)                                                   cAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCGAGaaacaa              <  1:2903308/69‑1 (MQ=255)                                                   cAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCGAGaaacaa              <  1:1293814/69‑1 (MQ=255)                                                   cAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCGAGaaacaa              <  1:2230392/69‑1 (MQ=255)                                                   cAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCGAGaaacaa              <  1:3405321/69‑1 (MQ=255)                                                    aGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCGAGaaacaa              <  1:3115844/68‑1 (MQ=255)                                                    aGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCGAGaaacaa              <  1:1716968/68‑1 (MQ=255)                                                    aGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCGAGaaacaa              <  1:1112500/68‑1 (MQ=255)                                                       ccAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCGAGAAACAAACCt          >  1:1950088/1‑69 (MQ=255)                                                       ccAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCGAGAAACAAACCt          >  1:1697217/1‑69 (MQ=255)                                                        cAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCGAGAAACAAACCTa         >  1:2291570/1‑69 (MQ=255)                                                        cAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCGAGAAACAAACCTa         >  1:1106002/1‑69 (MQ=255)                                                        cAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCGAGAAACAAACCTa         >  1:1592048/1‑69 (MQ=255)                                                               aCCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCGAGAAACAAACCTACGGCCg   >  1:2937438/1‑68 (MQ=255)                                                               aCCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCGAGAAACAAACCTACGGCCGc  >  1:902173/1‑69 (MQ=255)                                                               aCCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCGAGAAACAAACCTACGGCCGc  >  1:3414992/1‑69 (MQ=255)                                                               aCCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCGAGAAACAAACCTACGGCCGc  >  1:1511887/1‑69 (MQ=255)                                                               aCCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCGAGAAACAAACCTACGGCCGc  >  1:2724680/1‑69 (MQ=255)                                                               aCCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCGAGAAACAAACCTACGGCCGc  >  1:2116909/1‑69 (MQ=255)                                                               |                                                                     TAACCGTTCTCCAGCAAGGAGTTGATAAGCTTAGGCGAACCTGGCTGCGCCAGTCCAACATGACCTAACTCTTCATCAAGCTGGGTCACTTTGACGCTGTCGCCGTCACCGAGAAACAAACCTACGGCCGC  >  minE/2285953‑2286083

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