Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,813,536	C→A	66.6%	A547A (GCG→GCT)	alaS ←	alanyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,813,536	0	C	A	66.6%	29.8 / 28.4	42	A547A (GCG→GCT)	alaS	alanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (7/7);  new base A (10/18);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

CAGGTGCGTTGCGGAGTGATTCAGACGAATACGGGCGCGACGAGCCTCATCAACATCAGCCTGCACCGCGTCGCCCACTTTCAGAGAACCCGCAGCAAGTTTACCGATGTGACCAATCGCCTGGCCGTATT  >  minE/1813470‑1813600                                                                   |                                                                 cAGGTGCGTTGCGGAGTGATTCAGACGAATACGGGCGCGACGAGCCTCATCAACATCAGCCTGCACcgc                                                                >  1:2428848/1‑69 (MQ=255) cAGGTGCGTTGCGGAGTGATTCAGACGAATACGGGCGCGACGAGCCTCATCAACATCAGCCTGCACcgc                                                                >  1:1257541/1‑69 (MQ=255) cAGGTGCGTTGCGGAGTGATTCAGACGAATACGGGCGCGACGAGCCTCATCAACATCAGCCTGCACcgc                                                                >  1:268037/1‑69 (MQ=255) cAGGTGCGTTGCGGAGTGATTCAGACGAATACGGGCGCGACGAGCCTCATCAACATCAGCCTGCACcgc                                                                >  1:2744393/1‑69 (MQ=255) cAGGTGCGTTGCGGAGTGATTCAGACGAATACGGGCGCGACGAGCCTCATCAACATCAGCCTGCACcgc                                                                >  1:2801252/1‑69 (MQ=255) cAGGTGCGTTGCGGAGTGATTCAGACGAATACGGGCGCGACGAGCCTCATCAACATCAGCCTGCACcgc                                                                >  1:1910385/1‑69 (MQ=255) cAGGTGCGTTGCGGAGTGATTCAGACGAATACGGGCGCGACGAGCCTCATCAACATCAGCCTGCACcgc                                                                >  1:2036768/1‑69 (MQ=255)      gCGTTGCGGAGTGATTCAGACGAATACGGGCGCGACGAGCCTCATCAACATCAGCCTGCACCGCGTCg                                                            >  1:2392747/1‑68 (MQ=255)      gCGTTGCGGAGTGATTCAGACGAATACGGGCGCGACGAGCCTCATCAACATCAGCCTGCACCGCGTCg                                                            >  1:1750110/1‑68 (MQ=255)                          cGAATACGGGCGCGACGAGCCTCATCAACATCAGCCTGCACCGCGTCGCCCACTTTCAGAGAACCCgc                                        >  1:1404634/1‑68 (MQ=255)                            aaTACGGGCGCGACGAGCCTCATCAACATCAGCCTGCACCGCGTCGCCCACTTTCAGAGAACCCgcagc                                     <  1:3273664/69‑1 (MQ=255)                              tACGGGCGCGACGAGCCTCATCAACATCAGCCTGCACCGCGTCGCCCACTTTCAGAGAACCCgcag                                      >  1:3194097/1‑66 (MQ=255)                              tACGGGCGCGACGAGCCTCATCAACATCAGCCTGCACCGCGTCGCCCACTTTCAGAGAACCCGCAGCaa                                   >  1:2617456/1‑69 (MQ=255)                              tACGGGCGCGACGAGCCTCATCAACATCAGCCTGCACCGCGTCGCCCACTTTCAGAGAACCCGCAGCaa                                   >  1:789795/1‑69 (MQ=255)                              tACGGGCGCGACGAGCCTCATCAACATCAGCCTGCACCGCGTCGCCCACTTTCAGAGAACCCGCAGCaa                                   >  1:69129/1‑69 (MQ=255)                                     gcgACGAGCCTCATCAACATCAGCCTGCACCGCGTCGCCCACTTTCAGAGAACACGCAGCAAGTTTAcc                            <  1:1169021/69‑1 (MQ=255)                                     gcgACGAGCCTCATCAACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCg                                         <  1:341646/56‑1 (MQ=255)                                     gcgACGAGCCTCATCAACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCg                                         <  1:2779429/56‑1 (MQ=255)                                     gcgACGAGCCTCATCAACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCg                                         <  1:418370/56‑1 (MQ=255)                                     gcgACGAGCCTCATCAACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCg                                         <  1:641314/56‑1 (MQ=255)                                     gcgACGAGCCTCATCAACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCg                                         <  1:810181/56‑1 (MQ=255)                                     gcgACGAGCCTCATCAACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCg                                         <  1:857682/56‑1 (MQ=255)                                     gcgACGAGCCTCATCAACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCg                                         <  1:928467/56‑1 (MQ=255)                                     gcgACGAGCCTCATCAACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCg                                         <  1:255201/56‑1 (MQ=255)                                     gcgACGAGCCTCATCAACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCg                                         <  1:2470152/56‑1 (MQ=255)                                     gcgACGAGCCTCATCAACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCg                                         <  1:243321/56‑1 (MQ=255)                                     gcgACGAGCCTCATCAACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCg                                         <  1:2353467/56‑1 (MQ=255)                                     gcgACGAGCCTCATCAACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCg                                         <  1:2257522/56‑1 (MQ=255)                                     gcgACGAGCCTCATCAACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCg                                         <  1:1201185/56‑1 (MQ=255)                                     gcgACGAGCCTCATCAACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCg                                         <  1:1354247/56‑1 (MQ=255)                                     gcgACGAGCCTCATCAACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCg                                         <  1:1445479/56‑1 (MQ=255)                                     gcgACGAGCCTCATCAACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCg                                         <  1:1716232/56‑1 (MQ=255)                                     gcgACGAGCCTCATCAACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACACg                                         <  1:2580648/56‑1 (MQ=255)                                       gacgaGCCTCATCAACATCAGCCTTAACCGCGTCGCCCACTTTCAGAGAACCCgcagc                                     <  1:3261298/58‑1 (MQ=255)                                               tcatcaACATCAGCCTGCACCGCGTCGCCCACTTTCAGAGAACCCGCAGCAAGTTTACCGATGTGACCa                  <  1:2106596/69‑1 (MQ=255)                                               tcatcaACATCAGCCTGCACCGCGTCGCCCACTTTCAGAGAACCCGCAGCAAGTTTACCGATGTGACCa                  <  1:2974442/69‑1 (MQ=255)                                               tcatcaACATCAGCCTGCACAGCTTCGCCCACTTTCAGAGAAc                                            >  1:627428/1‑43 (MQ=38)                                               tcatcaACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCGCAGCAAGTTTACCGATGTGAcc                   >  1:2435171/1‑68 (MQ=255)                                               tcatcaACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCGCAGCAAGTTTACCGATGTGAcc                   >  1:1276651/1‑68 (MQ=255)                                               tcatcaACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCGCAGCAAGTTTACCGATGTGAcc                   >  1:735518/1‑68 (MQ=255)                                               tcatcaACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCGCAGCAAGTTTACCGATGTGAcc                   >  1:40017/1‑68 (MQ=255)                                               tcatcaACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCGCAGCAAGTTTACCGATGTGAcc                   >  1:1850125/1‑68 (MQ=255)                                               tcatcaACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCGCAGCAAGTTTACCGATGTGAcc                   >  1:3346308/1‑68 (MQ=255)                                               tcatcaACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCGCAGCAAGTTTACCGATGTGAcc                   >  1:1880341/1‑68 (MQ=255)                                               tcatcaACATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCGCAGCAAGTTTACCGATGTGAcc                   >  1:2212360/1‑68 (MQ=255)                                               tcatcaACAACAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCGCAGCAAGTTTACCGATGTGAc                    >  1:2018279/1‑67 (MQ=255)                                                     aCATCAGCCTGCACAGCGTCGCCCACTTTCAGAGAACCCg                                         <  1:2894774/40‑1 (MQ=255)                                                               gCACCGCGTCGCCCACTTTCAGAGAACCCGCAGCAAGTTTACCGATGTGACCAATCGCCTGGCCGTAtt  <  1:1372798/69‑1 (MQ=255)                                                                 aCCGCGTCGCCCACTTTCAGAGAACCCGCAGCAAGTTTACCGATGTGACCAATCGCCTGGCCGTAtt  <  1:1018976/67‑1 (MQ=255)                                                                   |                                                                 CAGGTGCGTTGCGGAGTGATTCAGACGAATACGGGCGCGACGAGCCTCATCAACATCAGCCTGCACCGCGTCGCCCACTTTCAGAGAACCCGCAGCAAGTTTACCGATGTGACCAATCGCCTGGCCGTATT  >  minE/1813470‑1813600

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