Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	349,124	C→T	47.5%	A474A (GCC→GCT)	gcl →	glyoxylate carboligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	349,124	0	C	T	47.5%	12.5 / 105.3	80	A474A (GCC→GCT)	gcl	glyoxylate carboligase
Reads supporting (aligned to +/- strand):  ref base C (21/21);  new base T (19/19);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGA  >  minE/349055‑349194                                                                      |                                                                       ccTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGcc                                                                        >  1:592674/1‑70 (MQ=255)   tGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTAt                                                                     >  1:2461500/1‑71 (MQ=255)     aTTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTAt                                                                     <  1:3513274/69‑1 (MQ=255)     aTTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTAt                                                                     <  1:3508938/69‑1 (MQ=255)     aTTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCt                                                                   >  1:1400423/1‑71 (MQ=255)      ttGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTg                                                                  >  1:1972569/1‑71 (MQ=255)      ttGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTg                                                                  >  1:3124504/1‑71 (MQ=255)      ttGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTg                                                                  >  1:3026600/1‑71 (MQ=255)      ttGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTg                                                                  >  1:653184/1‑71 (MQ=255)      ttGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATc                                                                    >  1:2984719/1‑69 (MQ=255)      ttGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTg                                                                  <  1:536380/71‑1 (MQ=255)       tGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTgg                                                                 >  1:1645846/1‑71 (MQ=255)       tGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTgg                                                                 >  1:3710021/1‑71 (MQ=255)         aaGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTgggg                                                               <  1:1668162/71‑1 (MQ=255)           gagTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCt                                                             <  1:1701628/71‑1 (MQ=255)           gagTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCt                                                             <  1:3363104/71‑1 (MQ=255)            agTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTAt                                                                     >  1:3210716/1‑62 (MQ=255)                aGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTc                                                        >  1:781262/1‑71 (MQ=255)                aGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCt                                                             >  1:717472/1‑66 (MQ=255)                aGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTc                                                        >  1:686891/1‑71 (MQ=255)                aGCTGTAGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTc                                                        >  1:1861724/1‑71 (MQ=255)                 gCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTc                                                        <  1:3653614/70‑1 (MQ=255)                 gCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTc                                                        <  1:2119516/70‑1 (MQ=255)                 gCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCg                                                       <  1:2891967/71‑1 (MQ=255)                    gTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgt                                                      >  1:1889077/1‑69 (MQ=255)                    gTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgt                                                      >  1:2842844/1‑69 (MQ=255)                     ttGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTc                                                        <  1:2174968/66‑1 (MQ=255)                     ttGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgtca                                                    <  1:3382189/70‑1 (MQ=255)                     ttGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCgt                                                      <  1:461249/68‑1 (MQ=255)                      tGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCgtcagt                                                  >  1:349141/1‑71 (MQ=255)                       ggCTCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCgtcagt                                                  <  1:1274314/70‑1 (MQ=255)                           cgcAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACaa                                             >  1:3628889/1‑71 (MQ=255)                            gcAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCt                                                                   >  1:3804572/1‑48 (MQ=255)                             cAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTcaca                                              <  1:2362117/68‑1 (MQ=255)                              aGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAAcgc                                          <  1:1526673/71‑1 (MQ=255)                                ttCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACgcgc                                        <  1:2675962/71‑1 (MQ=255)                                ttCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACgcgc                                        >  1:1259330/1‑71 (MQ=255)                                ttCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCt                                                             <  1:2980824/50‑1 (MQ=255)                                  cAACATTCCGTACATCCATGTGCTGGTCAACAACGCCt                                                                       >  1:470856/1‑38 (MQ=255)                                  cAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCtt                                      >  1:2153611/1‑71 (MQ=255)                                  cAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCtt                                      >  1:1474106/1‑71 (MQ=255)                                  cAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCt                                       >  1:2158031/1‑70 (MQ=255)                                   aaCATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgtcag                                                   >  1:3591573/1‑57 (MQ=255)                                   aaCATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCtt                                      <  1:2390789/70‑1 (MQ=255)                                    aCATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgt                                                      <  1:1484061/53‑1 (MQ=255)                                    aCATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCtttt                                    <  1:3722785/71‑1 (MQ=255)                                     cATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCgtcagt                                                  <  1:3290951/56‑1 (MQ=255)                                     cATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCgtcag                                                   <  1:3238978/55‑1 (MQ=255)                                      aTTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCg                                                       >  1:3524008/1‑50 (MQ=255)                                      aTTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGa                                  <  1:2100092/71‑1 (MQ=255)                                      aTTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGa                                  <  1:1758483/71‑1 (MQ=255)                                       ttCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCtt                                      <  1:1044787/66‑1 (MQ=255)                                       ttCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCtt                                      <  1:3470517/66‑1 (MQ=255)                                       ttCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGAc                                 <  1:2581270/71‑1 (MQ=255)                                        tCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACa                                >  1:3170895/1‑71 (MQ=255)                                        tCCGTACATCCATGTGCTGGTCAACAACGCCTATCt                                                                   <  1:3756450/36‑1 (MQ=255)                                        tCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACa                                >  1:1230859/1‑71 (MQ=255)                                         ccGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTcac                                               >  1:2138521/1‑55 (MQ=255)                                          cGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACAt                               <  1:1026065/70‑1 (MQ=255)                                          cGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATg                              <  1:3722309/71‑1 (MQ=255)                                          cGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACAt                               <  1:1039802/70‑1 (MQ=255)                                            tACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGAc                                 <  1:3040581/66‑1 (MQ=255)                                            tACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGAc                                 >  1:2072466/1‑66 (MQ=255)                                                  cATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCg                                                       <  1:1703767/38‑1 (MQ=255)                                                   aTGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGc                     <  1:3634760/71‑1 (MQ=255)                                                    tgtgCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGactac                        >  1:875291/1‑67 (MQ=255)                                                    tgtgCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGactac                        >  1:2743526/1‑67 (MQ=255)                                                     gtgCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGt                   >  1:3125509/1‑71 (MQ=255)                                                     gtgCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGt                   >  1:2217545/1‑71 (MQ=255)                                                      tgCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGAc                                 >  1:889735/1‑56 (MQ=255)                                                      tgCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGt                   <  1:450060/70‑1 (MQ=255)                                                      tgCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTg                  >  1:1851907/1‑71 (MQ=255)                                                           gTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCaa               <  1:1560225/69‑1 (MQ=255)                                                           gTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACt             >  1:3698661/1‑71 (MQ=255)                                                           gTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGt                   <  1:2758116/65‑1 (MQ=255)                                                           gTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACt             <  1:3167552/71‑1 (MQ=255)                                                              aacaacGCTTATCTGGGGCTGATTCGTCAGTCACAAc                                            >  1:3036686/1‑37 (MQ=255)                                                              aacaacGCTTATCTGGGGCTGATTCGTCAGTCACAAc                                            >  1:3236261/1‑37 (MQ=255)                                                                caacGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGc          <  1:1383317/69‑1 (MQ=255)                                                                  acGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCttt       <  1:2104215/70‑1 (MQ=255)                                                                  acGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTc      >  1:3334725/1‑71 (MQ=255)                                                                     ccTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGa                                  <  1:518108/40‑1 (MQ=255)                                                                      cTATCTGGGGCTGATTCGTCAGTCACAACGCGCtttt                                    >  1:1911451/1‑37 (MQ=255)                                                                      cTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCgaga  <  1:2076159/71‑1 (MQ=255)                                                                      |                                                                       CCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGA  >  minE/349055‑349194

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