Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,664,511	A→G	44.8%	G237G (GGT→GGC)	guaA ←	GMP synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,664,511	0	A	G	44.8%	6.1 / 45.6	58	G237G (GGT→GGC)	guaA	GMP synthetase
Reads supporting (aligned to +/- strand):  ref base A (18/14);  new base G (12/14);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

TACGCAAGTCAGGTTTTTACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATC  >  minE/1664448‑1664576                                                                |                                                                  tACGCAAGTCAGGTTTTTACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCAcacc                                                                 >  1:2135595/1‑66 (MQ=255) tACGCAAGTCAGGTTTTTACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCAcacc                                                                 >  1:3086760/1‑66 (MQ=255)  aCGCAAGTCAGGTTTTTACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACAccacca                                                             <  1:981703/69‑1 (MQ=255)  aCGCAAGTCAGGTTTTTACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACAccacca                                                             <  1:962511/69‑1 (MQ=255)  aCGCAAGTCAGGTTTTTACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACAccacca                                                             <  1:2403898/69‑1 (MQ=255)  aCGCAAGTCAGGTTTTTACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACAccacca                                                             <  1:886519/69‑1 (MQ=255)  aCGCAAGTCAGGTTTTTACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACAccacca                                                             <  1:1728454/69‑1 (MQ=255)              tttttACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGATAGGCCGAgg                                                 >  1:392890/1‑69 (MQ=255)              tttttACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAgg                                                 >  1:2450413/1‑69 (MQ=255)              tttttACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAgg                                                 >  1:843684/1‑69 (MQ=255)              tttttACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAgg                                                 >  1:2888129/1‑69 (MQ=255)                  tACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATg                                              <  1:1945038/68‑1 (MQ=255)                  tACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATg                                              <  1:3251797/68‑1 (MQ=255)                  tACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATg                                              <  1:3308472/68‑1 (MQ=255)                   aCCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGa                                             <  1:852185/68‑1 (MQ=255)                   aCCGATAGCGCGGTGCAACAGCATTGCGGTAACGGAGGAATCCACGCTACCAGAGAGGCCGAGGATGa                                             <  1:2329940/68‑1 (MQ=255)                   aCCGATAGCGCGGTGCAACAGCATTGCGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGa                                             <  1:2545695/68‑1 (MQ=255)                   aCCGATAGCGCGGTGCAACAGCATTGCGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGa                                             <  1:3009534/68‑1 (MQ=255)                   aCCGATAGCGCGGTGCAACAGCATTGCGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGa                                             <  1:2999162/68‑1 (MQ=255)                   aCCGATAGCGCGGTGCAACAGCATTGCGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGa                                             <  1:368756/68‑1 (MQ=255)                   aCCGATAGCGCGGTGCAACAGCATTGCGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGa                                             <  1:1075695/68‑1 (MQ=255)                   aCCGATAGCGCGGTGCAACAGCATTGCGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGa                                             <  1:2483128/68‑1 (MQ=255)                   aCCGATAGCGCGGTGCAACAGCATTGCGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGa                                             <  1:89895/68‑1 (MQ=255)                   aCCGATAGCGCGGTGCAACAGCATTGCGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGa                                             <  1:2306919/68‑1 (MQ=255)                   aCCGATAGCGCGGTGCAACAGCATTGCGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGa                                             <  1:2130448/68‑1 (MQ=255)                   aCCGATAGCGCGGTGCAACAGCATTGCGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGa                                             <  1:1785628/68‑1 (MQ=255)                   aCCGATAGCGCGGTGCAACAGCATTGCGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGa                                             <  1:1356068/68‑1 (MQ=255)                   aCCGATAGCGCGGTGCAAAAGCATTGCGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGa                                             <  1:3166251/68‑1 (MQ=255)                                       gcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCt                        >  1:2483079/1‑69 (MQ=255)                                           tGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTcgcg                    <  1:626416/69‑1 (MQ=255)                                           tGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTcgcg                    <  1:2001344/69‑1 (MQ=255)                                           tGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTcgcg                    <  1:3172429/69‑1 (MQ=255)                                           tGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTcgcg                    <  1:2845973/69‑1 (MQ=255)                                           tGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTcgcg                    <  1:2181118/69‑1 (MQ=255)                                           tGCGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGa                                             <  1:2349616/44‑1 (MQ=38)                                            gCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGAt                                               >  1:1686081/1‑41 (MQ=255)                                            gCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGAt                                               >  1:2182260/1‑41 (MQ=255)                                            gCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGAt                                               >  1:2238589/1‑41 (MQ=255)                                            gCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGAt                                               >  1:118422/1‑41 (MQ=255)                                            gCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGAt                                               >  1:1990306/1‑41 (MQ=255)                                            gCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGAt                                               >  1:2412355/1‑41 (MQ=255)                                            gCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGAt                                               >  1:2985053/1‑41 (MQ=255)                                            gCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGAt                                               >  1:207481/1‑41 (MQ=255)                                            gCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGAt                                               >  1:314260/1‑41 (MQ=255)                                            gCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCgg                   >  1:2802559/1‑69 (MQ=255)                                            gCGGTTACGGAGGAATCCACACCACAAGAGAGGCCGAGGAt                                               >  1:751400/1‑41 (MQ=255)                                            gCGGTTACGGAGGAATCCACACCA‑CAGAGAGGCCGAGGa                                                >  1:2257580/1‑39 (MQ=38)                                             cGGTAACGGAGGACTCCACGCCACCAGAGCGGCCGAGGATGACTTTATCGTCGCCCACTTGCTCGCgg                   >  1:2295416/1‑68 (MQ=255)                                             cGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGACTTTATCGTCGcc                                >  1:1225185/1‑55 (MQ=39)                                             cGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCCACTTGCTCGCgg                   >  1:1449993/1‑68 (MQ=255)                                             cGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCCACTTGCTCGCGGa                  >  1:2859770/1‑69 (MQ=255)                                             cGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCCACTTGCTCGCGGa                  >  1:134153/1‑69 (MQ=255)                                             cGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCCACTTGCTCGCGGa                  >  1:2834563/1‑69 (MQ=255)                                             cGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCCACTTGCTCGCGGa                  >  1:2490042/1‑69 (MQ=255)                                             cGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCCACTTGCTCGCGGa                  >  1:1430416/1‑69 (MQ=255)                                             cGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCCACTTGCTCGCGGa                  >  1:2788829/1‑69 (MQ=255)                                             cGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCCACTTGCTCGCGGa                  >  1:2775751/1‑69 (MQ=255)                                             cGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCCACTTGCTCGCGGa                  >  1:79854/1‑69 (MQ=255)                                             cGGTAACGGAGGAATCCACGCCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCAACTTGCTCGCGGa                  >  1:3177291/1‑69 (MQ=255)                                                          aTCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGc     >  1:2932894/1‑69 (MQ=255)                                                             cacaCCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCAtc  >  1:527215/1‑69 (MQ=255)                                                             cacaCCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCAtc  >  1:2434967/1‑69 (MQ=255)                                                             cacaCCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCAtc  >  1:3115431/1‑69 (MQ=255)                                                             cacaCCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCAtc  >  1:509471/1‑69 (MQ=255)                                                             cacaCCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCAtc  >  1:171667/1‑69 (MQ=255)                                                             cacaCCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCAtc  >  1:203508/1‑69 (MQ=255)                                                             cacaCCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCAt   >  1:2972835/1‑68 (MQ=255)                                                                |                                                                  TACGCAAGTCAGGTTTTTACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATC  >  minE/1664448‑1664576

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