Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,656,097	A→G	37.8%	A528A (GCA→GCG)	glgP →	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,656,097	0	A	G	37.8%	36.8 / 47.0	74	A528A (GCA→GCG)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (24/22);  new base G (21/7);  total (45/29)
Fisher's exact test for biased strand distribution p-value = 8.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGGTGGTGAATCCAAAGGCGTTGTTCGATGTACAAATCAAA  >  minE/2656033‑2656163                                                                 |                                                                   gtgATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCagag                                                                 >  1:386900/1‑68 (MQ=255)    aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGtata                                                             >  1:1981686/1‑69 (MQ=255)    aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGtata                                                             >  1:1210292/1‑69 (MQ=255)    aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGtata                                                             >  1:2922185/1‑69 (MQ=255)    aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGtata                                                             >  1:2849053/1‑69 (MQ=255)    aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGtata                                                             >  1:1484349/1‑69 (MQ=255)    aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGtat                                                              <  1:2789049/68‑1 (MQ=255)    aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGtat                                                              <  1:127145/68‑1 (MQ=255)    aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGtat                                                              <  1:38041/68‑1 (MQ=255)    aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGtat                                                              <  1:1482276/68‑1 (MQ=255)    aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGtat                                                              <  1:2899598/68‑1 (MQ=255)    aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGta                                                               <  1:2081792/67‑1 (MQ=255)    aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGta                                                               <  1:2838915/67‑1 (MQ=255)    aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGta                                                               <  1:1469446/67‑1 (MQ=255)           aaTGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCGGAGTATATCGCCca                                                      <  1:1863381/69‑1 (MQ=255)           aaTGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCGGAGTATATCGCCca                                                      <  1:2481983/69‑1 (MQ=255)           aaTGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCGGAGTATATCGCCca                                                      <  1:695706/69‑1 (MQ=255)           aaTGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCGGAGTATATCGCCca                                                      <  1:1238273/69‑1 (MQ=255)           aaTGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCGGAGTATATCGCCca                                                      <  1:2660311/69‑1 (MQ=255)           aaTGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCGGAGTATATCGCCca                                                      <  1:2799730/69‑1 (MQ=255)                ttAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCcagcag                                                  >  1:3434232/1‑68 (MQ=255)                ttAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCcagcag                                                  >  1:3421938/1‑68 (MQ=255)                      aCGCTGTGCATCAGGCGAAGCTGGTGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAAtg                                           >  1:3279228/1‑69 (MQ=255)                      aCGCTGTGCATCAGGCGAAGCTGGAGAACAATAAGCGTCTGGCAGAGTatat                                                            >  1:3038749/1‑52 (MQ=255)                      aCGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAAtg                                           >  1:95883/1‑69 (MQ=255)                      aCGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAAtg                                           >  1:2828716/1‑69 (MQ=255)                      aCGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAAtg                                           >  1:704057/1‑69 (MQ=255)                      aCGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAAtg                                           >  1:617160/1‑69 (MQ=255)                      aCGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAAtg                                           >  1:2325916/1‑69 (MQ=255)                      aCGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAAtg                                           >  1:2321349/1‑69 (MQ=255)                      aCGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAAtg                                           >  1:2307182/1‑69 (MQ=255)                      aCGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAAtg                                           >  1:533136/1‑69 (MQ=255)                      aCGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAAtg                                           >  1:2071701/1‑69 (MQ=255)                         cTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATgtgg                                        >  1:1101506/1‑69 (MQ=255)                              cATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATgtggtggtg                                   >  1:2046804/1‑69 (MQ=255)                                tCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGTTGGTGa                                  <  1:6729/68‑1 (MQ=255)                                tCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGGTGGTGa                                  <  1:3416923/68‑1 (MQ=255)                                tCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGGTGGTGa                                  <  1:3315385/68‑1 (MQ=255)                                tCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGGTGGTGa                                  <  1:1088664/68‑1 (MQ=255)                                tCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGGTGGTGa                                  <  1:726247/68‑1 (MQ=255)                                tCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGGTGGTGa                                  <  1:2940489/68‑1 (MQ=255)                                tCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGGTGGTGa                                  <  1:3192565/68‑1 (MQ=255)                                tCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGGTGGTGa                                  <  1:2358420/68‑1 (MQ=255)                                  aGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGGTGGTGAAt                                >  1:2241196/1‑68 (MQ=255)                                  aGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGGTGGTGAAt                                >  1:696091/1‑68 (MQ=255)                                  aGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGGTGGTGAAt                                >  1:1838106/1‑68 (MQ=255)                                  aGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGGTGGTGAAt                                >  1:3224737/1‑68 (MQ=255)                                        aGCTGGAGAACAAAAAGCGTCTGGCGGAGTATATCGCCca                                                      <  1:2636298/40‑1 (MQ=255)                                                        gCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGGTTGTGAATCCAAAGGCCTTGTTCGATGTACa         <  1:1558252/69‑1 (MQ=255)                                                        gCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGGTGGTGAATCCAAAGGCGTTGTTCGATGTACa         <  1:1453212/69‑1 (MQ=255)                                                        gCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGGTGGTGAATCCAAAGGCGTTGTTCGATGTACa         <  1:640190/69‑1 (MQ=255)                                                        gCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGGTGGTGAATCCAAAGGCGTTGTTCGATGTACa         <  1:1963316/69‑1 (MQ=255)                                                        gCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGGTGGTGAATCCAAAGGCGTTGTTCGATGTACa         <  1:2303322/69‑1 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCATCTGAATGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:2743777/1‑69 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCAGCTGATTGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:2797217/1‑69 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCAGCTGAATGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:777698/1‑69 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCAGCTGAATGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:810044/1‑69 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCAGCTGAATGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:816873/1‑69 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCAGCTGAATGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:1179001/1‑69 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCAGCTGAATGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:718403/1‑69 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCAGCTGAATGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:1288009/1‑69 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCAGCTGAATGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:1371023/1‑69 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCAGCTGAATGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:646432/1‑69 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCAGCTGAATGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:1019726/1‑69 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCAGCTGAATGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:549290/1‑69 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCAGCTGAATGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:1507586/1‑69 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCAGCTGAATGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:356017/1‑69 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCAGCTGAATGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:1685423/1‑69 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCAGCTGAATGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:1762229/1‑69 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCAGCTGAATGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:1814207/1‑69 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCAGCTGAATGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:1821021/1‑69 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCAGCTGAATGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:3024364/1‑69 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCAGCTGAATGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:3024235/1‑69 (MQ=255)                                                             tGGCGGAGTATATCGCCCAGCAGCTGAATGTGGTGGTAAATCCGAAAGCGCTGTTCGATGTGCAAATCa    >  1:1968227/1‑69 (MQ=255)                                                             tGGCAGAGTATATCGCCCAGCAGCTGAATGTGGTGGTGAATCCAAAGGCGTTGTTCGATGTACAAATCa    <  1:1669602/69‑1 (MQ=255)                                                                 agagTATATCGCCCAGCAGCTGAATGTGGTGGTGAATCCAAAGGCGTTGTTCGATGTACAAATCaaa  <  1:505165/67‑1 (MQ=255)                                                                 agagTATATCGCCCAGCAGCTGAATGTGGTGGTGAATCCAAAGGCGTTGTTCGATGTACAAATCaaa  <  1:1445452/67‑1 (MQ=255)                                                                 agagTATATCGCCCAGCAGCTGAATGTGGTGGTGAATCCAAAGGCGTTGTTCGATGTACAAATCaaa  <  1:1965852/67‑1 (MQ=255)                                                                 |                                                                   GTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGGTGGTGAATCCAAAGGCGTTGTTCGATGTACAAATCAAA  >  minE/2656033‑2656163

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