Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,078,055	A→G	67.6%	T154T (ACT→ACC)	uxaA ←	altronate hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,078,055	0	A	G	67.6%	38.8 / 22.8	50	T154T (ACT→ACC)	uxaA	altronate hydrolase
Reads supporting (aligned to +/- strand):  ref base A (8/8);  new base G (17/16);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

GCTGTGAGCAGCCGTAGGTGTGGCTGAAGAGGAACACGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGA  >  minE/2078012‑2078119                                            |                                                                 gCTGTGAGCAGCCGTAGGTGTGGCTGAAGAGGAACACGCCGTCAGTACCTTCGGCGTTGTTGGTCTCtt                                         <  1:3106581/69‑1 (MQ=255)    gtgAGCAGCCGTAGGTGTGGCTGAAGAGGAACACGCCGTCGGTACCTTcag                                                        >  1:699088/1‑49 (MQ=255)    gtgAGCAGCCGTAGGTGTGGCTGAAGAGGAACACGCCGTCGGTACCTTcag                                                        >  1:1023103/1‑49 (MQ=255)    gtgAGCAGCCGTAGGTGTGGCTGAAGAGGAACACGCCGTCGGTACCTTcag                                                        >  1:1372951/1‑49 (MQ=255)    gtgAGCAGCCGTAGGTGTGGCTGAAGAGGAACACGCCGTCGGTACCTTcag                                                        >  1:2629746/1‑49 (MQ=255)    gtgAGCAGCCGTAGGTGTGGCTGAAGAGGAACACGCCGTCGGTACCTTcag                                                        >  1:3159340/1‑49 (MQ=255)    gtgAGCAGCCGTAGGTGTGGCTGAAGAGGAACACGCCGTCGGTACCTTcag                                                        >  1:950493/1‑49 (MQ=255)    gtgAGCAGCCGTAGGTGTGGCTGAAGAGGAACACGCCGTCGGTACCTTcag                                                        >  1:1658686/1‑49 (MQ=255)    gtgAGCAGCCGTAGGTGTGGCTGAAGAGGAACACGCCGTCGGTACCTTcag                                                        >  1:2964967/1‑49 (MQ=255)    gtgAGCAGCCGTAGGTGTGGCTGAAGAGGAACACGCCGTCGGTACCTTcag                                                        >  1:2572529/1‑49 (MQ=255)    gtgAGCAGCCGTAGGTGTGGCTGAAGAGGAACACGCCGTCGGTACCTTcag                                                        >  1:2786354/1‑49 (MQ=255)    gtgAGCAGCCGTAGGTGTGGCTGAAGAGGAACACGCCGTCGGTACCTTcag                                                        >  1:2228404/1‑49 (MQ=255)    gtgAGCAGCCGTAGGTGTGGCTGAAGAGGAACACGCCGTCGGTACCTTcag                                                        >  1:2069650/1‑49 (MQ=255)                 ggTGTGGCTGAAGAGGAACACGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTg                         >  1:2567235/1‑69 (MQ=255)                   tgtgGCTGAAGAGGAACACGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTgg                        <  1:2947421/68‑1 (MQ=255)                    gtgGCTGAAGAGGAACACGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTgg                        <  1:709345/67‑1 (MQ=255)                               ggAACACGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGa           >  1:2084722/1‑69 (MQ=255)                               ggAACACGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGa           >  1:1852938/1‑69 (MQ=255)                               ggAACACGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGa           >  1:1685572/1‑69 (MQ=255)                               ggAACACGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGa           >  1:3159527/1‑69 (MQ=255)                               ggAACACGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGa           >  1:3329225/1‑69 (MQ=255)                                gAACACGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTgg                        <  1:1841839/55‑1 (MQ=255)                                  acacGCGGTCTGTACCTTCAGCGTTGTTGGTCTCTTTCTGGAAACGATTTTGGATCTGCCGCGCGATGc        <  1:3044091/69‑1 (MQ=255)                                  acacGCCGTCGGTGCCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc        <  1:1798262/69‑1 (MQ=255)                                  acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc        <  1:2997281/69‑1 (MQ=255)                                  acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc        <  1:728283/69‑1 (MQ=255)                                  acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc        <  1:2901926/69‑1 (MQ=255)                                  acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc        <  1:3058572/69‑1 (MQ=255)                                  acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc        <  1:2737380/69‑1 (MQ=255)                                  acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc        <  1:248017/69‑1 (MQ=255)                                  acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc        <  1:2067543/69‑1 (MQ=255)                                  acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc        <  1:818916/69‑1 (MQ=255)                                  acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc        <  1:1862837/69‑1 (MQ=255)                                  acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc        <  1:1727676/69‑1 (MQ=255)                                  acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc        <  1:867023/69‑1 (MQ=255)                                  acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc        <  1:1422785/69‑1 (MQ=255)                                  acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc        <  1:1039965/69‑1 (MQ=255)                                  acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTCGATCTGCCGCGCGATGc        <  1:180504/69‑1 (MQ=255)                                   cacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc        <  1:604527/68‑1 (MQ=255)                                   cacGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGcc       >  1:3428281/1‑69 (MQ=255)                                      gCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTg                   >  1:955038/1‑54 (MQ=255)                                      gCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCg      >  1:998354/1‑67 (MQ=255)                                      gCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGtt    >  1:1663247/1‑69 (MQ=255)                                      gCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGtt    >  1:1506129/1‑69 (MQ=255)                                      gCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGt     >  1:1911962/1‑68 (MQ=255)                                      gCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGt     >  1:846910/1‑68 (MQ=255)                                        cGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGa  <  1:3360240/69‑1 (MQ=255)                                        cGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGa  <  1:2281934/69‑1 (MQ=255)                                        cGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGa  <  1:835379/69‑1 (MQ=255)                                        cGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGa  <  1:1505307/69‑1 (MQ=255)                                            aGTACGTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc             <  1:278673/54‑1 (MQ=255)                                            aGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc             <  1:3164099/54‑1 (MQ=255)                                            aGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc             <  1:407657/54‑1 (MQ=255)                                            aGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc             <  1:274742/54‑1 (MQ=255)                                            aGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc             <  1:2696541/54‑1 (MQ=255)                                            aGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc             <  1:2612671/54‑1 (MQ=255)                                            aGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc             <  1:1114633/54‑1 (MQ=255)                                            aGTACCTTCGGCGTT‑TTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc             <  1:759826/53‑1 (MQ=255)                                            |                                                                 GCTGTGAGCAGCCGTAGGTGTGGCTGAAGAGGAACACGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGA  >  minE/2078012‑2078119

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