Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,656,058	T→G	38.8%	A515A (GCT→GCG)	glgP →	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,656,058	0	T	G	38.8%	41.9 / 78.9	85	A515A (GCT→GCG)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (38/14);  new base G (27/6);  total (65/20)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTAGCCTGCTTAATGAGCTGCAACAACACTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGG  >  minE/2656002‑2656125                                                         |                                                                    cTTAGCCTGCTTAATGAGCTGCAACAACACTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCg                                                         <  1:1117047/69‑1 (MQ=255) cTTAGCCTGCTTAATGAGCTGCAACAACACTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCg                                                         <  1:897246/69‑1 (MQ=255)   tAGCCTGCTTAATGAGCTGCAACAACACTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGaa                                                       >  1:2640486/1‑69 (MQ=255)   tAGCCTGCTTAATGAGCTGCAACAACACTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGaa                                                       >  1:2946854/1‑69 (MQ=255)   tAGCCTGCTTAATGAGCTGCAACAACACTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGaa                                                       >  1:304815/1‑69 (MQ=255)   tAGCCTGCTTAATGAGCTGCAACAACACTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGaa                                                       >  1:522236/1‑69 (MQ=255)   tAGCCTGCTTAATGAGCTGCAACAACACTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGaa                                                       >  1:69371/1‑69 (MQ=255)                     gCAACAACACTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCg                                     >  1:2301432/1‑69 (MQ=255)                     gCAACAACACTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCg                                     >  1:3081208/1‑69 (MQ=255)                     gCAACAACACTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCg                                     <  1:3267131/69‑1 (MQ=255)                     gCAACAACACTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCg                                     <  1:1773999/69‑1 (MQ=255)                     gCAACAACACTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCg                                     <  1:1628184/69‑1 (MQ=255)                     gCAACAACACTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCg                                     >  1:413880/1‑69 (MQ=255)                        acaacaCTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTc                                   <  1:2758789/68‑1 (MQ=255)                        acaacaCTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCt                                  >  1:1077200/1‑69 (MQ=255)                        acaacaCTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCt                                  >  1:2150451/1‑69 (MQ=255)                        acaacaCTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCt                                  >  1:2097584/1‑69 (MQ=255)                        acaacaCTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCt                                  >  1:1154055/1‑69 (MQ=255)                        acaacaCTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCt                                  >  1:924233/1‑69 (MQ=255)                        acaacaCTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTc                                   >  1:2839264/1‑68 (MQ=255)                        acaacaCTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTc                                   >  1:2945234/1‑68 (MQ=255)                        acaacaCTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTc                                   >  1:2926478/1‑68 (MQ=255)                        acaacaCTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTc                                   >  1:87175/1‑68 (MQ=255)                        acaacaCTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCt                                  >  1:1704849/1‑69 (MQ=255)                        acaacaCTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCt                                  >  1:1862064/1‑69 (MQ=255)                        acaacaCTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCt                                  >  1:1143840/1‑69 (MQ=255)                        acaacaCTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCt                                  >  1:3329217/1‑69 (MQ=255)                        acaacaCTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCt                                  >  1:1031084/1‑69 (MQ=255)                        acaacaCTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCt                                  >  1:2151987/1‑69 (MQ=255)                        acaacaCTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCt                                  >  1:962803/1‑69 (MQ=255)                        acaacaCTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCt                                  >  1:750352/1‑69 (MQ=255)                          aacaCTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGa                                                        >  1:1778122/1‑45 (MQ=255)                          aacaCTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGa                                                        >  1:1392490/1‑45 (MQ=255)                          aacaCTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGa                                                        >  1:3380287/1‑45 (MQ=255)                          aacaCTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGa                                                        >  1:674746/1‑45 (MQ=255)                           acacTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTggg                               >  1:3211693/1‑68 (MQ=255)                           acacTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTgg                                >  1:93880/1‑68 (MQ=255)                           acacTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTgg                                >  1:2875227/1‑68 (MQ=255)                           acacTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTgg                                >  1:331139/1‑68 (MQ=255)                           acacTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTgg                                >  1:792298/1‑68 (MQ=255)                           acacTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTgg                                >  1:1055992/1‑68 (MQ=255)                           acacTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTgg                                >  1:1111228/1‑68 (MQ=255)                           acacTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTgg                                >  1:154464/1‑68 (MQ=255)                           acacTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTgg                                >  1:2083761/1‑68 (MQ=255)                           acacTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTgg                                >  1:2392014/1‑68 (MQ=255)                           acacTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTgg                                >  1:133441/1‑68 (MQ=255)                           acacTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTgg                                >  1:1350037/1‑68 (MQ=255)                           acacTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTgg                                >  1:1419110/1‑68 (MQ=255)                           acacTGTGATTTCCCAATGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTg                                 >  1:128445/1‑67 (MQ=255)                           acacTGTGATTTCCCAATGGTTAATCACGCGGTACATCAGGCGAAGCTGGAGAACAAAAAGCGTCTg                                 >  1:881966/1‑67 (MQ=255)                                gtgATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCagag                           >  1:386900/1‑68 (MQ=255)                                   aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGtata                       >  1:1981686/1‑69 (MQ=255)                                   aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGtata                       >  1:1484349/1‑69 (MQ=255)                                   aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGtata                       >  1:2922185/1‑69 (MQ=255)                                   aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGtata                       >  1:2849053/1‑69 (MQ=255)                                   aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGtata                       >  1:1210292/1‑69 (MQ=255)                                   aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGtat                        <  1:1482276/68‑1 (MQ=255)                                   aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGtat                        <  1:38041/68‑1 (MQ=255)                                   aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGtat                        <  1:127145/68‑1 (MQ=255)                                   aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGtat                        <  1:2789049/68‑1 (MQ=255)                                   aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGtat                        <  1:2899598/68‑1 (MQ=255)                                   aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGta                         <  1:2838915/67‑1 (MQ=255)                                   aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGta                         <  1:2081792/67‑1 (MQ=255)                                   aTTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGta                         <  1:1469446/67‑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:2799730/69‑1 (MQ=255)                                          aaTGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCGGAGTATATCGCCca                <  1:2660311/69‑1 (MQ=255)                                          aaTGGTTAATCACGCGGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCGGAGTATATCGCCca                <  1:1863381/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)                                                     aCGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGc                                      >  1:2962633/1‑36 (MQ=255)                                                     aCGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAAtg     >  1:2828716/1‑69 (MQ=255)                                                     aCGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAAtg     >  1:95883/1‑69 (MQ=255)                                                     aCGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAAtg     >  1:2321349/1‑69 (MQ=255)                                                     aCGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAAtg     >  1:704057/1‑69 (MQ=255)                                                     aCGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAAtg     >  1:2325916/1‑69 (MQ=255)                                                     aCGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAAtg     >  1:617160/1‑69 (MQ=255)                                                     aCGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAAtg     >  1:533136/1‑69 (MQ=255)                                                     aCGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAAtg     >  1:2307182/1‑69 (MQ=255)                                                     aCGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAAtg     >  1:2071701/1‑69 (MQ=255)                                                        cTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATgtgg  >  1:1101506/1‑69 (MQ=255)                                                         |                                                                    CTTAGCCTGCTTAATGAGCTGCAACAACACTGTGATTTCCCAATGGTTAATCACGCTGTGCATCAGGCGAAGCTGGAGAACAAAAAGCGTCTGGCAGAGTATATCGCCCAGCAGCTGAATGTGG  >  minE/2656002‑2656125

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