Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	163,666	C→T	29.5%	E77K (GAA→AAA)	yaeI ←	predicted phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	163,666	0	C	T	29.5%	92.0 / 47.6	61	E77K (GAA→AAA)	yaeI	predicted phosphatase
Reads supporting (aligned to +/- strand):  ref base C (31/12);  new base T (7/11);  total (38/23)
Fisher's exact test for biased strand distribution p-value = 2.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATGGAG  >  minE/163596‑163726                                                                       |                                                             gCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTAGTTc                                                              <  1:2585603/71‑1 (MQ=255)       aaaTTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGc                                                         >  1:1101008/1‑70 (MQ=255)       aaaTTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGc                                                         >  1:1734749/1‑70 (MQ=255)            cAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAagag                                                   >  1:1184085/1‑71 (MQ=255)            cAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAagag                                                   >  1:3713947/1‑71 (MQ=255)            cAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAagag                                                   >  1:1111912/1‑71 (MQ=255)             aGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAagag                                                   <  1:1152206/70‑1 (MQ=255)             aGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAagag                                                   <  1:2927298/70‑1 (MQ=255)             aGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAaga                                                    <  1:3968234/69‑1 (MQ=255)                              taCGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCaa                                                      <  1:2590745/50‑1 (MQ=255)                              taCGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCaa                                                      <  1:473169/50‑1 (MQ=255)                               aCGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAg                                >  1:2759506/1‑71 (MQ=255)                               aCGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAg                                >  1:2678093/1‑71 (MQ=255)                               aCGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAg                                >  1:3495684/1‑71 (MQ=255)                               aCGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAg                                >  1:1685732/1‑71 (MQ=255)                                cGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAgg                               >  1:2180138/1‑71 (MQ=255)                                cGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAg                                >  1:3638736/1‑70 (MQ=255)                                   aaTCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCa                                                 <  1:1249998/50‑1 (MQ=255)                                       gccgccCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAg                        >  1:3952539/1‑71 (MQ=255)                                               gTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAgt                       <  1:1990896/64‑2 (MQ=255)                                               gTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAgg                       <  1:3951739/64‑1 (MQ=255)                                               gTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAgg                       <  1:626415/64‑1 (MQ=255)                                               gTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAgg                       <  1:3331700/64‑1 (MQ=255)                                               gTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAgg                       <  1:1498881/64‑1 (MQ=255)                                                tAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAAc               <  1:3475199/71‑1 (MQ=255)                                                tAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAAc               <  1:2582739/71‑1 (MQ=255)                                                tAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAAc               <  1:1632991/71‑1 (MQ=255)                                                 aaTATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGaa                     >  1:162510/1‑64 (MQ=255)                                                 aaTATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGaa                     >  1:3980006/1‑64 (MQ=255)                                                  atatCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTaaa                         >  1:3550537/1‑59 (MQ=255)                                                  atatCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTaaa                         >  1:290834/1‑59 (MQ=255)                                                  atatCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAAc               >  1:149855/1‑69 (MQ=255)                                                  atatCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAAc               >  1:3157280/1‑69 (MQ=255)                                                  atatCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACga             <  1:631577/71‑1 (MQ=255)                                                    atCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACga             >  1:3091541/1‑69 (MQ=255)                                                    atCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACga             >  1:3707690/1‑69 (MQ=255)                                                     tCAAATCGGTCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGa                      <  1:1383118/59‑1 (MQ=255)                                                     tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGa                      <  1:378359/59‑1 (MQ=255)                                                     tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGa                      <  1:3578750/59‑1 (MQ=255)                                                     tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGa                      <  1:2366021/59‑1 (MQ=255)                                                       aaaTCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaaaa                >  1:918438/1‑63 (MQ=255)                                                       aaaTCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaaaa                >  1:729908/1‑63 (MQ=255)                                                       aaaTCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaaaa                >  1:2401552/1‑63 (MQ=255)                                                       aaaTCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaaaa                >  1:2679120/1‑63 (MQ=255)                                                           cGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAAt                                   >  1:2718051/1‑40 (MQ=255)                                                           cGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAAt                                   >  1:326061/1‑40 (MQ=255)                                                           cGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAAt                                   >  1:2475682/1‑40 (MQ=255)                                                           cGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACga             >  1:2275987/1‑62 (MQ=255)                                                           cGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACga             >  1:3816006/1‑62 (MQ=255)                                                           cGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATg     >  1:761612/1‑70 (MQ=255)                                                              gCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATGgag  <  1:1150470/70‑1 (MQ=255)                                                               cTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAGACGAGAGTAATg     >  1:1846652/1‑66 (MQ=255)                                                               cTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATg     >  1:3676559/1‑66 (MQ=255)                                                                ttttGTTTTATGCCAAGAGCAATCGCGTCAGAAATc                                  >  1:739290/1‑36 (MQ=255)                                                                ttttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATg     >  1:3393171/1‑65 (MQ=255)                                                                ttttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATg     >  1:1699542/1‑65 (MQ=255)                                                                ttttGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaaaa                >  1:2145173/1‑54 (MQ=255)                                                                ttttGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaaaa                >  1:3882408/1‑54 (MQ=255)                                                                ttttGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaaaa                >  1:2214239/1‑54 (MQ=255)                                                                 tttGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt         <  1:3914370/60‑1 (MQ=255)                                                                 tttGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt         <  1:2163464/60‑1 (MQ=255)                                                                  ttGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt         <  1:2474623/59‑1 (MQ=255)                                                                       cTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACga             >  1:3085418/1‑50 (MQ=255)                                                                       cTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACga             >  1:451579/1‑50 (MQ=255)                                                                       cTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACga             >  1:3196175/1‑50 (MQ=255)                                                                       cTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACga             >  1:2409499/1‑50 (MQ=255)                                                                       |                                                             GCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATGGAG  >  minE/163596‑163726

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