Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	125,498	A→T	36.9%	P156P (CCT→CCA)	gcd ←	glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	125,498	0	A	T	36.9%	32.7 / 37.3	46	P156P (CCT→CCA)	gcd	glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (23/6);  new base T (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.75e-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: ‑