Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,910,340	G→A	54.9%	P444P (CCG→CCA)	mgtA →	magnesium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,910,340	0	G	A	54.9%	9.1 / 57.9	71	P444P (CCG→CCA)	mgtA	magnesium transporter
Reads supporting (aligned to +/- strand):  ref base G (22/10);  new base A (15/24);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACC  >  minE/2910272‑2910395                                                                     |                                                        cTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCg                                                         >  1:2628139/1‑69 (MQ=255) cTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCg                                                         >  1:2967982/1‑69 (MQ=255) cTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCg                                                         >  1:901407/1‑69 (MQ=255) cTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCg                                                         >  1:1911409/1‑69 (MQ=255) cTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCg                                                         >  1:185285/1‑69 (MQ=255) cTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCg                                                         >  1:62049/1‑69 (MQ=255)                 aaGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTc                                         >  1:907743/1‑69 (MQ=255)                   gagTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGa                                                    >  1:104408/1‑56 (MQ=255)                   gagTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGa                                                    >  1:3250561/1‑56 (MQ=255)                   gagTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGa                                                    >  1:480978/1‑56 (MQ=255)                   gagTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGa                                                    >  1:606981/1‑56 (MQ=255)                   gagTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGa                                                    >  1:60761/1‑56 (MQ=255)                   gagTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGa                                                    >  1:1783693/1‑56 (MQ=255)                   gagTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGa                                                    >  1:629842/1‑56 (MQ=255)                   gagTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGa                                                    >  1:663465/1‑56 (MQ=255)                   gagTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCAGTTTGa                                                    >  1:2829887/1‑56 (MQ=255)                        aGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGa                                                    >  1:2179063/1‑51 (MQ=255)                        aGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGAt                                  >  1:1823875/1‑69 (MQ=255)                        aGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGAt                                  >  1:2806331/1‑69 (MQ=255)                        aGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGAt                                  >  1:949883/1‑69 (MQ=255)                        aGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGAt                                  >  1:940956/1‑69 (MQ=255)                        aGCGCGCTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGAt                                  <  1:2713539/69‑1 (MQ=255)                        aGCGCGCTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGAt                                  <  1:2604718/69‑1 (MQ=255)                        aGCGCGCTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGAt                                  <  1:3053384/69‑1 (MQ=255)                        aGCGCGCTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGAt                                  <  1:355498/69‑1 (MQ=255)                        aGCGCGCTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGAt                                  <  1:1366962/69‑1 (MQ=255)                        aGCGCGCTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGAt                                  <  1:1637166/69‑1 (MQ=255)                         gcgcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGAt                                  <  1:32111/68‑1 (MQ=255)                          cgcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTT‑ATTTCGAGCGTCGTCGGATGTc                               <  1:1199823/69‑1 (MQ=255)                           gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                               <  1:3282355/69‑1 (MQ=255)                           gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                               <  1:3035407/69‑1 (MQ=255)                           gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                               <  1:283418/69‑1 (MQ=255)                           gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                               <  1:935915/69‑1 (MQ=255)                           gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                               <  1:1107283/69‑1 (MQ=255)                           gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                               <  1:636429/69‑1 (MQ=255)                           gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                               <  1:2730255/69‑1 (MQ=255)                           gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                               <  1:2688809/69‑1 (MQ=255)                           gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                               <  1:1377263/69‑1 (MQ=255)                           gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                               <  1:1667235/69‑1 (MQ=255)                           gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                               <  1:2113485/69‑1 (MQ=255)                           gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                               <  1:2181902/69‑1 (MQ=255)                           gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                               <  1:2262728/69‑1 (MQ=255)                           gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                               <  1:2595114/69‑1 (MQ=255)                                           cGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGt                      >  1:3066839/1‑62 (MQ=255)                                           cGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTg                     >  1:2880474/1‑63 (MQ=255)                                           cGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTg                     >  1:2910121/1‑63 (MQ=255)                                           cGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTg                     >  1:1291749/1‑63 (MQ=255)                                           cGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTg                     >  1:1489081/1‑63 (MQ=255)                                           cGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTg                     >  1:550551/1‑63 (MQ=255)                                           cGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTg                     >  1:410001/1‑63 (MQ=255)                                              tGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAAt            <  1:2561078/69‑1 (MQ=255)                                              tGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAAt            <  1:904131/69‑1 (MQ=255)                                              tGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAAt            <  1:889771/69‑1 (MQ=255)                                              tGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAAt            <  1:2205570/69‑1 (MQ=255)                                              tGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAAt            <  1:3328332/69‑1 (MQ=255)                                              tGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAAt            <  1:2385753/69‑1 (MQ=255)                                              tGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAAt            <  1:2504214/69‑1 (MQ=255)                                                  aaaaaaTTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCg        <  1:824636/69‑1 (MQ=255)                                                  aaaaaaTTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCg        <  1:1639146/69‑1 (MQ=255)                                                  aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                         >  1:1128311/1‑52 (MQ=255)                                                  aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                         >  1:2380716/1‑52 (MQ=255)                                                  aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                         >  1:1155714/1‑52 (MQ=255)                                                  aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                         >  1:2937491/1‑52 (MQ=255)                                                  aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                         >  1:2694299/1‑52 (MQ=255)                                                  aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                         >  1:2962636/1‑52 (MQ=255)                                                  aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                         >  1:1465577/1‑52 (MQ=255)                                                  aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                         >  1:2991938/1‑52 (MQ=255)                                                  aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                         >  1:2720282/1‑52 (MQ=255)                                                  aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                         >  1:2143238/1‑52 (MQ=255)                                                  aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                         >  1:2381738/1‑52 (MQ=255)                                                  aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                         >  1:263906/1‑52 (MQ=255)                                                  aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtgg                          >  1:3139104/1‑51 (MQ=255)                                                  aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCgg                             >  1:1610016/1‑48 (MQ=255)                                                  aaaaaaTTGATGAGATTCCATTTGATTT‑‑AGCGTCGTCGGATGTCggtggt                         >  1:2666274/1‑50 (MQ=255)                                                       aTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGAt                                  <  1:3171727/38‑1 (MQ=38)                                                         tgatgaGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCAcc  <  1:2261261/68‑1 (MQ=255)                                                          gatgaGATTCCATTTGATTTCGAGCGTCGTCGGAt                                  <  1:2741647/35‑1 (MQ=255)                                                                     |                                                        CTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACC  >  minE/2910272‑2910395

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑