Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,771,782	A→G	36.6%	P369P (CCA→CCG)	aceB →	malate synthase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,771,782	0	A	G	36.6%	52.0 / 56.2	90	P369P (CCA→CCG)	aceB	malate synthase A
Reads supporting (aligned to +/- strand):  ref base A (33/24);  new base G (8/25);  total (41/49)
Fisher's exact test for biased strand distribution p-value = 2.28e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

AAAGCGGATAAATCGCTGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACATTCTCGGCTCCCGTAAAAATCAGCTTGAA  >  minE/2771723‑2771848                                                            |                                                                   aaaGCGGATAAATCGCTGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCg                                                           >  1:404921/1‑69 (MQ=255) aaaGCGGATAAATCGCTGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCg                                                           >  1:1781721/1‑69 (MQ=255) aaaGCGGATAAATCGCTGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCg                                                           >  1:3116218/1‑69 (MQ=255)   aGCGGATAAATCGCTGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGa                                                         <  1:1924446/69‑1 (MQ=255)   aGCGGATAAATCGCTGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGa                                                         <  1:1566440/69‑1 (MQ=255)   aGCGGATAAATCGCTGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGa                                                         <  1:2112266/69‑1 (MQ=255)   aGCGGATAAATCGCTGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGa                                                         <  1:1275485/69‑1 (MQ=255)         tAAATCGCTGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGc                                                   <  1:1444754/69‑1 (MQ=255)         tAAATCGCTGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGc                                                   <  1:2347376/69‑1 (MQ=255)         tAAATCGCTGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGc                                                   <  1:1143453/69‑1 (MQ=255)               gCTGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGa                                                         <  1:994758/57‑1 (MQ=255)                 tGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCgg                                           <  1:2403996/69‑1 (MQ=255)                 tGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCgg                                           <  1:657401/69‑1 (MQ=255)                 tGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCgg                                           <  1:74631/69‑1 (MQ=255)                 tGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCgg                                           <  1:2602285/69‑1 (MQ=255)                 tGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCgg                                           <  1:3312033/69‑1 (MQ=255)                 tGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCgg                                           <  1:1236432/69‑1 (MQ=255)                 tGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCgg                                           <  1:1515541/69‑1 (MQ=255)                 tGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCgg                                           <  1:954464/69‑1 (MQ=255)                           tAACTGTCTCGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGa                                 <  1:3209891/69‑1 (MQ=255)                           tAACGGTCTCGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGa                                 <  1:3230379/69‑1 (MQ=255)                           tAACGGTCTCGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGa                                 <  1:760652/69‑1 (MQ=255)                           tAACGGTCTCGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGa                                 <  1:1950288/69‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCtt                                                              >  1:3155108/1‑39 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTg                                                             >  1:3307080/1‑40 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTg                                                             >  1:2751636/1‑40 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTg                                                             >  1:201934/1‑40 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTg                                                             >  1:1784732/1‑40 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTg                                                             >  1:883691/1‑40 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTg                                                             >  1:3132117/1‑40 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTg                                                             >  1:319827/1‑40 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:2606574/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:2710344/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:1991932/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:2811075/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:296773/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:3061773/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:358645/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:2347834/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:2306647/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:1200202/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:2155402/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:1272674/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:488774/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:13227/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:1419557/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:1975842/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:1830969/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:1819967/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:795188/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:147389/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:1585960/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:1568827/44‑1 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGAc                                >  1:37623/1‑69 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGAc                                >  1:298475/1‑69 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGAc                                >  1:1970250/1‑69 (MQ=255)                            aaCGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGAc                                >  1:332270/1‑69 (MQ=255)                            aaCGGTCACGATGG‑‑CATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACAt                              >  1:996563/1‑69 (MQ=255)                            aaCGGTCACCATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:472684/44‑1 (MQ=255)                             aCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:3234792/43‑1 (MQ=255)                             aCGGTCACGATGGCACATGGATCGCTCACCCGGGCCTTGCTGa                                                         <  1:1057287/43‑1 (MQ=255)                             aCGGTCACGATGGCACATGGATCGCTCACCCAGGc                                                                 >  1:2133296/1‑35 (MQ=255)                             aCGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCa                                                  <  1:2472894/50‑1 (MQ=255)                             aCGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACa                               >  1:2771966/1‑69 (MQ=255)                             aCGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACa                               >  1:2600212/1‑69 (MQ=255)                             aCGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACa                               >  1:2538855/1‑69 (MQ=255)                             aCGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACa                               >  1:1911627/1‑69 (MQ=255)                             aCGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACa                               >  1:2306914/1‑69 (MQ=255)                             aCGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACa                               >  1:522088/1‑69 (MQ=255)                                        ggCACATGGATCGCTCACCCAGGCCTTGCGGACACGGc                                                   <  1:1051860/38‑1 (MQ=255)                                         gaacaTGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACAtt                             >  1:2314901/3‑59 (MQ=255)                                         gCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACAtt                             <  1:320556/59‑1 (MQ=255)                                         gCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACAtt                             >  1:1018202/1‑59 (MQ=255)                                         gCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACAtt                             >  1:2238309/1‑59 (MQ=255)                                          cacaTGGATCGCTCACCCAGGCCTTGCGGACACGGc                                                   >  1:3385341/1‑36 (MQ=255)                                          cacaTGGATCGCTCACCCAGGCCTTGCGGACACGGc                                                   >  1:1632564/1‑36 (MQ=255)                                          cacaTGGATCGCTCACCCAGGCCTTGCGGACACGGc                                                   >  1:2355267/1‑36 (MQ=255)                                                   cGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACATTCTCGGCTCCCGTAAAAATCa         >  1:2989542/1‑69 (MQ=255)                                                   cGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACATTCTCGGCTCCCGTAAAAATCa         >  1:1992678/1‑69 (MQ=255)                                                   cGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACATTCTCGGCTCCCGTAAAAATCa         >  1:351140/1‑69 (MQ=255)                                                        aCCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGa                                 <  1:530173/40‑1 (MQ=255)                                                           cAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACATTCTCGGCTCCCGTAAAAATCAGCTTGaa  >  1:367041/1‑68 (MQ=255)                                                           cAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACATTCTCGGCTCCCGTAAAAATCAGCTTGaa  >  1:190230/1‑68 (MQ=255)                                                           cAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACATTCTCGGCTCCCGTAAAAATCAGCTTGaa  >  1:1586058/1‑68 (MQ=255)                                                           cAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACATTCTCGGCTCCCGTAAAAATCAGCTTGaa  >  1:1339228/1‑68 (MQ=255)                                                           cAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACATTCTCGGCTCCCGTAAAAATCAGCTTGaa  >  1:90220/1‑68 (MQ=255)                                                           cAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACATTCTCGGCTCCCGTAAAAATCAGCTTGaa  >  1:924383/1‑68 (MQ=255)                                                           cAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACATTCTCGGCTCCCGTAAAAATCAGCTTGaa  >  1:123382/1‑68 (MQ=255)                                                           cAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACATTCTCGGCTCCCGTAAAAATCAGCTTGaa  >  1:2980405/1‑68 (MQ=255)                                                           cAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACATTCTCGGCTCCCGTAAAAATCAGCTTGaa  >  1:1048966/1‑68 (MQ=255)                                                            |                                                                   AAAGCGGATAAATCGCTGGAAGCCAATAACGGTCACGATGGCACATGGATCGCTCACCCAGGCCTTGCGGACACGGCAATGGCGGTATTCAACGACATTCTCGGCTCCCGTAAAAATCAGCTTGAA  >  minE/2771723‑2771848

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