Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,093,769	A→G	8.3%	S361S (AGT→AGC)	tdcG ←	L‑serine dehydratase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,093,769	0	A	G	8.3%	128.6 / 2.8	48	S361S (AGT→AGC)	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base A (23/21);  new base G (2/2);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

ATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAG  >  minE/2093706‑2093838                                                                |                                                                      aTCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTgccgc                                                                <  1:3127635/71‑1 (MQ=255)    gCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCa                                                             <  1:3822503/71‑1 (MQ=255)     cAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCag                                                            <  1:3440753/71‑1 (MQ=255)         tCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCagtagt                                                        <  1:3609381/71‑1 (MQ=255)           aGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGtt                                                       <  1:577449/70‑1 (MQ=255)               cAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGt                                                   >  1:556013/1‑70 (MQ=255)                aaGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTa                                                 <  1:1068627/71‑1 (MQ=255)                aaGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTa                                                 <  1:1061863/71‑1 (MQ=255)                 aGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTaa                                                <  1:3242207/71‑1 (MQ=255)                  ggTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAg                                                    >  1:1757215/1‑66 (MQ=255)                  ggTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAg                                                    >  1:80479/1‑66 (MQ=255)                  ggTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAg                                                    >  1:2639012/1‑66 (MQ=255)                  ggTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAg                                                    >  1:2795890/1‑66 (MQ=255)                   gTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                              <  1:4757352/71‑1 (MQ=255)                   gTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                              <  1:829348/71‑1 (MQ=255)                            cATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgccg                                     >  1:5469700/1‑71 (MQ=255)                             aTCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgccgc                                    >  1:3856781/1‑71 (MQ=255)                              tCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgccgc                                    >  1:2843138/1‑70 (MQ=255)                               cgcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCagta                                                          >  1:4412515/1‑47 (MQ=255)                                gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGtt                                                  >  1:682396/1‑54 (MQ=255)                                gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                              >  1:1375060/1‑58 (MQ=255)                                gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                              >  1:5362224/1‑58 (MQ=255)                                gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                              >  1:2489800/1‑58 (MQ=255)                                gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCAt                                 >  1:849638/1‑71 (MQ=255)                                  gATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCAtt                                >  1:4173286/1‑70 (MQ=255)                                  gATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTg                               <  1:4048256/71‑1 (MQ=255)                                   aTTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTc                                                      <  1:1469940/47‑1 (MQ=255)                                    tttCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgccgc                                    <  1:2904729/64‑1 (MQ=255)                                       ccgccgCGATGCACACCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCATTGagc                            <  1:4323993/69‑3 (MQ=255)                                       ccgccgCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACa                           <  1:161852/70‑1 (MQ=255)                                        cgccgcGATGCACACCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCATTGAGCacg                         >  1:4570746/1‑69 (MQ=255)                                             cATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCa                                  >  1:2495466/1‑57 (MQ=255)                                               ttGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGc                   >  1:4963297/1‑70 (MQ=255)                                                tGCACACCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCATTGAGCAcgccacgcca                 <  1:543722/71‑1 (MQ=255)                                                 gCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAg                                                    >  1:2572045/1‑35 (MQ=255)                                                 gCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAAc                            >  1:1229003/1‑59 (MQ=255)                                                 gCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCa                 <  1:4668797/70‑1 (MQ=255)                                                 gCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCa                 <  1:4730031/70‑1 (MQ=255)                                                     aCCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCATTGAGCACGCCACGCCaactt             >  1:3027797/1‑67 (MQ=255)                                                       cTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTccc          >  1:4067880/1‑71 (MQ=255)                                                          cgcCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAAtc              >  1:5441428/1‑64 (MQ=255)                                                          cgcCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCt       <  1:756842/71‑1 (MQ=255)                                                          cgcCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCt       <  1:4555033/71‑1 (MQ=255)                                                            ccGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGa     <  1:663082/71‑1 (MQ=255)                                                            ccGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGa     <  1:344475/71‑1 (MQ=255)                                                              ggACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGaa                             >  1:3686379/1‑45 (MQ=255)                                                              ggACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACa   <  1:1509830/71‑1 (MQ=255)                                                               gACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAg  >  1:258567/1‑71 (MQ=255)                                                                |                                                                      ATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAG  >  minE/2093706‑2093838

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