Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	125,489	A→G	36.9%	A159A (GCT→GCC)	gcd ←	glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	125,489	0	A	G	36.9%	22.7 / 27.0	46	A159A (GCT→GCC)	gcd	glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (23/6);  new base G (6/11);  total (29/17)
Fisher's exact test for biased strand distribution p-value = 4.50e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

TGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGATCTCCTGCGGATCGTTAAATCCGGCCC  >  minE/125433‑125556                                                         |                                                                    tGACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggc                                                         >  1:2395529/1‑68 (MQ=255) tGACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggc                                                         >  1:1891082/1‑68 (MQ=255) tGACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggc                                                         >  1:776173/1‑68 (MQ=255) tGACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggc                                                         >  1:1884412/1‑68 (MQ=255) tGACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggc                                                         >  1:2802582/1‑68 (MQ=255) tGACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTgg                                                          >  1:1444255/1‑68 (MQ=255)  gACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggcg                                                        <  1:949902/69‑3 (MQ=255)  gACCTTCCTGATTGCGACCATAGGCAGGCCAGTCCTGATCGGCAACGGGGGAGATGGCTTCTGCTggcg                                                        <  1:3258584/69‑3 (MQ=255)              tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTgg                                                     >  1:1693892/1‑60 (MQ=255)              tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg                                             >  1:1917902/1‑68 (MQ=255)              tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg                                             >  1:2446655/1‑68 (MQ=255)              tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg                                             >  1:1532830/1‑68 (MQ=255)              tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg                                             >  1:631306/1‑68 (MQ=255)              tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg                                             >  1:1333462/1‑68 (MQ=255)              tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg                                             >  1:703967/1‑68 (MQ=255)              tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg                                             >  1:749228/1‑68 (MQ=255)              tACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcg                                             >  1:1146608/1‑68 (MQ=255)               aCGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGgcgc                                            >  1:1745596/1‑68 (MQ=255)                 gACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAgg                                                          <  1:1627714/52‑1 (MQ=255)                               cAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGAtc                           <  1:623050/69‑1 (MQ=255)                               cAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGAtc                           <  1:1099748/69‑1 (MQ=255)                                aGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGAtc                           <  1:1144812/68‑1 (MQ=255)                                   ccTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGATCTCCt                       >  1:1915954/1‑69 (MQ=255)                                    cTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGATCTCCTg                      >  1:122985/1‑69 (MQ=255)                                       aTCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTaa                                        >  1:1870494/1‑48 (MQ=255)                                       aTCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTaa                                        >  1:401159/1‑48 (MQ=255)                                       aTCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTaa                                        >  1:1500874/1‑48 (MQ=255)                                       aTCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTaa                                        >  1:2948590/1‑48 (MQ=255)                                         cGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGATCTCCTGCGGAt                 >  1:2094718/1‑69 (MQ=255)                                               cGGGGGAGATGGCTTCTGCTGGCGTGGTATCGGCGCTTAAGGTGCCGCTGATTTCCTGCGGATCGTTaa           <  1:3254646/69‑1 (MQ=255)                                               cGGGGGAGATGGCTTCTGCTGGCGTGGTATCGGCGCTTAAGGTGCCGCTGATTTCCTGCGGATCGTTaa           <  1:1150068/69‑1 (MQ=255)                                               cGGGGGAGATGGCTTCTGCTGGCGTGGTATCGGCGCTTAAGGTGCCGCTGATTTCCTGCGGATCGTTaa           <  1:1194373/69‑1 (MQ=255)                                               cGGGGGAGATGGCTTCTGCTGGCGTGGTATCGGCGCTTAAGGTGCCGCTGATTTCCTGCGGATCGTTaa           <  1:715675/69‑1 (MQ=255)                                               cGGGGGAGATGGCTTCTGCTGGCGTGGTATCGGCGCTTAAGGTGCCGCTGATTTCCTGCGGATCGTTaa           <  1:1987975/69‑1 (MQ=255)                                               cGGGGGAGATGGCTTCTGCTGGCGTGGTATCGGCGCTTAAGGTGCCGCTGATTTCCTGCGGATCGTTaa           <  1:70235/69‑1 (MQ=255)                                               cGGGGGAGATGGCTTCTGCTGGCGTGGTATCGGCGCTTAAGGTGCCGCTGATTTCCTGCGGATCGTTaa           <  1:1466602/69‑1 (MQ=255)                                               cGGGGGAGATGGCTTCTGCTGGCGTGGTATCGGCGCTTAAGGTGCCGCTGATTTCCTGCGGATCGTTaa           <  1:542094/69‑1 (MQ=255)                                               cGGGGGAGATGGCTTCTGCTGGCGTGGTATCGGCGCTTAAGGTGCCGCTGATTTCCTGCGGATCGTTaa           <  1:3340407/69‑1 (MQ=255)                                                 ggggAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGATCTCCTGCGGATCGTTaaa          <  1:1445541/68‑1 (MQ=255)                                                 ggggAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGATCTCCTGCGGATCGTTaaa          <  1:882617/68‑1 (MQ=255)                                                        tAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGcc                                  >  1:832360/1‑37 (MQ=255)                                                        tAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGATCTCCTGCGGATCGTTAAATCCGGccc  >  1:258166/1‑69 (MQ=255)                                                        tAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGATCTCCTGCGGATCGTTAAATCCGGccc  >  1:3373635/1‑69 (MQ=255)                                                        tAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGATCTCCTGCGGATCGTTAAATCCGGccc  >  1:1182706/1‑69 (MQ=255)                                                        tAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGATCTCCTGCGGATCGTTAAATCCGGccc  >  1:2871683/1‑69 (MQ=255)                                                        tAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGAACTCCTGCGGATCGTTAAATCCGGccc  >  1:1619373/1‑69 (MQ=255)                                                         |                                                                    TGACCTTCCTGATTACGACCATAGGCAGGCCAGTCCTGATCGGCTACGGGGGAGATAGCTTCAGCAGGTGTGGCATCGGCGCTTAAGGTGCCGTTGATCTCCTGCGGATCGTTAAATCCGGCCC  >  minE/125433‑125556

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