Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,910,358	C→T	44.3%	R450R (CGC→CGT)	mgtA →	magnesium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,910,358	0	C	T	44.3%	23.2 / 68.0	88	R450R (CGC→CGT)	mgtA	magnesium transporter
Reads supporting (aligned to +/- strand):  ref base C (29/20);  new base T (15/24);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 8.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGC  >  minE/2910295‑2910422                                                                |                                                                 aGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGAt                                                             >  1:949883/1‑69 (MQ=255) aGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGAt                                                             >  1:940956/1‑69 (MQ=255) aGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGAt                                                             >  1:1823875/1‑69 (MQ=255) aGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGAt                                                             >  1:2806331/1‑69 (MQ=255) aGCGCGCTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGAt                                                             <  1:2604718/69‑1 (MQ=255) aGCGCGCTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGAt                                                             <  1:355498/69‑1 (MQ=255) aGCGCGCTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGAt                                                             <  1:1637166/69‑1 (MQ=255) aGCGCGCTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGAt                                                             <  1:3053384/69‑1 (MQ=255) aGCGCGCTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGAt                                                             <  1:1366962/69‑1 (MQ=255) aGCGCGCTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGAt                                                             <  1:2713539/69‑1 (MQ=255)  gcgcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGAt                                                             <  1:32111/68‑1 (MQ=255)   cgcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTT‑ATTTCGAGCGTCGTCGGATGTc                                                          <  1:1199823/69‑1 (MQ=255)    gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                                                          <  1:283418/69‑1 (MQ=255)    gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                                                          <  1:2730255/69‑1 (MQ=255)    gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                                                          <  1:2595114/69‑1 (MQ=255)    gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                                                          <  1:3035407/69‑1 (MQ=255)    gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                                                          <  1:2688809/69‑1 (MQ=255)    gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                                                          <  1:2262728/69‑1 (MQ=255)    gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                                                          <  1:1107283/69‑1 (MQ=255)    gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                                                          <  1:2181902/69‑1 (MQ=255)    gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                                                          <  1:3282355/69‑1 (MQ=255)    gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                                                          <  1:2113485/69‑1 (MQ=255)    gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                                                          <  1:935915/69‑1 (MQ=255)    gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                                                          <  1:1667235/69‑1 (MQ=255)    gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                                                          <  1:636429/69‑1 (MQ=255)    gcgcTCGCTGGCCAGCCGTTGGCAAAAAATTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTc                                                          <  1:1377263/69‑1 (MQ=255)                    cGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGt                                                 >  1:3066839/1‑62 (MQ=255)                    cGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTg                                                >  1:410001/1‑63 (MQ=255)                    cGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTg                                                >  1:550551/1‑63 (MQ=255)                    cGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTg                                                >  1:1291749/1‑63 (MQ=255)                    cGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTg                                                >  1:2880474/1‑63 (MQ=255)                    cGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTg                                                >  1:2910121/1‑63 (MQ=255)                    cGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTg                                                >  1:1489081/1‑63 (MQ=255)                       tGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAAt                                       <  1:2504214/69‑1 (MQ=255)                       tGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAAt                                       <  1:904131/69‑1 (MQ=255)                       tGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAAt                                       <  1:889771/69‑1 (MQ=255)                       tGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAAt                                       <  1:2561078/69‑1 (MQ=255)                       tGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAAt                                       <  1:3328332/69‑1 (MQ=255)                       tGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAAt                                       <  1:2385753/69‑1 (MQ=255)                       tGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAAt                                       <  1:2205570/69‑1 (MQ=255)                           aaaaaaTTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCg                                   <  1:824636/69‑1 (MQ=255)                           aaaaaaTTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCg                                   <  1:1639146/69‑1 (MQ=255)                           aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                                                    >  1:2937491/1‑52 (MQ=255)                           aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                                                    >  1:2962636/1‑52 (MQ=255)                           aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                                                    >  1:2991938/1‑52 (MQ=255)                           aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                                                    >  1:2720282/1‑52 (MQ=255)                           aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                                                    >  1:1128311/1‑52 (MQ=255)                           aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                                                    >  1:2694299/1‑52 (MQ=255)                           aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                                                    >  1:263906/1‑52 (MQ=255)                           aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                                                    >  1:1465577/1‑52 (MQ=255)                           aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                                                    >  1:2143238/1‑52 (MQ=255)                           aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                                                    >  1:1155714/1‑52 (MQ=255)                           aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                                                    >  1:2381738/1‑52 (MQ=255)                           aaaaaaTTGATGAGATTCCATTTGATTTCGAGCGTCGTCGGATGTCggtggt                                                    >  1:2380716/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)                                                     ttCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAgg          <  1:793049/68‑1 (MQ=255)                                                     ttCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAgg          <  1:761401/68‑1 (MQ=255)                                                     ttCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAgg          <  1:2627048/68‑1 (MQ=255)                                                     ttCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAgg          <  1:3071541/68‑1 (MQ=255)                                                           cGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATc                          >  1:244788/1‑46 (MQ=255)                                                           cGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTTGTTTGCAAAGGTg        >  1:942780/1‑64 (MQ=255)                                                           cGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTg        >  1:3205066/1‑64 (MQ=255)                                                           cGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTg        <  1:2258652/64‑1 (MQ=255)                                                           cGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTg        >  1:127577/1‑64 (MQ=255)                                                           cGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTg        <  1:2633952/64‑1 (MQ=255)                                                           cGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTg        >  1:2885670/1‑64 (MQ=255)                                                           cGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTg        <  1:2337479/64‑1 (MQ=255)                                                            gTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAg                                 >  1:3331224/1‑38 (MQ=255)                                                            gTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGc  >  1:3312041/1‑69 (MQ=255)                                                            gTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGc  >  1:933357/1‑69 (MQ=255)                                                            gTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGc  >  1:1278135/1‑69 (MQ=255)                                                            gTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGc  >  1:798593/1‑69 (MQ=255)                                                            gTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGc  >  1:944393/1‑69 (MQ=255)                                                            gTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGc  >  1:1632006/1‑69 (MQ=255)                                                            gTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGc  >  1:187354/1‑69 (MQ=255)                                                            gTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGc  >  1:1999997/1‑69 (MQ=255)                                                            gTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGc  <  1:3317776/69‑1 (MQ=255)                                                            gTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGc  >  1:3316103/1‑69 (MQ=255)                                                            gTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGc  <  1:2156684/69‑1 (MQ=255)                                                            gTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGc  >  1:3280236/1‑69 (MQ=255)                                                            gTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGc  >  1:2792115/1‑69 (MQ=255)                                                            gTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGc  <  1:2780422/69‑1 (MQ=255)                                                            gTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCAACATCAGCTGGTTTGCAAAGGTGCATTGc  >  1:1277393/1‑69 (MQ=255)                                                                |                                                                 AGCGCGCTCGCTGGCCAGTCGTTGGCAAAAAATTGATGAGATTCCGTTTGATTTCGAGCGTCGCCGGATGTCGGTGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGC  >  minE/2910295‑2910422

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