Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,125,208	G→A	61.9%	A52A (GCG→GCA)	yhbQ →	predicted endonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,125,208	0	G	A	61.9%	37.1 / 49.1	76	A52A (GCG→GCA)	yhbQ	predicted endonuclease
Reads supporting (aligned to +/- strand):  ref base G (8/21);  new base A (17/30);  total (25/51)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

GCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCGTTCGCTGGCGTTACGGGCGGAATATCGCG  >  minE/2125142‑2125251                                                                   |                                            gcaCCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGcc                                           <  1:2476492/69‑1 (MQ=255) gcaCCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGcc                                           <  1:2194983/69‑1 (MQ=255) gcaCCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGcc                                           <  1:2824390/69‑1 (MQ=255) gcaCCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGcc                                           <  1:1644014/69‑1 (MQ=255)  caCCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGcc                                           <  1:14030/68‑1 (MQ=255)         gCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGCg                                   <  1:2827989/69‑1 (MQ=255)         gCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGCg                                   <  1:1500355/69‑1 (MQ=255)         gCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGCg                                   <  1:3291631/69‑1 (MQ=255)           ggCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGcc                                           <  1:2518876/59‑1 (MQ=255)             cAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCg                               <  1:3299810/69‑1 (MQ=255)             cAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCg                               <  1:1015388/69‑1 (MQ=255)             cAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCg                               <  1:2595519/69‑1 (MQ=255)             cAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCg                               <  1:2192851/69‑1 (MQ=255)             cAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCg                               <  1:2866096/69‑1 (MQ=255)             cAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCg                               <  1:2074028/69‑1 (MQ=255)             cAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCg                               <  1:996641/69‑1 (MQ=255)             cAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCg                               <  1:1335378/69‑1 (MQ=255)             cAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCg                               <  1:140205/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCa                                          >  1:8778/1‑46 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:581764/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   >  1:975997/1‑69 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   >  1:905326/1‑69 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:2557751/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:2579985/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:2776617/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:739373/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:716493/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:383327/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:289391/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   >  1:2968209/1‑69 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:3071744/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   >  1:3103140/1‑69 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:422974/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   >  1:408656/1‑69 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   >  1:2517033/1‑69 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:336920/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:3385887/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:3411335/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:2410044/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   >  1:1136730/1‑69 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:123437/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   >  1:1284761/1‑69 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   >  1:1467095/1‑69 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   >  1:1539620/1‑69 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:1539906/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   >  1:1569341/1‑69 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:1652695/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:1846510/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:2026424/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:2083461/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   >  1:2171894/1‑69 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:2510464/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:2219726/69‑1 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   >  1:2134749/1‑69 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   >  1:2158392/1‑69 (MQ=255)                         aGCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGt                    >  1:2534507/1‑68 (MQ=255)                         aGCACTGCGAGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   >  1:1459933/1‑69 (MQ=255)                          gCACTGCGCGGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:79587/68‑1 (MQ=255)                                  cgGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGTTACGGGca            <  1:3348938/67‑2 (MQ=255)                                  cgGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGTTACGGGca            <  1:789619/67‑2 (MQ=255)                                  cgGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGTTACGGGca            <  1:2091692/67‑2 (MQ=255)                                  cgGGAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGAGCGTTCGCTGGCGTTACGGGca            <  1:1719608/67‑2 (MQ=255)                                   gggAAGGGAGAACTAACGCTGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGTTACGGGca            <  1:3159536/66‑2 (MQ=255)                                   gggAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCGTTCGCTGGCGTTACGGGCGGaa         <  1:3310126/69‑1 (MQ=255)                                   gggAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCGTTCGCTGGCGTTACGGGCGGaa         <  1:2670418/69‑1 (MQ=255)                                    ggAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCGTTCGCTGGCGTTACGGGCGGaa         <  1:428741/68‑1 (MQ=255)                                     gAAGGGAGAACTAACGCAGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:2147224/57‑1 (MQ=255)                                          gagaACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCGTTCGCTGGCGTTACGGGCGGAATATcgc   >  1:796888/1‑68 (MQ=255)                                          gagaACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCGTTCGCTGGCGTTACGGGCGGAATATcgc   >  1:1712180/1‑68 (MQ=255)                                           agaACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCGTTCGCTGGCGTTACGGGCGGAATATcgcg  >  1:2874244/1‑68 (MQ=255)                                           agaACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCGTTCGCTGGCGTTACGGGCGGAATATcgcg  >  1:1411817/1‑68 (MQ=255)                                           agaACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCGTTCGCTGGCGTTACGGGCGGAATATcgcg  >  1:3399992/1‑68 (MQ=255)                                           agaACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCGTTCGCTGGCGTTACGGGCGGAATATcgcg  >  1:232490/1‑68 (MQ=255)                                           agaACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCGTTCGCTGGCGTTACGGGCGGAATATcgcg  >  1:1177361/1‑68 (MQ=255)                                           agaACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCGTTCGCTGGCGTTACGGGCGGAATATcgcg  >  1:2119987/1‑68 (MQ=255)                                                      tGGCGTTTTCCGCACCAGTGGGCGATCGTTCGCTGGCGtt                   <  1:228266/40‑1 (MQ=38)                                                                cgcgCCAGTGGGCGATCGTTCGCTGGCGTTACGGGCGGaa         <  1:216125/40‑1 (MQ=255)                                                                   |                                            GCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGATCGTTCGCTGGCGTTACGGGCGGAATATCGCG  >  minE/2125142‑2125251

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